Residual computation isn't working. Trying the MATLAB implementation of LSQR, even though it isn't good at managing memory, to see if that works.