MyOp< OP > Struct Template Reference

Operator used to check consistency between local coordinates and global cooridnates for integrated points and evaluated points. More...

Inheritance diagram for MyOp< OP >:

Collaboration diagram for MyOp< OP >:

Public Member Functions
	MyOp (std::array< double, 12 > &eval_points)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	MyOp (const char type, const char face_type)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)

Public Attributes
MatrixShallowArrayAdaptor< double >	evalPoints
const char	faceType

Detailed Description

template<typename OP>
struct MyOp< OP >

Operator used to check consistency between local coordinates and global cooridnates for integrated points and evaluated points.

Examples

field_evaluator.cpp, forces_and_sources_testing_edge_element.cpp, and forces_and_sources_testing_flat_prism_element.cpp.

Definition at line 21 of file field_evaluator.cpp.

Constructor & Destructor Documentation

MyOp() [1/2]

template<typename OP >

MyOp< OP >::MyOp ( std::array< double, 12 > &  eval_points )

inline

Definition at line 25 of file field_evaluator.cpp.

      : OP("FIELD1", OP::OPROW),
        evalPoints(
            getMatrixAdaptor(eval_points.data(), eval_points.size() / 3, 3)) {}

MyOp() [2/2]

template<typename OP >

MyOp< OP >::MyOp ( const char  type, const char  face_type  )

inline

Definition at line 19 of file forces_and_sources_testing_contact_prism_element.cpp.

      : ContactPrismElementForcesAndSourcesCore::UserDataOperator(
            "FIELD1", "FIELD1", type, face_type),
        faceType(face_type) {}

Member Function Documentation

doWork() [1/3]

template<typename OP >

MoFEMErrorCode MyOp< OP >::doWork ( int  row_side, int  col_side, EntityType  row_type, EntityType  col_type, EntitiesFieldData::EntData &  row_data, EntitiesFieldData::EntData &  col_data  )

inline

Definition at line 66 of file forces_and_sources_testing_contact_prism_element.cpp.

                                                            {
    MoFEMFunctionBeginHot;
 
    if (row_data.getFieldData().empty())
      MoFEMFunctionReturnHot(0);
 
    MOFEM_LOG("ATOM_TEST", Sev::inform) << "NH1NH1";
    MOFEM_LOG("ATOM_TEST", Sev::inform)
        << "row side: " << row_side << " row_type: " << row_type;
    MOFEM_LOG("ATOM_TEST", Sev::inform) << row_data;
    MOFEM_LOG("ATOM_TEST", Sev::inform) << "NH1NH1";
    MOFEM_LOG("ATOM_TEST", Sev::inform)
        << "col side: " << col_side << " col_type: " << col_type;
    MOFEM_LOG("ATOM_TEST", Sev::inform) << col_data;
 
    MoFEMFunctionReturnHot(0);
  }

doWork() [2/3]

template<typename OP >

MoFEMErrorCode MyOp< OP >::doWork ( int  side, EntityType  type, EntitiesFieldData::EntData &  data  )

inline

Definition at line 24 of file forces_and_sources_testing_contact_prism_element.cpp.

                                                        {
    MoFEMFunctionBeginHot;
 
    if (data.getFieldData().empty())
      MoFEMFunctionReturnHot(0);
 
    MOFEM_LOG("ATOM_TEST", Sev::inform) << "NH1";
    MOFEM_LOG("ATOM_TEST", Sev::inform)
        << "side: " << side << " type: " << type;
    MOFEM_LOG("ATOM_TEST", Sev::inform) << data;
 
    if (faceType == FACEMASTER) {
      MOFEM_LOG("ATOM_TEST", Sev::inform)
          << "coords Master " << std::fixed << std::setprecision(2)
          << getCoordsMaster();
      MOFEM_LOG("ATOM_TEST", Sev::inform)
          << "area Master " << std::fixed << std::setprecision(2)
          << getAreaMaster();
      MOFEM_LOG("ATOM_TEST", Sev::inform)
          << "normal Master " << std::fixed << std::setprecision(2)
          << getNormalMaster();
      MOFEM_LOG("ATOM_TEST", Sev::inform)
          << "coords at Gauss Pts Master " << std::fixed
          << std::setprecision(2) << getCoordsAtGaussPtsMaster();
    } else {
      MOFEM_LOG("ATOM_TEST", Sev::inform)
          << "coords Slave " << std::fixed << std::setprecision(2)
          << getCoordsSlave();
      MOFEM_LOG("ATOM_TEST", Sev::inform)
          << "area Slave " << std::fixed << std::setprecision(2)
          << getAreaSlave();
      MOFEM_LOG("ATOM_TEST", Sev::inform)
          << "normal Slave " << std::fixed << std::setprecision(2)
          << getNormalSlave();
      MOFEM_LOG("ATOM_TEST", Sev::inform)
          << "coords at Gauss Pts Slave " << std::fixed
          << std::setprecision(2) << getCoordsAtGaussPtsSlave();
    }
    MoFEMFunctionReturnHot(0);
  }

doWork() [3/3]

template<typename OP >

MoFEMErrorCode MyOp< OP >::doWork ( int  side, EntityType  type, EntitiesFieldData::EntData &  data  )

inline

Definition at line 30 of file field_evaluator.cpp.

                                                        {
    MoFEMFunctionBegin;
    if (type == MBVERTEX) {
 
      MOFEM_LOG("SELF", Sev::inform) << "FE " << OP::getFEEntityHandle();
 
      MOFEM_LOG("SELF", Sev::inform) << "Integration pts" << std::endl;
      MOFEM_LOG("SELF", Sev::inform) << OP::getGaussPts() << endl;
 
      MOFEM_LOG("SELF", Sev::inform) << "Global coordinates " << endl;
      MOFEM_LOG("SELF", Sev::inform) << OP::getCoordsAtGaussPts() << std::endl;
 
      for (int gg = 0; gg != OP::getCoordsAtGaussPts().size1(); ++gg) {
        int pt_number = OP::getGaussPts()(OP::getGaussPts().size1() - 1, gg);
 
        MOFEM_LOG("SELF", Sev::inform) << "gg " << gg << std::endl;
        MOFEM_LOG("SELF", Sev::inform) << "pt " << pt_number << std::endl;
 
        ublas::matrix_row<MatrixDouble> coord_at_gauss_pt(
            OP::getCoordsAtGaussPts(), gg);
        ublas::matrix_row<MatrixShallowArrayAdaptor<double>> eval_coord(
            evalPoints, pt_number);
 
        MOFEM_LOG("SELF", Sev::inform) << "coord_at_gauss_pt ";
        MOFEM_LOG("SELF", Sev::inform) << coord_at_gauss_pt << std::endl;
 
        MOFEM_LOG("SELF", Sev::inform) << "eval_coord ";
        MOFEM_LOG("SELF", Sev::inform) << eval_coord << std::endl;
 
        double error = norm_2(coord_at_gauss_pt - eval_coord);
 
        if (error > 1e-12)
          SETERRQ2(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
                   "Difference at %d error = %3.4e", pt_number, error);
      }
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

evalPoints

template<typename OP >

MatrixShallowArrayAdaptor<double> MyOp< OP >::evalPoints

Definition at line 23 of file field_evaluator.cpp.

faceType

template<typename OP >

const char MyOp< OP >::faceType

Definition at line 18 of file forces_and_sources_testing_contact_prism_element.cpp.

---

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

• atom_tests/field_evaluator.cpp
• atom_tests/forces_and_sources_testing_contact_prism_element.cpp