ParametricFace.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 PARAMETRICFACE_H_
#define PARAMETRICFACE_H_

#include "CMAFwd.h"

class ParametricFace {
 public:
  ParametricFace(std::vector<Point3d> in_ordered_nodes);

  virtual ~ParametricFace() {
  }

  inline void eval(Point2d in_xi,
                  Point3d & out_xyz) {
      v_eval(in_xi, out_xyz);
  }
  
  inline int order(void) {
    return v_order();
  }

  inline void deriv1(Point2d in_xi,
                    Point3d & out_dxyz_dxi1,
                    Point3d & out_dxyz_dxi2) {
      v_deriv1(in_xi, out_dxyz_dxi1, out_dxyz_dxi2);
  }

  inline void deriv1_fd(const Point2d in_xi,
                        Point3d & dxyz_dxi1,
                        Point3d & dxyz_dxi2,
                        double step = 1.0e-3) {
    Point3d pt3, pt2, pt1;
    Point2d xi = in_xi;
    // evaluate the starting point
    this->eval(xi, pt1);
    // evaluate after increment in xi1 direction
    this->eval(xi + Point2d(step, 0.0), pt2);
    // evaluate after increment in xi2 direction
    this->eval(xi + Point2d(0.0, step), pt3);
    dxyz_dxi1 = (pt2 - pt1) / step;
    dxyz_dxi2 = (pt3 - pt1) / step;
  }


 protected:
  virtual void v_eval(Point2d in_xi, Point3d & out_xyz) = 0;
  virtual int v_order(void) = 0;
  virtual void v_deriv1(Point2d in_xi,
                        Point3d & out_xyz1,
                        Point3d & out_xyz2) = 0;

  std::vector<Point3d> m_ordered_nodes;
};

class CADTri {
 public:
  CADTri(pMeshEnt in_vtx0, pMeshEnt in_vtx1, pMeshEnt in_vtx2);

  virtual ~CADTri() {
  }

  inline int eval(Point2d in_xi,
                  Point3d & out_xyz) {
    return v_eval(in_xi, out_xyz);
  }
  
  inline int order(void) {
    return v_order();
  }

  inline int deriv1(Point2d in_xi,
                    Point3d & out_dxyz_dxi1,
                    Point3d & out_dxyz_dxi2) {
    return v_deriv1(in_xi, out_dxyz_dxi1, out_dxyz_dxi2);
  }

  inline void deriv1_fd(const Point2d in_xi,
                        Point3d & dxyz_dxi1,
                        Point3d & dxyz_dxi2,
                        double step = 1.0e-3) {
    Point3d pt3, pt2, pt1;
    Point2d xi = in_xi;
    // evaluate the starting point
    this->eval(xi, pt1);
    // evaluate after increment in xi1 direction
    this->eval(xi + Point2d(step, 0.0), pt2);
    // evaluate after increment in xi2 direction
    this->eval(xi + Point2d(0.0, step), pt3);
    dxyz_dxi1 = (pt2 - pt1) / step;
    dxyz_dxi2 = (pt3 - pt1) / step;
  }

 protected:
  int v_eval(Point2d in_xi, Point3d & out_xyz);
  int v_order();
  int v_deriv1(Point2d in_xi,
               Point3d & out_xyz1,
               Point3d & out_xyz2);
 private:
  pMeshEnt m_vtx0;  // (1, 0, 0)
  pMeshEnt m_vtx1;  // (0, 1, 0)
  pMeshEnt m_vtx2;  // (0, 0, 1)
  pMeshEnt m_face;
  Point3d m_vert0;
  Point3d m_vert1;
  Point3d m_vert2;
};

#endif  // PARAMETRICFACE_H_