SmoothCurve.hpp

Go to the documentation of this file.

/*
 * SmoothCurve.hpp
 *
 */

#ifndef SMOOTHCURVE_HPP_
#define SMOOTHCURVE_HPP_

#include "moab/Interface.hpp"
#include "moab/Forward.hpp"
#include "moab/CartVect.hpp"

//#include "RefEdge.hpp"
//#include "SmoothFace.hpp"

#include <map>
#include <vector>

namespace moab {
class SmoothFace;
class GeomTopoTool;
//class SmoothVertex;

//class CMLEdgeMesher;
// evaluator for Camal Edge Mesher
// this is really copying what Cubit is doing

class SmoothCurve
{
public:
    //SmoothCurve(RefEdge * edge, SmoothFace * smoothFaceEval, int loopIndex);
    SmoothCurve(Interface * mb, EntityHandle curve,  GeomTopoTool * gTool); // the new constructor, which will use
    // sense entities to establish the control points on feature edges (geo edges, sets of mesh edges)
    virtual ~SmoothCurve();

    virtual double arc_length() ;

      virtual bool is_parametric() ;

      virtual bool is_periodic(double& period) ;

      virtual void get_param_range(double& u_start, double& u_end) ;

      virtual double u_from_arc_length(double u_root, double arc_length);


       // ! \param tg, if not null, return the tangent too at u
      virtual bool position_from_u(double u,
                                   double& x, double& y, double& z,
                                   double * tg = NULL) ;

      virtual void move_to_curve(double& x, double& y, double& z);

      virtual double u_from_position(double x, double y, double z, EntityHandle & v,
          int & indexEdge) ;

      virtual void start_coordinates(double& x, double& y, double& z) ;

      virtual void end_coordinates(double& x, double& y, double& z) ;

      // this will recompute the 2 tangents for each edge, considering the geo edge they are into
      void compute_tangents_for_each_edge();

      void compute_control_points_on_boundary_edges(double min_dot,
              std::map<EntityHandle, SmoothFace*> & mapSurfaces,
              Tag controlPointsTag, Tag markTag);

      ErrorCode evaluate_smooth_edge(EntityHandle eh, double &tt,
            CartVect & outv, CartVect & out_tangent) ;

private:

      std::vector<EntityHandle> _entities;// the mesh edges are stored here for fast access
      double _leng;
      std::vector<double> _fractions;// they are increasing from 0. to 1., do we need these?
      // this will be decided apriori, and eventually reset for paver
      // fractions will be from 0.0.. to 1.0, they will be decided upon the length of the geo edge

      Tag _edgeTag;

      Interface * _mb;
      EntityHandle _set;
      GeomTopoTool * _gtt;



};

} // namespace moab
#endif /* SMOOTHCURVE_HPP_ */