Margianal likelihood of Gaussian process by using DensityDist

You are not calling the tt.nlinalg.Det correctly, it should be tt.nlinalg.Det()(cov) - you initialized the theano_op first by doing tt.nlinalg.Det().