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/11139933/efficient-way-of-computing-matrix-product-axa