Projection10NodeCoordsOnField.hpp

Go to the documentation of this file.

/** \file Projection10NodeCoordsOnField.hpp
 
FIXME: Move code to cpp file.
 
Project displacements/coordinates from 10 node tetrahedra on hierarchical
approximation base.
 
This is example how to use MoFEM::DofMethod when some operator for each node
need to be applied.
 
*/
 
 
 
#ifndef __PROJECTION10NODECOORDSONFIELD_HPP__
#define __PROJECTION10NODECOORDSONFIELD_HPP__
 
using namespace boost::numeric;
 
namespace MoFEM {
 
/** \brief Projection of edge entities with one mid-node on hierarchical basis
 */
struct Projection10NodeCoordsOnField : public DofMethod {
 
  Projection10NodeCoordsOnField(Interface &m_field, std::string field_name,
                                int verb = 0);
 
  MoFEMErrorCode preProcess();
 
  MoFEMErrorCode operator()();
 
  MoFEMErrorCode postProcess();
 
protected:
 
  Interface &mField;
  std::string fieldName;
  int vErbose;
 
  VectorDouble coords;
  VectorDouble3 aveMidCoord;
  VectorDouble3 midNodeCoord;
  VectorDouble3 diffNodeCoord;
  VectorDouble3 dOf;
 
};
 
/** \brief Project edge mid-nodes onto a fine surface mesh using closest points
 */
struct Projection10NodeCoordsFromFineSurfaceOnField : public EntityMethod {
 
  Projection10NodeCoordsFromFineSurfaceOnField(Interface &m_field,
                                               moab::Interface &fine_moab,
                                               const Range &fine_tris,
                                               std::string field_name,
                                               double planar_angle_deg = 5.0,
                                               double max_disp_factor = 0.5,
                                               bool debug_surfaces = false,
                                               int verb = 0);
 
  MoFEMErrorCode preProcess();
 
  MoFEMErrorCode operator()();
 
  MoFEMErrorCode postProcess();
 
protected:
  FTENSOR_INDEX(3, i);
  FTENSOR_INDEX(3, j);
  Interface &mField;
  moab::Interface &fineMoab;
  Range fineTris;
  std::string fieldName;
  int vErbose;
 
  AdaptiveKDTree fineTree;
  EntityHandle fineRoot;
  AdaptiveKDTree fineFeatureTree;
  EntityHandle fineFeatureRoot;
  Range fineFeatureTris;
  double planarAngleCos;
  double maxDispFactor;
  bool debugSurfaces;
  Range coarseSkinTris;
  Range coarseSkinEdges;
 
  VectorDouble coords;
  VectorDouble3 aveMidCoord;
  VectorDouble3 midNodeCoord;
  VectorDouble3 diffNodeCoord;
  VectorDouble3 dOf;
};
 
struct ProjectionFieldOn10NodeTet : public Projection10NodeCoordsOnField {
 
 
  ProjectionFieldOn10NodeTet(Interface &m_field, std::string _fieldName,
                             bool set_nodes, bool on_coords,
                             std::string on_tag = "NoNE");
 
  MoFEMErrorCode preProcess();
 
  MoFEMErrorCode operator()(); 
 
 
  bool setNodes;
  bool onCoords;
  std::string onTag;
 
  const int maxApproximationOrder;
 
  Tag th;
 
protected:
 
  Field_multiIndex::index<FieldName_mi_tag>::type::iterator field_it;
  VectorDouble L;
  VectorDouble K;
 
};
 
} // namespace MoFEM
 
#endif // __PROJECTION10NODECOORDSONFIELD_HPP__