10#ifndef __FIELD_EVALUATOR_HPP__
11 #define __FIELD_EVALUATOR_HPP__
42 SetPtsData(boost::shared_ptr<MoFEM::ForcesAndSourcesCore> fe_method_ptr,
43 const double *eval_points,
const int nb_eval_points,
51 inline void setEvalPoints(
const double *ptr,
const int nb_eval_points) {
69 boost::scoped_ptr<AdaptiveKDTree>
treePtr;
80 int order_col,
int order_data);
105 template <
typename VE,
typename SPD = SetPtsData,
typename SP = SetPts>
106 boost::shared_ptr<SPD>
107 getData(
const double *ptr =
nullptr,
const int nb_eval_points = 0,
110 boost::scoped_ptr<VE> elePtr;
111 boost::scoped_ptr<SPD> setPtsDataPtr;
112 boost::scoped_ptr<SP> setPtsPtr;
114 boost::shared_ptr<PackData> pack_data(
new PackData());
116 pack_data->elePtr.reset(
new VE(m_field));
118 pack_data->setPtsDataPtr.reset(
119 new SPD(boost::shared_ptr<VE>(pack_data, pack_data->elePtr.get()), ptr,
120 nb_eval_points,
eps, verb));
121 pack_data->setPtsPtr.reset(
new SP(
122 boost::shared_ptr<SPD>(pack_data, pack_data->setPtsDataPtr.get())));
123 pack_data->elePtr->setRuleHook = boost::ref(*pack_data->setPtsPtr);
124 boost::shared_ptr<SPD> data(pack_data, pack_data->setPtsDataPtr.get());
131 const std::string finite_element);
136 const std::string finite_element);
140 const double *
const point,
const double distance,
141 const std::string problem,
const std::string finite_element,
142 boost::shared_ptr<SetPtsData> data_ptr,
int lower_rank,
int upper_rank,
148 const double *
const point,
const double distance,
149 const std::string problem,
const std::string finite_element,
150 boost::shared_ptr<SetPtsData> data_ptr,
int lower_rank,
int upper_rank,
162 const std::string finite_element);
210 const std::string problem,
const std::string finite_element,
211 boost::shared_ptr<SetPtsData> data_ptr,
int lower_rank,
212 int upper_rank, boost::shared_ptr<CacheTuple> cache_ptr,
218 const std::string finite_element);
223 const std::string problem,
224 const std::string finite_element,
225 boost::shared_ptr<SetPtsData> data_ptr,
int lower_rank,
226 int upper_rank, boost::shared_ptr<CacheTuple> cache_ptr,
233 const std::string finite_element);
238 const std::string finite_element);
242 const double *
const point,
const double distance,
const std::string problem,
243 const std::string finite_element, boost::shared_ptr<SetPtsData> data_ptr,
244 int lower_rank,
int upper_rank, boost::shared_ptr<CacheTuple> cache_ptr,
249 const double *
const point,
const double distance,
const std::string problem,
250 const std::string finite_element, boost::shared_ptr<SetPtsData> data_ptr,
251 int lower_rank,
int upper_rank, boost::shared_ptr<CacheTuple> cache_ptr,
VERBOSITY_LEVELS
Verbosity levels.
MoFEMTypes
Those types control how functions respond on arguments, f.e. error handling.
PetscErrorCode MoFEMErrorCode
MoFEM/PETSc error code.
implementation of Data Operators for Forces and Sources
Deprecated interface functions.
MatrixDouble shapeFunctions
boost::weak_ptr< MoFEM::ForcesAndSourcesCore > feMethodPtr
void setEvalPoints(const double *ptr, const int nb_eval_points)
boost::scoped_ptr< AdaptiveKDTree > treePtr
SetPtsData(boost::shared_ptr< MoFEM::ForcesAndSourcesCore > fe_method_ptr, const double *eval_points, const int nb_eval_points, const double eps, VERBOSITY_LEVELS verb=QUIET)
Set the Gauss Pts data.
std::vector< EntityHandle > evalPointEntityHandle
const double * evalPoints
Default evaluator for setting integration points.
MoFEMErrorCode operator()(ForcesAndSourcesCore *fe_raw_ptr, int order_row, int order_col, int order_data)
SetPts(boost::shared_ptr< SetPtsData > data_ptr)
boost::weak_ptr< SetPtsData > dataPtr
Field evaluator interface.
DEPRECATED MoFEMErrorCode buildTree3D(boost::shared_ptr< SetPtsData > spd_ptr, const std::string finite_element)
MoFEMErrorCode query_interface(boost::typeindex::type_index type_index, UnknownInterface **iface) const
DEPRECATED MoFEMErrorCode buildTree2D(boost::shared_ptr< SetPtsData > spd_ptr, const std::string finite_element)
MoFEMErrorCode evalFEAtThePointImpl(const double *const point, const double distance, const std::string problem, const std::string finite_element, boost::shared_ptr< SetPtsData > data_ptr, int lower_rank, int upper_rank, boost::shared_ptr< CacheTuple > cache_ptr, MoFEMTypes bh=MF_EXIST, VERBOSITY_LEVELS verb=QUIET)
DEPRECATED MoFEMErrorCode evalFEAtThePoint2D(const double *const point, const double distance, const std::string problem, const std::string finite_element, boost::shared_ptr< SetPtsData > data_ptr, int lower_rank, int upper_rank, boost::shared_ptr< CacheTuple > cache_ptr, MoFEMTypes bh=MF_EXIST, VERBOSITY_LEVELS verb=QUIET)
boost::shared_ptr< SPD > getData(const double *ptr=nullptr, const int nb_eval_points=0, const double eps=1e-12, VERBOSITY_LEVELS verb=QUIET)
Get the Data object.
FieldEvaluatorInterface(const MoFEM::Core &core)
MoFEMErrorCode buildTreeImpl(boost::shared_ptr< SetPtsData > spd_ptr, const std::string finite_element)
DEPRECATED MoFEMErrorCode evalFEAtThePoint3D(const double *const point, const double distance, const std::string problem, const std::string finite_element, boost::shared_ptr< SetPtsData > data_ptr, int lower_rank, int upper_rank, boost::shared_ptr< CacheTuple > cache_ptr, MoFEMTypes bh=MF_EXIST, VERBOSITY_LEVELS verb=QUIET)
MoFEMErrorCode buildTree(boost::shared_ptr< SetPtsData > spd_ptr, const std::string finite_element)
Build spatial tree.
MoFEMErrorCode evalFEAtThePoint(const double *const point, const double distance, const std::string problem, const std::string finite_element, boost::shared_ptr< SetPtsData > data_ptr, int lower_rank, int upper_rank, boost::shared_ptr< CacheTuple > cache_ptr, MoFEMTypes bh=MF_EXIST, VERBOSITY_LEVELS verb=QUIET)
Evaluate field at artbitray position.
structure to get information form mofem into EntitiesFieldData
base class for all interface classes