#include <assembly_lmh.hh>

Inheritance diagram for MHMatrixAssemblyLMH< dim >:

[legend]

Collaboration diagram for MHMatrixAssemblyLMH< dim >:

[legend]

Public Types
typedef DarcyLMH::EqFields	EqFields

typedef DarcyLMH::EqData	EqData

Public Types inherited from AssemblyBase< dim >
typedef GenericAssemblyBase::BulkIntegralData	BulkIntegralData

typedef GenericAssemblyBase::EdgeIntegralData	EdgeIntegralData

typedef GenericAssemblyBase::CouplingIntegralData	CouplingIntegralData

typedef GenericAssemblyBase::BoundaryIntegralData	BoundaryIntegralData

Public Member Functions
	DECLARE_EXCEPTION (ExcBCNotSupported,<< "BC type not supported.\n")

	MHMatrixAssemblyLMH (EqFields eq_fields, EqData eq_data)
	Constructor. More...

virtual	~MHMatrixAssemblyLMH ()
	Destructor. More...

void	initialize (ElementCacheMap *element_cache_map)
	Initialize auxiliary vectors and other data members. More...

void	cell_integral (DHCellAccessor cell, unsigned int element_patch_idx)

void	boundary_side_integral (DHCellSide cell_side)

void	dimjoin_intergral (DHCellAccessor cell_lower_dim, DHCellSide neighb_side)

void	begin () override
	Implements `AssemblyBase::begin`. More...

void	end () override
	Implements `AssemblyBase::end`. More...

Public Member Functions inherited from AssemblyBase< dim >
	AssemblyBase (unsigned int quad_order)
	Constructor. More...

virtual	~AssemblyBase ()
	Destructor. More...

virtual void	cell_integral (FMT_UNUSED DHCellAccessor cell, FMT_UNUSED unsigned int element_patch_idx)
	Assembles the volume integrals on cell. More...

virtual void	boundary_side_integral (FMT_UNUSED DHCellSide cell_side)
	Assembles the fluxes on the boundary. More...

virtual void	edge_integral (FMT_UNUSED RangeConvert< DHEdgeSide, DHCellSide > edge_side_range)
	Assembles the fluxes between sides on the edge. More...

virtual void	dimjoin_intergral (FMT_UNUSED DHCellAccessor cell_lower_dim, FMT_UNUSED DHCellSide neighb_side)
	Assembles the fluxes between elements of different dimensions. More...

virtual void	patch_reinit (FMT_UNUSED std::array< PatchElementsList, 4 > &patch_elements)

int	n_active_integrals () const
	Getter of active_integrals. More...

void	create_integrals (std::shared_ptr< EvalPoints > eval_points, AssemblyIntegrals &integrals)
	Create integrals according to dim of assembly object. More...

Range< BulkPoint >	bulk_points (unsigned int element_patch_idx) const
	Return BulkPoint range of appropriate dimension. More...

Range< EdgePoint >	edge_points (const DHCellSide &cell_side) const
	Return EdgePoint range of appropriate dimension. More...

Range< CouplingPoint >	coupling_points (const DHCellSide &cell_side) const
	Return CouplingPoint range of appropriate dimension. More...

Range< BoundaryPoint >	boundary_points (const DHCellSide &cell_side) const
	Return BoundaryPoint range of appropriate dimension. More...

virtual void	assemble_cell_integrals (const RevertableList< BulkIntegralData > &bulk_integral_data)
	Assembles the cell integrals for the given dimension. More...

void	assemble_boundary_side_integrals (const RevertableList< BoundaryIntegralData > &boundary_integral_data)
	Assembles the boundary side integrals for the given dimension. More...

void	assemble_edge_integrals (const RevertableList< EdgeIntegralData > &edge_integral_data)
	Assembles the edge integrals for the given dimension. More...

void	assemble_neighbour_integrals (const RevertableList< CouplingIntegralData > &coupling_integral_data)
	Assembles the neighbours integrals for the given dimension. More...

virtual void	add_patch_bulk_points (FMT_UNUSED const RevertableList< BulkIntegralData > &bulk_integral_data)
	Register cell points of volume integral. More...

virtual void	add_patch_bdr_side_points (FMT_UNUSED const RevertableList< BoundaryIntegralData > &boundary_integral_data)
	Register side points of boundary side integral. More...

virtual void	add_patch_edge_points (FMT_UNUSED const RevertableList< EdgeIntegralData > &edge_integral_data)
	Register side points of edge integral. More...

virtual void	add_patch_coupling_integrals (FMT_UNUSED const RevertableList< CouplingIntegralData > &coupling_integral_data)
	Register bulk and side points of coupling integral. More...

Static Public Member Functions
static constexpr const char *	name ()

Protected Member Functions
void	begin_mh_matrix ()
	Common code of begin method of MH matrix assembly (Darcy and Richards) More...

void	end_mh_matrix ()
	Common code of end method of MH matrix assembly (Darcy and Richards) More...

void	begin_reconstruct_schur ()
	Common code of begin method of Reconstruct Schur assembly (Darcy and Richards) More...

void	end_reconstruct_schur ()
	Common code of end method of Reconstruct Schur assembly (Darcy and Richards) More...

void	asm_sides (const DHCellAccessor &cell, BulkPoint &p, double conductivity)
	Part of cell_integral method, common in all descendants. More...

void	asm_element ()
	Part of cell_integral method, common in all descendants. More...

void	asm_source_term_darcy (const DHCellAccessor &cell, BulkPoint &p)
	Part of cell_integral method, specialized in Darcy equation. More...

void	precompute_boundary_side (DHCellSide &cell_side, BoundaryPoint &p_side, BulkPoint &p_bdr)
	Precompute values used in boundary side integral on given DHCellSide. More...

void	use_dirichlet_switch (DHCellSide &cell_side, const ElementAccessor< 3 > &b_ele, BulkPoint &p_bdr)
	Update BC switch dirichlet in MH matrix assembly if BC type is seepage. More...

void	boundary_side_integral_in (DHCellSide cell_side, const ElementAccessor< 3 > &b_ele, BulkPoint &p_bdr)

void	set_dofs ()
	Precompute loc_system and loc_schur data members. More...

void	set_loc_schur (const DHCellAccessor dh_cr_cell)
	Precompute loc_schur data member of given cell. More...

unsigned int	ngh_idx (DHCellAccessor &dh_cell, DHCellSide &neighb_side)
	Temporary method find neighbour index in higher-dim cell. More...

arma::vec	load_local_system (const DHCellAccessor &dh_cell)

void	postprocess_velocity (const DHCellAccessor &dh_cell, BulkPoint &p)

void	postprocess_velocity_darcy (const DHCellAccessor &dh_cell, BulkPoint &p, arma::vec &solution)
	Postprocess velocity during loading of local system and after calculating of cell integral. More...

Protected Member Functions inherited from AssemblyBase< dim >
	AssemblyBase ()

std::string	print_update_flags (UpdateFlags u) const
	Print update flags to string format. More...

Protected Attributes
FieldSet	used_fields_
	Sub field set contains fields used in calculation. More...

EqFields *	eq_fields_
	Data objects shared with DarcyFlow. More...

EqData *	eq_data_

FE_RT0< dim >	fe_rt_
	Assembly volume integrals. More...

QGauss	quad_rt_

FEValues< 3 >	fe_values_

shared_ptr< FiniteElement< dim > >	fe_
	Finite element for the solution of the advection-diffusion equation. More...

FEValues< 3 >	fe_values_side_
	FEValues of object (of P disc finite element type) More...

arma::vec	reconstructed_solution_
	Vector for reconstructed solution (velocity and pressure on element) from Schur complement. More...

double	coef_

double	source_term_
	Variables used in compute lumped source. More...

double	time_term_diag_

double	time_term_

double	time_term_rhs_
	Variables used in compute time term (unsteady) More...

double	cross_section_
	Precomputed cross-section value. More...

unsigned int	bulk_local_idx_
	Local idx of bulk element. More...

unsigned int	sidx_

unsigned int	side_row_

unsigned int	edge_row_
	Helper indices in boundary assembly. More...

DarcyLMH::EqFields::BC_Type	type_
	Type of boundary condition. More...

arma::vec3	nv_
	Normal vector. More...

double	ngh_value_
	Precomputed ngh value. More...

double	edge_scale_

double	edge_source_term_
	Precomputed values in postprocess_velocity. More...

LocalSystem	loc_schur_

Protected Attributes inherited from AssemblyBase< dim >
Quadrature *	quad_
	Quadrature used in assembling methods. More...

Quadrature *	quad_low_
	Quadrature used in assembling methods (dim-1). More...

int	active_integrals_
	Holds mask of active integrals. More...

DimIntegrals	integrals_
	Set of used integrals. More...

ElementCacheMap *	element_cache_map_
	ElementCacheMap shared with GenericAssembly object. More...

Friends
template<template< IntDim... > class DimAssembly>
class	GenericAssembly

Detailed Description

template<unsigned int dim>
class MHMatrixAssemblyLMH< dim >

Definition at line 116 of file assembly_lmh.hh.

Member Typedef Documentation

◆ EqData

template<unsigned int dim>

typedef DarcyLMH::EqData MHMatrixAssemblyLMH< dim >::EqData

Definition at line 120 of file assembly_lmh.hh.

◆ EqFields

template<unsigned int dim>

typedef DarcyLMH::EqFields MHMatrixAssemblyLMH< dim >::EqFields

Definition at line 119 of file assembly_lmh.hh.

Constructor & Destructor Documentation

◆ MHMatrixAssemblyLMH()

template<unsigned int dim>

MHMatrixAssemblyLMH< dim >::MHMatrixAssemblyLMH	(	EqFields *	eq_fields,
		EqData *	eq_data
	)

inline

Constructor.

Definition at line 127 of file assembly_lmh.hh.

◆ ~MHMatrixAssemblyLMH()

template<unsigned int dim>

virtual MHMatrixAssemblyLMH< dim >::~MHMatrixAssemblyLMH ( )

inlinevirtual

Destructor.

Definition at line 146 of file assembly_lmh.hh.

Member Function Documentation

◆ asm_element()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::asm_element ( )

inlineprotected

Part of cell_integral method, common in all descendants.

Definition at line 410 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ asm_sides()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::asm_sides	(	const DHCellAccessor &	cell,
		BulkPoint &	p,
		double	conductivity
	)

inlineprotected

Part of cell_integral method, common in all descendants.

Definition at line 361 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ asm_source_term_darcy()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::asm_source_term_darcy	(	const DHCellAccessor &	cell,
		BulkPoint &	p
	)

inlineprotected

Part of cell_integral method, specialized in Darcy equation.

Definition at line 427 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ begin()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::begin ( )

inlineoverridevirtual

Implements AssemblyBase::begin.

Reimplemented from AssemblyBase< dim >.

Reimplemented in ReconstructSchurAssemblyRichards< dim >, MHMatrixAssemblyRichards< dim >, and ReconstructSchurAssemblyLMH< dim >.

Definition at line 261 of file assembly_lmh.hh.

◆ begin_mh_matrix()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::begin_mh_matrix ( )

inlineprotected

Common code of begin method of MH matrix assembly (Darcy and Richards)

Definition at line 275 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ begin_reconstruct_schur()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::begin_reconstruct_schur ( )

inlineprotected

Common code of begin method of Reconstruct Schur assembly (Darcy and Richards)

Definition at line 327 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ boundary_side_integral()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::boundary_side_integral ( DHCellSide cell_side )

inline

Assembles between boundary element and corresponding side on bulk element.

Override in descendants.

Definition at line 199 of file assembly_lmh.hh.

◆ boundary_side_integral_in()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::boundary_side_integral_in	(	DHCellSide	cell_side,
		const ElementAccessor< 3 > &	b_ele,
		BulkPoint &	p_bdr
	)

inlineprotected

Common code of boundary_side_integral.

Definition at line 530 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ cell_integral()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::cell_integral	(	DHCellAccessor	cell,
		unsigned int	element_patch_idx
	)

inline

Integral over element.

Override in descendants.

Definition at line 180 of file assembly_lmh.hh.

◆ DECLARE_EXCEPTION()

template<unsigned int dim>

MHMatrixAssemblyLMH< dim >::DECLARE_EXCEPTION	(	ExcBCNotSupported	,
		<< "BC type not supported.\n"
	)

◆ dimjoin_intergral()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::dimjoin_intergral	(	DHCellAccessor	cell_lower_dim,
		DHCellSide	neighb_side
	)

inline

Assembles the fluxes between elements of different dimensions.

Common in all descendants.

Definition at line 223 of file assembly_lmh.hh.

◆ end()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::end ( )

inlineoverridevirtual

Implements AssemblyBase::end.

Reimplemented from AssemblyBase< dim >.

Reimplemented in ReconstructSchurAssemblyRichards< dim >, MHMatrixAssemblyRichards< dim >, and ReconstructSchurAssemblyLMH< dim >.

Definition at line 268 of file assembly_lmh.hh.

◆ end_mh_matrix()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::end_mh_matrix ( )

inlineprotected

Common code of end method of MH matrix assembly (Darcy and Richards)

Definition at line 292 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ end_reconstruct_schur()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::end_reconstruct_schur ( )

inlineprotected

Common code of end method of Reconstruct Schur assembly (Darcy and Richards)

Definition at line 336 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ initialize()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::initialize ( ElementCacheMap * element_cache_map )

inline

Initialize auxiliary vectors and other data members.

Definition at line 149 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ load_local_system()

template<unsigned int dim>

arma::vec MHMatrixAssemblyLMH< dim >::load_local_system ( const DHCellAccessor & dh_cell )

inlineprotected

Loads the local system from a map: element index -> LocalSystem, if it exits, or if the full solution is not yet reconstructed, and reconstructs the full solution on the element. Currently used only for seepage BC.

Definition at line 695 of file assembly_lmh.hh.

◆ name()

template<unsigned int dim>

static constexpr const char* MHMatrixAssemblyLMH< dim >::name ( )

inlinestaticconstexpr

Definition at line 124 of file assembly_lmh.hh.

◆ ngh_idx()

template<unsigned int dim>

unsigned int MHMatrixAssemblyLMH< dim >::ngh_idx	(	DHCellAccessor &	dh_cell,
		DHCellSide &	neighb_side
	)

inlineprotected

Temporary method find neighbour index in higher-dim cell.

Definition at line 680 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ postprocess_velocity()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::postprocess_velocity	(	const DHCellAccessor &	dh_cell,
		BulkPoint &	p
	)

inlineprotected

Precompute edge_scale and edge_source_term.

This method must be calls in methods postprocess_velocity_darcy and postprocess_velocity_richards

Definition at line 723 of file assembly_lmh.hh.

◆ postprocess_velocity_darcy()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::postprocess_velocity_darcy	(	const DHCellAccessor &	dh_cell,
		BulkPoint &	p,
		arma::vec &	solution
	)

inlineprotected

Postprocess velocity during loading of local system and after calculating of cell integral.

Definition at line 736 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ precompute_boundary_side()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::precompute_boundary_side	(	DHCellSide &	cell_side,
		BoundaryPoint &	p_side,
		BulkPoint &	p_bdr
	)

inlineprotected

Precompute values used in boundary side integral on given DHCellSide.

Definition at line 475 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ set_dofs()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::set_dofs ( )

inlineprotected

Precompute loc_system and loc_schur data members.

Definition at line 597 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ set_loc_schur()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::set_loc_schur ( const DHCellAccessor dh_cr_cell )

inlineprotected

Precompute loc_schur data member of given cell.

Definition at line 647 of file assembly_lmh.hh.

Here is the caller graph for this function:

◆ use_dirichlet_switch()

template<unsigned int dim>

void MHMatrixAssemblyLMH< dim >::use_dirichlet_switch	(	DHCellSide &	cell_side,
		const ElementAccessor< 3 > &	b_ele,
		BulkPoint &	p_bdr
	)

inlineprotected

Update BC switch dirichlet in MH matrix assembly if BC type is seepage.

Definition at line 488 of file assembly_lmh.hh.

Here is the caller graph for this function:

Friends And Related Function Documentation

◆ GenericAssembly

template<unsigned int dim>

double MHMatrixAssemblyLMH< dim >::edge_scale_

protected

Definition at line 784 of file assembly_lmh.hh.

◆ edge_source_term_

template<unsigned int dim>

double MHMatrixAssemblyLMH< dim >::edge_source_term_

protected

Precomputed values in postprocess_velocity.

Definition at line 784 of file assembly_lmh.hh.

template<unsigned int dim>

FEValues<3> MHMatrixAssemblyLMH< dim >::fe_values_

protected

Definition at line 768 of file assembly_lmh.hh.

◆ fe_values_side_

template<unsigned int dim>

FEValues<3> MHMatrixAssemblyLMH< dim >::fe_values_side_

protected

FEValues of object (of P disc finite element type)

Definition at line 771 of file assembly_lmh.hh.

template<unsigned int dim>

QGauss MHMatrixAssemblyLMH< dim >::quad_rt_

protected

Definition at line 767 of file assembly_lmh.hh.

◆ reconstructed_solution_

template<unsigned int dim>

arma::vec MHMatrixAssemblyLMH< dim >::reconstructed_solution_

protected

Vector for reconstructed solution (velocity and pressure on element) from Schur complement.

Definition at line 774 of file assembly_lmh.hh.

◆ side_row_

template<unsigned int dim>

unsigned int MHMatrixAssemblyLMH< dim >::side_row_

protected

Definition at line 780 of file assembly_lmh.hh.

◆ sidx_

template<unsigned int dim>

unsigned int MHMatrixAssemblyLMH< dim >::sidx_

protected

Definition at line 780 of file assembly_lmh.hh.

◆ source_term_

template<unsigned int dim>

double MHMatrixAssemblyLMH< dim >::source_term_

protected

Variables used in compute lumped source.

Definition at line 776 of file assembly_lmh.hh.

◆ time_term_

template<unsigned int dim>

double MHMatrixAssemblyLMH< dim >::time_term_

protected

Definition at line 777 of file assembly_lmh.hh.

template<unsigned int dim>

FieldSet MHMatrixAssemblyLMH< dim >::used_fields_

protected

Sub field set contains fields used in calculation.

Definition at line 759 of file assembly_lmh.hh.

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

/home/runner/work/flow123d/flow123d/src/flow/assembly_lmh.hh

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Friends

Detailed Description

template<unsigned int dim> class MHMatrixAssemblyLMH< dim >

Member Typedef Documentation

◆ EqData

◆ EqFields

Constructor & Destructor Documentation

◆ MHMatrixAssemblyLMH()

◆ ~MHMatrixAssemblyLMH()

Member Function Documentation

◆ asm_element()

◆ asm_sides()

◆ asm_source_term_darcy()

◆ begin()

◆ begin_mh_matrix()

◆ begin_reconstruct_schur()

◆ boundary_side_integral()

◆ boundary_side_integral_in()

◆ cell_integral()

◆ DECLARE_EXCEPTION()

◆ dimjoin_intergral()

◆ end()

◆ end_mh_matrix()

◆ end_reconstruct_schur()

◆ initialize()

◆ load_local_system()

◆ name()

◆ ngh_idx()

◆ postprocess_velocity()

◆ postprocess_velocity_darcy()

◆ precompute_boundary_side()

◆ set_dofs()

◆ set_loc_schur()

◆ use_dirichlet_switch()

Friends And Related Function Documentation

◆ GenericAssembly

Member Data Documentation

◆ bulk_local_idx_

◆ coef_

◆ cross_section_

◆ edge_row_

◆ edge_scale_

◆ edge_source_term_

◆ eq_data_

◆ eq_fields_

◆ fe_

◆ fe_rt_

◆ fe_values_

◆ fe_values_side_

◆ loc_schur_

◆ ngh_value_

◆ nv_

◆ quad_rt_

◆ reconstructed_solution_

◆ side_row_

◆ sidx_

◆ source_term_

◆ time_term_

◆ time_term_diag_

◆ time_term_rhs_

◆ type_

◆ used_fields_

template<unsigned int dim>
class MHMatrixAssemblyLMH< dim >