CrvTetBlending.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_CRVTETBLENDING_H
#define MESHADAPT_CRVTETBLENDING_H

#include "CrvTet.h"
#include <stdlib.h>
#include <vector>

typedef class mEntity * pMeshEnt;
class ParametricCurve;

class CrvTetBlending : public CrvTet {
 public:
    CrvTetBlending(VtxPtrVec in_vert_vec);

    // the destructor
    ~CrvTetBlending();

 private:
   
    CrvEdgePtrVec m_crv_edges;

    CrvFacePtrVec m_crv_faces;

//    std::vector<pTriGeom> m_face_shape_vector;

    // evaluate the mapping (by blending)
    virtual int v_eval(Point3d in_xi, Point3d & out_pos);
    // evaluate first derivative
    virtual int v_eval_deriv1(Point3d in_xi,
                              Mat3x3 & out_jac);
    // evaluate first derivative
    virtual int v_eval_deriv1_old(Point3d in_xi,
                                Mat3x3 & out_jac);

    void setup_edges();

    void setup_faces();

    void face_blending_eval(int index,
                                   double in_xi_1,
                                   double in_xi_2,
                                   double in_xi_3,
                                   double & out_fb);

    void face_blending_eval_deriv1(int index,
                                   double in_xi_1,
                                   double in_xi_2,
                                   double in_xi_3,
                                   Point3d & dfb_dxi);

    void projected_face_coords(int face_index,
                                   double in_xi1,
                                   double in_xi2,
                                   double in_xi3,
                                   Point2d & out_xip);

    void projected_face_coords_deriv1(int face_index,
                                      double in_xi1,
                                      double in_xi2,
                                      double in_xi3,
                                      Point3d & dxip1_dxi,
                                      Point3d & dxip2_dxi);
                                   
    void edge_blending_eval(int index,
                                   double in_xi_1,
                                   double in_xi_2,
                                   double in_xi_3,
                                   double & out_eb);

    void edge_blending_eval_deriv1(int index,
                                   double in_xi_1,
                                   double in_xi_2,
                                   double in_xi_3,
                                   Point3d & deb_dxi);

    void projected_edge_coords(int edge_index,
                                   double in_xi1,
                                   double in_xi2,
                                   double in_xi3,
                                   Point1d & out_xip);

    void projected_edge_coords_deriv1(int edge_index,
                                      double in_xi1,
                                      double in_xi2,
                                      double in_xi3,
                                      Point3d & dxip1_dxi);
};
#endif//MESHADAPT_CRVTETBLENDING_H