How to print out the variance of the inferred variational parameter?

for the example, I modified the data generation procedure by setting an arbitrary value for w and b (instead of using the np.tanh function).

After HMC inferencing, I can print out the mean by print(qw.eval(session=sess)). How do I print out the variance of qw and qb?

You can fetch the attributes of the underlying distribution: