f529deccde4a185a6c0f98f0968154f922abb86e,Lib/sandbox/pyem/densities.py,,_diag_gauss_den,#Any#Any#Any#Any#,132

Before Change

    if not log:
        inva = 1/va[0, 0]
        fac = (2*N.pi) ** (-d/2.0) * N.sqrt(inva)
        y =  (x[:, 0] - mu[0, 0]) ** 2 * inva * -0.5
        for i in range(1, d):
            inva = 1/va[0, i]
            fac *= N.sqrt(inva)
            y += (x[:, i] - mu[0, i]) ** 2 * inva * -0.5
        y = fac * N.exp(y)
    else:
        y = _scalar_gauss_den(x[:, 0], mu[0, 0], va[0, 0], log)
        for i in range(1, d):

After Change

        inva *= -0.5
        x = x - mu
        x **= 2
        y = fac * N.exp(N.dot(x, inva))
    else:
        // XXX optimize log case as non log case above
        y = _scalar_gauss_den(x[:, 0], mu[0, 0], va[0, 0], log)

Instances