d8c83117426d14a996d233ad0b43a1f3041b54bb,astroML/stats/_point_statistics.py,,fit_bivariate_normal,#,192

Before Change

    // There is a degeneracy.  Given our input, we set sigma1 > sigma2,
    //  so alpha must be positive
    alpha = abs(alpha)

    sigma1 = np.sqrt((0.5 * (sigma_x ** 2 + sigma_y ** 2)
                      + np.sqrt(0.25 * (sigma_x ** 2 - sigma_y ** 2) ** 2
                                + sigma_xy ** 4)))

After Change

    // There is a degeneracy.  Given our input, we set sigma1 > sigma2,
    //  so alpha must be positive
    if alpha < 0:
        alpha += 0.5 * np.pi

    sigma1 = np.sqrt((0.5 * (sigma_x ** 2 + sigma_y ** 2)
                      + np.sqrt(0.25 * (sigma_x ** 2 - sigma_y ** 2) ** 2
                                + sigma_xy ** 4)))
    sigma2 = np.sqrt((0.5 * (sigma_x ** 2 + sigma_y ** 2)

Instances