MoFEM::DataOperator Struct Reference

base operator to do operations at Gauss Pt. level More...

#include <src/finite_elements/DataOperators.hpp>

Inheritance diagram for MoFEM::DataOperator:

Collaboration diagram for MoFEM::DataOperator:

Public Types
using	DoWorkLhsHookFunType = boost::function< MoFEMErrorCode(DataOperator *op_ptr, int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)>
using	DoWorkRhsHookFunType = boost::function< MoFEMErrorCode(DataOperator *op_ptr, int side, EntityType type, EntitiesFieldData::EntData &data)>

Public Member Functions
	DataOperator (const bool symm=true)
virtual	~DataOperator ()=default
virtual MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
	Operator for bi-linear form, usually to calculate values on left hand side. More...
virtual MoFEMErrorCode	opLhs (EntitiesFieldData &row_data, EntitiesFieldData &col_data)
virtual MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	Operator for linear form, usually to calculate values on right hand side. More...
virtual MoFEMErrorCode	opRhs (EntitiesFieldData &data, const bool error_if_no_base=false)
bool	getSymm () const
	Get if operator uses symmetry of DOFs or not. More...
void	setSymm ()
	set if operator is executed taking in account symmetry More...
void	unSetSymm ()
	unset if operator is executed for non symmetric problem More...

Public Attributes
DoWorkLhsHookFunType	doWorkLhsHook
DoWorkRhsHookFunType	doWorkRhsHook
bool	sYmm
	If true assume that matrix is symmetric structure. More...
std::array< bool, MBMAXTYPE >	doEntities
	If true operator is executed for entity. More...
bool &	doVertices
	\deprectaed If false skip vertices More...
bool &	doEdges
	\deprectaed If false skip edges More...
bool &	doQuads
	\deprectaed More...
bool &	doTris
	\deprectaed More...
bool &	doTets
	\deprectaed More...
bool &	doPrisms
	\deprectaed More...

Private Member Functions
template<bool Symm>
MoFEMErrorCode	opLhs (EntitiesFieldData &row_data, EntitiesFieldData &col_data)
template<bool ErrorIfNoBase>
MoFEMErrorCode	opRhs (EntitiesFieldData &data, const std::array< bool, MBMAXTYPE > &do_entities)

Detailed Description

base operator to do operations at Gauss Pt. level

Examples

child_and_parent.cpp, free_surface.cpp, hanging_node_approx.cpp, level_set.cpp, poisson_2d_dis_galerkin.cpp, and simple_l2_only.cpp.

Definition at line 24 of file DataOperators.hpp.

Member Typedef Documentation

DoWorkLhsHookFunType

using MoFEM::DataOperator::DoWorkLhsHookFunType = boost::function<MoFEMErrorCode( DataOperator *op_ptr, int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)>

Definition at line 33 of file DataOperators.hpp.

DoWorkRhsHookFunType

using MoFEM::DataOperator::DoWorkRhsHookFunType = boost::function<MoFEMErrorCode( DataOperator *op_ptr, int side, EntityType type, EntitiesFieldData::EntData &data)>

Definition at line 60 of file DataOperators.hpp.

Constructor & Destructor Documentation

DataOperator()

MoFEM::DataOperator::DataOperator

(

const bool

symm = true

)

This not yet implemented, switch off.

Definition at line 21 of file DataOperators.cpp.

    :
 
      sYmm(symm),
 
      doEntities{true, true, true, true, true, true,
                 true, true, true, true, true, true},
 
      doVertices(doEntities[MBVERTEX]), doEdges(doEntities[MBEDGE]),
      doQuads(doEntities[MBQUAD]), doTris(doEntities[MBTRI]),
      doTets(doEntities[MBTET]), doPrisms(doEntities[MBPRISM]) {
 
  /// This not yet implemented, switch off.
  doEntities[MBPOLYGON] = false;
  doEntities[MBPYRAMID] = false;
  doEntities[MBKNIFE] = false;
  doEntities[MBPOLYHEDRON] = false;
}

~DataOperator()

virtual MoFEM::DataOperator::~DataOperator ( )

virtualdefault

Member Function Documentation

doWork() [1/2]

virtual MoFEMErrorCode MoFEM::DataOperator::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inlinevirtual

Operator for bi-linear form, usually to calculate values on left hand side.

Reimplemented in OpCalculateMassMatrix, MyOp2, OpSimpleRodK, CallingOp, OpS, PrismOpLhs, OpFace, OpK, QuadOpLhs, AnalyticalDirichletBC::ApproxField::OpLhs, OpRowCol, OpVolume, MetaSpringBC::OpSpringALEMaterialLhs_dX_dX, MetaSpringBC::OpSpringALEMaterialLhs_dX_dx, MetaSpringBC::OpSpringALEMaterialLhs, OpAssemble, MetaSpringBC::SpringALEMaterialVolOnSideLhs, MetaSpringBC::OpSpringKs_dX, NavierStokesElement::OpAssembleLhs, NonlinearPoissonOps::OpBoundaryTangentMatrix, MetaSpringBC::OpSpringKs, ReactionDiffusionEquation::OpAssembleStiffLhs< DIM >, Poisson2DLagrangeMultiplierOperators::OpBoundaryLhsC, Poisson2DNonhomogeneousOperators::OpBoundaryLhs, ReactionDiffusionEquation::OpAssembleMass, NonlinearPoissonOps::OpDomainTangentMatrix, Poisson2DLagrangeMultiplierOperators::OpDomainLhsK, and Poisson2DNonhomogeneousOperators::OpDomainLhs.

Examples

free_surface.cpp, and MagneticElement.hpp.

Definition at line 40 of file DataOperators.hpp.

                                                                    {
    MoFEMFunctionBeginHot;
    if (doWorkLhsHook) {
      CHKERR doWorkLhsHook(this, row_side, col_side, row_type, col_type,
                           row_data, col_data);
    } else {
      SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
              "doWork function is not implemented for this operator");
    }
    MoFEMFunctionReturnHot(0);
  }

doWork() [2/2]

virtual MoFEMErrorCode MoFEM::DataOperator::doWork ( int  side, EntityType  type, EntitiesFieldData::EntData &  data  )

inlinevirtual

Operator for linear form, usually to calculate values on right hand side.

Reimplemented in MoFEM::OpInvertMatrix< DIM >, MoFEM::OpSetInvJacH1ForFlatPrism, MoFEM::OpCalculateInvJacForFlatPrism, MoFEM::OpSetInvJacH1ForFatPrism, MoFEM::OpCalculateInvJacForFatPrism, MoFEM::OpMultiplyDeterminantOfJacobianAndWeightsForFatPrisms, MoFEM::OpSetContravariantPiolaTransformOnEdge2D, MoFEM::OpSetContravariantPiolaTransformOnFace2DImpl< 3 >, MoFEM::OpSetContravariantPiolaTransformOnFace2DImpl< 2 >, MoFEM::OpSetCovariantPiolaTransformOnFace2DImpl< 2 >, MoFEM::OpMakeHdivFromHcurl, MoFEM::OpSetInvJacHcurlFaceImpl< 3 >, MoFEM::OpSetInvJacHcurlFaceImpl< 2 >, MoFEM::OpSetInvJacSpaceForFaceImpl< 2, 2 >, MoFEM::OpSetInvJacSpaceForFaceImpl< 3, 1 >, MoFEM::OpSetInvJacSpaceForFaceImpl< 2, 1 >, MoFEM::OpCalculateHVecTensorDivergence< Tensor_Dim0, Tensor_Dim1, CoordSys >, MoFEM::OpCalculateHTensorTensorField< Tensor_Dim0, Tensor_Dim1 >, MoFEM::OpCalculateHVecTensorField< Tensor_Dim0, Tensor_Dim1 >, MoFEM::OpCalculateHcurlVectorCurl< 1, 3 >, MoFEM::OpCalculateHcurlVectorCurl< 1, 2 >, MoFEM::OpCalculateHcurlVectorCurl< 3, 3 >, MoFEM::OpCalculateHdivVectorDivergenceDot< Tensor_Dim1, Tensor_Dim2 >, MoFEM::OpCalculateHVecVectorHessian< BASE_DIM, SPACE_DIM >, MoFEM::OpCalculateHVecVectorGradient< BASE_DIM, SPACE_DIM >, MoFEM::OpCalculateHdivVectorDivergence< BASE_DIM, SPACE_DIM >, MoFEM::OpCalculateHVecVectorFieldDot< 3, Field_Dim >, MoFEM::OpCalculateHVecVectorField_General< 3, Field_Dim, double, ublas::row_major, DoubleAllocator >, MoFEM::OpCalculateVectorFieldHessian< Tensor_Dim0, Tensor_Dim1 >, MoFEM::OpCalculateTensor2SymmetricFieldGradient_General< Tensor_Dim0, Tensor_Dim1, double, ublas::row_major, DoubleAllocator >, MoFEM::OpCalculateVectorFieldGradientDot< Tensor_Dim0, Tensor_Dim1 >, MoFEM::OpCalculateVectorFieldGradient_General< Tensor_Dim0, Tensor_Dim1, double, ublas::row_major, DoubleAllocator >, MoFEM::OpCalculateScalarFieldHessian< Tensor_Dim >, MoFEM::OpLoopSide< E >, MoFEM::OpCalculateScalarFieldGradient_General< Tensor_Dim, double, ublas::row_major, DoubleAllocator >, MoFEM::OpCalculateTensor2SymmetricFieldValuesDot< Tensor_Dim >, MoFEM::OpCalculateTensor2SymmetricFieldValues< Tensor_Dim >, MoFEM::OpCalculateTensor2FieldValuesDot< Tensor_Dim0, Tensor_Dim1 >, MoFEM::OpCalculateTensor2FieldValues_General< Tensor_Dim0, Tensor_Dim1, double, ublas::row_major, DoubleAllocator >, MoFEM::OpCalculateTensor2FieldValues_General< Tensor_Dim0, Tensor_Dim1, T, L, A >, MoFEM::OpPostProcMapInMoab< DIM1, DIM2 >, MoFEM::OpCalculateVectorFieldValuesFromPetscVecImpl< Tensor_Dim, CTX >, MoFEM::OpSetCovariantPiolaTransformOnEdge, MoFEM::OpCalculateDivergenceVectorFieldValues< Tensor_Dim, COORDINATE_SYSTEM >, MoFEM::OpGetHOTangentOnEdge, MoFEM::OpSetCovariantPiolaTransformOnFace, LevelSet::OpLhsSkeleton, MoFEM::OpSetContravariantPiolaTransformOnFace, LevelSet::OpRhsSkeleton, MoFEM::OpGetCoordsAndNormalsOnPrism, MoFEM::OpGetHOTangentsOnEdge< FIELD_DIM >, MoFEM::OpCopyGeomDataToE< 2 >, MoFEM::OpCalculateVectorFieldValues_General< Tensor_Dim, double, ublas::row_major, DoubleAllocator >, MoFEM::OpHOSetCovariantPiolaTransformOnFace3D, MoFEM::OpHOSetContravariantPiolaTransformOnEdge3D, MoFEM::OpCalculateVectorFieldValues_General< Tensor_Dim, T, L, A >, MoFEM::OpHOSetContravariantPiolaTransformOnFace3D, MoFEM::OpGetDataAndGradient< RANK, DIM >, MoFEM::OpGetHONormalsOnFace< FIELD_DIM >, MoFEM::OpCalculateHOJacForFaceImpl< 3 >, MoFEM::OpCalculateHOJacForFaceImpl< 2 >, OpPressure, MoFEM::OpSetCovariantPiolaTransform, OpCalculateLameStress< DIM >, OpSimpleRodPreStress, MoFEM::OpSetHOCovariantPiolaTransform, MoFEM::OpSetContravariantPiolaTransform, MoFEM::OpSetHOContravariantPiolaTransform, OpCheckValsDiffVals, OpCheckValsDiffVals, MoFEM::OpSetHOWeights, MoFEM::OpCalculateScalarFieldValuesFromPetscVecImpl< CTX >, MoFEM::OpSetInvJacHdivAndHcurl, OpVolumeAssemble, MoFEM::OpDGProjectionEvaluation, MoFEM::OpSetInvJacH1, MoFEM::OpSetHOWeightsOnEdge, MyOp2, MoFEM::OpDGProjectionCoefficients, MoFEM::OpSetHOWeightsOnFace, OpVolumeSide, MoFEM::OpSchurAssembleEnd< SCHUR_DGESV >, Example::OpCalcSurfaceAverageTemperature, MoFEM::OpDGProjectionMassMatrix, MoFEM::OpSchurAssembleEnd< SCHUR_DSYSV >, MoFEM::OpSetHOInvJacVectorBase, MoFEM::OpBaseDerivativesNext< 3 >, OpVolumeTest, OpFaceSide, MoFEM::OpCalculateScalarFieldValues, SkeletonFE, CallingOp, SkeletonFE, SkeletonFE, MoFEM::OpLoopThis< E >, OpFace, MoFEM::OpBaseDerivativesNext< 1 >, MoFEM::OpCalcNormL2Tensor2< DIM_1, DIM_2 >, AnalyticalDirichletBC::ApproxField::OpRhs< FUNEVAL >, OpVolume, PrismOpRhs, QuadOpRhs, OpFlux, MoFEM::OpBaseDerivativesSetHOInvJacobian< 2 >, MoFEM::OpSetHOInvJacToScalarBasesImpl, PrismOpCheck, OpFace, MoFEM::OpCalcNormL2Tensor1< DIM >, OpFlux, PrismOp, QuadOpCheck, MoFEM::OpCalculateScalarFieldValues_General< T, A >, MoFEM::OpCalculateScalarFieldValues_General< double, DoubleAllocator >, MoFEM::OpUnSetBc, MoFEM::OpCalculateHOCoords< FIELD_DIM >, MoFEM::OpRunParent, OpFacesRot, MoFEM::OpBaseDerivativesMass< 1 >, MoFEM::OpSetBc, OpFacesFluxes, OpVolumeSet, SkeletonFE::OpFaceSide, SkeletonFE::OpFaceSide, SkeletonFE::OpFaceSide, OpDivergence, OpVolCurl, MoFEM::OpCalcNormL2Tensor0, OpDivergence, OpVolDivergence, MoFEM::OpSchurAssembleBegin, OpFace, MoFEM::OpCalculateHOJacForVolume, OpRow, OpCheck, OpVolume, MetaSpringBC::OpGetNormalSpEle, MetaSpringBC::OpGetTangentSpEle, MetaSpringBC::OpSpringFsMaterial, MoFEM::OpScaleMatrix, MoFEM::OpSymmetrizeTensor< DIM >, MoFEM::OpTensorTimesSymmetricTensor< DIM_01, DIM_23, S >, NavierStokesElement::OpCalcVolumeFlux, NavierStokesElement::OpPostProcVorticity, NavierStokesElement::OpPostProcDrag, NavierStokesElement::OpCalcDragTraction, NavierStokesElement::OpCalcDragForce, NonlinearPoissonOps::OpBoundaryResidualVector, AtomTest::OpErrorSkel< 1 >, Poisson2DNonhomogeneousOperators::OpBoundaryRhs, Poisson2DLagrangeMultiplierOperators::OpBoundaryRhsG, Example::BoundaryOp, MetaSpringBC::OpSpringFs, ReactionDiffusionEquation::OpAssembleStiffRhs< DIM >, NonlinearPoissonOps::OpDomainResidualVector, PhotonDiffusion::OpCameraInteg, Example::OpSecond, Poisson2DLagrangeMultiplierOperators::OpDomainRhsF, Poisson2DNonhomogeneousOperators::OpDomainRhs, ReactionDiffusionEquation::OpAssembleSlowRhs, Example::OpFirst, Example::OpZero, ElasticExample::OpCalculateTraction, OpVolumeCalculation, OpAleLhsPre_dX_dx< S >, NavierStokesElement::OpAssembleRhs, OpAssemble, OpCalculateStiffnessScaledByDensityField, OpCalculateHomogeneousStiffness< S >, OpCalculateEshelbyStress, OpCalculateEnergy, OpCalculateStress< S >, HookeInternalStressElement::OpSaveStress, OpCalculateStrainAle, OpCalculateStrain< D >, HookeInternalStressElement::OpGetAnalyticalInternalStress< S >, and HookeInternalStressElement::OpGetInternalStress.

Definition at line 67 of file DataOperators.hpp.

                                                                {
    MoFEMFunctionBeginHot;
    if (doWorkRhsHook) {
      CHKERR doWorkRhsHook(this, side, type, data);
    } else {
      SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
              "doWork function is not implemented for this operator");
    }
    MoFEMFunctionReturnHot(0);
  }

getSymm()

bool MoFEM::DataOperator::getSymm ( ) const

inline

Get if operator uses symmetry of DOFs or not.

If symmetry is used, only not repeating combinations of entities are looped. For an example pair of (Vertex, Edge_0) and (Edge_0, Vertex) will calculate the same matrices only transposed. Implementing that this can be exploited by integrating only one pair.

Returns

true if symmetry

Definition at line 107 of file DataOperators.hpp.

{ return sYmm; }

opLhs() [1/2]

MoFEMErrorCode MoFEM::DataOperator::opLhs ( EntitiesFieldData &  row_data, EntitiesFieldData &  col_data  )

virtual

Definition at line 90 of file DataOperators.cpp.

                                                                {
  if (getSymm())
    return opLhs<true>(row_data, col_data);
  else
    return opLhs<false>(row_data, col_data);
}

opLhs() [2/2]

template<bool Symm>

MoFEMErrorCode MoFEM::DataOperator::opLhs ( EntitiesFieldData &  row_data, EntitiesFieldData &  col_data  )

inlineprivate

Definition at line 41 of file DataOperators.cpp.

                                                                {
  MoFEMFunctionBeginHot;
 
  auto do_col_entity =
      [&](boost::ptr_vector<EntitiesFieldData::EntData> &row_ent_data,
          const int ss, const EntityType row_type, const EntityType low_type,
          const EntityType hi_type) {
        MoFEMFunctionBeginHot;
        for (EntityType col_type = low_type; col_type != hi_type; ++col_type) {
          auto &col_ent_data = col_data.dataOnEntities[col_type];
          for (size_t SS = 0; SS != col_ent_data.size(); SS++) {
            if (col_ent_data[SS].getFieldData().size())
              CHKERR doWork(ss, SS, row_type, col_type, row_ent_data[ss],
                            col_ent_data[SS]);
          }
        }
        MoFEMFunctionReturnHot(0);
      };
 
  auto do_row_entity = [&](const EntityType type) {
    MoFEMFunctionBeginHot;
    auto &row_ent_data = row_data.dataOnEntities[type];
    for (size_t ss = 0; ss != row_ent_data.size(); ++ss) {
      if constexpr (!Symm)
        CHKERR do_col_entity(row_ent_data, ss, type, MBVERTEX, type);
      size_t SS = 0;
      if constexpr (Symm)
        SS = ss;
      for (; SS < col_data.dataOnEntities[type].size(); ++SS) {
        CHKERR doWork(ss, SS, type, type, row_ent_data[ss],
                      col_data.dataOnEntities[type][SS]);
      }
      CHKERR do_col_entity(row_ent_data, ss, type,
                           static_cast<EntityType>(type + 1), MBMAXTYPE);
    }
    MoFEMFunctionReturnHot(0);
  };
 
  for (EntityType row_type = MBVERTEX; row_type != MBMAXTYPE; ++row_type) {
    if (doEntities[row_type]) {
      CHKERR do_row_entity(row_type);
    }
  }
 
 
  MoFEMFunctionReturnHot(0);
}

opRhs() [1/2]

MoFEMErrorCode MoFEM::DataOperator::opRhs ( EntitiesFieldData &  data, const bool  error_if_no_base = false  )

virtual

Reimplemented in MoFEM::OpAddParentEntData.

Definition at line 140 of file DataOperators.cpp.

                                                                {
  if (error_if_no_base)
    return opRhs<true>(data, doEntities);
  else
    return opRhs<false>(data, doEntities);
}

opRhs() [2/2]

template<bool ErrorIfNoBase>

MoFEMErrorCode MoFEM::DataOperator::opRhs ( EntitiesFieldData &  data, const std::array< bool, MBMAXTYPE > &  do_entities  )

inlineprivate

Definition at line 100 of file DataOperators.cpp.

                                                                  {
  MoFEMFunctionBeginHot;
 
  auto do_entity = [&](auto type) {
    MoFEMFunctionBeginHot;
 
    auto &ent_data = data.dataOnEntities[type];
    const size_t size = ent_data.size();
    for (size_t ss = 0; ss != size; ++ss) {
 
      auto &side_data = ent_data[ss];
 
      if (ErrorIfNoBase) {
        if (side_data.getFieldData().size() &&
            (side_data.getBase() == NOBASE ||
             side_data.getBase() == LASTBASE)) {
          for (VectorDofs::iterator it = side_data.getFieldDofs().begin();
               it != side_data.getFieldDofs().end(); it++)
            if ((*it) && (*it)->getActive())
              SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "No base on");
        }
      }
 
      CHKERR doWork(ss, type, side_data);
    }
 
    MoFEMFunctionReturnHot(0);
  };
 
  for (EntityType row_type = MBVERTEX; row_type != MBMAXTYPE; ++row_type) {
    if (do_entities[row_type]) {
      CHKERR do_entity(row_type);
    }
  }
 
 
  MoFEMFunctionReturnHot(0);
}

setSymm()

void MoFEM::DataOperator::setSymm ( )

inline

set if operator is executed taking in account symmetry

Definition at line 110 of file DataOperators.hpp.

{ sYmm = true; }

unSetSymm()

void MoFEM::DataOperator::unSetSymm ( )

inline

unset if operator is executed for non symmetric problem

Definition at line 113 of file DataOperators.hpp.

{ sYmm = false; }

Member Data Documentation

doEdges

bool& MoFEM::DataOperator::doEdges

\deprectaed If false skip edges

Examples

ElasticityMixedFormulation.hpp.

Definition at line 91 of file DataOperators.hpp.

doEntities

std::array<bool, MBMAXTYPE> MoFEM::DataOperator::doEntities

If true operator is executed for entity.

Examples

HookeInternalStressElement.hpp.

Definition at line 85 of file DataOperators.hpp.

doPrisms

bool& MoFEM::DataOperator::doPrisms

\deprectaed

Examples

ElasticityMixedFormulation.hpp.

Definition at line 95 of file DataOperators.hpp.

doQuads

bool& MoFEM::DataOperator::doQuads

\deprectaed

Examples

ElasticityMixedFormulation.hpp.

Definition at line 92 of file DataOperators.hpp.

doTets

bool& MoFEM::DataOperator::doTets

\deprectaed

Examples

ElasticityMixedFormulation.hpp.

Definition at line 94 of file DataOperators.hpp.

doTris

bool& MoFEM::DataOperator::doTris

\deprectaed

Examples

ElasticityMixedFormulation.hpp.

Definition at line 93 of file DataOperators.hpp.

doVertices

bool& MoFEM::DataOperator::doVertices

\deprectaed If false skip vertices

Examples

ElasticityMixedFormulation.hpp.

Definition at line 90 of file DataOperators.hpp.

doWorkLhsHook

DoWorkLhsHookFunType MoFEM::DataOperator::doWorkLhsHook

Definition at line 35 of file DataOperators.hpp.

doWorkRhsHook

DoWorkRhsHookFunType MoFEM::DataOperator::doWorkRhsHook

Definition at line 62 of file DataOperators.hpp.

sYmm

bool MoFEM::DataOperator::sYmm

If true assume that matrix is symmetric structure.

Examples

UnsaturatedFlow.hpp.

Definition at line 82 of file DataOperators.hpp.

---

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

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