Implement deterministic nontrivial function of random variables

Questions
4

Looking at the function torusmodel itself it doesnt look too difficult to rewrite it into theano - have you try import theano.tensor as tt and just replace np.* with tt.*?

Implementing something like tt.trapz() for numerical integration over a set domain?