22#ifndef INCLUDE_EXADG_TIME_INTEGRATION_TIME_STEP_CALCULATION_H_
23#define INCLUDE_EXADG_TIME_INTEGRATION_TIME_STEP_CALCULATION_H_
31limit_time_step_change(
double & new_time_step,
double const & last_time_step,
double const & fac)
33 if(new_time_step >= fac * last_time_step)
35 new_time_step = fac * last_time_step;
37 else if(new_time_step <= last_time_step / fac)
39 new_time_step = last_time_step / fac;
47adjust_time_step_to_hit_end_time(
double const start_time,
48 double const end_time,
49 double const time_step)
51 return (end_time - start_time) / (1 + int((end_time - start_time) / time_step));
59calculate_const_time_step(
double const dt,
unsigned int const n_refine_time)
61 double const time_step = dt / std::pow(2., n_refine_time);