Efficient way of computing matrix product AXA'?

问题

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 AXA^T, where A is an arbitrary n×n matrix and X is a symmetric n×n matrix?

来源：https://stackoverflow.com/questions/13292325/m-x-s-x-tm-multiplication-with-blas