solve_lsqr works (but is too slow)