Hi all,

I could use some help implementing a cross product in Theano. I’m porting over a model from Jackman and am stuck on finding the pymc3 equivalent of:

```
Sigma = crossprod(diag_pre_multiply(Sigma_scale, Sigma_corr_L))
```

I can implement in numpy easily enough but Theano doesn’t seem to have the equivalent. Many thanks!