IntersectionUtils.hpp

Go to the documentation of this file.

/*
 * IntersectionUtils.hpp
 *
 *  Created on: 6 Apr 2023
 *      Author: Valerio Formato
 */
 
#ifndef HERD_INTERSECTION_UTILS
#define HERD_INTERSECTION_UTILS
 
namespace Herd::Geom {
struct IntersectionResult {
  size_t volID;
  std::array<Point, 2> points;
};
 
template <class GeoParams>
[[nodiscard]] std::vector<IntersectionResult> Intersect(const Line &line, const GeoParams &geoParams) {
  std::vector<IntersectionResult> result;
 
  // First, check if the line intersects the geoparams bounding box
  if (const auto &bb_intersections = geoParams.BoundingBox().Intersections(line); !bb_intersections) {
    return result;
  }
 
  for (size_t id = 0; id < geoParams.NTotElements(); ++id) {
    if (auto intersections = geoParams.ElementShape(id).Intersections(line); intersections) {
      result.push_back(IntersectionResult{id, intersections.value()});
    }
  }
 
  return result;
}
} // namespace Herd::Geom
 
#endif // HERD_INTERSECTION_UTILS