如何只计算八度矩阵乘积的对角线?

有没有一种用八度来计算和存储矩阵积对角线的方法?

基本上喜欢做: vector = diag(A*B);

我不在乎 A*B 除了对角线上的那些。矩阵大小大约是 80k x 12 和 12 x 80k ,所以即使我不在乎速度/额外的内存,它也不能放在RAM中。

奇怪的是,因为八度是一个巨大数据集的包,对角线非常重要,所以应该是可能的。