MiCell.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 _MiCell_H
#define _MiCell_H
 
#include <MeshVizXLM/mesh/geometry/MiGeometryI.h>
#include <MeshVizXLM/MxMeshVizXLM.h>
#include <MeshVizXLM/MiAbstractMethodError.h>
 
#include <MeshVizXLM/MbVec3.h>
#include <Inventor/STL/vector>
#include <Inventor/STL/ostream>
 
class MESHVIZXLM_EXTR_API MiCell
{
public:
 
  virtual ~MiCell() {}
 
  virtual size_t getNumNodes() const = 0; 
 
  virtual size_t getNodeIndex(size_t node) const = 0;
 
  virtual double getRelativeSize(const MiGeometryI* SO_UNUSED_PARAM(meshGeometry)) const 
  {
    return 1;
  }
 
  virtual void getWeight(const MbVec3d& SO_UNUSED_PARAM(ipcoord), std::vector<double>& SO_UNUSED_PARAM(weight)) const 
  {
    throw MiAbstractMethodError("MiCell::getWeight(const MbVec3d &, std::vector<double>&)");
  }
 
  virtual MbVec3d getIsoParametricCoord(size_t SO_UNUSED_PARAM(nodeIndex)) const 
  {
    throw MiAbstractMethodError("MiCell::getIsoParametricCoord(size_t nodeIndex)");
  }
 
  virtual MbVec3d getCenter(const MiGeometryI& geometry) const;
 
 
  virtual bool isPointInsideCell(const MiGeometryI& SO_UNUSED_PARAM(meshGeometry), 
                                 const MbVec3d& SO_UNUSED_PARAM(point), 
                                 std::vector<double>& SO_UNUSED_PARAM(weights)) const
  {
    throw MiAbstractMethodError("MiCell::isPointInsideCell(const MbVec3d &point)");
  }
 
  friend std::ostream& operator<< (std::ostream& s, const MiCell& cell);
 
private: protected:
  virtual std::ostream& toStream(std::ostream& s) const;
 
};
 
//-----------------------------------------------------------------------------
inline MbVec3d
MiCell::getCenter(const MiGeometryI& geometry) const
{
  MbVec3d center(0);
  for (size_t i=0; i<getNumNodes(); ++i) 
    center += geometry.getCoord(getNodeIndex(i));
  return ( center / (double) getNumNodes() );
}
 
//-----------------------------------------------------------------------------
inline std::ostream& 
MiCell::toStream(std::ostream& s) const
{
  s << "# num cell's nodes" << std::endl;
  s << getNumNodes() << std::endl;
  s << "[";
  for (size_t i=0; i<getNumNodes(); ++i)
    s << getNodeIndex(i) << " ";
  s << "]";
 
  return s;
}
 
//-----------------------------------------------------------------------------
inline std::ostream& 
operator << (std::ostream& s, const MiCell& cell)
{
  return cell.toStream(s);
}
 
 
#endif