MiVolumeCell.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-2023 BY FEI S.A.S,                        ***
 ***                        BORDEAUX, FRANCE                                        ***
 ***                      ALL RIGHTS RESERVED                                       ***
**=======================================================================*/
/*=======================================================================
** Author      : VSG (MMM YYYY)
**=======================================================================*/
#ifndef _MiVolumeCell_H
#define _MiVolumeCell_H
 
#ifdef _WIN32
#  pragma warning( push )
#  pragma warning(disable:4250)
#endif
 
#include <MeshVizXLM/mesh/cell/MiCell.h>
#include <MeshVizXLM/mesh/data/MiDataSet.h>
 
#include <Inventor/STL/vector>
#include <Inventor/STL/utility>
 
class MESHVIZXLM_EXTR_API MiVolumeCell : public MiCell
{
public:
 
  virtual size_t getNumFacets() const = 0; 
 
  virtual size_t appendNodesIndexOfFacet (size_t facet, std::vector<size_t>& facetNodes) const = 0;
 
  virtual size_t getSubTriangleNodesIndex(size_t SO_UNUSED_PARAM(facet), std::vector<size_t>& SO_UNUSED_PARAM(triangleNodeIds)) const 
  {
    return 0;
  }
 
  virtual size_t getSubTetrahedronNodesIndex (std::vector<size_t>& SO_UNUSED_PARAM(tetrahedronNodeIds)) const 
  {
    return 0;
  }
 
  virtual size_t appendLocalNodesIndexOfFacet (size_t SO_UNUSED_PARAM(facet), std::vector<size_t>& SO_UNUSED_PARAM(facetNodes)) const
  {
    throw MiAbstractMethodError("MiVolumeCell::appendLocalNodesIndexOfFacet(size_t facet, std::vector<size_t>& facetNodes)");
  }
 
  virtual void getIsosurfTopology(unsigned char caseId, std::vector< std::pair<size_t,size_t> >& edgeList ) const = 0;
 
  virtual unsigned char getMarchingCaseId(std::vector<bool> &nodesSign, size_t beginNodeId=0) const;
 
  virtual MbVec3d getFacetCenter(size_t SO_UNUSED_PARAM(facet), const MiGeometryI& SO_UNUSED_PARAM(geometry)) const;
 
private: protected:
  virtual std::ostream& toStream(std::ostream& s) const;
 
};
 
//-----------------------------------------------------------------------------
inline unsigned char 
MiVolumeCell::getMarchingCaseId(std::vector<bool> &nodesSign, size_t beginNodeId) const
{
  unsigned char caseId = 0;
  unsigned char maxId = 0;
  size_t i;
  for (i=0; i<getNumNodes(); ++i)
  {
    if (nodesSign[getNodeIndex(i)-beginNodeId]) caseId |= (1<<i);
    maxId |= (1<<i);
  }
  if (caseId == maxId) 
    caseId = 0;
  return caseId;
}
 
//-----------------------------------------------------------------------------
inline MbVec3d
MiVolumeCell::getFacetCenter(size_t facet, const MiGeometryI& geometry) const
{
  std::vector<size_t> facetNodes;
  MbVec3d facetCenter(0);
  size_t numFacetNodes = this->appendNodesIndexOfFacet(facet,facetNodes);
  for ( size_t n = 0; n < numFacetNodes; ++n) 
    facetCenter += geometry.getCoord(facetNodes[n]);
  return ( facetCenter / (double) numFacetNodes );
}
 
//-----------------------------------------------------------------------------
inline std::ostream& 
MiVolumeCell::toStream(std::ostream& s) const
{
  s << "# num cell's facets" << std::endl;
  s << getNumFacets() << std::endl;
 
  s << "# num cell's nodes" << std::endl;
  s << getNumNodes();
 
  std::vector<size_t> facetNodes;
  for (size_t i=0; i<getNumFacets(); ++i)
  {
    facetNodes.clear();
    s << std::endl << "[";
    appendNodesIndexOfFacet(i,facetNodes);
    for (size_t j=0; j<facetNodes.size(); ++j)
      s << facetNodes[j] << " ";
    s << "]";
  }
 
  return s;
}
 
 
 
#ifdef _WIN32
#  pragma warning( pop )
#endif
 
#endif