MxWedge18CellExtract.h

/*=======================================================================
 *** THE CONTENT OF THIS WORK IS PROPRIETARY TO FEI S.A.S, (FEI S.A.S.),            ***
 ***              AND IS DISTRIBUTED UNDER A LICENSE AGREEMENT.                     ***
 ***                                                                                ***
 ***  REPRODUCTION, DISCLOSURE,  OR USE,  IN WHOLE OR IN PART,  OTHER THAN AS       ***
 ***  SPECIFIED  IN THE LICENSE ARE  NOT TO BE  UNDERTAKEN  EXCEPT WITH PRIOR       ***
 ***  WRITTEN AUTHORIZATION OF FEI S.A.S.                                           ***
 ***                                                                                ***
 ***                        RESTRICTED RIGHTS LEGEND                                ***
 ***  USE, DUPLICATION, OR DISCLOSURE BY THE GOVERNMENT OF THE CONTENT OF THIS      ***
 ***  WORK OR RELATED DOCUMENTATION IS SUBJECT TO RESTRICTIONS AS SET FORTH IN      ***
 ***  SUBPARAGRAPH (C)(1) OF THE COMMERCIAL COMPUTER SOFTWARE RESTRICTED RIGHT      ***
 ***  CLAUSE  AT FAR 52.227-19  OR SUBPARAGRAPH  (C)(1)(II)  OF  THE RIGHTS IN      ***
 ***  TECHNICAL DATA AND COMPUTER SOFTWARE CLAUSE AT DFARS 52.227-7013.             ***
 ***                                                                                ***
 ***                   COPYRIGHT (C) 1996-2017 BY FEI S.A.S,                        ***
 ***                        BORDEAUX, FRANCE                                        ***
 ***                      ALL RIGHTS RESERVED                                       ***
**=======================================================================*/
/*=======================================================================
** Author      : VSG (MMM YYYY)
**=======================================================================*/
#ifndef  _MxWedge18CellExtract_h
#define  _MxWedge18CellExtract_h
 
#include <MeshVizXLM/mesh/geometry/MiGeometryI.h>
 
#include <Inventor/STL/vector>
#include <Inventor/STL/utility>
 
class MiVolumeCell;
 
class MESHVIZXLM_EXTR_API MxWedge18CellExtract 
{
public:
 
  static MbVec3d getIsoParametricCoord(const MiGeometryI& meshGeometry, const MiVolumeCell* wedgeCell, const MbVec3d &point);
  
  static MbVec3d getIsoParametricCoord(size_t nodeIndex);
 
  static void getWeight(const MiGeometryI& meshGeometry, const MiVolumeCell* wedgeCell, const MbVec3d &point, std::vector<double>& weights)
  {
    MbVec3d ipcoord = getIsoParametricCoord(meshGeometry,wedgeCell,point);
    getWeight(ipcoord,weights);
  }
 
  static void getWeight(const MbVec3d &ipcoord, std::vector<double>& weights);
 
  static void getDerivs(const MbVec3d &ipcoord, std::vector<double>& derivs);
 
  static bool isPointInsideCell(const MiGeometryI& meshGeometry, const MiVolumeCell* wedgeCell, const MbVec3d &point, std::vector<double>& weights)
  {
    MbVec3d pcoord = getIsoParametricCoord(meshGeometry,wedgeCell,point);
    bool inside =  ( -1.E-5 < pcoord[0] &&  pcoord[0] < 1 + 1.E-5 &&
      -1.E-5 < pcoord[1] &&  pcoord[1] < 1 + 1.E-5 &&
      -1.E-5 < pcoord[2] &&  pcoord[2] < 1 + 1.E-5 &&
      -1.E-5 < (pcoord[0]+pcoord[1]) &&  (pcoord[0]+pcoord[1]) < 1 + 1.E-5);
    if (inside)
      getWeight(pcoord,weights);
    return inside;
  }
 
  static size_t getNodesIndexOfFacet (size_t facet, std::vector<size_t>& facetNodes)
  {
    facetNodes.assign(s_linearNodesPerFacet[facet],s_linearNodesPerFacet[facet]+s_numLinearNodesPerFacet[facet]);
    return s_numLinearNodesPerFacet[facet];
  }
 
  static size_t getSubTriangleNodesIndex(size_t face, std::vector<size_t>& triangleNodeIds)  
  {
    triangleNodeIds.assign(s_subTrianglePerFacet[face],s_subTrianglePerFacet[face]+(3*s_numSubTrianglePerFacet[face]));
    return s_numSubTrianglePerFacet[face];
  }
 
  static size_t getSubTetrahedronNodesIndex (std::vector<size_t>& tetrahedronNodeIds)
  {
    tetrahedronNodeIds.assign(s_subTetrahedronNodeIds,s_subTetrahedronNodeIds+(24*4));
    return 24;
  }
 
 
private:
  static void get18Derivs(double r, double s, double t, size_t numNodes, std::vector<double>& derivs);
 
  friend std::ostream& operator << (std::ostream& s, const MxWedge18CellExtract& cell);
private:
 
  static bool initClass();
  static void exitClass();
 
  static void computeSubTetrahedronNodesIndex();
  static void addSubWedge(size_t n0,size_t n1,size_t n2, size_t n3,size_t n4,size_t n5, size_t orient, std::vector<size_t>& wedgeNodeIds);
 
  template <typename _F>
  static void addSymetricWedgesByPlane(_F& symetricByPlane, std::vector<size_t>& wedgeNodeIds);
  static void decomposeWedge(size_t orient, size_t n0,size_t n1,size_t n2, size_t n3,size_t n4,size_t n5, std::vector<size_t>& tetraNodeIds);
 
  static void computeSubTriangleNodesIndex();
  template <typename _F>
  static std::vector<size_t> getSymetricTriangles(_F& symetricFunctor, const std::vector<size_t>& triangleNodeIds);
  template <typename _F>
  static std::vector<size_t> getRotateTriangles(_F& rotateFunctor, const std::vector<size_t>& triangleNodeIds);
  static void addTriangle(size_t n0, size_t n1, size_t n2, std::vector<size_t>& triangleNodeIds);
 
  static bool            s_isInit;
 
  static double s_nodesIsoParametricCoords[18][3];
 
  static size_t s_linearNodesPerFacet[5][4];
  static size_t s_numLinearNodesPerFacet[5];
 
  static size_t s_numSubTrianglePerFacet[5];
  static size_t s_subTrianglePerFacet[5][24]; 
 
  static size_t s_subTetrahedronNodeIds[24*4]; 
 
  static int s_imat0[36];
  static int s_imat1[24];
  static int s_imat2[24];
  static int s_imat3[18];
  static int s_iperm[24];
};
 
#endif