meshOctreeLibOL.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 MESH_OCTREE_LIB_OL
#define MESH_OCTREE_LIB_OL
 
#include <vector>
#include <array>
#include <string>
#include "GPoint.h"
#include "GFace.h"
 
class MTriangle;
 
class SurfaceProjector {
public:
  SurfaceProjector() : gf(NULL), OctIdx(0) {}
  SurfaceProjector(GFace *gf); /* read triangles and quads from GFace */
  SurfaceProjector(SurfaceProjector const &) = delete;
  SurfaceProjector &operator=(SurfaceProjector const &) = delete;
  ~SurfaceProjector();
 
  bool initialize(GFace *gf, const std::vector<MTriangle *> &triangles, bool useCADStl = false);
 
  void clear();
 
  bool setAnalyticalProjection(GFace *gf);
 
  GPoint closestPoint(const double query[3], bool evalOnCAD = false,
                      bool projectOnCAD = false) const;
 
public:
  GFace *gf;
 
protected:
  std::vector<std::array<double, 3> > points;
  std::vector<std::array<int32_t, 3> > triangles;
  std::vector<std::array<std::array<double, 2>, 3> > triangle_uvs;
  std::vector<bool> triangle_no_uv_eval; /* no eval. at param. singularity */
  int64_t OctIdx; /* pointer to libOL octree (C structure) */
 
  /* For simple CAD shapes, we have analytical formula for projection */
  bool useAnalyticalFormula = false;
  GFace::GeomType analyticalShape = GFace::GeomType::Unknown;
  std::array<double, 10> analyticalParameters;
};
 
 
 
enum class libOLTypTag {  LolTypVer=1, LolTypEdg, LolTypTri, LolTypQad,
               LolTypTet,   LolTypPyr, LolTypPri, LolTypHex, LolNmbTyp };
 
class libOLwrapper {
 
  public:
    libOLwrapper(
        const std::vector<std::array<double,3> >& points,
        const std::vector<std::array<int32_t,2> >& edges,
        const std::vector<std::array<int32_t,3> >& triangles,
        const std::vector<std::array<int32_t,4> >& quads,
        const std::vector<std::array<int32_t,4> >& tetrahedra,
        const std::vector<std::array<int32_t,5> >& pyramids,
        const std::vector<std::array<int32_t,6> >& prisms,
        const std::vector<std::array<int32_t,8> >& hexahedra);
 
    ~libOLwrapper();
 
    int elementsInsideBoundingBox(libOLTypTag elementType,
        double* bboxMin, double* bboxMax,
        std::vector<int32_t>& elements);
 
  protected:
    std::vector<std::array<double,3> > points;
    std::vector<std::array<int32_t,2> > edges;
    std::vector<std::array<int32_t,3> > triangles;
    std::vector<std::array<int32_t,4> > quads;
    std::vector<std::array<int32_t,4> > tetrahedra;
    std::vector<std::array<int32_t,5> > pyramids;
    std::vector<std::array<int32_t,6> > prisms;
    std::vector<std::array<int32_t,8> > hexahedra;
    int64_t OctIdx; /* pointer to libOL octree (C structure) */
 
};
 
#endif