D对角线的XDX ^ T的快速矩阵乘法

Question

考虑XDX ^ T对于X和n矩阵的快速矩阵乘法，以及D和m乘以m对角矩阵。这里m >> n（假设n大约1000，m大约100000）。在我的应用程序中，X是一个固定矩阵，D的值可以在每次迭代时改变。

什么是快速的计算方法？目前我只是在C ++中进行简单的乘法运算。

编辑：我应该澄清我目前的程序，它不是“简单的乘法”。特别是，我将列乘以X乘以D的对角线条目的平方根得到A：= XD ^ {1/2}。然后我直接计算A * t（A）（这是n乘m矩阵与其转置的乘法）。

谢谢。

Answer 1

如果您知道D是对角线，那么您可以进行简单的乘法运算。希望你没有乘以零。