From 4d9c49d0d366efc4d6d821ce63c4d743b8f749c8 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 14 Mar 2012 15:10:34 +0100
Subject: [PATCH] Only enable FPI if the state shall evolve
Use size_t
---
src/one-body-sample.cc | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 93bad145..3efbab56 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -236,8 +236,11 @@ int main(int argc, char *argv[]) {
u4_diff[i][0] /= refinement_factor;
}
- for (int state_fpi = 0;
- state_fpi < parset.get<int>("solver.tnnmg.fixed_point_iterations");
+ auto const fixed_point_iterations =
+ parset.get<bool>("boundary.friction.state.evolve")
+ ? parset.get<size_t>("solver.tnnmg.fixed_point_iterations")
+ : 1;
+ for (size_t state_fpi = 0; state_fpi < fixed_point_iterations;
++state_fpi) {
auto myGlobalNonlinearity =
assemble_nonlinearity<VectorType, OperatorType>(
--
GitLab