Flow123d  master-f44eb46
Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | Friends | List of all members
MHMatrixAssemblyLMH< dim > Class Template Reference

#include <assembly_lmh.hh>

Inheritance diagram for MHMatrixAssemblyLMH< dim >:
Inheritance graph
[legend]
Collaboration diagram for MHMatrixAssemblyLMH< dim >:
Collaboration graph
[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...
 
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< BulkPointbulk_points (unsigned int element_patch_idx) const
 Return BulkPoint range of appropriate dimension. More...
 
Range< EdgePointedge_points (const DHCellSide &cell_side) const
 Return EdgePoint range of appropriate dimension. More...
 
Range< CouplingPointcoupling_points (const DHCellSide &cell_side) const
 Return CouplingPoint range of appropriate dimension. More...
 
Range< BoundaryPointboundary_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...
 

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...
 
EqFieldseq_fields_
 Data objects shared with DarcyFlow. More...
 
EqDataeq_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 >
Quadraturequad_
 Quadrature used in assembling methods. More...
 
Quadraturequad_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...
 
ElementCacheMapelement_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

◆ 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

◆ 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>
template<template< IntDim... > class DimAssembly>
friend class GenericAssembly
friend

Definition at line 789 of file assembly_lmh.hh.

Member Data Documentation

◆ bulk_local_idx_

template<unsigned int dim>
unsigned int MHMatrixAssemblyLMH< dim >::bulk_local_idx_
protected

Local idx of bulk element.

Definition at line 779 of file assembly_lmh.hh.

◆ coef_

template<unsigned int dim>
double MHMatrixAssemblyLMH< dim >::coef_
protected

Definition at line 776 of file assembly_lmh.hh.

◆ cross_section_

template<unsigned int dim>
double MHMatrixAssemblyLMH< dim >::cross_section_
protected

Precomputed cross-section value.

Definition at line 778 of file assembly_lmh.hh.

◆ edge_row_

template<unsigned int dim>
unsigned int MHMatrixAssemblyLMH< dim >::edge_row_
protected

Helper indices in boundary assembly.

Definition at line 780 of file assembly_lmh.hh.

◆ edge_scale_

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.

◆ eq_data_

template<unsigned int dim>
EqData* MHMatrixAssemblyLMH< dim >::eq_data_
protected

Definition at line 763 of file assembly_lmh.hh.

◆ eq_fields_

template<unsigned int dim>
EqFields* MHMatrixAssemblyLMH< dim >::eq_fields_
protected

Data objects shared with DarcyFlow.

Definition at line 762 of file assembly_lmh.hh.

◆ fe_

template<unsigned int dim>
shared_ptr<FiniteElement<dim> > MHMatrixAssemblyLMH< dim >::fe_
protected

Finite element for the solution of the advection-diffusion equation.

Definition at line 770 of file assembly_lmh.hh.

◆ fe_rt_

template<unsigned int dim>
FE_RT0<dim> MHMatrixAssemblyLMH< dim >::fe_rt_
protected

Assembly volume integrals.

Definition at line 766 of file assembly_lmh.hh.

◆ fe_values_

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.

◆ loc_schur_

template<unsigned int dim>
LocalSystem MHMatrixAssemblyLMH< dim >::loc_schur_
protected

Definition at line 786 of file assembly_lmh.hh.

◆ ngh_value_

template<unsigned int dim>
double MHMatrixAssemblyLMH< dim >::ngh_value_
protected

Precomputed ngh value.

Definition at line 783 of file assembly_lmh.hh.

◆ nv_

template<unsigned int dim>
arma::vec3 MHMatrixAssemblyLMH< dim >::nv_
protected

Normal vector.

Definition at line 782 of file assembly_lmh.hh.

◆ quad_rt_

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.

◆ time_term_diag_

template<unsigned int dim>
double MHMatrixAssemblyLMH< dim >::time_term_diag_
protected

Definition at line 777 of file assembly_lmh.hh.

◆ time_term_rhs_

template<unsigned int dim>
double MHMatrixAssemblyLMH< dim >::time_term_rhs_
protected

Variables used in compute time term (unsteady)

Definition at line 777 of file assembly_lmh.hh.

◆ type_

template<unsigned int dim>
DarcyLMH::EqFields::BC_Type MHMatrixAssemblyLMH< dim >::type_
protected

Type of boundary condition.

Definition at line 781 of file assembly_lmh.hh.

◆ used_fields_

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: