GregorySurface.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 CURVEMESHADAPT_GREGORY_H
#define CURVEMESHADAPT_GREGORY_H

#include <map>

#include "CMAFwd.h"

class EdgeCtrlPt; 

class GregorySurface {
  public:
    GregorySurface();

    GregorySurface(std::vector<Point3d> in_ordered_nodes,
                   std::vector<Point3d> in_ordered_normals);

    ~GregorySurface();

    int eval(double u, double v, double *retval);

    void deriv1(Point2d in_xi,
                Point3d & out_dxyz_dxi1,
                Point3d & out_dxyz_dxi2);

    void InputPos(double *q0, double *q1, double *q2);
    void InputNormal(double *n0, double *n1, double *n2);

    void SetEdgeCtrlPoints(int i, Point3d *ctrl_pts);

    void Setup();

    void SetupBdryCrv(int i);

    void print_var();

    Point3d m_G01_, m_G02_, m_G11_, m_G12_, m_G21_, m_G22_;

  protected:
    int snap_to_model();

    int m_order_;
    int m_dim_;

    Point3d m_q0_;
    Point3d m_q1_;
    Point3d m_q2_;

    Point3d m_n0_;
    Point3d m_n1_;
    Point3d m_n2_;

    // P -- edge control points
    Point3d P004;
    Point3d P013, P103;
    Point3d P022,       P202;
    Point3d P031,             P301;
    Point3d P040, P130, P220, P310, P400;

    // W, A -- aux matrices
    Point3d W[3][3];   
    Point3d A[3][3];

    std::vector<EdgeCtrlPt*> edge_ctrl_points;

};
#endif//CURVEMESHADAPT_GREGORY_H