Assembler< Discretization_, Model_ > Class Template Reference#

DiFfRG: DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ > Class Template Reference
DiFfRG
Discretization Framework for functional Renormalization Group flows
DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ > Class Template Reference

#include <KurganovTadmor.hh>

Inheritance diagram for DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >:
DiFfRG::AbstractAssembler< Discretization_::VectorType, Discretization_::SparseMatrixType, Discretization_::dim >

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::AbstractAssembler< Discretization_::VectorType, Discretization_::SparseMatrixType, Discretization_::dim >
using NumberType = typename get_type::NumberType< Discretization_::VectorType >
 

Public Member Functions

 Assembler (Discretization &discretization, Model &model, const JSONValue &json)
 
virtual void reinit_vector (VectorType &vec) const override
 
virtual IndexSet get_differential_indices () const override
 Obtain the dofs which contain time derivatives.
 
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())
 
virtual void reinit () override
 Reinitialize the assembler. This is necessary if the mesh has changed, e.g. after a mesh refinement.
 
virtual void set_time (double t) override
 Set the current time. The assembler should usually just forward this to the numerical model.
 
virtual const SparsityPattern & get_sparsity_pattern_jacobian () const override
 Obtain the sparsity pattern of the jacobian matrix.
 
virtual const SparseMatrix< NumberType > & get_mass_matrix () const override
 Obtain the mass matrix.
 
virtual void residual_variables (VectorType &residual, const VectorType &variables, const VectorType &) override
 
virtual void jacobian_variables (FullMatrix< NumberType > &jacobian, const VectorType &variables, const VectorType &) override
 
void readouts (DataOutput< dim, VectorType > &data_out, const VectorType &, const VectorType &variables) const
 
virtual void mass (VectorType &mass, const VectorType &solution_global, const VectorType &solution_global_dot, NumberType weight) override
 
virtual void residual (VectorType &residual, const VectorType &solution_global, NumberType weight, const VectorType &solution_global_dot, NumberType weight_mass, const VectorType &variables=VectorType()) override
 
virtual void jacobian_mass (SparseMatrix< NumberType > &jacobian, const VectorType &solution_global, const VectorType &solution_global_dot, NumberType alpha=1., NumberType beta=1.) override
 
virtual void jacobian (SparseMatrix< NumberType > &jacobian, const VectorType &solution_global, NumberType weight, const VectorType &solution_global_dot, NumberType alpha, NumberType beta, const VectorType &variables=VectorType()) override
 
void build_sparsity (SparsityPattern &sparsity_pattern, const DoFHandler< dim > &to_dofh, const DoFHandler< dim > &from_dofh, const int stencil=1, bool add_extractor_dofs=false) const
 
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::AbstractAssembler< Discretization_::VectorType, Discretization_::SparseMatrixType, Discretization_::dim >
virtual void attach_data_output (DataOutput< dim, Discretization_::VectorType > &data_out, const Discretization_::VectorType &solution, const Discretization_::VectorType &variables=Discretization_::VectorType(), const Discretization_::VectorType &dt_solution=Discretization_::VectorType(), const Discretization_::VectorType &residual=Discretization_::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 (Discretization_::VectorType &vector) const=0
 Reinitialize an arbitrary vector so that it has the correct size and structure.
 
virtual void residual_variables (Discretization_::VectorType &residual, const Discretization_::VectorType &variables, const Discretization_::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 Discretization_::VectorType &variables, const Discretization_::VectorType &spatial_solution)
 When coupling the spatial discretization to additional variables, this function should calculate the jacobian for the additional variables.
 
void mass (Discretization_::VectorType &mass, const Discretization_::VectorType &solution_global, NumberType weight)
 Calculates the mass \(m_i(u)\) for an ODE.
 
virtual void mass (Discretization_::VectorType &mass, const Discretization_::VectorType &solution_global, const Discretization_::VectorType &solution_global_dot, NumberType weight)=0
 Calculates the mass \(m(u, \partial_t u)\) for a DAE.
 
void residual (Discretization_::VectorType &residual, const Discretization_::VectorType &solution_global, NumberType weight, NumberType weight_mass, const Discretization_::VectorType &variables=Discretization_::VectorType())
 Calculates the residual for an ODE.
 
virtual void residual (Discretization_::VectorType &residual, const Discretization_::VectorType &solution_global, NumberType weight, const Discretization_::VectorType &solution_global_dot, NumberType weight_mass, const Discretization_::VectorType &variables=Discretization_::VectorType())=0
 Calculates the residual for a DAE.
 
void jacobian_mass (Discretization_::SparseMatrixType &jacobian, const Discretization_::VectorType &solution_global, NumberType mass_weight=1.)
 Calculates the jacobian of the mass function for an ODE.
 
virtual void jacobian_mass (Discretization_::SparseMatrixType &jacobian, const Discretization_::VectorType &solution_global, const Discretization_::VectorType &solution_global_dot, NumberType alpha=1., NumberType beta=1.)=0
 Calculates the jacobian of the mass function for a DAE.
 
void jacobian (Discretization_::SparseMatrixType &jacobian, const Discretization_::VectorType &solution_global, NumberType weight, NumberType mass_weight, const Discretization_::VectorType &variables=Discretization_::VectorType())
 Calculates the jacobian of the residual function for an ODE.
 
virtual void jacobian (Discretization_::SparseMatrixType &jacobian, const Discretization_::VectorType &solution_global, NumberType weight, const Discretization_::VectorType &solution_global_dot, NumberType alpha, NumberType beta, const Discretization_::VectorType &variables=Discretization_::VectorType())=0
 Calculates the jacobian of the residual function for a DAE.
 

Static Public Attributes

static constexpr uint dim = Discretization::dim
 
static constexpr uint n_components = Components::count_fe_functions(0)
 

Protected Member Functions

template<typename... T>
auto fv_tie (T &&...t)
 

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

Discretizationdiscretization
 
Modelmodel
 
const DoFHandler< dim > & dof_handler
 
const Mapping< dim > & mapping
 
const Triangulation< dim > & triangulation
 
const JSONValuejson
 
uint threads
 
const uint batch_size
 
Triangulation< dim >::cell_iterator EoM_cell
 
Triangulation< dim >::cell_iterator old_EoM_cell
 
const double EoM_abs_tol
 
const uint EoM_max_iter
 
const QGauss< dimquadrature
 
SparsityPattern sparsity_pattern_mass
 
SparsityPattern sparsity_pattern_jacobian
 
SparseMatrix< NumberTypemass_matrix
 
std::vector< doubletimings_reinit
 
std::vector< doubletimings_residual
 
std::vector< doubletimings_jacobian
 

Member Typedef Documentation

◆ Components

◆ Discretization

◆ Model

◆ NumberType

◆ VectorType

Constructor & Destructor Documentation

◆ Assembler()

DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::Assembler ( Discretization discretization,
Model model,
const JSONValue json 
)
inline

Member Function Documentation

◆ attach_data_output()

virtual void DiFfRG::FV::KurganovTadmor::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() 
)
inlinevirtual

◆ average_time_jacobian_assembly()

double DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::average_time_jacobian_assembly ( ) const
inline

◆ average_time_reinit()

◆ average_time_residual_assembly()

double DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::average_time_residual_assembly ( ) const
inline

◆ build_sparsity()

void DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::build_sparsity ( SparsityPattern &  sparsity_pattern,
const DoFHandler< dim > &  to_dofh,
const DoFHandler< dim > &  from_dofh,
const int  stencil = 1,
bool  add_extractor_dofs = false 
) const
inline

◆ e_tie()

template<typename... T>
static constexpr auto DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::e_tie ( T &&...  t)
inlinestaticconstexprprotected

◆ fv_tie()

template<typename... T>
auto DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::fv_tie ( T &&...  t)
inlineprotected

◆ get_differential_indices()

virtual IndexSet DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::get_differential_indices ( ) const
inlineoverridevirtual

Obtain the dofs which contain time derivatives.

Returns
IndexSet The indices of the dofs which contain time derivatives

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

◆ get_mass_matrix()

Obtain the mass matrix.

Returns
const SparseMatrixType& The mass matrix

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

◆ get_sparsity_pattern_jacobian()

virtual const SparsityPattern & DiFfRG::FV::KurganovTadmor::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()

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

◆ jacobian_mass()

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

◆ jacobian_variables()

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

◆ log()

◆ mass()

virtual void DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::mass ( VectorType mass,
const VectorType solution_global,
const VectorType solution_global_dot,
NumberType  weight 
)
inlineoverridevirtual

◆ num_jacobians()

◆ num_reinits()

◆ num_residuals()

◆ readouts()

◆ reinit()

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

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

◆ reinit_vector()

virtual void DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::reinit_vector ( VectorType vec) const
inlineoverridevirtual

◆ residual()

virtual void DiFfRG::FV::KurganovTadmor::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

◆ residual_variables()

virtual void DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::residual_variables ( VectorType residual,
const VectorType variables,
const VectorType  
)
inlineoverridevirtual

◆ set_time()

Set the current time. The assembler should usually just forward this to the numerical model.

Parameters
tThe current time

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

◆ v_tie()

template<typename... T>
static constexpr auto DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::v_tie ( T &&...  t)
inlinestaticconstexprprotected

Member Data Documentation

◆ batch_size

◆ dim

◆ discretization

◆ dof_handler

◆ EoM_abs_tol

◆ EoM_cell

◆ EoM_max_iter

◆ json

◆ mapping

◆ mass_matrix

◆ model

◆ n_components

constexpr uint DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::n_components = Components::count_fe_functions(0)
staticconstexpr

◆ old_EoM_cell

◆ quadrature

◆ sparsity_pattern_jacobian

SparsityPattern DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::sparsity_pattern_jacobian
protected

◆ sparsity_pattern_mass

SparsityPattern DiFfRG::FV::KurganovTadmor::Assembler< Discretization_, Model_ >::sparsity_pattern_mass
protected

◆ threads

◆ timings_jacobian

◆ timings_reinit

◆ timings_residual

◆ triangulation


The documentation for this class was generated from the following file:
  • /__w/DiFfRG_current/DiFfRG_current/DiFfRG/include/DiFfRG/discretization/FV/assembler/KurganovTadmor.hh