Detailed Description

template<typename LocalAssemblerInterface, template< typename, typename, int > class LocalAssemblerDataMatrix, template< typename, typename, int > class LocalAssemblerDataMatrixNearFracture, template< typename, typename, int > class LocalAssemblerDataFracture, int DisplacementDim, typename... ConstructorArgs>
class ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >

The LocalDataInitializer is a functor creating a local assembler data with corresponding to the mesh element type shape functions and calling initialization of the new local assembler data. For example for MeshLib::Quad a local assembler data with template argument NumLib::ShapeQuad4 is created.

Attention: This is modified version of the ProcessLib::LocalDataInitializer class which does not include line elements, allows only shapefunction of order 2, and provides additional template argument DisplacementDim.

Definition at line 44 of file LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h.

#include <LocalDataInitializer.h>

Collaboration diagram for ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >:

[legend]

Classes
struct	IsElementEnabled

Public Types
using	LADataIntfPtr = std::unique_ptr<LocalAssemblerInterface>

Public Member Functions
	LocalDataInitializer (NumLib::LocalToGlobalIndexMap const &dof_table, NumLib::IntegrationOrder const integration_order)
LADataIntfPtr	operator() (std::size_t const id, MeshLib::Element const &mesh_item, ConstructorArgs &&... args) const

Private Types
using	LADataBuilder

Static Private Member Functions
template<typename ShapeFunctionDisplacement, typename ShapeFunctionPressure>
static LADataBuilder	makeLocalAssemblerBuilder (NumLib::IntegrationOrder const integration_order)

Private Attributes
std::unordered_map< std::type_index, LADataBuilder >	_builder
	Mapping of element types to local assembler constructors.
NumLib::LocalToGlobalIndexMap const &	_dof_table

Member Typedef Documentation

◆ LADataBuilder

template<typename LocalAssemblerInterface, template< typename, typename, int > class LocalAssemblerDataMatrix, template< typename, typename, int > class LocalAssemblerDataMatrixNearFracture, template< typename, typename, int > class LocalAssemblerDataFracture, int DisplacementDim, typename... ConstructorArgs>

using ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >::LADataBuilder

private

Initial value:

 std::function<LADataIntfPtr(
        MeshLib::Element const& e,
        std::size_t const n_variables,
        std::size_t const local_matrix_size,
        std::vector<unsigned> const& dofIndex_to_localIndex,
        ConstructorArgs&&...)>

Definition at line 174 of file LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h.

◆ LADataIntfPtr

template<typename LocalAssemblerInterface, template< typename, typename, int > class LocalAssemblerDataMatrix, template< typename, typename, int > class LocalAssemblerDataMatrixNearFracture, template< typename, typename, int > class LocalAssemblerDataFracture, int DisplacementDim, typename... ConstructorArgs>

using ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >::LADataIntfPtr = std::unique_ptr<LocalAssemblerInterface>

Definition at line 63 of file LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h.

Constructor & Destructor Documentation

◆ LocalDataInitializer()

template<typename LocalAssemblerInterface, template< typename, typename, int > class LocalAssemblerDataMatrix, template< typename, typename, int > class LocalAssemblerDataMatrixNearFracture, template< typename, typename, int > class LocalAssemblerDataFracture, int DisplacementDim, typename... ConstructorArgs>

ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >::LocalDataInitializer	(	NumLib::LocalToGlobalIndexMap const &	dof_table,
		NumLib::IntegrationOrder const	integration_order )

inlineexplicit

Definition at line 65 of file LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h.

        : _dof_table(dof_table)
    {
        using EnabledElementTraits =
            decltype(BaseLib::TMP::filter<EnabledElementTraitsLagrange>(
                std::declval<IsElementEnabled>()));
 
        BaseLib::TMP::foreach<EnabledElementTraits>(
            [this, integration_order]<typename ET>(ET*)
            {
                using MeshElement = typename ET::Element;
                using ShapeFunction = typename ET::ShapeFunction;
                using LowerOrderShapeFunction =
                    typename ET::LowerOrderShapeFunction;
 
                _builder[std::type_index(typeid(MeshElement))] =
                    makeLocalAssemblerBuilder<ShapeFunction,
                                              LowerOrderShapeFunction>(
                        integration_order);
            });
    }

References _builder, _dof_table, BaseLib::TMP::filter(), BaseLib::TMP::foreach(), and makeLocalAssemblerBuilder().

Member Function Documentation

◆ makeLocalAssemblerBuilder()

template<typename LocalAssemblerInterface, template< typename, typename, int > class LocalAssemblerDataMatrix, template< typename, typename, int > class LocalAssemblerDataMatrixNearFracture, template< typename, typename, int > class LocalAssemblerDataFracture, int DisplacementDim, typename... ConstructorArgs>

template<typename ShapeFunctionDisplacement, typename ShapeFunctionPressure>

LADataBuilder ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >::makeLocalAssemblerBuilder ( NumLib::IntegrationOrder const integration_order )

inlinestaticprivate

Generates a function that creates a new LocalAssembler of type LAData<ShapeFunctionDisplacement>. Only functions with shape function's dimension less or equal to the global dimension are instantiated, e.g. following combinations of shape functions and global dimensions: (Line2, 1), (Line2, 2), (Line2, 3), (Hex20, 3) but not (Hex20, 2) or (Hex20, 1).

Definition at line 189 of file LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h.

    {
        return [integration_order](
                   MeshLib::Element const& e,
                   std::size_t const n_variables,
                   std::size_t const local_matrix_size,
                   std::vector<unsigned> const& dofIndex_to_localIndex,
                   ConstructorArgs&&... args)
        {
            auto const& integration_method = NumLib::IntegrationMethodRegistry::
                template getIntegrationMethod<
                    typename ShapeFunctionDisplacement::MeshElement>(
                    integration_order);
 
            if (e.getDimension() == DisplacementDim)
            {
                if (n_variables == 2)
                {
                    return LADataIntfPtr{
                        new LocalAssemblerDataMatrix<ShapeFunctionDisplacement,
                                                     ShapeFunctionPressure,
                                                     DisplacementDim>{
                            e, n_variables, local_matrix_size,
                            dofIndex_to_localIndex, integration_method,
                            std::forward<ConstructorArgs>(args)...}};
                }
                return LADataIntfPtr{new LocalAssemblerDataMatrixNearFracture<
                    ShapeFunctionDisplacement, ShapeFunctionPressure,
                    DisplacementDim>{e, n_variables, local_matrix_size,
                                     dofIndex_to_localIndex, integration_method,
                                     std::forward<ConstructorArgs>(args)...}};
            }
            return LADataIntfPtr{new LocalAssemblerDataFracture<
                ShapeFunctionDisplacement, ShapeFunctionPressure,
                DisplacementDim>{e, local_matrix_size, dofIndex_to_localIndex,
                                 integration_method,
                                 std::forward<ConstructorArgs>(args)...}};
        };
    }

References MeshLib::Element::getDimension().

Referenced by LocalDataInitializer().

◆ operator()()

template<typename LocalAssemblerInterface, template< typename, typename, int > class LocalAssemblerDataMatrix, template< typename, typename, int > class LocalAssemblerDataMatrixNearFracture, template< typename, typename, int > class LocalAssemblerDataFracture, int DisplacementDim, typename... ConstructorArgs>

LADataIntfPtr ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >::operator()	(	std::size_t const	id,
		MeshLib::Element const &	mesh_item,
		ConstructorArgs &&...	args ) const

inline

Returns data pointer to the newly created local assembler data.

Attention: The index id is not necessarily the mesh item's id. Especially when having multiple meshes it will differ from the latter.

Definition at line 94 of file LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h.

    {
        auto const type_idx = std::type_index(typeid(mesh_item));
        auto const it = _builder.find(type_idx);
 
        if (it == _builder.end())
        {
            OGS_FATAL(
                "You are trying to build a local assembler for an unknown mesh "
                "element type ({:s})."
                " Maybe you have disabled this mesh element type in your build "
                "configuration, or a mesh element order does not match shape "
                "function order given in the project file.",
                type_idx.name());
        }
 
        auto const n_local_dof = _dof_table.getNumberOfElementDOF(id);
        auto const varIDs = _dof_table.getElementVariableIDs(id);
        bool const isPressureDeactivated = (varIDs.front() != 0);
        std::vector<int> involved_varIDs;  // including deactivated elements
        involved_varIDs.reserve(varIDs.size() + 1);
        if (isPressureDeactivated)
        {
            involved_varIDs.push_back(0);  // always pressure come in
        }
        involved_varIDs.insert(involved_varIDs.end(), varIDs.begin(),
                               varIDs.end());
 
        std::vector<unsigned> dofIndex_to_localIndex;
 
        // matrix and fracture assemblers with enrichments
        dofIndex_to_localIndex.resize(n_local_dof);
        std::vector<unsigned> vec_n_element_nodes;
        // TODO how to get the shape function order for each variable?
        vec_n_element_nodes.push_back(
            mesh_item.getNumberOfBaseNodes());  // pressure
        auto const max_varID = *std::max_element(varIDs.begin(), varIDs.end());
        for (int i = 1; i < max_varID + 1; i++)
        {
            vec_n_element_nodes.push_back(
                mesh_item.getNumberOfNodes());  // displacements
        }
 
        unsigned local_id = 0;
        unsigned dof_id = 0;
        for (unsigned i = 0; i < involved_varIDs.size(); i++)
        {
            auto const var_id = involved_varIDs[i];
            auto const n_var_comp =
                _dof_table.getNumberOfVariableComponents(var_id);
            auto const n_var_element_nodes = vec_n_element_nodes[i];
            for (int var_comp_id = 0; var_comp_id < n_var_comp; var_comp_id++)
            {
                auto const& ms = _dof_table.getMeshSubset(var_id, var_comp_id);
                auto const mesh_id = ms.getMeshID();
                for (unsigned k = 0; k < n_var_element_nodes; k++)
                {
                    MeshLib::Location l(mesh_id,
                                        MeshLib::MeshItemType::Node,
                                        getNodeIndex(mesh_item, k));
                    auto global_index =
                        _dof_table.getGlobalIndex(l, var_id, var_comp_id);
                    if (global_index != NumLib::MeshComponentMap::nop &&
                        dof_id < n_local_dof)
                    {
                        dofIndex_to_localIndex[dof_id++] = local_id;
                    }
                    local_id++;
                }
            }
        }
 
        return it->second(mesh_item, involved_varIDs.size(), n_local_dof,
                          dofIndex_to_localIndex,
                          std::forward<ConstructorArgs>(args)...);
    }

References _builder, _dof_table, MeshLib::Element::getNumberOfBaseNodes(), MeshLib::Element::getNumberOfNodes(), MeshLib::Node, NumLib::MeshComponentMap::nop, and OGS_FATAL.

Member Data Documentation

◆ _builder

template<typename LocalAssemblerInterface, template< typename, typename, int > class LocalAssemblerDataMatrix, template< typename, typename, int > class LocalAssemblerDataMatrixNearFracture, template< typename, typename, int > class LocalAssemblerDataFracture, int DisplacementDim, typename... ConstructorArgs>

std::unordered_map<std::type_index, LADataBuilder> ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >::_builder

private

Mapping of element types to local assembler constructors.

Definition at line 231 of file LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h.

Referenced by LocalDataInitializer(), and operator()().

◆ _dof_table

template<typename LocalAssemblerInterface, template< typename, typename, int > class LocalAssemblerDataMatrix, template< typename, typename, int > class LocalAssemblerDataMatrixNearFracture, template< typename, typename, int > class LocalAssemblerDataFracture, int DisplacementDim, typename... ConstructorArgs>

NumLib::LocalToGlobalIndexMap const& ProcessLib::LIE::HydroMechanics::LocalDataInitializer< LocalAssemblerInterface, LocalAssemblerDataMatrix, LocalAssemblerDataMatrixNearFracture, LocalAssemblerDataFracture, DisplacementDim, ConstructorArgs >::_dof_table

private

Definition at line 233 of file LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h.

Referenced by LocalDataInitializer(), and operator()().

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

LIE/HydroMechanics/LocalAssembler/LocalDataInitializer.h

Detailed Description

Classes

Public Types

Public Member Functions

Private Types

Static Private Member Functions

Private Attributes

Member Typedef Documentation

◆ LADataBuilder

◆ LADataIntfPtr

Constructor & Destructor Documentation

◆ LocalDataInitializer()

Member Function Documentation

◆ makeLocalAssemblerBuilder()

◆ operator()()

Member Data Documentation

◆ _builder

◆ _dof_table