#include <assembly_richards.hh>

Inheritance diagram for AssemblyRichards< dim >:

Collaboration diagram for AssemblyRichards< dim >:

[legend]

Public Types
typedef std::shared_ptr< RichardsLMH::EqData >	AssemblyDataPtrRichards

Public Types inherited from AssemblyLMH< dim >
typedef std::shared_ptr< DarcyLMH::EqData >	AssemblyDataPtrLMH

Public Member Functions
	AssemblyRichards (AssemblyDataPtrRichards data)

Public Member Functions inherited from AssemblyLMH< dim >
	AssemblyLMH (AssemblyDataPtrLMH data)

	~AssemblyLMH () override

void	fix_velocity (const DHCellAccessor &) override

void	assemble_reconstruct (const DHCellAccessor &dh_cell) override

void	assemble (const DHCellAccessor &dh_cell) override

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

virtual	~AssemblyBase ()

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

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

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

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	begin ()
	Method prepares object before assemblation (e.g. balance, ...). More...

virtual void	end ()
	Method finishes object after assemblation (e.g. balance, ...). 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< 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...

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

virtual	~AssemblyBase ()

Protected Member Functions
void	reset_soil_model (const DHCellAccessor &dh_cell)

double	compute_conductivity (ElementAccessor< 3 > ele)

void	update_water_content (const DHCellAccessor &dh_cell)
	Updates water content in Richards. More...

void	assemble_sides (const DHCellAccessor &dh_cell) override

void	assemble_source_term (const DHCellAccessor &dh_cell) override

void	update_dofs (const DHCellAccessor &dh_cell)

void	postprocess_velocity_specific (const DHCellAccessor &dh_cell, arma::vec &solution, double edge_scale, double edge_source_term) override

Protected Member Functions inherited from AssemblyLMH< dim >
void	assemble_local_system (const DHCellAccessor &dh_cell, bool use_dirichlet_switch)

void	load_local_system (const DHCellAccessor &dh_cell)

void	save_local_system (const DHCellAccessor &dh_cell)

void	set_dofs (const DHCellAccessor &dh_cell)

void	assemble_bc (const DHCellAccessor &dh_cell, bool use_dirichlet_switch)

void	assemble_side_bc (const DHCellSide &side, double cross_section, bool use_dirichlet_switch)

void	assemble_sides_scale (const DHCellAccessor &dh_cell, double scale)

void	assemble_element (FMT_UNUSED const DHCellAccessor &dh_cell)

void	assembly_dim_connections (const DHCellAccessor &dh_cell)

void	postprocess_velocity (const DHCellAccessor &dh_cell, arma::vec &solution)

Protected Member Functions inherited from AssemblyBase< dim >
std::string	print_update_flags (UpdateFlags u) const
	Print update flags to string format. More...

Protected Attributes
AssemblyDataPtrRichards	ad_

bool	genuchten_on

double	cross_section

LocDofVec	cr_disc_dofs
	Dofs of discontinuous fields on element edges. More...

LocDofVec	edge_indices_
	Dofs of discontinuous fields on element edges. More...

Protected Attributes inherited from AssemblyLMH< dim >
FE_RT0< dim >	fe_rt_

QGauss	quad_

FEValues< 3 >	fe_values_

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...

Additional Inherited Members
Static Public Member Functions inherited from AssemblyBase< dim >
template<template< int dim > class Impl, class Data >
static MultidimAssembly	create (Data data)

Static Protected Member Functions inherited from AssemblyLMH< dim >
static unsigned int	size ()

Detailed Description

template<int dim>
class AssemblyRichards< dim >

Prove of concept for general assembly scheme. Ideas:

Local assembly class for given DIM provides assembly methods for:
- local assembly on elements of dimension DIM
- local assembly on edges/faces of neighbouring elemnts of dimension DIM
- local assembly of interactiong elements DIM and OTHER_DIM with OTHER_DIM < DIM
Such class template would be necessary to create particular instance of Assembler class, that takes job of pass through the mesh, precomputing necessary fields, etc. This allows to manage order of assembled entities in arbitrary way in order to use caches and vectorisation efficiently.
Separate local assembly template would be necessary for every pass through the mesh, but it may delegate its actions to a single general local assembly class.
The local assembly class gets an Accessor object for particular domain of the assembly, i.e. ElementAccessor for element assembly, FaceAccessor for interface intergrals assembly, and InteractionAccessor for assembly of interaction of two elements, e.g. for XFEM there may be an interaction between every 3D element within the enrichment with appropriate 1D elemnt on the well. These accessors provides methods to access fields values as well as the DOF indices on the element(s).
Results of the assembly are passed to the global linear algebra objects collected in the RichardSystem class, global indices (which are still local indicies of the single process) provided by accessors are used.

TODO:

finish functional Richards model
move whole internals of the assembly_mh_matrix into local assembly classes
mean while: imporve accessors, imporve mesh and DOF handler interface; possibly create new mesh implementation used in Darcy first and then apply it to other equations

Definition at line 48 of file assembly_richards.hh.

Member Typedef Documentation

template<int dim>

typedef std::shared_ptr<RichardsLMH::EqData> AssemblyRichards< dim >::AssemblyDataPtrRichards

Definition at line 51 of file assembly_richards.hh.

Constructor & Destructor Documentation

template<int dim>

AssemblyRichards< dim >::AssemblyRichards ( AssemblyDataPtrRichards data )

inline

Definition at line 53 of file assembly_richards.hh.

Member Function Documentation

template<int dim>

void AssemblyRichards< dim >::assemble_sides ( const DHCellAccessor & dh_cell )

inlineoverrideprotectedvirtual

Reimplemented from AssemblyLMH< dim >.

Definition at line 124 of file assembly_richards.hh.

template<int dim>

void AssemblyRichards< dim >::assemble_source_term ( const DHCellAccessor & dh_cell )

inlineoverrideprotectedvirtual

Reimplemented from AssemblyLMH< dim >.

Definition at line 139 of file assembly_richards.hh.

template<int dim>

double AssemblyRichards< dim >::compute_conductivity ( ElementAccessor< 3 > ele )

inlineprotected

Definition at line 78 of file assembly_richards.hh.

Here is the caller graph for this function:

template<int dim>

void AssemblyRichards< dim >::postprocess_velocity_specific	(	const DHCellAccessor &	dh_cell,
		arma::vec &	solution,
		double	edge_scale,
		double	edge_source_term
	)

inlineoverrideprotectedvirtual

Reimplemented from AssemblyLMH< dim >.

Definition at line 203 of file assembly_richards.hh.

template<int dim>

void AssemblyRichards< dim >::reset_soil_model ( const DHCellAccessor & dh_cell )

inlineprotected

Definition at line 62 of file assembly_richards.hh.

Here is the caller graph for this function:

template<int dim>

void AssemblyRichards< dim >::update_dofs ( const DHCellAccessor & dh_cell )

inlineprotected

Updates DoFs for edge pressure vector (dh CR) and for water content vector (dh CR_disc) Be sure to call it before update_water_content().

Definition at line 197 of file assembly_richards.hh.

Here is the caller graph for this function:

template<int dim>

void AssemblyRichards< dim >::update_water_content ( const DHCellAccessor & dh_cell )

inlineprotectedvirtual

Updates water content in Richards.

Reimplemented from AssemblyLMH< dim >.

Definition at line 95 of file assembly_richards.hh.

Here is the caller graph for this function:

Member Data Documentation

template<int dim>

AssemblyDataPtrRichards AssemblyRichards< dim >::ad_

protected

Definition at line 225 of file assembly_richards.hh.

template<int dim>

LocDofVec AssemblyRichards< dim >::cr_disc_dofs

protected

Dofs of discontinuous fields on element edges.

Definition at line 229 of file assembly_richards.hh.

template<int dim>

double AssemblyRichards< dim >::cross_section

protected

Definition at line 228 of file assembly_richards.hh.

template<int dim>

LocDofVec AssemblyRichards< dim >::edge_indices_

protected

Dofs of discontinuous fields on element edges.

Definition at line 230 of file assembly_richards.hh.

template<int dim>

bool AssemblyRichards< dim >::genuchten_on

protected

Definition at line 227 of file assembly_richards.hh.

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

/opt/flow123d/flow123d/src/flow/assembly_richards.hh

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

template<int dim> class AssemblyRichards< dim >

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

template<int dim>
class AssemblyRichards< dim >