R：应对"double sum of three number product"的策略

Question

假设矩阵 A 有 5000 行和 6000 列，其中每个元素a_ij，还有另一个矩阵 S，它有 6000

行和 6000 列，它是一个对称矩阵，其中每个元素都是s_lk的，我想计算m_ij = sum(k=1..6000, l=1..6000)(a_il * a_jk * s_lk)，（即这个公式），应该是 5000 * 5000 矩阵。做这样的问题（多数乘积的双和）有什么好技巧？我们应该先简化它吗？

谢谢。

Answer 1

@nicola，尼古拉给了正确的安瑟

。 谢谢。

A %*% S %*% t(A)

在这里，我只是发布尼古拉的答案。如果不合适，请将其删除。