X
#(',,i->ii', s, s, d)
sum((#('ij,i->j', X, w)-y)^2)
log(trace(X))
exp(trace(X))
trace(#('ia,ij,jb->ab', X, B, X))
sum((A-#('ij,kj->ik', U, V))^2)
sum((X - #('ijk,ia,jb,kc ->abc', K, A, B, C))^2)
||(X - #('i,j,k->ijk', u, v, w))||^2
log(sum(exp(X)))
log(det(X))
#('ab,bc,cd,da->', X, B, X, B)