Open
Description
There's a scan-based implementation in the docs, but it's super cryptic. Scan is not necessary to implement this, we can just call the function in a loop. Still might be a nice place to start:
https://pytensor.readthedocs.io/en/latest/library/scan.html#using-shared-variables-gibbs-sampling