v0.14.0
Loading...
Searching...
No Matches
Public Member Functions | Static Public Member Functions | List of all members
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 545 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 896 of file HODataOperators.cpp.

900 {
902
903 if (!jac)
904 jac = boost::make_shared<MatrixDouble>();
905 if (!det)
906 det = boost::make_shared<VectorDouble>();
907 if (!inv_jac)
908 inv_jac = boost::make_shared<MatrixDouble>();
909
910 if (geom_field_name.empty()) {
911
912 pipeline.push_back(new OpCalculateHOJac<3>(jac));
913
914 } else {
915
916 pipeline.push_back(new OpCalculateHOCoords<3>(geom_field_name));
917 pipeline.push_back(
918 new OpCalculateVectorFieldGradient<3, 3>(geom_field_name, jac));
919 }
920
921 pipeline.push_back(new OpInvertMatrix<3>(jac, det, inv_jac));
922 pipeline.push_back(new OpSetHOWeights(det));
923
924 for (auto s : spaces) {
925 switch (s) {
926 case NOSPACE:
927 break;
928 case H1:
929 pipeline.push_back(new OpSetHOInvJacToScalarBases<3>(H1, inv_jac));
930 break;
931 case HCURL:
932 pipeline.push_back(new OpSetHOCovariantPiolaTransform(HCURL, inv_jac));
933 pipeline.push_back(new OpSetHOInvJacVectorBase(HCURL, inv_jac));
934 break;
935 case HDIV:
936 pipeline.push_back(
937 new OpSetHOContravariantPiolaTransform(HDIV, det, jac));
938 pipeline.push_back(new OpSetHOInvJacVectorBase(HDIV, inv_jac));
939 break;
940 case L2:
941 pipeline.push_back(new OpSetHOInvJacToScalarBases<3>(L2, inv_jac));
942 break;
943 default:
944 SETERRQ1(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
945 "Space %s not yet implemented", FieldSpaceNames[s]);
946 break;
947 }
948 }
949
951}
@ L2
field with C-1 continuity
Definition: definitions.h:88
@ H1
continuous field
Definition: definitions.h:85
@ NOSPACE
Definition: definitions.h:83
@ HCURL
field with continuous tangents
Definition: definitions.h:86
@ HDIV
field with continuous normal traction
Definition: definitions.h:87
static const char *const FieldSpaceNames[]
Definition: definitions.h:92
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:346
@ MOFEM_NOT_IMPLEMENTED
Definition: definitions.h:32
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:416

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