tlunttil

Hi Elise, have you noticed that you are using two-dimensional models while the papers use 3D? In Jacques (2014) the geometry appears to be genuinely three-dimensional. In Thompson et al....

PRNG initialisation uses `time()`, which has granularity of a second. A better way would be to use random_device (available since C++11). Correlation between runs started a couple of seconds apart...

Similarly, there doesn't seem to be any checks to prevent divisions by zero when calculating the propagation length in case there's a mus=0 medium somewhere.