The LDG assembler that can be used for any LDG scheme, with as many levels as one wants. More...

#include <ldg.hh>

Inheritance diagram for DiFfRG::LDG::Assembler< Discretization_, Model_ >:

Public Types
using	Discretization = Discretization_

using	Model = Model_

using	NumberType = typename Discretization::NumberType

using	VectorType = typename Discretization::VectorType

using	Components = typename Discretization::Components

Public Types inherited from DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >
using	Discretization = Discretization_

using	Model = Model_

using	NumberType = typename Discretization::NumberType

using	VectorType = typename Discretization::VectorType

using	Components = typename Discretization::Components

Public Types inherited from DiFfRG::AbstractAssembler< Discretization_::VectorType, Discretization_::SparseMatrixType, Discretization_::dim >
using	NumberType

Public Member Functions
	Assembler (Discretization &discretization, Model &model, const JSONValue &json)

virtual void	reinit_vector (VectorType &vec) const override

virtual void	attach_data_output (DataOutput< dim, VectorType > &data_out, const VectorType &solution, const VectorType &variables, const VectorType &dt_solution=VectorType(), const VectorType &residual=VectorType()) override
	Attach all intermediate (ldg) vectors to the data output.

virtual void	reinit () override
	Reinitialize the assembler. This is necessary if the mesh has changed, e.g. after a mesh refinement.

virtual void	rebuild_jacobian_sparsity () override

virtual void	refinement_indicator (Vector< double > &indicator, const VectorType &solution_global) override

virtual const BlockSparsityPattern &	get_sparsity_pattern_jacobian () const override
	Obtain the sparsity pattern of the jacobian matrix.

virtual const BlockSparseMatrix< NumberType > &	get_mass_matrix () const override
	Obtain the mass matrix.

virtual void	mass (VectorType &residual, const VectorType &solution_global, const VectorType &solution_global_dot, NumberType weight) override
	Construct the mass.

virtual void	residual (VectorType &residual, const VectorType &solution_global, NumberType weight, const VectorType &solution_global_dot, NumberType weight_mass, const VectorType &variables=VectorType()) override
	Construct the system residual, i.e. Res = grad(flux) - source.

virtual void	jacobian_mass (BlockSparseMatrix< NumberType > &jacobian, const VectorType &solution_global, const VectorType &solution_global_dot, NumberType alpha=1., NumberType beta=1.) override

virtual void	jacobian (BlockSparseMatrix< NumberType > &jacobian, const VectorType &solution_global, NumberType weight, const VectorType &solution_global_dot, NumberType alpha, NumberType beta, const VectorType &variables=VectorType()) override
	Construct the system jacobian, i.e. dRes/du.

void	log (const std::string logger)

double	average_time_reinit () const

uint	num_reinits () const

double	average_time_residual_assembly () const

uint	num_residuals () const

double	average_time_jacobian_assembly () const

uint	num_jacobians () const

Public Member Functions inherited from DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >
	LDGAssemblerBase (Discretization &discretization, Model &model, const JSONValue &json)

virtual IndexSet	get_differential_indices () const override
	Obtain the dofs which contain time derivatives.

const auto &	get_discretization () const

auto &	get_discretization ()

virtual void	set_time (double t) override
	Set the current time. The assembler should usually just forward this to the numerical model.

double	average_time_variable_residual_assembly ()

uint	num_variable_residuals () const

double	average_time_variable_jacobian_assembly ()

uint	num_variable_jacobians () const

Public Member Functions inherited from DiFfRG::AbstractAssembler< Discretization_::VectorType, Discretization_::SparseMatrixType, Discretization_::dim >
virtual void	attach_data_output (DataOutput< dim, VectorType > &data_out, const VectorType &solution, const VectorType &variables=VectorType(), const VectorType &dt_solution=VectorType(), const VectorType &residual=VectorType())=0
	Attach any data output to the DataOutput object provided. This can be used to extract additional data from the solution and write it to the output file. This includes both derivatives and other spatial functions, as well as single values that can be appended to the .csv file.

virtual void	reinit_vector (VectorType &vector) const=0
	Reinitialize an arbitrary vector so that it has the correct size and structure.

virtual void	residual_variables (VectorType &residual, const VectorType &variables, const VectorType &spatial_solution)
	When coupling the spatial discretization to additional variables, this function should calculate the residual for the additional variables.

virtual void	jacobian_variables (FullMatrix< NumberType > &jacobian, const VectorType &variables, const VectorType &spatial_solution)
	When coupling the spatial discretization to additional variables, this function should calculate the jacobian for the additional variables.

void	mass (VectorType &mass, const VectorType &solution_global, NumberType weight)
	Calculates the mass \(m_i(u)\) for an ODE.

virtual void	mass (VectorType &mass, const VectorType &solution_global, const VectorType &solution_global_dot, NumberType weight)=0
	Calculates the mass \(m(u, \partial_t u)\) for a DAE.

void	residual (VectorType &residual, const VectorType &solution_global, NumberType weight, NumberType weight_mass, const VectorType &variables=VectorType())
	Calculates the residual for an ODE.

virtual void	residual (VectorType &residual, const VectorType &solution_global, NumberType weight, const VectorType &solution_global_dot, NumberType weight_mass, const VectorType &variables=VectorType())=0
	Calculates the residual for a DAE.

void	jacobian_mass (SparseMatrixType &jacobian, const VectorType &solution_global, NumberType mass_weight=1.)
	Calculates the jacobian of the mass function for an ODE.

virtual void	jacobian_mass (SparseMatrixType &jacobian, const VectorType &solution_global, const VectorType &solution_global_dot, NumberType alpha=1., NumberType beta=1.)=0
	Calculates the jacobian of the mass function for a DAE.

void	jacobian (SparseMatrixType &jacobian, const VectorType &solution_global, NumberType weight, NumberType mass_weight, const VectorType &variables=VectorType())
	Calculates the jacobian of the residual function for an ODE.

virtual void	jacobian (SparseMatrixType &jacobian, const VectorType &solution_global, NumberType weight, const VectorType &solution_global_dot, NumberType alpha, NumberType beta, const VectorType &variables=VectorType())=0
	Calculates the jacobian of the residual function for a DAE.

Static Public Attributes
static constexpr uint	dim = Discretization::dim

static constexpr uint	stencil = Components::count_fe_subsystems()

Static Public Attributes inherited from DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >
static constexpr uint	dim = Discretization::dim

Protected Member Functions
void	rebuild_ldg_vectors (const VectorType &sol) const

template<int k>
void	rebuild_ldg_jacobian (const VectorType &sol) const

template<int from, int to, typename VectorType , typename VectorTypeldg >
void	build_ldg_vector (const VectorType &solution_global, VectorTypeldg &ldg_vector, VectorTypeldg &ldg_vector_tmp) const
	Build the LDG vector at level 'to', which takes information from level 'from'.

template<int from, int to, typename VectorType >
void	build_ldg_jacobian (const VectorType &solution_global, BlockSparseMatrix< NumberType > &ldg_jacobian, BlockSparseMatrix< NumberType > &ldg_jacobian_tmp) const
	Build the LDG jacobian at level 'to', which takes information from level 'from'.

void	build_ldg_sparsity (BlockSparsityPattern &sparsity_pattern, const DoFHandler< dim > &to_dofh, const DoFHandler< dim > &from_dofh, const int stencil=1, bool add_extractor_dofs=false) const
	Create a sparsity pattern for matrices between the DoFs of two DoFHandlers, with given stencil size.

void	build_inverse (const SparseMatrix< NumberType > &in, SparseMatrix< NumberType > &out) const
	Build the inverse of matrix in and save the result to out.

void	readouts (DataOutput< dim, VectorType > &data_out, const VectorType &solution_global, const VectorType &variables) const

void	extract (std::array< NumberType, Components::count_extractors()> &data, const VectorType &solution_global, const VectorType &variables, bool search_EoM, bool set_EoM, bool postprocess) const

bool	jacobian_extractors (FullMatrix< NumberType > &extractor_jacobian, const VectorType &solution_global, const VectorType &variables)

virtual void	residual_variables (VectorType &residual, const VectorType &variables, const VectorType &spatial_solution) override

virtual void	jacobian_variables (FullMatrix< NumberType > &jacobian, const VectorType &variables, const VectorType &spatial_solution) override

Static Protected Member Functions
template<typename... T>
static constexpr auto	v_tie (T &&...t)

template<typename... T>
static constexpr auto	e_tie (T &&...t)

Protected Attributes
QGauss< dim >	quadrature

QGauss< dim - 1 >	quadrature_face

std::vector< const DoFHandler< dim > * >	dof_handler_list

array< BlockVector< NumberType >, Components::count_fe_subsystems()>	sol_vector

array< BlockVector< NumberType >, Components::count_fe_subsystems()>	sol_vector_tmp

array< Vector< NumberType >, Components::count_fe_subsystems()>	sol_vector_vec_tmp

BlockSparsityPattern	sparsity_pattern_jacobian

BlockSparsityPattern	sparsity_pattern_mass

array< BlockSparsityPattern, Components::count_fe_subsystems()>	sparsity_pattern_ug

array< BlockSparsityPattern, Components::count_fe_subsystems()>	sparsity_pattern_gu

array< BlockSparsityPattern, Components::count_fe_subsystems()>	sparsity_pattern_wg

array< BlockSparseMatrix< NumberType >, Components::count_fe_subsystems()>	jacobian_tmp

BlockSparseMatrix< NumberType >	mass_matrix

SparseMatrix< NumberType >	component_mass_matrix_inverse

array< BlockSparseMatrix< NumberType >, Components::count_fe_subsystems()>	j_ug

array< BlockSparseMatrix< NumberType >, Components::count_fe_subsystems()>	j_gu

array< BlockSparseMatrix< NumberType >, Components::count_fe_subsystems()>	j_wg

array< BlockSparseMatrix< NumberType >, Components::count_fe_subsystems()>	j_wg_tmp

std::vector< double >	timings_reinit

std::vector< double >	timings_residual

std::vector< double >	timings_jacobian

array< bool, Components::count_fe_subsystems()>	ldg_matrix_built

array< bool, Components::count_fe_subsystems()>	jacobian_tmp_built

Discretization &	discretization

const DoFHandler< dim > &	dof_handler

const FiniteElement< dim > &	fe

const Mapping< dim > &	mapping

Model &	model

uint	batch_size

uint	threads

const double	EoM_abs_tol

const uint	EoM_max_iter

std::vector< types::global_dof_index >	extractor_dof_indices

Point< dim >	EoM

DoFHandler< dim >::cell_iterator	EoM_cell

FullMatrix< NumberType >	extractor_jacobian_u

DoFHandler< dim >::cell_iterator	old_EoM_cell

std::vector< double >	timings_variable_jacobian

std::vector< double >	timings_variable_residual

Protected Attributes inherited from DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >
Discretization &	discretization

Model &	model

const FiniteElement< dim > &	fe

const DoFHandler< dim > &	dof_handler

const Mapping< dim > &	mapping

uint	threads

uint	batch_size

DoFHandler< dim >::cell_iterator	EoM_cell

DoFHandler< dim >::cell_iterator	old_EoM_cell

const double	EoM_abs_tol

const uint	EoM_max_iter

Point< dim >	EoM

FullMatrix< NumberType >	extractor_jacobian

FullMatrix< NumberType >	extractor_jacobian_u

FullMatrix< NumberType >	extractor_jacobian_du

FullMatrix< NumberType >	extractor_jacobian_ddu

std::vector< types::global_dof_index >	extractor_dof_indices

std::vector< double >	timings_variable_residual

std::vector< double >	timings_variable_jacobian

Static Protected Attributes
static constexpr int	nothing = 0

Private Types
using	Base = LDGAssemblerBase<Discretization_, Model_>

Private Member Functions
template<typename... T>
auto	fe_conv (std::tuple< T &... > &t) const

template<typename... T>
auto	fe_more_conv (std::tuple< T &... > &t) const

template<typename... T>
auto	ref_conv (std::tuple< T &... > &t) const

Detailed Description

template<typename Discretization_, typename Model_>
class DiFfRG::LDG::Assembler< Discretization_, Model_ >

The LDG assembler that can be used for any LDG scheme, with as many levels as one wants.

Template Parameters

Discretization	Discretization on which to assemble
Model	The model class which contains the physical equations.

Member Typedef Documentation

◆ Base

template<typename Discretization_ , typename Model_ >

using DiFfRG::LDG::Assembler< Discretization_, Model_ >::Base = LDGAssemblerBase<Discretization_, Model_>

private

◆ Components

template<typename Discretization_ , typename Model_ >

using DiFfRG::LDG::Assembler< Discretization_, Model_ >::Components = typename Discretization::Components

◆ Discretization

template<typename Discretization_ , typename Model_ >

using DiFfRG::LDG::Assembler< Discretization_, Model_ >::Discretization = Discretization_

◆ Model

template<typename Discretization_ , typename Model_ >

using DiFfRG::LDG::Assembler< Discretization_, Model_ >::Model = Model_

◆ NumberType

template<typename Discretization_ , typename Model_ >

using DiFfRG::LDG::Assembler< Discretization_, Model_ >::NumberType = typename Discretization::NumberType

◆ VectorType

template<typename Discretization_ , typename Model_ >

using DiFfRG::LDG::Assembler< Discretization_, Model_ >::VectorType = typename Discretization::VectorType

Constructor & Destructor Documentation

◆ Assembler()

template<typename Discretization_ , typename Model_ >

DiFfRG::LDG::Assembler< Discretization_, Model_ >::Assembler	(	Discretization &	discretization,
		Model &	model,
		const JSONValue &	json )

inline

Member Function Documentation

◆ attach_data_output()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::attach_data_output	(	DataOutput< dim, VectorType > &	data_out,
		const VectorType &	solution,
		const VectorType &	variables,
		const VectorType &	dt_solution = VectorType(),
		const VectorType &	residual = VectorType() )

inlineoverridevirtual

Attach all intermediate (ldg) vectors to the data output.

Parameters

data_out	A DataOutput to which we attach data
sol	The current global solution

◆ average_time_jacobian_assembly()

template<typename Discretization_ , typename Model_ >

double DiFfRG::LDG::Assembler< Discretization_, Model_ >::average_time_jacobian_assembly ( ) const

inline

◆ average_time_reinit()

template<typename Discretization_ , typename Model_ >

double DiFfRG::LDG::Assembler< Discretization_, Model_ >::average_time_reinit ( ) const

inline

◆ average_time_residual_assembly()

template<typename Discretization_ , typename Model_ >

double DiFfRG::LDG::Assembler< Discretization_, Model_ >::average_time_residual_assembly ( ) const

inline

◆ build_inverse()

template<typename Discretization_ , typename Model_ >

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::build_inverse	(	const SparseMatrix< NumberType > &	in,
		SparseMatrix< NumberType > &	out ) const

inlineprotected

Build the inverse of matrix in and save the result to out.

Parameters

in	The matrix to invert. Must have a valid sparsity pattern.
out	The matrix to store the result. Must have a valid sparsity pattern.

◆ build_ldg_jacobian()

template<typename Discretization_ , typename Model_ >

template<int from, int to, typename VectorType >

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::build_ldg_jacobian	(	const VectorType &	solution_global,
		BlockSparseMatrix< NumberType > &	ldg_jacobian,
		BlockSparseMatrix< NumberType > &	ldg_jacobian_tmp ) const

inlineprotected

Build the LDG jacobian at level 'to', which takes information from level 'from'.

Template Parameters

from	Level to take information from.
to	Level to write information to.

Parameters

solution_global	The solution at level 'from'
ldg_jacobian	Where to store the result.
ldg_jacobian_tmp	A temporary of the same size as ldg_jacobian.

◆ build_ldg_sparsity()

template<typename Discretization_ , typename Model_ >

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::build_ldg_sparsity	(	BlockSparsityPattern &	sparsity_pattern,
		const DoFHandler< dim > &	to_dofh,
		const DoFHandler< dim > &	from_dofh,
		const int	stencil = 1,
		bool	add_extractor_dofs = false ) const

inlineprotected

Create a sparsity pattern for matrices between the DoFs of two DoFHandlers, with given stencil size.

Parameters

sparsity_pattern	The pattern to store the result in.
dofh_to	DoFHandler giving the row DoFs.
dofh_from	DoFHandler giving the column DoFs.
stencil	Stencil size of the resulting pattern.

◆ build_ldg_vector()

template<typename Discretization_ , typename Model_ >

template<int from, int to, typename VectorType , typename VectorTypeldg >

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::build_ldg_vector	(	const VectorType &	solution_global,
		VectorTypeldg &	ldg_vector,
		VectorTypeldg &	ldg_vector_tmp ) const

inlineprotected

Build the LDG vector at level 'to', which takes information from level 'from'.

Template Parameters

from	Level to take information from.
to	Level to write information to.

Parameters

solution_global	The solution at level 'from'
ldg_vector	Where to store the result.
ldg_vector_tmp	A temporary of the same size as ldg_vector.

◆ e_tie()

template<typename Discretization_ , typename Model_ >

template<typename... T>

static constexpr auto DiFfRG::LDG::Assembler< Discretization_, Model_ >::e_tie ( T &&... t )

inlinestaticconstexprprotected

◆ extract()

template<typename Discretization_ , typename Model_ >

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::extract	(	std::array< NumberType, Components::count_extractors()> &	data,
		const VectorType &	solution_global,
		const VectorType &	variables,
		bool	search_EoM,
		bool	set_EoM,
		bool	postprocess ) const

inlineprotected

◆ fe_conv()

template<typename Discretization_ , typename Model_ >

template<typename... T>

auto DiFfRG::LDG::Assembler< Discretization_, Model_ >::fe_conv ( std::tuple< T &... > & t ) const

inlineprivate

◆ fe_more_conv()

template<typename Discretization_ , typename Model_ >

template<typename... T>

auto DiFfRG::LDG::Assembler< Discretization_, Model_ >::fe_more_conv ( std::tuple< T &... > & t ) const

inlineprivate

◆ get_mass_matrix()

template<typename Discretization_ , typename Model_ >

virtual const BlockSparseMatrix< NumberType > & DiFfRG::LDG::Assembler< Discretization_, Model_ >::get_mass_matrix ( ) const

inlineoverridevirtual

Obtain the mass matrix.

Returns: const SparseMatrixType& The mass matrix

Implements DiFfRG::AbstractAssembler< Discretization_::VectorType, Discretization_::SparseMatrixType, Discretization_::dim >.

◆ get_sparsity_pattern_jacobian()

template<typename Discretization_ , typename Model_ >

virtual const BlockSparsityPattern & DiFfRG::LDG::Assembler< Discretization_, Model_ >::get_sparsity_pattern_jacobian ( ) const

inlineoverridevirtual

Obtain the sparsity pattern of the jacobian matrix.

Returns: const SparsityPattern<VectorType>& The sparsity pattern of the jacobian matrix

Implements DiFfRG::AbstractAssembler< Discretization_::VectorType, Discretization_::SparseMatrixType, Discretization_::dim >.

◆ jacobian()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::jacobian	(	BlockSparseMatrix< NumberType > &	jacobian,
		const VectorType &	solution_global,
		NumberType	weight,
		const VectorType &	solution_global_dot,
		NumberType	alpha,
		NumberType	beta,
		const VectorType &	variables = VectorType() )

inlineoverridevirtual

Construct the system jacobian, i.e. dRes/du.

Parameters

jacobian	The result is stored here.
solution_global	The current global solution.
weight	A factor to multiply the whole jacobian with.

◆ jacobian_extractors()

template<typename Discretization_ , typename Model_ >

bool DiFfRG::LDG::Assembler< Discretization_, Model_ >::jacobian_extractors	(	FullMatrix< NumberType > &	extractor_jacobian,
		const VectorType &	solution_global,
		const VectorType &	variables )

inlineprotected

◆ jacobian_mass()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::jacobian_mass	(	BlockSparseMatrix< NumberType > &	jacobian,
		const VectorType &	solution_global,
		const VectorType &	solution_global_dot,
		NumberType	alpha = 1.,
		NumberType	beta = 1. )

inlineoverridevirtual

◆ jacobian_variables()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::jacobian_variables	(	FullMatrix< NumberType > &	jacobian,
		const VectorType &	variables,
		const VectorType &	spatial_solution )

inlineoverrideprotectedvirtual

◆ log()

template<typename Discretization_ , typename Model_ >

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::log ( const std::string logger )

inline

◆ mass()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::mass	(	VectorType &	residual,
		const VectorType &	solution_global,
		const VectorType &	solution_global_dot,
		NumberType	weight )

inlineoverridevirtual

Construct the mass.

Parameters

residual	The result is stored here.
solution_global	The current global solution.
weight	A factor to multiply the whole residual with.

◆ num_jacobians()

template<typename Discretization_ , typename Model_ >

uint DiFfRG::LDG::Assembler< Discretization_, Model_ >::num_jacobians ( ) const

inline

◆ num_reinits()

template<typename Discretization_ , typename Model_ >

uint DiFfRG::LDG::Assembler< Discretization_, Model_ >::num_reinits ( ) const

inline

◆ num_residuals()

template<typename Discretization_ , typename Model_ >

uint DiFfRG::LDG::Assembler< Discretization_, Model_ >::num_residuals ( ) const

inline

◆ readouts()

template<typename Discretization_ , typename Model_ >

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::readouts	(	DataOutput< dim, VectorType > &	data_out,
		const VectorType &	solution_global,
		const VectorType &	variables ) const

inlineprotected

◆ rebuild_jacobian_sparsity()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::rebuild_jacobian_sparsity ( )

inlineoverridevirtual

Implements DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >.

◆ rebuild_ldg_jacobian()

template<typename Discretization_ , typename Model_ >

template<int k>

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::rebuild_ldg_jacobian ( const VectorType & sol ) const

inlineprotected

◆ rebuild_ldg_vectors()

template<typename Discretization_ , typename Model_ >

void DiFfRG::LDG::Assembler< Discretization_, Model_ >::rebuild_ldg_vectors ( const VectorType & sol ) const

inlineprotected

◆ ref_conv()

template<typename Discretization_ , typename Model_ >

template<typename... T>

auto DiFfRG::LDG::Assembler< Discretization_, Model_ >::ref_conv ( std::tuple< T &... > & t ) const

inlineprivate

◆ refinement_indicator()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::refinement_indicator	(	Vector< double > &	indicator,
		const VectorType &	solution_global )

inlineoverridevirtual

Implements DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >.

◆ reinit()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::reinit ( )

inlineoverridevirtual

Reinitialize the assembler. This is necessary if the mesh has changed, e.g. after a mesh refinement.

Reimplemented from DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >.

◆ reinit_vector()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::reinit_vector ( VectorType & vec ) const

inlineoverridevirtual

◆ residual()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::residual	(	VectorType &	residual,
		const VectorType &	solution_global,
		NumberType	weight,
		const VectorType &	solution_global_dot,
		NumberType	weight_mass,
		const VectorType &	variables = VectorType() )

inlineoverridevirtual

Construct the system residual, i.e. Res = grad(flux) - source.

Parameters

residual	The result is stored here.
solution_global	The current global solution.
weight	A factor to multiply the whole residual with.

◆ residual_variables()

template<typename Discretization_ , typename Model_ >

virtual void DiFfRG::LDG::Assembler< Discretization_, Model_ >::residual_variables	(	VectorType &	residual,
		const VectorType &	variables,
		const VectorType &	spatial_solution )

inlineoverrideprotectedvirtual

◆ v_tie()

template<typename Discretization_ , typename Model_ >

template<typename... T>

static constexpr auto DiFfRG::LDG::Assembler< Discretization_, Model_ >::v_tie ( T &&... t )

inlinestaticconstexprprotected

Member Data Documentation

◆ batch_size

template<typename Discretization_ , typename Model_ >

uint DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::batch_size

protected

◆ component_mass_matrix_inverse

template<typename Discretization_ , typename Model_ >

SparseMatrix<NumberType> DiFfRG::LDG::Assembler< Discretization_, Model_ >::component_mass_matrix_inverse

protected

◆ dim

template<typename Discretization_ , typename Model_ >

uint DiFfRG::LDG::Assembler< Discretization_, Model_ >::dim = Discretization::dim

staticconstexpr

◆ discretization

template<typename Discretization_ , typename Model_ >

Discretization& DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::discretization

protected

◆ dof_handler

template<typename Discretization_ , typename Model_ >

const DoFHandler<dim>& DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::dof_handler

protected

◆ dof_handler_list

template<typename Discretization_ , typename Model_ >

std::vector<const DoFHandler<dim> *> DiFfRG::LDG::Assembler< Discretization_, Model_ >::dof_handler_list

protected

◆ EoM

template<typename Discretization_ , typename Model_ >

Point<dim> DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::EoM

mutableprotected

◆ EoM_abs_tol

template<typename Discretization_ , typename Model_ >

const double DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::EoM_abs_tol

protected

◆ EoM_cell

template<typename Discretization_ , typename Model_ >

DoFHandler<dim>::cell_iterator DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::EoM_cell

mutableprotected

◆ EoM_max_iter

template<typename Discretization_ , typename Model_ >

const uint DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::EoM_max_iter

protected

◆ extractor_dof_indices

template<typename Discretization_ , typename Model_ >

std::vector<types::global_dof_index> DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::extractor_dof_indices

protected

◆ extractor_jacobian_u

template<typename Discretization_ , typename Model_ >

FullMatrix<NumberType> DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::extractor_jacobian_u

protected

◆ fe

template<typename Discretization_ , typename Model_ >

const FiniteElement<dim>& DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::fe

protected

◆ j_gu

template<typename Discretization_ , typename Model_ >

array<BlockSparseMatrix<NumberType>, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::j_gu

mutableprotected

◆ j_ug

template<typename Discretization_ , typename Model_ >

array<BlockSparseMatrix<NumberType>, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::j_ug

protected

◆ j_wg

template<typename Discretization_ , typename Model_ >

array<BlockSparseMatrix<NumberType>, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::j_wg

mutableprotected

◆ j_wg_tmp

template<typename Discretization_ , typename Model_ >

array<BlockSparseMatrix<NumberType>, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::j_wg_tmp

mutableprotected

◆ jacobian_tmp

template<typename Discretization_ , typename Model_ >

array<BlockSparseMatrix<NumberType>, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::jacobian_tmp

protected

◆ jacobian_tmp_built

template<typename Discretization_ , typename Model_ >

array<bool, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::jacobian_tmp_built

mutableprotected

◆ ldg_matrix_built

template<typename Discretization_ , typename Model_ >

array<bool, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::ldg_matrix_built

mutableprotected

◆ mapping

template<typename Discretization_ , typename Model_ >

const Mapping<dim>& DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::mapping

protected

◆ mass_matrix

template<typename Discretization_ , typename Model_ >

BlockSparseMatrix<NumberType> DiFfRG::LDG::Assembler< Discretization_, Model_ >::mass_matrix

protected

◆ model

template<typename Discretization_ , typename Model_ >

Model& DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::model

protected

◆ nothing

template<typename Discretization_ , typename Model_ >

int DiFfRG::LDG::Assembler< Discretization_, Model_ >::nothing = 0

staticconstexprprotected

◆ old_EoM_cell

template<typename Discretization_ , typename Model_ >

DoFHandler<dim>::cell_iterator DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::old_EoM_cell

protected

◆ quadrature

template<typename Discretization_ , typename Model_ >

QGauss<dim> DiFfRG::LDG::Assembler< Discretization_, Model_ >::quadrature

protected

◆ quadrature_face

template<typename Discretization_ , typename Model_ >

QGauss<dim - 1> DiFfRG::LDG::Assembler< Discretization_, Model_ >::quadrature_face

protected

◆ sol_vector

template<typename Discretization_ , typename Model_ >

array<BlockVector<NumberType>, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::sol_vector

mutableprotected

◆ sol_vector_tmp

template<typename Discretization_ , typename Model_ >

array<BlockVector<NumberType>, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::sol_vector_tmp

mutableprotected

◆ sol_vector_vec_tmp

template<typename Discretization_ , typename Model_ >

array<Vector<NumberType>, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::sol_vector_vec_tmp

mutableprotected

◆ sparsity_pattern_gu

template<typename Discretization_ , typename Model_ >

array<BlockSparsityPattern, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::sparsity_pattern_gu

protected

◆ sparsity_pattern_jacobian

template<typename Discretization_ , typename Model_ >

BlockSparsityPattern DiFfRG::LDG::Assembler< Discretization_, Model_ >::sparsity_pattern_jacobian

protected

◆ sparsity_pattern_mass

template<typename Discretization_ , typename Model_ >

BlockSparsityPattern DiFfRG::LDG::Assembler< Discretization_, Model_ >::sparsity_pattern_mass

protected

◆ sparsity_pattern_ug

template<typename Discretization_ , typename Model_ >

array<BlockSparsityPattern, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::sparsity_pattern_ug

protected

◆ sparsity_pattern_wg

template<typename Discretization_ , typename Model_ >

array<BlockSparsityPattern, Components::count_fe_subsystems()> DiFfRG::LDG::Assembler< Discretization_, Model_ >::sparsity_pattern_wg

protected

◆ stencil

template<typename Discretization_ , typename Model_ >

uint DiFfRG::LDG::Assembler< Discretization_, Model_ >::stencil = Components::count_fe_subsystems()

staticconstexpr

◆ threads

template<typename Discretization_ , typename Model_ >

uint DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::threads

protected

◆ timings_jacobian

template<typename Discretization_ , typename Model_ >

std::vector<double> DiFfRG::LDG::Assembler< Discretization_, Model_ >::timings_jacobian

protected

◆ timings_reinit

template<typename Discretization_ , typename Model_ >

std::vector<double> DiFfRG::LDG::Assembler< Discretization_, Model_ >::timings_reinit

protected

◆ timings_residual

template<typename Discretization_ , typename Model_ >

std::vector<double> DiFfRG::LDG::Assembler< Discretization_, Model_ >::timings_residual

protected

◆ timings_variable_jacobian

template<typename Discretization_ , typename Model_ >

std::vector<double> DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::timings_variable_jacobian

protected

◆ timings_variable_residual

template<typename Discretization_ , typename Model_ >

std::vector<double> DiFfRG::LDG::LDGAssemblerBase< Discretization_, Model_ >::timings_variable_residual

protected

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

/home/runner/work/DiFfRG/DiFfRG/DiFfRG/include/DiFfRG/discretization/FEM/assembler/ldg.hh

Public Types

Public Member Functions

Static Public Attributes

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Types

Private Member Functions

Detailed Description

Member Typedef Documentation

◆ Base

◆ Components

◆ Discretization

◆ Model

◆ NumberType

◆ VectorType

Constructor & Destructor Documentation

◆ Assembler()

Member Function Documentation

◆ attach_data_output()

◆ average_time_jacobian_assembly()

◆ average_time_reinit()

◆ average_time_residual_assembly()

◆ build_inverse()

◆ build_ldg_jacobian()

◆ build_ldg_sparsity()

◆ build_ldg_vector()

◆ e_tie()

◆ extract()

◆ fe_conv()

◆ fe_more_conv()

◆ get_mass_matrix()

◆ get_sparsity_pattern_jacobian()

◆ jacobian()

◆ jacobian_extractors()

◆ jacobian_mass()

◆ jacobian_variables()

◆ log()

◆ mass()

◆ num_jacobians()

◆ num_reinits()

◆ num_residuals()

◆ readouts()

◆ rebuild_jacobian_sparsity()

◆ rebuild_ldg_jacobian()

◆ rebuild_ldg_vectors()

◆ ref_conv()

◆ refinement_indicator()

◆ reinit()

◆ reinit_vector()

◆ residual()

◆ residual_variables()

◆ v_tie()

Member Data Documentation

◆ batch_size

◆ component_mass_matrix_inverse

◆ dim

◆ discretization

◆ dof_handler

◆ dof_handler_list

◆ EoM

◆ EoM_abs_tol

◆ EoM_cell

◆ EoM_max_iter

◆ extractor_dof_indices

◆ extractor_jacobian_u

◆ fe

◆ j_gu

◆ j_ug

◆ j_wg

◆ j_wg_tmp

◆ jacobian_tmp

◆ jacobian_tmp_built

◆ ldg_matrix_built

◆ mapping

◆ mass_matrix

◆ model

◆ nothing

◆ old_EoM_cell