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

Before Change

                v_new = v_new.copy()

            //     ++ orthogonalize
            hcur = []
            for v in vs:
                alpha = dot(v, v_new)
                hcur.append(alpha)
                v_new = axpy(v, v_new, v.shape[0], -alpha)  // v_new -= alpha*v

After Change

                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