Skip to content
Snippets Groups Projects
Commit 9e5d2428 authored by Sander, Oliver's avatar Sander, Oliver
Browse files

Merge branch 'smallRTRFix' into 'master'

Small change/fix in Mixed Riemanian Trust region solver

See merge request !69
parents 051f34ea a82357d3
No related branches found
No related tags found
1 merge request!69Small change/fix in Mixed Riemanian Trust region solver
Pipeline #5779 passed
......@@ -414,8 +414,10 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis,Basis0,TargetSpace0,Basis1,
double oldEnergy = 0;
Dune::Timer solutionTimer;
int ii = 0;
CorrectionType diff{corr_global};
try {
for (; ii<innerIterations_; ii++) {
diff=corr_global;
residual[_0] = rhs_global[_0];
stiffnessMatrix[_0][_1].mmv(corr_global[_1], residual[_0]);
mmgStep0->setRhs(residual[_0]);
......@@ -436,7 +438,8 @@ void MixedRiemannianTrustRegionSolver<GridType,Basis,Basis0,TargetSpace0,Basis1,
std::cout << "Warning: energy increase!" << std::endl;
oldEnergy = energy;
if (corr_global.infinity_norm() < innerTolerance_)
diff -= corr_global;
if (diff.infinity_norm() < innerTolerance_)
break;
}
} catch (Dune::Exception &e) {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment