diff --git a/NumLib/ODESolver/NonlinearSolver.cpp b/NumLib/ODESolver/NonlinearSolver.cpp index 9383e04f458a5969fa87241e34f72f3e7ce88e3d..50bd6dd53a3b104974cc707333ebbaa4feaf7fa0 100644 --- a/NumLib/ODESolver/NonlinearSolver.cpp +++ b/NumLib/ODESolver/NonlinearSolver.cpp @@ -34,6 +34,8 @@ void NonlinearSolver<NonlinearSolverTag::Picard>:: return; } + INFO("Calculate non-equilibrium initial residuum."); + auto& A = NumLib::GlobalMatrixProvider::provider.getMatrix(_A_id); auto& rhs = NumLib::GlobalVectorProvider::provider.getVector(_rhs_id); _equation_system->assemble(x, x_prev, process_id); @@ -235,6 +237,8 @@ void NonlinearSolver<NonlinearSolverTag::Newton>:: return; } + INFO("Calculate non-equilibrium initial residuum."); + _equation_system->assemble(x, x_prev, process_id); _r_neq = &NumLib::GlobalVectorProvider::provider.getVector(_r_neq_id); _equation_system->getResidual(*x[process_id], *x_prev[process_id], *_r_neq); diff --git a/ProcessLib/TimeLoop.cpp b/ProcessLib/TimeLoop.cpp index 228d6b29578937ea36b8aa870c8c4170cf21cf1a..20e7bb31e67f60ebc95409acd26f90b753e5be2a 100644 --- a/ProcessLib/TimeLoop.cpp +++ b/ProcessLib/TimeLoop.cpp @@ -209,7 +209,6 @@ void calculateNonEquilibriumInitialResiduum( process_solutions, std::vector<GlobalVector*> const& process_solutions_prev) { - INFO("Calculate non-equilibrium initial residuum."); for (auto& process_data : per_process_data) { auto& time_disc = *process_data->time_disc;