MxTetrahedronCellExtract.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  _MxTetrahedronCellExtract_h
#define  _MxTetrahedronCellExtract_h
 
#include <MeshVizXLM/mesh/geometry/MiGeometryI.h>
 
#include <Inventor/STL/vector>
#include <Inventor/STL/utility>
 
class MxTabOfMarchingCase;
class MiVolumeCell;
 
class MESHVIZXLM_EXTR_API MxTetrahedronCellExtract 
{
public:
 
  static void getIsosurfTopology(unsigned char caseId, std::vector< std::pair<size_t,size_t> >& edgeList );
 
  static MbVec3d getIsoParametricCoord(const MiGeometryI& meshGeometry, const MiVolumeCell* tetrahedronCell, const MbVec3d &point);
  
  static MbVec3d getIsoParametricCoord(size_t nodeIndex)
  {
    return MbVec3d(s_nodesIsoParametricCoords[nodeIndex]);
  }
 
  static void getWeight(const MiGeometryI& meshGeometry, const MiVolumeCell* tetrahedronCell, const MbVec3d &point, std::vector<double>& weights)
  {
    MbVec3d ipcoord = getIsoParametricCoord(meshGeometry,tetrahedronCell,point);
    getWeight(ipcoord,weights);
  }
 
  static void getWeight(const MbVec3d &ipcoord, std::vector<double>& weights)
  {
    ipcoord.getValue(weights[1],weights[2],weights[3]);
    weights[0] = 1.0 - weights[1] - weights[2] - weights[3];
  }
 
  static void getDerivs(const MbVec3d &ipcoord, std::vector<double>& derivs);
 
  static bool isPointInsideCell(const MiGeometryI& meshGeometry, const MiVolumeCell* tetrahedronCell, const MbVec3d &point, std::vector<double>& weights)
  {
    getWeight(getIsoParametricCoord(meshGeometry,tetrahedronCell,point),weights);
    return ( -1.E-5 < weights[0] &&  weights[0] < 1 + 1.E-5 &&
             -1.E-5 < weights[1] &&  weights[1] < 1 + 1.E-5 &&
             -1.E-5 < weights[2] &&  weights[2] < 1 + 1.E-5 &&
             -1.E-5 < weights[3] &&  weights[3] < 1 + 1.E-5);      
  }
 
  static double getVolume(const MiGeometryI& meshGeometry, const MiVolumeCell* tetrahedronCell);
 
  static MbVec3d getFacetNormal(size_t facetId, const MiGeometryI& meshGeometry, const MiVolumeCell* tetrahedronCell);
 
  static double getLongestEdgeLength(const MiGeometryI& meshGeometry, const MiVolumeCell* cell);
 
  static double getShortestEdgeLength(const MiGeometryI& meshGeometry, const MiVolumeCell* cell);
 
private:
 
  static bool initClass();
  static void exitClass();
 
  static void initMarchingCase0();
  static void initMarchingCase1();
  static void initMarchingCase2();
 
  static bool                 s_isInit;
  static MxTabOfMarchingCase* s_tabOfMarchingTetra;
  static double               s_nodesIsoParametricCoords[4][3];
};
 
#endif