def gauss_log_prob(mu, logstd, x): var = tf.exp(2*logstd) gp = -tf.square(x - mu)/(2*var) - .5*tf.log(tf.constant(2*np.pi)) - logstd return tf.reduce_sum(gp, [1])