start tackling the diagonal covariance problem.
* add better starting value for xi
* add a lower bound to precision tauh
* print three convergence metrics