在 C++ 中存储大矩阵（犰狳）

Question

我在 C++ 中使用 Armadillo 库来存储/计算大型矩阵。我的理解是应该动态地（在堆上）存储大型数组/矩阵。

假设我声明一个矩阵

mat X;

并将大小设置为（例如）500 行，500 列，随机条目：

X.randn(500,500);

尽管没有使用new 或delete，犰狳是否会动态存储X（即在堆上）？我问的原因是，似乎犰狳允许我将变量声明为：

mat::fixed<n_rows, n_cols> 

其中，我quote：“一般比动态内存分配要快，但是以后不能（直接或间接）改变矩阵的大小”。

不管上述情况——我应该使用这个吗：

mat A;
A.set_size(n-1,n-1);  

或者这个：

mat *A = new mat;
(*A).set_size(n-1,n-1);

其中 n 介于 1000 或 100000 之间，并且事先不知道。

Answer 1

犰狳是否动态存储 X（即在堆上），尽管没有 使用new 或delete.?

是的。库代码中会有某种形式的new 或delete。你只是没有从外面注意到它。

我问的原因是因为它似乎是犰狳 允许我将变量声明为 (mat::fixed ...)

您必须查看源代码才能了解这里到底发生了什么。我的猜测是它具有某种内部逻辑，可以根据大小决定如何处理事物。不过，对于小型矩阵，您通常会使用 mat::fixed。

接下来，你应该使用

mat A(n-1,n-1);

如果您当时已经知道尺寸。在某些情况下，

mat A;
A.set_size(n-1,n-1);  

也可能没问题。

我想不出将第二个选项与@​​987654331@ 指针一起使用的充分理由。首先，像犰狳这样的库在内部处理它们的内存分配，开发人员非常小心地把它弄好。此外，即使库中的内存代码被破坏，您的想法new mat 也无法修复它：您将为mat 对象分配内存，但该对象肯定相当小。很大一部分可能隐藏在 mat 类中的成员变量 T* data 之类的东西后面，您无法从外部影响它的分配方式。

我最初错过了您对n 大小的评论。正如 Mikhail 所说，处理 100000x100000 矩阵需要比简单地考虑实例化它们的方式更加小心。