#include <assembly_richards.hh>
Inheritance diagram for AssemblyRichards< dim >:
Collaboration diagram for AssemblyRichards< dim >:
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
◆ AssemblyDataPtrRichards
template<int dim>
| typedef std::shared_ptr<RichardsLMH::EqData> AssemblyRichards< dim >::AssemblyDataPtrRichards |
Definition at line 51 of file assembly_richards.hh.
Constructor & Destructor Documentation
◆ AssemblyRichards()
template<int dim>
|
inline |
Definition at line 53 of file assembly_richards.hh.
Member Function Documentation
◆ assemble_sides()
template<int dim>
|
inlineoverrideprotectedvirtual |
Reimplemented from AssemblyLMH< dim >.
Definition at line 124 of file assembly_richards.hh.
◆ assemble_source_term()
template<int dim>
|
inlineoverrideprotectedvirtual |
Reimplemented from AssemblyLMH< dim >.
Definition at line 139 of file assembly_richards.hh.
◆ compute_conductivity()
template<int dim>
|
inlineprotected |
◆ postprocess_velocity_specific()
template<int dim>
|
inlineoverrideprotectedvirtual |
Reimplemented from AssemblyLMH< dim >.
Definition at line 203 of file assembly_richards.hh.
◆ reset_soil_model()
template<int dim>
|
inlineprotected |
◆ update_dofs()
template<int dim>
|
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:
◆ update_water_content()
template<int dim>
|
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
◆ ad_
template<int dim>
|
protected |
Definition at line 225 of file assembly_richards.hh.
◆ cr_disc_dofs
template<int dim>
|
protected |
Dofs of discontinuous fields on element edges.
Definition at line 229 of file assembly_richards.hh.
◆ cross_section
template<int dim>
|
protected |
Definition at line 228 of file assembly_richards.hh.
◆ edge_indices_
template<int dim>
|
protected |
Dofs of discontinuous fields on element edges.
Definition at line 230 of file assembly_richards.hh.
◆ genuchten_on
template<int dim>
|
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
