AtomTest::OpError< 1 > Struct Reference

Inheritance diagram for AtomTest::OpError< 1 >:

Collaboration diagram for AtomTest::OpError< 1 >:

Public Member Functions
	OpError (boost::shared_ptr< CommonData > &common_data_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
	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 88 of file child_and_parent.cpp.

Constructor & Destructor Documentation

OpError() [1/2]

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

inline

Definition at line 90 of file child_and_parent.cpp.

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

OpError() [2/2]

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

inline

Definition at line 226 of file hanging_node_approx.cpp.

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

Member Function Documentation

doWork() [1/2]

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

inline

Definition at line 92 of file child_and_parent.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 - AtomTest::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);
  }

doWork() [2/2]

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

inline

Definition at line 228 of file hanging_node_approx.cpp.

                                                                  {
    MoFEMFunctionBegin;
 
    if (const size_t nb_dofs = data.getIndices().size()) {
 
      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_grad_val =
          getFTensor1FromMat<SPACE_DIM>(*(commonDataPtr->divApproxVals));
      auto t_coords = getFTensor1CoordsAtGaussPts();
 
      VectorDouble nf(nb_dofs, false);
      nf.clear();
 
      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 - AtomTest::approxFunction(t_coords(0), t_coords(1),
                                                       t_coords(2));
 
        auto t_grad_diff =
            AtomTest::divApproxFunction(t_coords(0), t_coords(1), t_coords(2));
        t_grad_diff(i) -= t_grad_val(i);
 
        error += alpha * (pow(diff, 2) + t_grad_diff(i) * t_grad_diff(i));
 
        for (size_t r = 0; r != nb_dofs; ++r) {
          nf[r] += alpha * t_row_base * diff;
          ++t_row_base;
        }
 
        ++t_w;
        ++t_val;
        ++t_grad_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 > AtomTest::OpError< 1 >::commonDataPtr

Definition at line 89 of file child_and_parent.cpp.

---

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

• atom_tests/child_and_parent.cpp
• atom_tests/hanging_node_approx.cpp