Improved output of solvers.

Improved the information printed out at every iteration of the iterative solvers by using format() strings.
Fixed a typo.
Use direct loop in diagonal inversion

It's more efficient to use a direct par_loop to calculate the inverse of the diagonal matrix here.