v0.14.0
Public Member Functions | Static Public Member Functions | List of all members
MoFEM::AddHOOps< 2, 2, 2 > 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 > 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)
 

Detailed Description

Definition at line 417 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

Definition at line 801 of file HODataOperators.cpp.

806  {
808 
809  if (!jac_ptr)
810  jac_ptr = boost::make_shared<MatrixDouble>();
811  if (!det_ptr)
812  det_ptr = boost::make_shared<VectorDouble>();
813  if (!inv_jac_ptr)
814  inv_jac_ptr = boost::make_shared<MatrixDouble>();
815 
816  if (geom_field_name.empty()) {
817 
818  pipeline.push_back(new OpCalculateHOJac<2>(jac_ptr));
819 
820  } else {
821 
822  pipeline.push_back(new OpCalculateHOCoords<2>(geom_field_name));
823  pipeline.push_back(
824  new OpCalculateVectorFieldGradient<2, 2>(geom_field_name, jac_ptr));
825  pipeline.push_back(new OpGetHONormalsOnFace<2>(geom_field_name));
826  }
827 
828  pipeline.push_back(new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
829  pipeline.push_back(new OpSetHOWeightsOnFace());
830 
831  for (auto s : spaces) {
832  switch (s) {
833  case NOSPACE:
834  break;
835  case H1:
836  case L2:
837  pipeline.push_back(new OpSetHOInvJacToScalarBases<2>(s, inv_jac_ptr));
838  break;
839  case HCURL:
840  pipeline.push_back(new OpSetCovariantPiolaTransformOnFace2D(inv_jac_ptr));
841  pipeline.push_back(new OpSetInvJacHcurlFace(inv_jac_ptr));
842  break;
843  case HDIV:
844  pipeline.push_back(new OpMakeHdivFromHcurl());
845  pipeline.push_back(new OpSetContravariantPiolaTransformOnFace2D(jac_ptr));
846  pipeline.push_back(new OpSetInvJacHcurlFace(inv_jac_ptr));
847  break;
848  default:
849  SETERRQ1(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
850  "Space %s not yet implemented", FieldSpaceNames[s]);
851  }
852  }
853 
855 }

The documentation for this struct was generated from the following files:
NOSPACE
@ NOSPACE
Definition: definitions.h:83
MoFEM::OpSetInvJacHcurlFace
OpSetInvJacHcurlFaceImpl< 2 > OpSetInvJacHcurlFace
Definition: UserDataOperators.hpp:2978
H1
@ H1
continuous field
Definition: definitions.h:85
MoFEM::OpSetContravariantPiolaTransformOnFace2D
OpSetContravariantPiolaTransformOnFace2DImpl< 2 > OpSetContravariantPiolaTransformOnFace2D
Definition: UserDataOperators.hpp:3076
L2
@ L2
field with C-1 continuity
Definition: definitions.h:88
MoFEM::OpSetCovariantPiolaTransformOnFace2D
OpSetCovariantPiolaTransformOnFace2DImpl< 2 > OpSetCovariantPiolaTransformOnFace2D
Definition: UserDataOperators.hpp:3028
FieldSpaceNames
const static char *const FieldSpaceNames[]
Definition: definitions.h:92
HCURL
@ HCURL
field with continuous tangents
Definition: definitions.h:86
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
HDIV
@ HDIV
field with continuous normal traction
Definition: definitions.h:87
MOFEM_NOT_IMPLEMENTED
@ MOFEM_NOT_IMPLEMENTED
Definition: definitions.h:32
OpCalculateVectorFieldGradient
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359