BEAMnrc error: IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG

[walleludvig]

Hi EGSnrc community.

I am implementing an X-ray tube with BEAMnrc. My model works well most of the time, however as I increase incident particle number above ~2M I start getting a warning which kills my simulation. However, if I run the same model with incident particle number = 2M (or less), the simulation completes successfully every time without any warnings.

The warning that kills my simulations is:
Note: The following floating-point exceptions are signalling: IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG

I read in a previous forum that this problem should have been resolved already, but maybe it is recurring for a different reason?

I am running a recently installed EGSnrc package on MacBook M1 (Ventura 13.6), with GNU Fortran 13.1.0 (Homebrew GCC 13.1.0).

Any help to resolving this issue would be greatly appreciated!