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

Specialization for 2D face elements in 2D problems embedded in 3D 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 embedded in 3D.

Detailed Description

Specialization for 2D face elements in 2D problems embedded in 3D space.

This specialization adds high-order operators for 2D face elements operating in 2D problems but embedded within 3D coordinate space. Typically used for surface elements or membranes in 3D environments.

Definition at line 640 of file HODataOperators.hpp.

Constructor & Destructor Documentation

AddHOOps()

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

delete

Member Function Documentation

add()

MoFEMErrorCode AddHOOps< 2, 2, 3 >::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 embedded in 3D.

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 831 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()) {
 
  } else {
 
    pipeline.push_back(new OpCalculateHOCoords<3>(geom_field_name));
    pipeline.push_back(new OpGetHONormalsOnFace<3>(geom_field_name));
  }
 
  pipeline.push_back(new OpCalculateHOJacForFaceEmbeddedIn3DSpace(jac_ptr));
  pipeline.push_back(new OpInvertMatrix<3>(jac_ptr, det_ptr, inv_jac_ptr));
  pipeline.push_back(new OpSetHOWeightsOnFace());
 
  for (auto s : spaces) {
    switch (s) {
    case NOSPACE:
      break;
    case H1:
      pipeline.push_back(
          new OpSetInvJacH1ForFaceEmbeddedIn3DSpace(inv_jac_ptr));
      break;
    case HDIV:
      pipeline.push_back(new OpMakeHdivFromHcurl());
      pipeline.push_back(
          new OpSetContravariantPiolaTransformOnFace2DEmbeddedIn3DSpace(
              jac_ptr));
      pipeline.push_back(
          new OpSetInvJacHcurlFaceEmbeddedIn3DSpace(inv_jac_ptr));
      break;
    case L2:
      pipeline.push_back(
          new OpSetInvJacL2ForFaceEmbeddedIn3DSpace(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