7dc870df00c8fbb7fd8a73f70b18e24c41664188,osmnx/utils.py,,great_circle_vec,#Any#Any#Any#Any#Any#,272

Before Change

    theta1 = np.deg2rad(lng1)
    theta2 = np.deg2rad(lng2)

    cos = (np.sin(phi1) * np.sin(phi2) * np.cos(theta1 - theta2) + np.cos(phi1) * np.cos(phi2))

    // ignore warnings during this calculation because numpy warns it cannot
    // calculate arccos for self-loops since u==v
    with warnings.catch_warnings():
        warnings.simplefilter("ignore")
        arc = np.arccos(cos)

    // return distance in units of earth_radius
    distance = arc * earth_radius
    return distance

After Change

    h = np.sin(dphi/2)**2 + np.cos(phi1) * np.cos(phi2) * np.sin(dtheta/2)**2 
    h = np.minimum(1.0, h) // protection against floating point errors

    arc = 2 * np.arcsin(np.sqrt(h))

    // return distance in units of earth_radius
    distance = arc * earth_radius

Instances