AtomTest::OpErrorSkel< 1 > Struct Reference

Inheritance diagram for AtomTest::OpErrorSkel< 1 >:

Collaboration diagram for AtomTest::OpErrorSkel< 1 >:

Public Member Functions
	OpErrorSkel (boost::shared_ptr< CommonData > &common_data_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Public Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

Definition at line 281 of file hanging_node_approx.cpp.

Constructor & Destructor Documentation

OpErrorSkel()

AtomTest::OpErrorSkel< 1 >::OpErrorSkel ( boost::shared_ptr< CommonData > &  common_data_ptr )

inline

Definition at line 283 of file hanging_node_approx.cpp.

      : BoundaryEleOp(H1, OPSPACE), commonDataPtr(common_data_ptr) {}

Member Function Documentation

doWork()

MoFEMErrorCode AtomTest::OpErrorSkel< 1 >::doWork ( int  side, EntityType  type, EntData &  data  )

inline

Definition at line 285 of file hanging_node_approx.cpp.

                                                                  {
    MoFEMFunctionBegin;
 
    FTensor::Index<'i', SPACE_DIM> i;
 
    const int nb_integration_pts = getGaussPts().size2();
    auto t_w = getFTensor0IntegrationWeight();
    auto t_val = getFTensor0FromVec(*(commonDataPtr->approxVals));
    auto t_coords = getFTensor1CoordsAtGaussPts();
 
    const double volume = getMeasure();
 
    double error2 = 0;
    for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
      const double alpha = t_w * volume;
      double diff = t_val - AtomTest::approxFunction(t_coords(0), t_coords(1),
                                                     t_coords(2));
      error2 += alpha * (pow(diff, 2));
 
      ++t_w;
      ++t_val;
      ++t_coords;
    }
 
    MOFEM_LOG("SELF", Sev::verbose) << "Boundary error " << sqrt(error2);
 
    constexpr double eps = 1e-8;
    if (sqrt(error2) > eps)
      SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
               "Error on boundary = %6.4e", sqrt(error2));
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> AtomTest::OpErrorSkel< 1 >::commonDataPtr

Definition at line 282 of file hanging_node_approx.cpp.

---

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

• atom_tests/hanging_node_approx.cpp