MCDC icon indicating copy to clipboard operation
MCDC copied to clipboard

"A particle is lost" recursive loop

Open spasmann opened this issue 1 year ago • 0 comments

Running the 2D C5G7, I've run into a "Particle is lost" bug using both Monte Carlo and iQMC. A particle being lost could be due to an error in the input deck created by the user. So I'm not as concerned that a particle was lost, but that the "Particle is Lost" statement repeats infinitely. Shouldn't it print the error, kill the particle, and move on?

When get_particle_cell() fails to find a cell_ID we set:

  • P["alive"]=False,
  • return a cell_ID = -1,
  • and expect the while True loop in get_particle_material() to break?

However, it gets stuck in this recursive loop where it repeats the "Particle is Lost" statement, until the output file exceeds my disk quote and the run fails.

Sample of output:

A particle is lost at ( 16.612528973121755 -3.824586145564854 10000000000.0 )
A particle is lost at ( 16.612528973121755 -3.824586145564854 10000000000.0 )
A particle is lost at ( 16.612528973121755 -3.824586145564854 10000000000.0 )
A particle is lost at ( 16.612528973121755 -3.824586145564854 10000000000.0 )

spasmann avatar Apr 02 '24 19:04 spasmann