matrix - Armadillo C++ : Is there a specific way for creating efficiently triangular or symmetric matrix

Question

我主要使用 Armadillo 来处理对称矩阵和三角形矩阵。我希望在内存存储方面保持高效。然而，似乎没有其他方法，只能创建一个新的垫子并用零(对于三角形)或重复项(对于对称)填充矩阵的下/上部分。

是否有更有效的方法使用 Armadillo 来使用三角形/对称矩阵？

谢谢，安托万

Answer 1

Armadillo 中没有对三角形或带状矩阵的特定支持。然而，从 3.4 版本开始，逐渐添加了对稀疏矩阵的支持。根据您需要的 Armadillo 函数以及矩阵的稀疏性，您可能会从使用SpMat<type>中获益。它实现了 compressed sparse column (CSC) format 。对于矩阵中的每个非零值，CSC 格式将行索引与值一起存储，因此您可能不会为三角矩阵节省太多内存。然而，带状对角矩阵消耗的内存应该少得多。