ParametricCurve.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 PARAMETRICCURVE_H
#define PARAMETRICCURVE_H

#include "CMAFwd.h"

class ParametricCurve{
  public:
    virtual ~ParametricCurve() {
    }

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

    inline int deriv1(Point1d in_t, Point3d & out_xyz)
    {
      return v_deriv1(in_t, out_xyz);
    }

    inline void deriv1_fd(const Point1d in_t, Point3d & dxyz_dt, double step = 1.0e-3) {
      Point3d pt2, pt1;
      Point1d t = in_t;
      this->eval(t + step, pt2);
      this->eval(t,        pt1);
      dxyz_dt = (pt2 - pt1) / step;
    }

  protected:
    virtual int v_eval(Point1d in_t, Point3d & out_xyz) = 0;
    virtual int v_order(void) = 0;
    virtual int v_deriv1(Point1d in_t, Point3d & out_xyz) = 0;

};
#endif //PARAMETRICCURVE_H