I wonder whether this operation (in R)
t(A) %*% B %*% A
can be simplified via any algebraic trick (beyond crossprod(A,B)%*%B
crossprod(A,B)%*%B