MoFEM::FieldEvaluatorInterface::SetPts Struct Reference

Default evaluator for setting integration points. More...

#include <src/interfaces/FieldEvaluator.hpp>

Collaboration diagram for MoFEM::FieldEvaluatorInterface::SetPts:

Public Member Functions
	SetPts ()=delete
	SetPts (boost::shared_ptr< SetPtsData > data_ptr)
MoFEMErrorCode	operator() (ForcesAndSourcesCore *fe_raw_ptr, 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 76 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 )

inline

Definition at line 78 of file FieldEvaluator.hpp.

: dataPtr(data_ptr) {}

Member Function Documentation

operator()()

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

Definition at line 83 of file FieldEvaluator.cpp.

                                                            {
  MoFEMFunctionBegin;
 
  if (auto data_ptr = dataPtr.lock()) {
 
    const auto nb_eval_points = data_ptr->nbEvalPoints;
    const auto &shape_functions = data_ptr->shapeFunctions;
    const auto &eval_pointentity_handle = data_ptr->evalPointEntityHandle;
 
    if (auto fe_ptr = data_ptr->feMethodPtr.lock()) {
 
      auto &fe = static_cast<ForcesAndSourcesCore &>(*fe_ptr);
 
#ifndef NDEBUG
      if (fe_ptr.get() != fe_raw_ptr)
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Wrong FE ptr");
#endif
 
      const auto fe_ent = fe.numeredEntFiniteElementPtr->getEnt();
      const auto fe_type = type_from_handle(fe_ent);
      const auto fe_dim = moab::CN::Dimension(fe_type);
 
      auto &gauss_pts = fe.gaussPts;
      int nb_gauss_pts;
 
      if (fe_dim == 3) {
        gauss_pts.resize(4, nb_eval_points, false);
        FTensor::Tensor1<FTensor::PackPtr<const double *, 4>, 4> t_shape{
            &shape_functions(0, 0), &shape_functions(0, 1),
            &shape_functions(0, 2), &shape_functions(0, 3)};
        FTensor::Tensor1<FTensor::PackPtr<double *, 1>, 3> t_gauss_pts{
            &gauss_pts(0, 0), &gauss_pts(1, 0), &gauss_pts(2, 0)};
 
        nb_gauss_pts = 0;
        for (int nn = 0; nn != nb_eval_points; ++nn) {
          if (eval_pointentity_handle[nn] == fe_ent) {
            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;
        }
        gauss_pts.resize(4, nb_gauss_pts, true);
      } else if (fe_dim == 2) {
        gauss_pts.resize(3, nb_eval_points, false);
        FTensor::Tensor1<FTensor::PackPtr<const double *, 3>, 3> t_shape{
            &shape_functions(0, 0), &shape_functions(0, 1),
            &shape_functions(0, 2)};
        FTensor::Tensor1<FTensor::PackPtr<double *, 1>, 2> t_gauss_pts{
            &gauss_pts(0, 0), &gauss_pts(1, 0)};
        nb_gauss_pts = 0;
        for (int nn = 0; nn != nb_eval_points; ++nn) {
          if (eval_pointentity_handle[nn] == fe_ent) {
            for (const int i : {0, 1}) {
              t_gauss_pts(i) = t_shape(i + 1);
            }
            gauss_pts(2, nb_gauss_pts) = nn;
            ++t_gauss_pts;
            ++nb_gauss_pts;
          }
          ++t_shape;
        }
        gauss_pts.resize(3, nb_gauss_pts, true);
      } else {
        SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
                "Dimension not implemented");
      }
 
#ifndef NDEBUG
      MOFEM_LOG("SELF", Sev::noisy)
          << "nbEvalOPoints / nbGau_sspt_s: " << nb_eval_points << " / "
          << nb_gauss_pts;
      MOFEM_LOG("SELF", Sev::noisy) << "gauss pts: " << gauss_pts;
#endif
 
    } 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 83 of file FieldEvaluator.hpp.

---

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

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