StandardPoisson::OpError Struct Reference

Inheritance diagram for StandardPoisson::OpError:

Collaboration diagram for StandardPoisson::OpError:

Public Member Functions
	OpError (std::string domain_field, boost::shared_ptr< CommonData > &common_data_ptr, MoFEM::Interface &m_field)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)
	[OpError] More...

Public Attributes
std::string	domainField
boost::shared_ptr< CommonData >	commonDataPtr
MoFEM::Interface &	mField

Detailed Description

Definition at line 84 of file standard_poisson.cpp.

Constructor & Destructor Documentation

OpError()

StandardPoisson::OpError::OpError ( std::string  domain_field, boost::shared_ptr< CommonData > &  common_data_ptr, MoFEM::Interface &  m_field  )

inline

Definition at line 88 of file standard_poisson.cpp.

        : DomainEleOp(domain_field, OPROW), domainField(domain_field),
          commonDataPtr(common_data_ptr), mField(m_field) {
      std::fill(&doEntities[MBVERTEX], &doEntities[MBMAXTYPE], false);
      doEntities[MBTRI] = doEntities[MBQUAD] = true;
    }

Member Function Documentation

doWork()

MoFEMErrorCode StandardPoisson::OpError::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

[OpError]

Definition at line 366 of file standard_poisson.cpp.

                                                               {
  MoFEMFunctionBegin;
  const int nb_integration_pts = getGaussPts().size2();
  const double area = getMeasure();
  auto t_w = getFTensor0IntegrationWeight();
  auto t_val = getFTensor0FromVec(*(commonDataPtr->approxVals));
  auto t_val_grad = getFTensor1FromMat<2>(*(commonDataPtr->approxValsGrad));
  auto t_coords = getFTensor1CoordsAtGaussPts();
  FTensor::Tensor1<double, 2> t_diff;
  FTensor::Index<'i', 2> i;
 
  double error_l2 = 0;
  double error_h1 = 0;
 
  for (int gg = 0; gg != nb_integration_pts; ++gg) {
    const double alpha = t_w * area;
 
    double diff = t_val - StandardPoisson::analyticalFunction(
                              t_coords(0), t_coords(1), t_coords(2));
    error_l2 += alpha * sqr(diff);
 
    VectorDouble vec = StandardPoisson::analyticalFunctionGrad(
        t_coords(0), t_coords(1), t_coords(2));
    auto t_fun_grad = getFTensor1FromArray<2, 2>(vec);
    t_diff(i) = t_val_grad(i) - t_fun_grad(i);
 
    error_h1 += alpha * t_diff(i) * t_diff(i);
 
    ++t_w;
    ++t_val;
    ++t_val_grad;
    ++t_coords;
  }
 
  int index = CommonData::ERROR_L2_NORM;
  constexpr std::array<int, 2> indices = {CommonData::ERROR_L2_NORM,
                                          CommonData::ERROR_H1_SEMINORM};
  std::array<double, 2> values;
  values[0] = error_l2;
  values[1] = error_h1;
  CHKERR VecSetValues(commonDataPtr->petscVec, 2, indices.data(), values.data(),
                      ADD_VALUES);
  MoFEMFunctionReturn(0);
}

Member Data Documentation

commonDataPtr

boost::shared_ptr<CommonData> StandardPoisson::OpError::commonDataPtr

Definition at line 86 of file standard_poisson.cpp.

domainField

std::string StandardPoisson::OpError::domainField

Definition at line 85 of file standard_poisson.cpp.

mField

MoFEM::Interface& StandardPoisson::OpError::mField

Definition at line 87 of file standard_poisson.cpp.

---

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

• users_modules/tutorials/mix-0/standard_poisson.cpp