From fed4c2ae32b744a6da3e7d14723c74ea51243145 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Thu, 10 May 2012 15:43:44 +0200
Subject: [PATCH] Do not loop over parset.get needlessly
---
src/one-body-sample.cc | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index b5b824d2..adb3d148 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -226,6 +226,12 @@ int main(int argc, char *argv[]) {
double const h = 1.0 / timesteps;
auto const &neumannFunction = functions.get("neumannCondition");
auto const &dirichletFunction = functions.get("dirichletCondition");
+ auto const L = parset.get<double>("boundary.friction.ruina.L");
+ auto const a = parset.get<double>("boundary.friction.ruina.a");
+ auto const b = parset.get<double>("boundary.friction.ruina.b");
+ auto const eta = parset.get<double>("boundary.friction.eta");
+ auto const mu = parset.get<double>("boundary.friction.mu");
+
for (size_t run = 1; run <= timesteps; ++run) {
if (parset.get<bool>("printProgress")) {
std::cout << '*';
@@ -272,7 +278,6 @@ int main(int argc, char *argv[]) {
for (size_t i = 0; i < frictionalNodes.size(); ++i) {
if (frictionalNodes[i][0]) {
- auto const L = parset.get<double>("boundary.friction.ruina.L");
double const unorm = u4_diff[i].two_norm();
// the (logarithmic) steady state corresponding to the current
@@ -310,11 +315,6 @@ int main(int argc, char *argv[]) {
double const V = u4_diff[i].two_norm();
double const state = (*s4_new)[i];
- auto const L = parset.get<double>("boundary.friction.ruina.L");
- auto const a = parset.get<double>("boundary.friction.ruina.a");
- auto const b = parset.get<double>("boundary.friction.ruina.b");
- auto const eta = parset.get<double>("boundary.friction.eta");
- auto const mu = parset.get<double>("boundary.friction.mu");
coefficient_writer << mu + a *std::log(V * eta) +
b * (state - std::log(eta * L))
--
GitLab