Regularizing the posterior of a parameter with the KL divergence from a different posterior

Maybe one method could be to actually use pm.sample on the first step and then use the posterior samples obtained from the sampling as the prior in the next fit or sample or whatever you want to do? As far as I know, using multidimensional posteriors as priors is not completely straightforward but there is a nice discussion here: