to allow for sampling controls from a custom distribution (e.g. correlated noise) rather than the default: ```julia randn!(m.noise) lmul!(m.covar_ul, reshape(m.noise, (nu, :))) ```