Question about finding BFMI outputs from Stan results

Hi! I am working on diagnostics from stan results. I would like to check the sampling efficiency using Bayesian fraction missing information(BFMI). May I ask where to find the output of BFMI from stan results or I have to check it using arviz from python? Thanks!

GitHub - betanalpha/mcmc_diagnostics: Markov chain Monte Carlo general, and Hamiltonian Monte Carlo specific, diagnostics for Stan demonstrates how to access this information from RStan and PyStan.

is it also compatible with pystan 3?