stop rays when exiting model

openmc-dev · Dec 20, 2023 · f58d45e · f58d45e
1 parent 4adabf0
commit f58d45e
Showing 1 changed file with 7 additions and 0 deletions.
diff --git a/src/plot.cpp b/src/plot.cpp
@@ -1659,6 +1659,13 @@ void Ray::trace()
       if (stop_)
         return;
 
+      // There are no more intersections to process
+      // if we hit the edge of the model, so stop
+      // the particle in that case:
+      if (dist_.distance == INFTY || dist_.distance == INFINITY) {
+        return;
+      }
+
       // Advance particle, prepare for next intersection
       for (int lev = 0; lev < n_coord(); ++lev) {
         coord(lev).r += dist_.distance * coord(lev).u;