That would be great, thanks. Right now, I’m not able to tell for certain if this is an error in fast_sample_posterior_predictive (and you say you get a related error in sample_posterior_predictive, right?), in pm.summary() (which is really arviz.summary()), or in the way the two interact.
I’ll look forward to seeing your test case.
Best,
r