Fix conestack infinite loop

Fix an infinite loop that could occur in long simulations with a conestack with more than one layer. Previously, particles could get stuck on the surface between layers of the conestack.
nrc-cnrc · Sep 20, 2018 · f45ff00 · f45ff00
1 parent b4b6bc9
commit f45ff00
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/HEN_HOUSE/egs++/geometry/egs_cones/egs_cones.h b/HEN_HOUSE/egs++/geometry/egs_cones/egs_cones.h
@@ -1554,6 +1554,12 @@ class EGS_ConeStack : public EGS_BaseGeometry {
  t = tp; // set t = distance to plane
  }
 
+ // If we are right on a plane, return -1
+ if(tp < boundaryTolerance) {
+ t = 0;
+ return -1;
+ }
+
  // distance to outer cone
  int ir = ireg - il*nmax; // cone index in current layer
  bool hitc = false; // assume we don't hit the cone