DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec > Class Template Reference

A class to perform time stepping using the adaptive Boost Runge-Kutta method for the explicit part and SUNDIALS IDA for the implicit part. In this scheme, the IDA stepper is the controller and the Boost RK stepper solves the explicit part of the problem on-demand. More...

#include <sundials_ida_boost_rk.hh>

Inheritance diagram for DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >:

Public Types
using	NumberType = typename Base::NumberType
using	BlockVectorType = typename Base::BlockVectorType
using	error_stepper_type = typename stepperChoice<prec>::value

Public Member Functions
	TimeStepperSUNDIALS_IDA_BoostRK_impl (const ConfigTree &config, AbstractAssembler< VectorType, SparseMatrixType, dim > &assembler, OutputSession_impl< dim, VectorType > &data_out)
	TimeStepperSUNDIALS_IDA_BoostRK_impl (const ConfigTree &config, AbstractAssembler< VectorType, SparseMatrixType, dim > &assembler, OutputSession_impl< dim, VectorType > &data_out, AbstractAdaptor< VectorType > &adaptor)
virtual void	run (AbstractFlowingVariables< NumberType, VectorType > &initial_condition, const double t_start, const double t_stop) override
	Any derived class must implement this method to run the timestepping algorithm.
Public Member Functions inherited from DiFfRG::AbstractTimestepper< VectorType, SparseMatrixType, dim >
	AbstractTimestepper (const ConfigTree &config, AbstractAssembler< VectorType, SparseMatrixType, dim > &assembler, OutputSession_impl< dim, VectorType > &data_out, const bool implicit_stepper, const bool explicit_stepper)
	Construct a new Abstract Timestepper object.
	AbstractTimestepper (const ConfigTree &config, AbstractAssembler< VectorType, SparseMatrixType, dim > &assembler, OutputSession_impl< dim, VectorType > &data_out, AbstractAdaptor< VectorType > &adaptor, const bool implicit_stepper, const bool explicit_stepper)

Private Types
using	Base = AbstractTimestepper<VectorType, SparseMatrixType, dim>

Private Member Functions
void	run (BlockVectorType &initial_data, const double t_start, const double t_stop)
void	update_variables (VectorType &variable_y, const VectorType &spatial_y, const double t)

Additional Inherited Members
Protected Types inherited from DiFfRG::AbstractTimestepper< VectorType, SparseMatrixType, dim >
using	VectorType
using	NumberType
using	SparseMatrixType
using	BlockVectorType
Protected Member Functions inherited from DiFfRG::AbstractTimestepper< VectorType, SparseMatrixType, dim >
void	drain_output ()
void	finalize_output_after_failure (EmitFinalFrame &&emit_final_frame)
	Best-effort output cleanup while a timestepping failure is being reported.
bool	is_implicit () const
bool	is_explicit () const
Protected Attributes inherited from DiFfRG::AbstractTimestepper< VectorType, SparseMatrixType, dim >
const ConfigTree	config
NoAdaptivity< VectorType >	adaptor_default
AbstractAssembler< VectorType, SparseMatrixType, dim > &	assembler
OutputSession_impl< dim, VectorType > &	data_out
AbstractAdaptor< VectorType > &	adaptor
ReportPort	log
const bool	m_is_implicit
const bool	m_is_explicit
double	output_dt
struct DiFfRG::AbstractTimestepper::ImplicitParameters	impl
struct DiFfRG::AbstractTimestepper::ExplicitParameters	expl
std::size_t	next_jacobian_build_id
TimestepperJacobianDiagnosticsState	jacobian_diagnostics_state
Static Protected Attributes inherited from DiFfRG::AbstractTimestepper< VectorType, SparseMatrixType, dim >
static constexpr uint	dim

Detailed Description

template<typename VectorType, typename SparseMatrixType, uint dim, template< typename... > typename LinearSolver, int prec>
class DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >

A class to perform time stepping using the adaptive Boost Runge-Kutta method for the explicit part and SUNDIALS IDA for the implicit part. In this scheme, the IDA stepper is the controller and the Boost RK stepper solves the explicit part of the problem on-demand.

Template Parameters

VectorType	Type of the vector
dim	Dimension of the problem
prec	Algorithm choice: 0 for Cash-Karp54 (5th order), 1 for Fehlberg78 (8th order)

Member Typedef Documentation

Base

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

using DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::Base = AbstractTimestepper<VectorType, SparseMatrixType, dim>

private

BlockVectorType

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

using DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::BlockVectorType = typename Base::BlockVectorType

error_stepper_type

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

using DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::error_stepper_type = typename stepperChoice<prec>::value

NumberType

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

using DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::NumberType = typename Base::NumberType

Constructor & Destructor Documentation

TimeStepperSUNDIALS_IDA_BoostRK_impl() [1/2]

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::TimeStepperSUNDIALS_IDA_BoostRK_impl ( const ConfigTree & config, AbstractAssembler< VectorType, SparseMatrixType, dim > & assembler, OutputSession_impl< dim, VectorType > & data_out )

inline

Forwards to the base with this stepper's kind. Not using Base::Base;: the base needs to know which /timestepping/ section to read, and it cannot ask a virtual function for that from inside its own constructor.

TimeStepperSUNDIALS_IDA_BoostRK_impl() [2/2]

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::TimeStepperSUNDIALS_IDA_BoostRK_impl ( const ConfigTree & config, AbstractAssembler< VectorType, SparseMatrixType, dim > & assembler, OutputSession_impl< dim, VectorType > & data_out, AbstractAdaptor< VectorType > & adaptor )

inline

Member Function Documentation

run() [1/2]

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

virtual void DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::run ( AbstractFlowingVariables< NumberType, VectorType > & initial_condition, const double t_start, const double t_stop )

overridevirtual

Any derived class must implement this method to run the timestepping algorithm.

Parameters

initial_condition	The flowing variables object that contains the initial condition.
t_start	The start time of the simulation.
t_stop	The run method will evolve the system from t_start to t_stop.

Implements DiFfRG::AbstractTimestepper< VectorType, SparseMatrixType, dim >.

run() [2/2]

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

void DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::run ( BlockVectorType & initial_data, const double t_start, const double t_stop )

private

update_variables()

template<typename VectorType , typename SparseMatrixType , uint dim, template< typename... > typename LinearSolver, int prec>

void DiFfRG::TimeStepperSUNDIALS_IDA_BoostRK_impl< VectorType, SparseMatrixType, dim, LinearSolver, prec >::update_variables ( VectorType & variable_y, const VectorType & spatial_y, const double t )

private

---

The documentation for this class was generated from the following file:

• /home/runner/work/DiFfRG_current/DiFfRG_current/DiFfRG/include/DiFfRG/timestepping/sundials_ida_boost_rk.hh