curveUtil.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_CURVEUTIL_H_
#define MESHADAPT_CURVEUTIL_H_

#include "AdaptTypes.h"
#include "MeshTools.h"

#include <pumi.h>
#include <pumi_mesh.h>

#include <iostream>
#include <vector>
#include <string>

class CurvedRegionShapeInfo{
 public:
  int index;
  int index1[3];
  int isCorner;
  int processed;
  double shape;
  double qs;
  double qc;      //same as above
  //double minDetJ;  // minimum of the Jacobian Determinant evaluated over a curved tet
  //double maxDetJ;  // maximum of the above
};
typedef class CurvedRegionShapeInfo crShpInfo;

class ReShapeEdgeInfo{
public:
  pMeshEdge edge;
  double shape;
  double ptxyz[3];
  int newR;
};
typedef class ReShapeEdgeInfo rEdShpInfo;

int XYZ_RadiusRatio(double (*)[3], double&);

void FIXME_LATER();
void CRV_ERROR();

namespace curveUtil {

  void gauss_legendre_quadrature(int num_pts, double * pts, double * weights);

  int R_numEdgesOnGeomRegn(pMeshRgn pTheRgn);

  int R_numEdgesOnGeomFace(pMeshRgn pTheRgn);

  int R_numEdgesOnGeomEdge(pMeshRgn pTheRgn);

  int R_numFacesOnGeomRegn(pMeshRgn pTheRgn);

  int R_numFacesOnGeomFace(pMeshRgn pTheRgn);

  void M_checkSurfaceMeshByAngle(pPart pmesh, double angle, std::vector<pMeshFace> &vec_facesLargeAngle);

#ifdef MA_PARALLEL
  int CR_isOnPartBdry(pMeshRgn region);
#endif

  void F_bezierSplit(pMeshFace face, pMeshEdge edge, double *ptxyz);
  pMeshVtx F_opEdgeVert(pMeshFace face, pMeshEdge edge);
  void computePtForOppEdge(pMeshEdge oppEdge, pMeshFace *faces, double *xyz);
  void writeRgnDataForMaple(pMeshRgn, char*);
  int HO_XYZ_isValid(double (*)[3], crShpInfo*, pMeshRgn region=NULL);
  int CR_isValid(pMeshRgn, crShpInfo*);

  int XYZ_JacRatio(double (*)[3], double &shape); 

  pumi::pNode extractPt(pPart msh, pMeshFace mf, pMeshEdge edge, double u, double v);// for face
  void F_bezierXYZ(pMeshFace face, double (*fxyz)[3]);
  void F_bezierEval(double (*fxyz)[3], double u, double v, double *ptxyz);

  void writeOutCorner(pMeshRgn region, pMeshVtx vertex, char* name);
 
  void M_checkCurvedRegionShape(pPart, char*);
  void writeOutCurvedRegion(pMeshEnt, int, char *);
  void writeOutCurvedRList(pumi::pPList *, int, char*);
  int R_worstShape(pumi::pPList *, double &, double &,  pumi::pPList *);

  /* --- tagged data communication --- */

  void unifyTaggedEntitiesWithValuesDbl(pMeshMdl FmdbMesh, pTag FmdbTag, std::list<mEntity*> *LIST);
}

#endif //MESHADAPT_CURVEUTIL_H_