createLibmeshElements.h

Go to the documentation of this file.

#ifndef included_AMP_createLibmeshElements
#define included_AMP_createLibmeshElements
 
// AMP headers
#include "AMP/AMP_TPLs.h"
#include "AMP/mesh/MeshIterator.h"
 
 
#ifdef AMP_USE_LIBMESH
 
 
// Libmesh headers
DISABLE_WARNINGS
    #include "libmesh/libmesh_config.h"
    #undef LIBMESH_ENABLE_REFERENCE_COUNTING
    #include "libmesh/elem.h"
    #include "libmesh/enum_quadrature_type.h"
    #include "libmesh/fe_base.h"
    #include "libmesh/fe_type.h"
ENABLE_WARNINGS
 
 
namespace AMP::Discretization {
 
 
class createLibmeshElements
{
public:
    createLibmeshElements();
 
    virtual ~createLibmeshElements();
 
    void reinit( const AMP::Mesh::MeshIterator &iterator );
 
    void reinit( const AMP::Mesh::MeshIterator &iterator,
                 libMeshEnums::QuadratureType qtype,
                 libMeshEnums::Order qorder,
                 std::shared_ptr<const libMesh::FEType> type,
                 bool cache_fe = false );
 
    const libMesh::Elem *getElement( const AMP::Mesh::MeshElementID &id ) const;
 
    const libMesh::FEBase *getFEBase( const AMP::Mesh::MeshElementID &id ) const;
 
    const libMesh::QBase *getQBase( const AMP::Mesh::MeshElementID &id ) const;
 
    const libMesh::FEType *getFEType() const { return d_type.get(); }
 
    static libMesh::Elem *createElement( const AMP::Mesh::MeshElement &elem );
 
private:
    std::vector<AMP::Mesh::MeshElementID> d_ids;
    std::vector<size_t> d_index;
    libMeshEnums::QuadratureType d_qtype;
    libMeshEnums::Order d_qorder;
    std::shared_ptr<const libMesh::FEType> d_type;
    mutable AMP::Mesh::MeshElementID d_last_id;
    std::shared_ptr<libMesh::FEBase> d_base;
    std::shared_ptr<libMesh::QBase> d_rule;
    std::vector<libMesh::Elem *> d_elements;
    std::vector<libMesh::FEBase *> d_base_element;
    std::vector<libMesh::QBase *> d_rule_element;
};
} // namespace AMP::Discretization
 
#endif
#endif