SpECTRE
v2025.03.17
|
#include <ImexRungeKutta.hpp>
Classes | |
struct | ImplicitButcherTableau |
Implicit part of the Butcher tableau. Most parts of the tableau must be the same as the explicit part, and so are omitted. More... | |
Public Member Functions | |
virtual size_t | implicit_stage_order () const =0 |
virtual const ImplicitButcherTableau & | implicit_butcher_tableau () const =0 |
![]() | |
uint64_t | number_of_substeps () const override |
Number of substeps in this TimeStepper. More... | |
uint64_t | number_of_substeps_for_error () const override |
Number of substeps in this TimeStepper when providing an error measure for adaptive time-stepping. More... | |
size_t | number_of_past_steps () const override |
Number of past time steps needed for multi-step method. More... | |
bool | monotonic () const override |
Whether computational and temporal orderings of operations match. More... | |
TimeStepId | next_time_id (const TimeStepId ¤t_id, const TimeDelta &time_step) const override |
The TimeStepId after the current substep. More... | |
TimeStepId | next_time_id_for_error (const TimeStepId ¤t_id, const TimeDelta &time_step) const override |
The TimeStepId after the current substep when providing an error measure for adaptive time-stepping. More... | |
virtual const ButcherTableau & | butcher_tableau () const =0 |
![]() | |
WRAPPED_PUPable_abstract (TimeStepper) | |
template<typename Vars > | |
void | update_u (const gsl::not_null< Vars * > u, const TimeSteppers::History< Vars > &history, const TimeDelta &time_step) const |
Set u to the value at the end of the current substep. More... | |
template<typename Vars > | |
std::optional< StepperErrorEstimate > | update_u (const gsl::not_null< Vars * > u, const TimeSteppers::History< Vars > &history, const TimeDelta &time_step, const std::optional< StepperErrorTolerances > &tolerances) const |
Set u to the value at the end of the current substep; report the error measure when available. More... | |
template<typename Vars > | |
void | clean_history (const gsl::not_null< TimeSteppers::History< Vars > * > history) const |
Remove old entries from the history. More... | |
template<typename Vars > | |
bool | dense_update_u (const gsl::not_null< Vars * > u, const TimeSteppers::History< Vars > &history, const double time) const |
Compute the solution value at a time between steps. To evaluate at a time within a given step, call this method at the start of the step containing the time. The function returns true on success, otherwise the call should be retried after the next substep. More... | |
virtual size_t | order () const =0 |
The convergence order of the stepper. More... | |
virtual uint64_t | number_of_substeps () const =0 |
Number of substeps in this TimeStepper. More... | |
virtual uint64_t | number_of_substeps_for_error () const =0 |
Number of substeps in this TimeStepper when providing an error measure for adaptive time-stepping. More... | |
virtual size_t | number_of_past_steps () const =0 |
Number of past time steps needed for multi-step method. More... | |
virtual double | stable_step () const =0 |
Rough estimate of the maximum step size this method can take stably as a multiple of the step for Euler's method. More... | |
virtual bool | monotonic () const =0 |
Whether computational and temporal orderings of operations match. More... | |
virtual TimeStepId | next_time_id (const TimeStepId ¤t_id, const TimeDelta &time_step) const =0 |
The TimeStepId after the current substep. More... | |
virtual TimeStepId | next_time_id_for_error (const TimeStepId ¤t_id, const TimeDelta &time_step) const =0 |
The TimeStepId after the current substep when providing an error measure for adaptive time-stepping. More... | |
template<typename Vars > | |
bool | can_change_step_size (const TimeStepId &time_id, const TimeSteppers::History< Vars > &history) const |
Whether a change in the step size is allowed before taking a step. Step sizes can never be changed on a substep. More... | |
![]() | |
WRAPPED_PUPable_abstract (ImexTimeStepper) | |
virtual size_t | imex_order () const =0 |
Convergence order of the integrator when used in IMEX mode. More... | |
template<typename Vars > | |
void | add_inhomogeneous_implicit_terms (const gsl::not_null< Vars * > u, const TimeSteppers::History< Vars > &implicit_history, const TimeDelta &time_step) const |
Add the change for the current implicit substep, | |
template<typename Vars > | |
double | implicit_weight (const TimeSteppers::History< Vars > &implicit_history, const TimeDelta &time_step) const |
The coefficient | |
Additional Inherited Members | |
![]() | |
using | provided_time_stepper_interfaces = implementation defined |
![]() | |
using | provided_time_stepper_interfaces = implementation defined |
![]() | |
static constexpr bool | local_time_stepping = false |
static constexpr bool | imex = false |
![]() | |
static constexpr bool | imex = true |
Intermediate base class implementing a generic IMEX Runge-Kutta scheme.
Implements most of the virtual methods of ImexTimeStepper for a generic Runge-Kutta method. Derived classes must implement the requirements of the RungeKutta
base class, as well as imex_order()
, implicit_stage_order()
, and implicit_butcher_tableau()
.
|
pure virtual |
Smallest order of the intermediate result at any substep. For the methods supported by this class, this cannot exceed 2.
Implemented in TimeSteppers::Heun2, TimeSteppers::Rk3Kennedy, TimeSteppers::Rk3Pareschi, and TimeSteppers::Rk4Kennedy.