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

for the bayesian_logistic_regression.py 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:

ed.get_session().run(qw.variance())