def gaussian_nll(x, mu, log_sigma): sigma_squared = T.exp(2*log_sigma) return ( lib.floatX(0.5*numpy.log(2*numpy.pi)) + (2*log_sigma) + ( ((x-mu)**2) / (2*sigma_squared) ) )