Skip to content
Snippets Groups Projects
Commit fc5fd4c5 authored by Markus Blatt's avatar Markus Blatt
Browse files

[Bugfix] Corrects current residual after postsmoothing.

Previously the current residual was not updated after the coarse
grid correction. This patch fixes this by storing the only the
prolongated coarse grid correction in context.lhs and updating the
accumulated update. During postsmoothing the currently residual is
updated according to context.lhs.
parent 3e8db323
No related branches found
No related tags found
No related merge requests found
......@@ -402,9 +402,10 @@ public:
policy_->moveToCoarseLevel(*context.rhs);
InverseOperatorResult res;
coarseSolver_->apply(policy_->getCoarseLevelLhs(), policy_->getCoarseLevelRhs(), res);
policy_->moveToFineLevel(*context.update);
// Postsmoothing
*context.lhs=0;
policy_->moveToFineLevel(*context.lhs);
*context.update += *context.lhs;
// Postsmoothing
postsmooth(context, postSteps_);
}
......
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