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

#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)

Detailed Description

Definition at line 455 of file HODataOperators.hpp.

Constructor & Destructor Documentation

AddHOOps()

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

delete

Member Function Documentation

add()

MoFEMErrorCode MoFEM::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

Definition at line 1000 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));
      pipeline.push_back(new OpSetHOInvJacVectorBase(HDIV, inv_jac));
      break;
    case L2:
      pipeline.push_back(new OpSetHOInvJacToScalarBases<3>(L2, inv_jac));
      break;
    default:
      SETERRQ1(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