Homec4science
Diffusion SpecMiCP / ReactMiCP (cython) Browse

specmicp/src/reactmicp/solver/runner.cpp
a1bb082f067bcython

Last Change
Raw File

runner.cpp
View Options

#include "runner.hpp"
#include "utils/io/reactive_transport.hpp"
#include "utils/timer.hpp"
#include <fstream>
#include "staggers_base/variables_base.hpp"
namespace specmicp {
namespace reactmicp {
namespace solver {
void ReactiveTransportRunner::run_until(scalar_t target, VariablesBasePtr variables)
{
Timer total_time;
total_time.start();
m_timestepper.set_total_target(target);
m_output_target = m_timestepper.get_total()+m_simulinfo.output_step;
std::ofstream out_iter;
if (m_simulinfo.print_iter_info)
{
out_iter.open(m_simulinfo.complete_filepath("iter", "dat"));
io::print_reactmicp_header(&out_iter);
}
scalar_t dt = m_timestepper.get_options().restart_timestep;
ReactiveTransportReturnCode retcode = ReactiveTransportReturnCode::NotConvergedYet;
while(retcode >= ReactiveTransportReturnCode::NotConvergedYet
and m_timestepper.get_total() < m_timestepper.get_total_target())
{
Timer step_timer;
step_timer.start();
reactmicp::solver::ReactiveTransportReturnCode retcode = m_solver.solve_timestep(dt, variables);
step_timer.stop();
if (m_simulinfo.print_iter_info)
io::print_reactmicp_performance_long(&out_iter, m_cnt, m_timestepper.get_total()+dt, m_solver.get_perfs());
dt = m_timestepper.next_timestep(dt, retcode, m_solver.get_perfs().nb_iterations);
if (retcode <= reactmicp::solver::ReactiveTransportReturnCode::NotConvergedYet)
{
dt = m_timestepper.get_options().restart_timestep;
variables->reset_main_variables();
retcode = m_solver.solve_timestep(dt, variables);
if (m_simulinfo.print_iter_info)
io::print_reactmicp_performance_long(&out_iter, m_cnt, m_timestepper.get_total()+dt, m_solver.get_perfs());
dt = m_timestepper.next_timestep(dt, retcode, m_solver.get_perfs().nb_iterations);
}
++m_cnt;
// output
if (m_timestepper.get_total() > m_output_target)
{
m_output_function(m_timestepper.get_total(), variables);
m_output_target += m_simulinfo.output_step;
}
}
total_time.stop();
if (m_simulinfo.print_iter_info)
{
io::print_reactmicp_end(&out_iter, total_time, get_timer());
out_iter.close();
}
}
} // end namespace solver
} // end namespace reactmicp
} // end namespace specmicp
c4science · Help