Matrix 计算矩阵积AXA'；的有效方法；？

我目前正在使用BLAS函数DSYMM来计算

Y=AX

，然后使用DGEMM来计算

YA'

，但我想知道是否有更有效的方法来计算矩阵积AXAT，其中A是任意

n×n

矩阵，X是对称

n×n

矩阵？

计算中是否有重复？也就是说，这两个项是否至少在几个迭代中保持不变？什么尺寸的矩阵？