How to increase posterior sample count when fitting GP model with find_MAP()?

The docs need updating! I think what changed is now pm.sample_posterior_predictive will draw one sample for each element in the trace. What you can do now is repeat the result from find_MAP, and also subsample the results from NUTS.