Example::OpError< 1 > Struct Reference

Inheritance diagram for Example::OpError< 1 >:

Collaboration diagram for Example::OpError< 1 >:

Public Member Functions
	OpError (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 86 of file approximaton.cpp.

Constructor & Destructor Documentation

OpError()

Example::OpError< 1 >::OpError ( boost::shared_ptr< CommonData > &  common_data_ptr )

inline

Definition at line 88 of file approximaton.cpp.

      : DomainEleOp(FIELD_NAME, OPROW), commonDataPtr(common_data_ptr) {}

Member Function Documentation

doWork()

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

inline

Definition at line 90 of file approximaton.cpp.

                                                                  {
    MoFEMFunctionBegin;
 
    if (const size_t nb_dofs = data.getIndices().size()) {
 
      const int nb_integration_pts = getGaussPts().size2();
      auto t_w = getFTensor0IntegrationWeight();
      auto t_val = getFTensor0FromVec(*(commonDataPtr->approxVals));
      auto t_coords = getFTensor1CoordsAtGaussPts();
 
      VectorDouble nf(nb_dofs, false);
      nf.clear();
 
      FTensor::Index<'i', 3> i;
      const double volume = getMeasure();
 
      auto t_row_base = data.getFTensor0N();
      double error = 0;
      for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
        const double alpha = t_w * volume;
        double diff = t_val - Example::approxFunction(t_coords(0), t_coords(1),
                                                      t_coords(2));
        error += alpha * pow(diff, 2);
 
        for (size_t r = 0; r != nb_dofs; ++r) {
          nf[r] += alpha * t_row_base * diff;
          ++t_row_base;
        }
 
        ++t_w;
        ++t_val;
        ++t_coords;
      }
 
      const int index = 0;
      CHKERR VecSetValue(commonDataPtr->L2Vec, index, error, ADD_VALUES);
      CHKERR VecSetValues(commonDataPtr->resVec, data, &nf[0], ADD_VALUES);
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> Example::OpError< 1 >::commonDataPtr

Definition at line 87 of file approximaton.cpp.

---

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

• tutorials/scl-0/approximaton.cpp