BackfaceCullingStudyTest.C

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
 
#include "BackfaceCullingStudyTest.h"
 
// Local includes
#include "TraceRay.h"
 
registerMooseObject("RayTracingTestApp", BackfaceCullingStudyTest);
 
InputParameters
BackfaceCullingStudyTest::validParams()
{
  auto params = LotsOfRaysRayStudy::validParams();
  return params;
}
 
BackfaceCullingStudyTest::BackfaceCullingStudyTest(const InputParameters & parameters)
  : LotsOfRaysRayStudy(parameters)
{
  // Use backface culling
  for (const auto tid : make_range(libMesh::n_threads()))
    traceRay(tid).setBackfaceCulling(true);
}
 
void
BackfaceCullingStudyTest::initialSetup()
{
  LotsOfRaysRayStudy::initialSetup();
  generateNormals();
}
 
void
BackfaceCullingStudyTest::meshChanged()
{
  LotsOfRaysRayStudy::meshChanged();
  generateNormals();
}
 
void
BackfaceCullingStudyTest::generateNormals()
{
  _normals.clear();
 
  for (const Elem * elem : _mesh.getMesh().active_element_ptr_range())
  {
    std::vector<Point> normals(elem->n_sides());
    for (const auto s : elem->side_index_range())
      normals[s] = RayTracingStudy::getSideNormal(elem, s, 0);
    _normals.emplace(elem, std::move(normals));
  }
}
 
const Point *
BackfaceCullingStudyTest::getElemNormals(const Elem * elem, const THREAD_ID)
{
  const auto find = _normals.find(elem);
  if (find == _normals.end())
    mooseError("Failed to find normal for elem ", elem->id());
 
  const auto & normals = find->second;
  mooseAssert(normals.size() == elem->n_sides(), "Normals not sized properly");
  return normals.data();
}