CrvEdge.h

Go to the documentation of this file.

/****************************************************************************** 

  (c) 2004-2014 Scientific Computation Research Center, 
      Rensselaer Polytechnic Institute. All rights reserved.
  
  The LICENSE-SCOREC file included with this distribution describes the terms
  of the SCOREC Non-Commercial License this program is distributed under.
 
*******************************************************************************/

#ifndef CURVEMESH_CRVEDGE_H
#define CURVEMESH_CRVEDGE_H

#include "CrvEnt.h"
#include "ParametricCurve.h"

namespace pumi {
typedef class mEntity * pNode;
}

class CrvEdge : public CrvEnt {

public:
  CrvEdge(pMeshEdge pin_edge);

  CrvEdge(VtxPtrVec in_vert_vec);

  virtual ~CrvEdge();

  /* common interface */

  int display(std::ostream& file, int in_sample_pts);

  virtual int v_ent_type() const;

  void eval_deriv1_fd(const Point1d in_param, Point3d & dxyz_dxi, double step = 1.0e-3);

  double arc_length_between(double a, double b);

  double arc_length();


  void get_interp_point_xyz(int index, Point3d & out, int mode = 0);

  double interp_error();

  void get_edge_ctrl_pt(int index, Point3d & out_pt);

  /* pure virtual functions */

  int get_edge_node_coords(int in_node_index, double *xyz);

  virtual int v_data_size() const = 0;

  virtual std::string v_tag_name() const = 0;

  double dxi_deta(double eta) {
    return -0.5;
  }
  int get_vertex_coords(double *vxyz0, double *vxyz1);   

  int get_orient() {
    return m_orient;
  }

  int get_geom_order() {
    return m_geom_shape->order();
  }

  int LinearMidPt(double * dout_coords);

  int GetNumNode();

  int DelNode(int n);

  pumi::pNode GetNode(int n);

  void get_mid_edge_ctrl_pt(double * out_xyz);

  void eval_on_mesh_geom(double in_t, Point3d & out_xyz);

  int update_or_create_node(double *xyz, double *param);

  void save_ctrl_pt_data();
  
 protected:
  virtual void v_eval_by_xi (Point3d in, Point3d & out);

  virtual void v_eval_by_eta(Point3d in, Point3d & out);

  virtual void v_eval_by_coll_eta(Point3d in, Point3d & out);

  virtual void v_diff_by_xi(Point3d in, Mat3x3 & out);

  virtual void v_diff_by_xi(Point3d in, int j, Point3d & out);

  virtual void v_diff_by_eta(Point3d in, Mat3x3 & out);

  virtual void v_diff_by_eta(Point3d in, int j, Point3d & out);

  virtual int v_get_exp_dim() {
    return 1;
  }

  pCurveGeom get_geom_shape() {
    return m_geom_shape;
  }

  void set_geom_shape(pCurveGeom & in_geom) {
    m_geom_shape = in_geom;
  }

  pCurveGeom m_geom_shape;

  int m_orient;
 private:

  void get_interp_point_bisect(int index, Point3d & out);

  void get_interp_point_newton_iter(int index, Point3d & out);
};
#endif//CURVEMESH_CRVEDGE_H