Example::OpRhs Struct Reference

Inheritance diagram for Example::OpRhs:

Collaboration diagram for Example::OpRhs:

Public Member Functions
	OpRhs (boost::shared_ptr< MatrixDouble > u_grad_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)

Protected Attributes
boost::shared_ptr< MatrixDouble >	uGradPtr

Detailed Description

Use problem specific implementation for second stage of heat methid

Definition at line 69 of file heat_method.cpp.

Constructor & Destructor Documentation

OpRhs()

Example::OpRhs::OpRhs

(

boost::shared_ptr< MatrixDouble >

u_grad_ptr

)

Examples

heat_method.cpp.

Definition at line 422 of file heat_method.cpp.

    : DomainEleOp("U", DomainEleOp::OPROW), uGradPtr(u_grad_ptr) {}

Member Function Documentation

doWork()

MoFEMErrorCode Example::OpRhs::doWork	(	int	side,
		EntityType	type,
		EntitiesFieldData::EntData &	data
	)

Examples

heat_method.cpp.

Definition at line 425 of file heat_method.cpp.

                                                                      {
  MoFEMFunctionBegin;
 
  FTensor::Index<'i', 3> i;
 
  auto nb_dofs = data.getIndices().size();
  if (nb_dofs) {
 
    auto t_grad = getFTensor1FromMat<3>(*uGradPtr);
 
    auto nb_base_functions = data.getN().size2();
    auto nb_gauss_pts = getGaussPts().size2();
    std::array<double, MAX_DOFS_ON_ENTITY> nf;
    std::fill(nf.begin(), &nf[nb_dofs], 0);
 
    auto t_diff_base = data.getFTensor1DiffN<3>();
    auto t_w = getFTensor0IntegrationWeight();
    auto a = getMeasure();
 
    for (size_t gg = 0; gg != nb_gauss_pts; ++gg) {
      double alpha = t_w * a;
 
      const auto l2 = t_grad(i) * t_grad(i);
      FTensor::Tensor1<double, 3> t_one;
      if (l2 > std::numeric_limits<double>::epsilon())
        t_one(i) = t_grad(i) / std::sqrt(l2);
      else
        t_one(i) = 0;
 
      size_t bb = 0;
      for (; bb != nb_dofs; ++bb) {
        nf[bb] -= alpha * t_diff_base(i) * t_one(i);
        ++t_diff_base;
      }
 
      for (; bb < nb_base_functions; ++bb) {
        ++t_diff_base;
      }
 
      ++t_grad;
    }
 
    CHKERR VecSetValues<MoFEM::EssentialBcStorage>(getKSPf(), data, &nf[0],
                                                   ADD_VALUES);
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

uGradPtr

boost::shared_ptr<MatrixDouble> Example::OpRhs::uGradPtr

protected

Examples

heat_method.cpp.

Definition at line 77 of file heat_method.cpp.

---

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

• tutorials/scl-9/heat_method.cpp