TraceData.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

// Local includes
#include "Ray.h"

// MOOSE includes
#include "MooseError.h"

// libMesh includes
#include "libmesh/point.h"

struct TracePointData
{
  TracePointData(const libMesh::Point & point) : _point(point)
  {
    mooseAssert(_point != RayTracingCommon::invalid_point, "Invalid point");
  }

  libMesh::Point _point;
  std::vector<RayData> _data;
  std::vector<RayData> _aux_data;
};

struct TraceData
{
  TraceData(const std::shared_ptr<Ray> & ray)
    : _ray_id(ray->id()),
      _intersections(ray->intersections()),
      _processor_crossings(ray->processorCrossings()),
      _trajectory_changes(ray->trajectoryChanges()),
      _last(false)
  {
    mooseAssert(_ray_id != Ray::INVALID_RAY_ID, "Invalid Ray ID");
    addPoint(ray->currentPoint());
  }

  void addPoint(const libMesh::Point & point) { _point_data.emplace_back(point); }

  TracePointData & lastPoint() { return _point_data.back(); }

  bool stationary() const { return _last && _intersections == 0 && _point_data.size() == 1; }

  unsigned int numSegments() const
  {
    mooseAssert(!_point_data.empty(), "Should not be empty");
    return _point_data.size() - 1;
  }

  const RayID _ray_id;
  const unsigned long int _intersections;
  const unsigned int _processor_crossings;
  const unsigned int _trajectory_changes;
  bool _last;
  std::vector<TracePointData> _point_data;
};