Implemented block-tridiagonal solver, found out that anisotropic RW implementation is wrong
??