MoFEM::FieldEvaluatorInterface::SetPts Struct Reference

Default evaluator for setting integration points. More...

#include <src/interfaces/FieldEvaluator.hpp>

Public Member Functions
	SetPts ()=delete
	SetPts (boost::shared_ptr< SetPtsData > data_ptr)
MoFEMErrorCode	operator() (int order_row, int order_col, int order_data)

Private Attributes
boost::weak_ptr< SetPtsData >	dataPtr

Detailed Description

Default evaluator for setting integration points.

Definition at line 87 of file FieldEvaluator.hpp.

Constructor & Destructor Documentation

SetPts() [1/2]

MoFEM::FieldEvaluatorInterface::SetPts::SetPts ( )

delete

SetPts() [2/2]

MoFEM::FieldEvaluatorInterface::SetPts::SetPts

(

boost::shared_ptr< SetPtsData >

data_ptr

)

Definition at line 89 of file FieldEvaluator.hpp.

: dataPtr(data_ptr) {}

Member Function Documentation

operator()()

MoFEMErrorCode MoFEM::FieldEvaluatorInterface::SetPts::operator()	(	int	order_row,
		int	order_col,
		int	order_data
	)

Definition at line 76 of file FieldEvaluator.cpp.

                                                         {
  MoFEMFunctionBegin;
 
  if (auto data_ptr = dataPtr.lock()) {
 
    decltype(data_ptr->nbEvalPoints) nbEvalPoints = data_ptr->nbEvalPoints;
    decltype(data_ptr->verb) verb = data_ptr->verb;
 
    decltype(data_ptr->shapeFunctions) &shapeFunctions =
        data_ptr->shapeFunctions;
    decltype(data_ptr->evalPointEntityHandle) &evalPointEntityHandle =
        data_ptr->evalPointEntityHandle;
 
    if (auto fe_ptr = data_ptr->feMethodPtr.lock()) {
 
      VolumeElementForcesAndSourcesCore &fe =
          static_cast<VolumeElementForcesAndSourcesCore &>(*fe_ptr);
 
      MatrixDouble &gauss_pts = fe.gaussPts;
      gauss_pts.resize(4, nbEvalPoints, false);
      gauss_pts.clear();
 
      FTensor::Tensor1<FTensor::PackPtr<double *, 4>, 4> t_shape = {
          &shapeFunctions(0, 0), &shapeFunctions(0, 1), &shapeFunctions(0, 2),
          &shapeFunctions(0, 3)};
      FTensor::Tensor1<FTensor::PackPtr<double *, 1>, 3> t_gauss_pts = {
          &gauss_pts(0, 0), &gauss_pts(1, 0), &gauss_pts(2, 0)};
 
      int nb_gauss_pts = 0;
      for (int nn = 0; nn != nbEvalPoints; ++nn) {
        if (evalPointEntityHandle[nn] ==
            fe.numeredEntFiniteElementPtr->getEnt()) {
          for (const int i : {0, 1, 2}) {
            t_gauss_pts(i) = t_shape(i + 1);
            gauss_pts(3, nb_gauss_pts) = nn;
          }
          ++t_gauss_pts;
          ++nb_gauss_pts;
        }
        ++t_shape;
      }
 
      if (verb >= VERY_NOISY)
        std::cout << "nbEvalOPoints / nbGaussPts: " << nbEvalPoints << " / "
                  << nb_gauss_pts << std::endl;
      gauss_pts.resize(4, nb_gauss_pts, true);
 
      if (verb >= VERY_NOISY)
        std::cout << "gauss pts: " << gauss_pts << std::endl;
 
    } else
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Pointer to element does not exists");
 
  } else
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "Pointer to data does not exists");
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataPtr

boost::weak_ptr<SetPtsData> MoFEM::FieldEvaluatorInterface::SetPts::dataPtr

private

Definition at line 93 of file FieldEvaluator.hpp.

---

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

• src/interfaces/FieldEvaluator.hpp
• src/interfaces/impl/FieldEvaluator.cpp