I'm currently using BLAS function DSYMM to compute Y = AX and then DGEMM for YA', but I'm wondering is there some more efficient way of computing the matrix product AXAT, where A is an arbitrary n×n matrix and X is a symmetric n×n matrix?
Y = AX
YA'
n×n