def gauss_prob(mu, logstd, x): std = tf.exp(logstd) var = tf.square(std) gp = tf.exp(-(x - mu)/(2*var)) / ((2*np.pi)**.5 * std) return tf.reduce_prod(gp, [1])