SubdomainGrainIDGenerator.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "SBMSubdomainGeneratorBase.h"
#include "libmesh/bounding_box.h"
#include "SurfaceElement.h"
#include "SurfaceEdge2.h"
#include "SurfaceTri3.h"
#include "AdaptiveRayContainmentCheck.h"
 
class SubdomainGrainIDGenerator : public SBMSubdomainGeneratorBase
{
public:
  static InputParameters validParams();
  SubdomainGrainIDGenerator(const InputParameters & parameters);
 
  std::unique_ptr<libMesh::MeshBase> generate() override;
 
protected:
  struct SubdomainCheckerEntry
  {
    SubdomainID subdomain_id;
    BoundingBox bbox;
    AdaptiveRayContainmentCheck * checker;
  };
 
  std::unique_ptr<libMesh::MeshBase> & _boundary_mesh;
 
protected:
  std::unordered_map<SubdomainID, std::vector<Point>> _centroids_by_subdomain;
  std::unordered_map<SubdomainID, std::vector<std::unique_ptr<SurfaceElement>>>
      _boundary_elements_by_subdomain;
  std::unordered_map<SubdomainID, BoundingBox> _bboxes_by_subdomain;
 
  std::unordered_map<SubdomainID, std::unique_ptr<AdaptiveRayContainmentCheck>>
      _subdomain_id_checkers;
  std::vector<SubdomainCheckerEntry> _checker_entries;
 
  Point _ray_direction;
 
  Real _eps;
  unsigned int _leaf_max_size;
 
  std::unique_ptr<libMesh::MeshBase> _boundary_owned;
 
  void buildSubdomainGroupedData();
 
  void buildInOutTesters();
 
  static bool
  boundingBoxesIntersect(const BoundingBox & lhs, const BoundingBox & rhs, const Real tolerance);
};