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!