OpCheckGaussCoords< ELE_OP, PARENT_ELE > Struct Template Reference

Inheritance diagram for OpCheckGaussCoords< ELE_OP, PARENT_ELE >:

Collaboration diagram for OpCheckGaussCoords< ELE_OP, PARENT_ELE >:

Public Member Functions
	OpCheckGaussCoords ()
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)

Detailed Description

template<typename ELE_OP, typename PARENT_ELE>
struct OpCheckGaussCoords< ELE_OP, PARENT_ELE >

Examples

child_and_parent.cpp.

Definition at line 137 of file child_and_parent.cpp.

Constructor & Destructor Documentation

OpCheckGaussCoords()

template<typename ELE_OP , typename PARENT_ELE >

OpCheckGaussCoords< ELE_OP, PARENT_ELE >::OpCheckGaussCoords ( )

inline

Definition at line 138 of file child_and_parent.cpp.

: ELE_OP(NOSPACE, ELE_OP::OPSPACE) {}

Member Function Documentation

doWork()

template<typename ELE_OP , typename PARENT_ELE >

MoFEMErrorCode OpCheckGaussCoords< ELE_OP, PARENT_ELE >::doWork ( int  side, EntityType  type, EntitiesFieldData::EntData &  data  )

inline

Definition at line 140 of file child_and_parent.cpp.

                                                        {
    MoFEMFunctionBegin;
 
    MatrixDouble parent_coords;
 
    PARENT_ELE parent_fe(this->getPtrFE()->mField);
    auto op = new ELE_OP(NOSPACE, ELE_OP::OPSPACE);
    op->doWorkRhsHook = [&](DataOperator *op_ptr, int side,
                                   EntityType type,
                                   EntitiesFieldData::EntData &data) {
      MoFEMFunctionBegin;
 
      MOFEM_LOG("SELF", Sev::noisy)
          << "parent_coords in op "
          << static_cast<ELE_OP *>(op_ptr)->getCoordsAtGaussPts();
 
      parent_coords = static_cast<ELE_OP *>(op_ptr)->getCoordsAtGaussPts();
      MoFEMFunctionReturn(0);
    };
    parent_fe.getOpPtrVector().push_back(op);
 
    MOFEM_LOG("SELF", Sev::noisy) << "fe name " << this->getFEName();
    CHKERR this->loopParent(this->getFEName(), &parent_fe);
    MOFEM_LOG("SELF", Sev::noisy) << "parent_coords " << parent_coords;
 
    MatrixDouble child_coords = this->getCoordsAtGaussPts();
    MOFEM_LOG("SELF", Sev::noisy) << "child_coords " << child_coords;
 
    child_coords -= parent_coords;
 
    MOFEM_LOG("SELF", Sev::noisy) << "Corrds diffs" << child_coords;
 
    double n = 0;
    for (auto d : child_coords.data())
      n += d * d;
 
    if (sqrt(n) > 1e-12)
      SETERRQ1(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
               "Parent and child global coords at integration points are "
               "diffrent norm = %3.2e",
               sqrt(n));
 
    MoFEMFunctionReturn(0);
  }

---

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

• atom_tests/child_and_parent.cpp