MoFEM::AddHOOps< 2, 2, 2 > Struct Reference

Specialization for 2D face elements in 2D problems with 2D space. More...

#include "src/finite_elements/HODataOperators.hpp"

Public Member Functions
	AddHOOps ()=delete

Static Public Member Functions
static MoFEMErrorCode	add (boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &pipeline, std::vector< FieldSpace > spaces, std::string geom_field_name="", boost::shared_ptr< MatrixDouble > jac=nullptr, boost::shared_ptr< VectorDouble > det=nullptr, boost::shared_ptr< MatrixDouble > inv_jac=nullptr)
	Add high-order operators to the pipeline for 2D face elements.

Detailed Description

Specialization for 2D face elements in 2D problems with 2D space.

This specialization adds high-order operators for 2D face elements operating in 2D problems within 2D coordinate space. Commonly used for planar face elements in 2D finite element analysis.

Definition at line 608 of file HODataOperators.hpp.

Constructor & Destructor Documentation

AddHOOps()

MoFEM::AddHOOps< 2, 2, 2 >::AddHOOps ( )

delete

Member Function Documentation

add()

MoFEMErrorCode MoFEM::AddHOOps< 2, 2, 2 >::add ( boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &  pipeline, std::vector< FieldSpace >  spaces, std::string  geom_field_name = "", boost::shared_ptr< MatrixDouble >  jac = nullptr, boost::shared_ptr< VectorDouble >  det = nullptr, boost::shared_ptr< MatrixDouble >  inv_jac = nullptr  )

static

Add high-order operators to the pipeline for 2D face elements.

Parameters

pipeline	Pipeline to add operators to
spaces	Vector of field spaces to be handled
geom_field_name	Name of the geometry field (default: empty string)
jac	Shared pointer to Jacobian matrix. If provided, values will be set by this function. If null, Jacobian will be handled internally.
det	Shared pointer to determinant vector. If provided, values will be set by this function. If null, determinant will be handled internally.
inv_jac	Shared pointer to inverse Jacobian matrix. If provided, values will be set by this function. If null, inverse Jacobian will be handled internally.

Returns

MoFEMErrorCode Error code indicating success or failure

Definition at line 775 of file HODataOperators.cpp.

                                               {
  MoFEMFunctionBegin;
 
  if (!jac_ptr)
    jac_ptr = boost::make_shared<MatrixDouble>();
  if (!det_ptr)
    det_ptr = boost::make_shared<VectorDouble>();
  if (!inv_jac_ptr)
    inv_jac_ptr = boost::make_shared<MatrixDouble>();
 
  if (geom_field_name.empty()) {
 
    pipeline.push_back(new OpCalculateHOJac<2>(jac_ptr));
 
  } else {
 
    pipeline.push_back(new OpCalculateHOCoords<2>(geom_field_name));
    pipeline.push_back(
        new OpCalculateVectorFieldGradient<2, 2>(geom_field_name, jac_ptr));
    pipeline.push_back(new OpGetHONormalsOnFace<2>(geom_field_name));
  }
 
  pipeline.push_back(new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
  pipeline.push_back(new OpSetHOWeightsOnFace());
 
  for (auto s : spaces) {
    switch (s) {
    case NOSPACE:
      break;
    case H1:
    case L2:
      pipeline.push_back(new OpSetHOInvJacToScalarBases<2>(s, inv_jac_ptr));
      break;
    case HCURL:
      pipeline.push_back(new OpSetCovariantPiolaTransformOnFace2D(inv_jac_ptr));
      pipeline.push_back(new OpSetInvJacHcurlFace(inv_jac_ptr));
      break;
    case HDIV:
      pipeline.push_back(new OpMakeHdivFromHcurl());
      pipeline.push_back(new OpSetContravariantPiolaTransformOnFace2D(jac_ptr));
      pipeline.push_back(new OpSetInvJacHcurlFace(inv_jac_ptr));
      break;
    default:
      SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
               "Space %s not yet implemented", FieldSpaceNames[s]);
    }
  }
 
  MoFEMFunctionReturn(0);
}

---

The documentation for this struct was generated from the following files:

• src/finite_elements/HODataOperators.hpp
• src/finite_elements/impl/HODataOperators.cpp