How could we get log p(z,x) (joint) and log q(z) when using autoguide.AutoDiagonalNormal()? We could not find where to get these?

We would like to compute the ELBO manually. How does Pyro compute:

ELBO = -svi.step(x,y)

where svi = pyro.infer.SVI(model, guide, optimizer, loss=JitTraceGraph_ELBO() )?

Could you please help? Thank you!