Skip to content
Snippets Groups Projects
Commit 171fe064 authored by Elias Pipping's avatar Elias Pipping
Browse files

[Cleanup] Make logic easier to understand

parent b94c63ee
No related branches found
No related tags found
No related merge requests found
...@@ -74,16 +74,16 @@ bool AdaptiveTimeStepper<Factory, Updaters, ErrorNorm>::coarsen() { ...@@ -74,16 +74,16 @@ bool AdaptiveTimeStepper<Factory, Updaters, ErrorNorm>::coarsen() {
template <class Factory, class Updaters, class ErrorNorm> template <class Factory, class Updaters, class ErrorNorm>
void AdaptiveTimeStepper<Factory, Updaters, ErrorNorm>::refine() { void AdaptiveTimeStepper<Factory, Updaters, ErrorNorm>::refine() {
while (true) { while (true) {
Updaters F2 = current_.clone(); Updaters F1 = current_.clone();
Updaters F1; Updaters F2;
{ {
MyCoupledTimeStepper coupledTimeStepper(finalTime_, factory_, parset_, MyCoupledTimeStepper coupledTimeStepper(finalTime_, factory_, parset_,
globalFriction_, F2, errorNorm_, globalFriction_, F1, errorNorm_,
externalForces_); externalForces_);
stepAndReport("F1", coupledTimeStepper, relativeTime_, stepAndReport("F1", coupledTimeStepper, relativeTime_,
relativeTau_ / 2.0); relativeTau_ / 2.0);
F1 = F2.clone(); F2 = F1.clone();
stepAndReport("F2", coupledTimeStepper, stepAndReport("F2", coupledTimeStepper,
relativeTime_ + relativeTau_ / 2.0, relativeTau_ / 2.0); relativeTime_ + relativeTau_ / 2.0, relativeTau_ / 2.0);
} }
......
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