CrvTet.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_CRVTET_H
#define MESHADAPT_CRVTET_H

#include "CrvEnt.h"

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

  CrvTet(pMeshEnt in_tet);

  virtual ~CrvTet();

  void deriv1_fd(const Point3d in_xi, Mat3x3 & dxyz_dxi, double step = 1.0e-3);

  void deriv1_by_eta_fd(const Point3d in_eta, Mat3x3 & dxyz_deta, double step = 1.0e-3);

  void dxi_deta(const Point3d in_eta, Mat3x3 & out_dxi_deta);

  virtual int v_ent_type() const;

  virtual int v_data_size() const;

  virtual std::string v_tag_name() const;

// **** the two shape measure for the ST-frame **** //

  double RadiusRatio();

  double AspectRatio();

// **** curved mesh shape measure **** //

  double JacRatio();

 protected:
  void setup_vertices();

  virtual void v_eval_by_xi (Point3d in, Point3d & out);

  virtual void v_eval_by_eta(Point3d in, Point3d & out);

  virtual void v_eval_by_coll_eta(Point3d in, Point3d & out);

  virtual void v_diff_by_xi(Point3d in, Mat3x3 & out);

  virtual void v_diff_by_xi(Point3d in, int j, Point3d & out);

  virtual void v_diff_by_eta(Point3d in, Mat3x3 & out);

  virtual void v_diff_by_eta(Point3d in, int j, Point3d & out);

  virtual int v_get_exp_dim() {
    return 3;
  }

  EdgePtrVec m_edges;

  FacePtrVec m_faces;

 private:
  virtual int v_eval(Point3d in_xi, Point3d & out_pos) = 0;

  virtual int v_eval_deriv1(Point3d in_xi, Mat3x3 & out_dxyz_dxi) = 0;
};
#endif//MESHADAPT_CRVTET_H