CrvTetPolyMap.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 MESHADAPT_CRVTETPOLYMAP_H_
#define MESHADAPT_CRVTETPOLYMAP_H_

#include "CrvTet.h"

typedef class mMesh * pMeshMdl;
class curveMesh;

class CrvTetPolyMap : public CrvTet {
 public:
    // the only constructor
    CrvTetPolyMap(pMeshEnt in_rgn, curveMesh * in_crv_mesh);
    // the destructor
    ~CrvTetPolyMap();
    // evaluate the mapping (by blending)
    virtual int v_eval(Point3d in_xi, Point3d & pos);
    // evaluate first derivative
    virtual int v_eval_deriv1(Point3d in_xi,
                              Mat3x3 & dxyz_dxi);
 protected:
    // first derivative w.r.t. u
    int dxyz_du(double u, double v, double w, double * out);
    // first derivative w.r.t. v
    int dxyz_dv(double u, double v, double w, double * out);
    // first derivative w.r.t. w
    int dxyz_dw(double u, double v, double w, double * out);

 private:
    // vertex coordinates             
    Point3d m_xyz_3000_;  // V1, xi_1,  u: V0 -> V1 
    Point3d m_xyz_0300_;  // V2, xi_2,  v: V0 -> V2
    Point3d m_xyz_0030_;  // V3, xi_3,  w: V0 -> V3
    Point3d m_xyz_0003_;  // V0, xi_4,  1.0 - u -v - w

    // edge control pts
    Point3d m_ctrl_2100_;
    Point3d m_ctrl_1200_;
    
    Point3d m_ctrl_0210_;
    Point3d m_ctrl_0120_;

    Point3d m_ctrl_0021_;
    Point3d m_ctrl_0012_;

    Point3d m_ctrl_2010_;
    Point3d m_ctrl_1020_;

    Point3d m_ctrl_2001_;
    Point3d m_ctrl_1002_;

    Point3d m_ctrl_0201_;
    Point3d m_ctrl_0102_;

    // face control pts  
    Point3d m_ctrl_1110_;
    Point3d m_ctrl_1101_;
    Point3d m_ctrl_1011_;
    Point3d m_ctrl_0111_;
};

#endif  // MESHADAPT_CRVTETPOLYMAP_H_