In the documentation of the LSQR function, it says :
"lsqr(A,b,tol,maxit,M1) and lsqr(A,b,tol,maxit,M1,M2) use n-by-n preconditioner M or M = M1*M2 and effectively solve the system A*inv(M)*y= b for y, where x = M*y. "
Why x = M*y, whereas we are supposed to solve the system A*x = b, so could it be x = inv(M)*y ?
Best Answer