¿Manera eficiente de computar el producto matricial AXA '?
Actualmente estoy utilizando la función DSYMM de BLAS para calcularY = AX
y luego DGEMM paraYA'
, pero me pregunto si hay alguna forma más eficiente de calcular el producto matricial AXA.T, donde A es un arbitrarion×n
matriz y X es una simetrian×n
¿matriz?