afd38df1eff037f0d27168320616533dc1ab189c,GPy/util/univariate_Gaussian.py,,std_norm_cdf,#Any#,11

Before Change

    return_val = 0.5*(1.0 + sign*erf);
    
    x = float(x)
    return weave.inline(code,arg_names=["x"],support_code=support_code)

def inv_std_norm_cdf(x):
    

After Change

    Based on Abramowitz, M. and Stegun, I. (1970)
    
    //Generalize for many x
    x = np.asarray(x).copy()
    cdf_x = np.zeros_like(x)
    N = x.size
    support_code = "//include <math.h>"
    code = 

    double sign, t, erf;

Instances