[ ] implement AIES and ESS methods, which work with gradient free approaches. (Could be a lot faster with my model).
[ ] Can the ßmethod be parallelized? It shares information between chains, but the proposals are computed independently, so I can imagine that parallelization is possible (although there may be numpyro specific constraints.)