def GaussianCdfInverse(p, mu=<NUM_LIT:0>, sigma=<NUM_LIT:1>):
x = ROOT2 * erfinv(<NUM_LIT:2> * p - <NUM_LIT:1>)<EOL>return mu + x * sigma<EOL>
Evaluates the inverse CDF of the gaussian distribution.

    See http://en.wikipedia.org/wiki/Normal_distribution#Quantile_function  

    Args:
        p: float

        mu: mean parameter

        sigma: standard deviation parameter

    Returns:
        float