                                                                  {
    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

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ OpErrorSkel()

Member Function Documentation

◆ doWork()

Member Data Documentation

◆ commonDataPtr