CrvEnt.h

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

  (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_CRVENT_H
#define CURVEMESH_CRVENT_H

#include <string>
#include <vector>
#include <iostream>

#include "CMAFwd.h"

class CrvEnt {

public:
  CrvEnt(pMeshEnt in_mesh_ent);

  CrvEnt(VtxPtrVec in_vert_vec);

  virtual ~CrvEnt();

  /* coordinate system transformation */

  void get_xi_by_eta(Point3d in_eta, Point3d & out_xi);

  void get_eta_by_ceta(Point3d in_ceta, Point3d & out_eta);

  void get_xi_by_ceta(Point3d in_ceta, Point3d & out_xi);

  /* evaluation */

  void eval_by_xi(Point3d in, Point3d & out);

  void eval_by_eta(Point3d in, Point3d & out);

  void eval_by_coll_eta(Point3d in, Point3d & out);

  /* differentiate */

  void diff_by_xi(Point3d in, Mat3x3 & out);

  void diff_by_xi(Point3d in, int j, Point3d & out);

  void diff_by_eta(Point3d in, Mat3x3 & out);

  void diff_by_eta(Point3d in, int j, Point3d & out); 

  /* helper */

  void set_id(int i);

  int get_id();

  int get_exp_dim();

  int get_coord_dim();

  pMeshEnt get_vertex(int i);

  int ent_type();

  int data_size();

  std::string tag_name();

  /* common interface */

  int save_nodal_data_to_mesh(double * in_data);

  int load_nodal_data_from_mesh(double ** out_data);

  bool has_data();

  bool edge_on_model_bdry(pMeshEnt vtx0, pMeshEnt vtx1);

  bool face_on_model_bdry(pMeshEnt vtx0, pMeshEnt vtx1, pMeshEnt vtx2);

 protected:

  /* pure virtual functions */

  /* -- eval -- */
  virtual void v_eval_by_xi (Point3d in, Point3d & out) = 0;
  virtual void v_eval_by_eta(Point3d in, Point3d & out) = 0;
  virtual void v_eval_by_coll_eta(Point3d in, Point3d & out) = 0;

  /* -- diff -- */
  virtual void v_diff_by_xi(Point3d in, Mat3x3 & out) = 0;
  virtual void v_diff_by_xi(Point3d in, int j, Point3d & out) = 0;
  virtual void v_diff_by_eta(Point3d in, Mat3x3 & out) = 0;
  virtual void v_diff_by_eta(Point3d in, int j, Point3d & out) = 0;

  /* -- help -- */
  virtual int  v_ent_type() const = 0;
  virtual int  v_data_size() const = 0;
  virtual std::string v_tag_name() const = 0;
  virtual int  v_get_exp_dim() = 0;

  virtual int v_get_coord_dim() {
    return 3;
  }

  pMeshEnt get_mesh_ent();

  void set_mesh_ent(pMeshEnt in_ent);

  VtxPtrVec m_vert_vec;

  pMeshEnt m_mesh_ent;

  std::vector<Point3d> m_ctrl_pt_vec;

 private:
  int m_id;
};
#endif//CURVEMESH_CRVENT_H