diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 64916565afe242c8254b8c64bac059097975f29d..676c7317d66d4c243d65e63c31d06977aee4a0c2 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -374,10 +374,11 @@ int main(int argc, char *argv[]) { stateUpdater->setup(tau); timeSteppingScheme->setup(ell, tau, time, problem_rhs, problem_iterate, problem_A); - EnergyNorm<MatrixType, VectorType> energyNorm(problem_A); + EnergyNorm<MatrixType, VectorType> velocityEnergyNorm(problem_A); LoopSolver<VectorType> overallSolver( multigridStep, parset.get<size_t>("solver.tnnmg.maxiterations"), - solver_tolerance, &energyNorm, verbosity, false); // absolute error + solver_tolerance, &velocityEnergyNorm, verbosity, + false); // absolute error size_t iterationCounter; auto solveDisplacementProblem = [&](VectorType &_problem_iterate, @@ -435,7 +436,7 @@ int main(int argc, char *argv[]) { std::cerr.flush(); } if (state_fpi > 1) { - double const correctionNorm = energyNorm.diff(u_saved, u); + double const correctionNorm = velocityEnergyNorm.diff(u_saved, u); if (correctionNorm < fixedPointTolerance) { if (parset.get<bool>("printProgress")) { std::cerr << '#';