PointContainmentClassifier.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "SurfaceSide.h"
#include "MooseTypes.h"
 
#include "libmesh/bounding_box.h"
#include "libmesh/parallel.h"
 
#include <memory>
 
namespace libMesh
{
class MeshBase;
}
class SurfaceElementSet;
class AdaptiveRayContainmentCheck;
class TriangleManifold;
 
class PointContainmentClassifier
{
public:
  enum class Method
  {
    PCA_RAY,
    USER_SELECTED_RAY,
    FIXED_X_RAY
  };
 
  struct RayOptions
  {
    Point ray_direction;
    int leaf_max_size = 10;
    FileName obb_file_name = "";
    FileName ray_file_name = "";
    const libMesh::Parallel::Communicator * comm = nullptr; // debug output only
  };
 
  PointContainmentClassifier(libMesh::MeshBase & mesh,
                             const SurfaceElementSet * set,
                             Method method,
                             Real tolerance);
 
  PointContainmentClassifier(libMesh::MeshBase & mesh,
                             const SurfaceElementSet * set,
                             Method method,
                             Real tolerance,
                             const RayOptions & options);
 
  ~PointContainmentClassifier();
 
  SurfaceGeometry::SurfaceSide sideness(const Point & point) const;
 
  bool contains(const Point & point) const
  {
    return sideness(point) != SurfaceGeometry::SurfaceSide::OUTSIDE;
  }
 
  const libMesh::BoundingBox & boundingBox() const { return _bounding_box; }
 
  std::size_t numElements() const { return _num_elements; }
 
  Point rayDirection() const;
 
private:
  const Method _method;
 
  std::unique_ptr<AdaptiveRayContainmentCheck> _pca; 
  std::unique_ptr<TriangleManifold> _tri;            
 
  libMesh::BoundingBox _bounding_box;
  std::size_t _num_elements = 0;
};