c7308c4ab8870e6506ee13801da4fe88f1a43004,scipy/sparse/linalg/isolve/lgmres.py,,lgmres,#Any#Any#Any#Any#Any#Any#Any#Any#Any#Any#Any#,14

Before Change

        // -- Get the LSQ problem solution
        y, info = trtrs(R[:j,:j], Q[0,:j].conj())
        if info != 0:
            // Zero diagonal -> exact solution, but we handled that above
            raise RuntimeError("QR solution failed")
        y *= inner_res_0

        if not np.isfinite(y).all():
            // Floating point over/underflow, non-finite result from

After Change

            if inner_res <= tol * inner_res_0 or breakdown:
                break

        if not np.isfinite(R[j-1,j-1]):
            // nans encountered, bail out
            return postprocess(x), k_outer + 1

        // -- Get the LSQ problem solution
        //
        // The problem is triangular, but the condition number may be
        // bad (or in case of breakdown the last diagonal entry may be
        // zero), so use lstsq instead of trtrs.
        y, _, _, _, = lstsq(R[:j,:j], Q[0,:j].conj(), lapack_driver="gelsy")
        y *= inner_res_0

        if not np.isfinite(y).all():

Instances