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);
};