Skip to content

PUMAS v1.2

Latest
Compare
Choose a tag to compare
@niess niess released this 16 May 11:59

This release improves the sampling efficiency of backward collisions, for radiative processes. It also introduces a breaking change to the API.

Full details on PUMAS v1.2 physics are published in CPC (arXiv:2206.01457).

Physics update

  • The power-law exponent for sampling backward decays has been updated from 1.4 (previously) to 2.0. The new value fits radiative DCSs for fractional energy losses approximately in [1, 10]%. As a result, the Monte Carlo efficiency of backward collisions is improved.

  • For very long distances and small energy cutoffs, the current sampling method for backward collisions is numerically problematic. Therefore, for now PUMAS forbids using energy cutoffs values lower than 1% in backward mixed or straggled modes. The latter value has been checked to be numerically robust at least up to 10 km of standard rock. Note that in forward mode, arbitrary cutoffs values can still be used.

API update

  • The pumas_version function now has an additional patch argument. I.e., full version number uses the traditonal MAJOR.MINOR.PATCH semantic. However, GitHub releases are still limited to minor version changes.

Bug patches

  • Patch 1 corrects a memory alignment issue with some compilers / versions combinations when using agressive optimisation settings, typically -O3 with gcc. Note that this bug appeared in v1.1 already.

  • Patch 3 should slightly improve the accuracy for small steps, at high energy and in CSDA mode.