Advice for speeding up sampling for models with multiple evaluations of Cholesky factors

Maybe try the pm.GP with approximation? From a quick look these is nothing obvious to speed up :sweat_smile: