OCCFace.h

Go to the documentation of this file.

// Gmsh - Copyright (C) 1997-2022 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file in the Gmsh root directory for license information.
// Please report all issues on https://gitlab.onelab.info/gmsh/gmsh/issues.
 
#ifndef OCC_FACE_H
#define OCC_FACE_H
 
#include "GmshConfig.h"
#include "GFace.h"
#include "GModel.h"
#include "Range.h"
 
#if defined(HAVE_OCC)
 
#include <BRepLProp_SLProps.hxx>
#include <GeomAPI_ProjectPointOnSurf.hxx>
#include <TopoDS_Face.hxx>
 
class OCCFace : public GFace {
private:
  TopoDS_Face _s;
  Handle(Geom_Surface) _occface;
  const BRepAdaptor_Surface _sf;
  double _umin, _umax, _vmin, _vmax;
  bool _periodic[2];
  double _period[2];
  double _radius;
  SPoint3 _center;
  void _setup();
  mutable GeomAPI_ProjectPointOnSurf _projector;
  bool _project(const double p[3], double uv[2], double xyz[3]) const;
 
public:
  OCCFace(GModel *m, TopoDS_Face s, int num);
  virtual ~OCCFace() {}
  virtual SBoundingBox3d bounds(bool fast = false);
  virtual Range<double> parBounds(int i) const;
  virtual GPoint point(double par1, double par2) const;
  virtual GPoint closestPoint(const SPoint3 &queryPoint,
                              const double initialGuess[2]) const;
  virtual bool containsPoint(const SPoint3 &pt) const;
  virtual bool containsParam(const SPoint2 &pt);
  virtual SVector3 normal(const SPoint2 &param) const;
  virtual Pair<SVector3, SVector3> firstDer(const SPoint2 &param) const;
  virtual void secondDer(const SPoint2 &, SVector3 &, SVector3 &,
                         SVector3 &) const;
  virtual GEntity::GeomType geomType() const;
  virtual ModelType getNativeType() const { return OpenCascadeModel; }
  virtual void *getNativePtr() const { return (void *)&_s; }
  virtual SPoint2 parFromPoint(const SPoint3 &, bool onSurface = true,
                               bool convTestXYZ = false) const;
  virtual double curvatureMax(const SPoint2 &param) const;
  virtual double curvatures(const SPoint2 &param, SVector3 &dirMax,
                            SVector3 &dirMin, double &curvMax,
                            double &curvMin) const;
  // tells if it's a sphere, and if it is, returns parameters
  virtual bool isSphere(double &radius, SPoint3 &center) const;
  virtual bool periodic(int dim) const { return _periodic[dim]; }
  virtual double period(int dim) const { return _period[dim]; }
  // save itself in BREP format (for debug e.g.)
  void writeBREP(const char *filename);
  bool buildSTLTriangulation(bool force = false);
};
 
#endif
 
#endif