How to specify this potential?

Hi @junpenglao! Thanks for answering.

See my code below:

var_alpha = theano.shared(value = 1., borrow = False)

phi_e = pm.Uniform(‘phi_e’, lower = lb_phi_e, upper = ub_phi_e, shape = (ub_phi_e.size))

S1_phi_e = (phi_e/phi_e.sum()*np.log(phi_e/phi_e.sum())).sum()

pm.DensityDist(‘phi_e_potential’, logp = -var_alpha*S1_phi_e)


The error I receive with this approach is “TypeError: ‘TensorVariable’ object is not callable”