scale the normal distribution in update() to avoid an extra multiplication
use that ak and -ak are equivalent random numbers to simplify dphi