#include <source/euler_poisson/parabolic_module.h>

Inheritance diagram for ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >:

Collaboration diagram for ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >:

Public Member Functions
Constructor and setup
	ParabolicModule (const MPIEnsemble &mpi_ensemle, const OfflineData< dim, Number > &offline_data, const HyperbolicSystem &hyperbolic_system, const ParabolicSystem &parabolic_system, const InitialValues< HyperbolicDescription, dim, Number > &initial_values, const std::string &subsection="/ParabolicModule")

void	prepare ()

Functions for performing explicit time steps
void	reinit_state_vector (StateVector< Number > &state_vector) const

void	prepare_state_vector (StateVector< Number > &state_vector, Number t) const

template<int stages>
void	backward_euler_step (const StateVector< Number > &old_state_vector, const Number old_t, std::array< std::reference_wrapper< const StateVector< Number > >, stages > stage_state_vectors, const std::array< Number, stages > stage_weights, StateVector< Number > &new_state_vector, Number tau) const

void	crank_nicolson_step (const StateVector< Number > &old_state_vector, const Number old_t, StateVector< Number > &new_state_vector, Number tau) const

void	set_id_violation_strategy (const IDViolationStrategy &strategy) const

Information and statistics
void	print_solver_statistics (std::ostream &output) const

Typedefs and constexpr constants
using	HyperbolicSystem = typename HyperbolicDescription::HyperbolicSystem

using	View = typename HyperbolicSystem::template View< dim, Number >

using	ParabolicSystem = EulerPoisson::ParabolicSystem

using	BlockHostVector = dealii::LinearAlgebra::distributed::BlockVector< Number >

using	ScalarNumber = typename View::ScalarNumber

using	state_type = typename View::state_type

static constexpr auto	problem_dimension = View::problem_dimension

Detailed Description

template<typename HyperbolicDescription, int dim, typename Number>
class ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >

Implicit backward-Euler time stepping for the parabolic limiting equation for the Euler-Poisson system

Definition at line 133 of file parabolic_module.h.

Member Typedef Documentation

◆ HyperbolicSystem

template<typename HyperbolicDescription , int dim, typename Number >

using ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::HyperbolicSystem = typename HyperbolicDescription::HyperbolicSystem

Definition at line 141 of file parabolic_module.h.

◆ View

template<typename HyperbolicDescription , int dim, typename Number >

using ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::View = typename HyperbolicSystem::template View<dim, Number>

Definition at line 143 of file parabolic_module.h.

◆ ParabolicSystem

template<typename HyperbolicDescription , int dim, typename Number >

using ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::ParabolicSystem = EulerPoisson::ParabolicSystem

Definition at line 145 of file parabolic_module.h.

◆ BlockHostVector

template<typename HyperbolicDescription , int dim, typename Number >

using ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::BlockHostVector = dealii::LinearAlgebra::distributed::BlockVector<Number>

Definition at line 147 of file parabolic_module.h.

◆ ScalarNumber

template<typename HyperbolicDescription , int dim, typename Number >

using ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::ScalarNumber = typename View::ScalarNumber

Definition at line 150 of file parabolic_module.h.

◆ state_type

template<typename HyperbolicDescription , int dim, typename Number >

using ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::state_type = typename View::state_type

Definition at line 154 of file parabolic_module.h.

Constructor & Destructor Documentation

◆ ParabolicModule()

template<typename HyperbolicDescription , int dim, typename Number >

ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::ParabolicModule	(	const MPIEnsemble &	mpi_ensemle,
		const OfflineData< dim, Number > &	offline_data,
		const HyperbolicSystem &	hyperbolic_system,
		const ParabolicSystem &	parabolic_system,
		const InitialValues< HyperbolicDescription, dim, Number > &	initial_values,
		const std::string &	subsection = `"/ParabolicModule< HyperbolicDescription, dim, Number >"`
	)

Constructor.

Definition at line 32 of file parabolic_module.template.h.

References ryujin::EulerPoisson::no_restart.

Member Function Documentation

◆ prepare()

template<typename HyperbolicDescription , int dim, typename Number >

void ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::prepare ( )

Prepare time stepping. A call to prepare() allocates temporary storage and is necessary before any of the following time-stepping functions can be called.

Definition at line 135 of file parabolic_module.template.h.

References ryujin::cg_q1, and ryujin::dg_q1.

◆ reinit_state_vector()

template<typename HyperbolicDescription , int dim, typename Number >

void ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::reinit_state_vector ( StateVector< Number > & state_vector ) const

(Re)initialize the parabolic state vector component of the state vector.

Note: This routine does not modify the hyperbolic state vector or the precomputed vector component.

Definition at line 243 of file parabolic_module.template.h.

References ryujin::precomputed.

◆ prepare_state_vector()

template<typename HyperbolicDescription , int dim, typename Number >

void ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::prepare_state_vector	(	StateVector< Number > &	state_vector,
		Number	t
	)		const

This function preprocesses a given state vector U in preparation for a high order IMEX time step. This function exists because some time stepping variants have to precompute quantities before we can perform an IMEX step. In addition, this function is called whenever we perform a mesh transfer or output operation.

Definition at line 262 of file parabolic_module.template.h.

References ryujin::EulerPoisson::correction, ryujin::EulerPoisson::full_restart, and ryujin::EulerPoisson::static_full_restart.

◆ backward_euler_step()

template<typename HyperbolicDescription , int dim, typename Number >

template<int stages>

void ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::backward_euler_step	(	const StateVector< Number > &	old_state_vector,
		const Number	old_t,
		std::array< std::reference_wrapper< const StateVector< Number > >, stages >	stage_state_vectors,
		const std::array< Number, stages >	stage_weights,
		StateVector< Number > &	new_state_vector,
		Number	tau
	)		const

Given a reference to a previous state vector old_U at time old_t and a time-step size tau perform an implicit backward euler step (and store the result in new_U).

The function takes an optional array of states stage_U together with a an array of weights stage_weights to construct a modified high-order right-hand side / flux.

Definition at line 298 of file parabolic_module.template.h.

◆ crank_nicolson_step()

template<typename HyperbolicDescription , int dim, typename Number >

void ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::crank_nicolson_step	(	const StateVector< Number > &	old_state_vector,
		const Number	old_t,
		StateVector< Number > &	new_state_vector,
		Number	tau
	)		const

Given a reference to a previous state vector old_U at time old_t and a time-step size tau perform an implicit Crank-Nicolson step (and store the result in new_U).

This variant is used in the TimeIntegrator class for the Strang split variants.

Definition at line 317 of file parabolic_module.template.h.

◆ set_id_violation_strategy()

template<typename HyperbolicDescription , int dim, typename Number >

void ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::set_id_violation_strategy ( const IDViolationStrategy & strategy ) const

inline

Sets the invariant domain violation strategy.

Definition at line 240 of file parabolic_module.h.

◆ print_solver_statistics()

template<typename HyperbolicDescription , int dim, typename Number >

void ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::print_solver_statistics ( std::ostream & output ) const

Print a status line with solver statistics. This function is used for constructing the status message displayed periodically in the TimeLoop.

Definition at line 352 of file parabolic_module.template.h.

Member Data Documentation

◆ problem_dimension

template<typename HyperbolicDescription , int dim, typename Number >

constexpr auto ryujin::EulerPoisson::ParabolicModule< HyperbolicDescription, dim, Number >::problem_dimension = View::problem_dimension

staticconstexpr

Definition at line 152 of file parabolic_module.h.

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

source/euler_poisson/parabolic_module.h
source/euler_poisson/parabolic_module.template.h

Public Member Functions

Typedefs and constexpr constants

Detailed Description

Member Typedef Documentation

◆ HyperbolicSystem

◆ View

◆ ParabolicSystem

◆ BlockHostVector

◆ ScalarNumber

◆ state_type

Constructor & Destructor Documentation

◆ ParabolicModule()

Member Function Documentation

◆ prepare()

◆ reinit_state_vector()

◆ prepare_state_vector()

◆ backward_euler_step()

◆ crank_nicolson_step()

◆ set_id_violation_strategy()

◆ print_solver_statistics()

Member Data Documentation

◆ problem_dimension