Efficient AXA 'matrix product calculation?

Question

I am currently using the BLAS DSYMM function to calculate Y = AX and then DGEMM for YA' , but I wonder if there is an even more efficient way to calculate the matrix product AXA ^T where A is an arbitrary n×n matrix and X is a symmetric matrix n×n ?