Accord.Math.Decompositions.EigenvalueDecompositionF.cdiv C# (CSharp) Метод

        private static void cdiv(Single xr, Single xi, Single yr, Single yi,
            out Single cdivr, out Single cdivi)
        {
            // Complex scalar division.
            Single r;
            Single d;
            if (System.Math.Abs(yr) > System.Math.Abs(yi))
            {
                r = yi / yr;
                d = yr + r * yi;
                cdivr = (xr + r * xi) / d;
                cdivi = (xi - r * xr) / d;
            }
            else
            {
                r = yr / yi;
                d = yi + r * yr;
                cdivr = (r * xr + xi) / d;
                cdivi = (r * xi - xr) / d;
            }
        }