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

Before Change

                alpha = dot(v, v_new)
                hcur.append(alpha)
                v_new = axpy(v, v_new, v.shape[0], -alpha)  // v_new -= alpha*v
            hcur.append(nrm2(v_new))

            if hcur[-1] == 0:
                // Exact solution.
                break

After Change

                alpha = dot(v, v_new)
                hcur[i] = alpha
                v_new = axpy(v, v_new, v.shape[0], -alpha)  // v_new -= alpha*v
            hcur[-1] = nrm2(v_new)

            with np.errstate(over="ignore", divide="ignore"):
                // Careful with denormals
                alpha = 1/hcur[-1]

            if np.isfinite(alpha):
                v_new = scal(alpha, v_new)
            else:
                // v_new either zero (solution in span of previous

Instances