Fri, 24 Jun 2011 10:49:35 +0000, Pauli Virtanen wrote: [clip] >> print ' system too large trying iteration' result = >> scipy.sparse.linalg.gmres(matrix,b)[0] Also, you can specify the required tolerance with the tol= parameter: often getting the solution down to machine precision is not needed.