MCDC
MCDC copied to clipboard
"A particle is lost" recursive loop
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 Trueloop inget_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 )