MeshVizXLM/extractors/MxPyramidCellExtract.h

Go to the documentation of this file.

/*=======================================================================
 *** 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  _MxPyramidCellExtract_h
#define  _MxPyramidCellExtract_h

#include <MeshVizXLM/mesh/geometry/MiGeometryI.h>
#include <MeshVizXLM/extractors/MxHalfPyramidCell.h>

#include <Inventor/STL/vector>
#include <Inventor/STL/utility>

class MxTabOfMarchingCase;
class MiVolumeCell;

class MESHVIZXLM_EXTR_API MxPyramidCellExtract 
{
public:

  static void getIsosurfTopology(unsigned char caseId, std::vector< std::pair<size_t,size_t> >& edgeList );

  static MbVec3d getIsoParametricCoord(const MiGeometryI& meshGeometry, const MiVolumeCell* pyramidCell, const MbVec3d &point);
  
  static void getWeight(const MiGeometryI& meshGeometry, const MiVolumeCell* pyramidCell, const MbVec3d &point, std::vector<double>& weights)
  {
    getWeight(getIsoParametricCoord(meshGeometry,pyramidCell,point),weights);
  }

  static void getWeight(const MbVec3d &ipcoord, std::vector<double>& weights)
  {
    double rm, sm, tm;

    rm = 1. - ipcoord[0];
    sm = 1. - ipcoord[1];
    tm = 1. - ipcoord[2];

    weights[0] = rm * sm * tm;
    weights[1] = ipcoord[0] * sm * tm;
    weights[2] = ipcoord[0] * ipcoord[1] * tm;
    weights[3] = ipcoord[1] * rm * tm;
    weights[4] = ipcoord[2];
  }

  static void getDerivs(const MbVec3d &ipcoord, std::vector<double>& derivs);

  static bool isPointInsideCell(const MiGeometryI& meshGeometry, const MiVolumeCell* pyramidCell, const MbVec3d &point, std::vector<double>& weights)
  {
    size_t n0 = pyramidCell->getNodeIndex(0);
    size_t n1 = pyramidCell->getNodeIndex(1);
    size_t n2 = pyramidCell->getNodeIndex(2);
    size_t n3 = pyramidCell->getNodeIndex(3);
    size_t n4 = pyramidCell->getNodeIndex(4);
    s_halfPyramid.setNodeIds(n0,n1,n2,n4);
    bool inside = s_halfPyramid.isPointInsideCell(meshGeometry,point,weights);
    if (!inside)
    {
      s_halfPyramid.setNodeIds(n0,n2,n3,n4);
      inside = s_halfPyramid.isPointInsideCell(meshGeometry,point,weights);
    }
    if(inside)
      getWeight(meshGeometry,pyramidCell,point,weights);
    return inside;
  }

  static double getVolume(const MiGeometryI& meshGeometry, const MiVolumeCell* pyramidCell);

  static MbVec3d getFacetNormal(size_t facetId, const MiGeometryI& meshGeometry, const MiVolumeCell* pyramidCell);

  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 void initMarchingCase3();
  static void initMarchingCase4();
  static void initMarchingCase5();


  static bool                 s_isInit;
  static MxTabOfMarchingCase* s_tabOfMarchingPyramid;
  static MxHalfPyramidCell    s_halfPyramid;
};

#endif 

---

Open Inventor Toolkit reference manual, generated on 15 Mar 2023
Copyright © Thermo Fisher Scientific All rights reserved. http://www.openinventor.com/