How to do VI or ADVI in pystan?

How to do Variational Inference (VI) or Automatic Differentiation Variational Inference (ADVI) with pystan? I would greatly appreciate a simple code example.

They are not currently implemented with PyStan.

Old pystan has them (see pystan2).

I recommend using CmdStanPy for these.