Preconditioned CG solver works!

This release also has an important bug fix when creating matrices from grid resolutions. Only grids with ni == nj == nk were working, because at one point it was written nj*nj instead of ni*nj