SelfShadowSideUserObject.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 "SideUserObject.h"
#include <vector>
#include "libmesh/data_type.h"
 
class SelfShadowSideUserObject : public SideUserObject
{
public:
  using SideIDType = std::pair<dof_id_type, unsigned int>;
 
  static InputParameters validParams();
 
  SelfShadowSideUserObject(const InputParameters & parameters);
 
  virtual void initialize() override;
  virtual void execute() override;
  virtual void finalize() override;
 
  virtual void threadJoin(const UserObject & y) override;
 
  // API to check if a QP is illuminated
  unsigned int illumination(const SideIDType & id) const;
 
  bool useDisplacedMesh() const { return getParam<bool>("use_displaced_mesh"); }
 
  using Triangle = std::tuple<Point, Point, Point, SideIDType>;
  using LineSegment = std::tuple<Point, Point, SideIDType>;
 
protected:
  const unsigned int _dim;
 
  std::vector<const PostprocessorValue *> _raw_direction;
 
  RealTensorValue _rotation;
 
  std::map<SideIDType, unsigned int> & _illumination_status;
 
  std::vector<Triangle> _triangles;
 
  std::vector<LineSegment> _lines;
 
  std::vector<std::pair<SideIDType, MooseArray<Point>>> _local_qps;
 
private:
  void addLines(const SideIDType & id);
  void addTriangles(const SideIDType & id);
  bool check2DIllumination(const Point & qp, const SideIDType & id);
  bool check3DIllumination(const Point & qp, const SideIDType & id);
 
  void rotate(MooseArray<Point> & points);
  void rotate(Triangle & triangle);
  void rotate(LineSegment & line);
};