Commit 2c48fb54 authored by GILLES Sebastien's avatar GILLES Sebastien
Browse files

#820 Put back the Newton FP loop (only one call was done previously...).

parent e5c426d5
......@@ -71,7 +71,19 @@ namespace HappyHeart
auto& newton_fp_index = variable_holder.newton_fp_index;
newton_fp_index = 0u;
current_norm = inner_varf.Perform();
double absolute_error = 1.e20;
double relative_error = 1.e20;
while (absolute_error > 1.e-10 && relative_error > 1.e-6 && newton_fp_index < 2u) // \todo #820 Should not remain hardcoded!
// \todo #820 iteration index sup value is obviously temporary for dev purposes...
// \todo Absolute and relative error should be recomputed!
{
std::cout << "NEWTON FP LOOP" << std::endl;
current_norm = inner_varf.Perform();
++newton_fp_index;
}
// Update fluidmass, velocity fluid and pressure.
Wrappers::Petsc::MatMult(inner_varf.GetMonolithic2Mass().GetInterpolationMatrix(),
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment