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

Specialization for 3D volume elements in 3D problems with 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 > space, 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 3D volume elements.

Detailed Description

Specialization for 3D volume elements in 3D problems with 3D space.

This specialization adds high-order operators for 3D volume elements operating in 3D problems within 3D coordinate space. This is the most common case for solid mechanics, heat transfer, and other 3D finite element applications involving volumetric elements like tetrahedra, hexahedra, and prisms.

Definition at line 793 of file HODataOperators.hpp.

Constructor & Destructor Documentation

AddHOOps()

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

delete

Member Function Documentation

add()

MoFEMErrorCode AddHOOps< 3, 3, 3 >::add ( boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator > &  pipeline, std::vector< FieldSpace >  space, 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 3D volume elements.

Parameters

pipeline	Pipeline to add operators to
space	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 1067 of file HODataOperators.cpp.

                                           {
  MoFEMFunctionBegin;
 
  if (!jac)
    jac = boost::make_shared<MatrixDouble>();
  if (!det)
    det = boost::make_shared<VectorDouble>();
  if (!inv_jac)
    inv_jac = boost::make_shared<MatrixDouble>();
 
  if (geom_field_name.empty()) {
 
    pipeline.push_back(new OpCalculateHOJac<3>(jac));
 
  } else {
 
    pipeline.push_back(new OpCalculateHOCoords<3>(geom_field_name));
    pipeline.push_back(
        new OpCalculateVectorFieldGradient<3, 3>(geom_field_name, jac));
  }
 
  pipeline.push_back(new OpInvertMatrix<3>(jac, det, inv_jac));
  pipeline.push_back(new OpSetHOWeights(det));
 
  for (auto s : spaces) {
    switch (s) {
    case NOSPACE:
      break;
    case H1:
      pipeline.push_back(new OpSetHOInvJacToScalarBases<3>(H1, inv_jac));
      break;
    case HCURL:
      pipeline.push_back(new OpSetHOCovariantPiolaTransform(HCURL, inv_jac));
      pipeline.push_back(new OpSetHOInvJacVectorBase(HCURL, inv_jac));
      break;
    case HDIV:
      pipeline.push_back(
          new OpSetHOContravariantPiolaTransform(HDIV, det, jac));
      break;
    case L2:
      pipeline.push_back(new OpSetHOInvJacToScalarBases<3>(L2, inv_jac));
      break;
    default:
      SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
               "Space %s not yet implemented", FieldSpaceNames[s]);
      break;
    }
  }
 
  MoFEMFunctionReturn(0);
}

---

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

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