---- the correct implementation is in pm_new.c, which makes use of the Rayleigh coefficient ----