在 C++ 中存储大矩阵（Armadillo）

Question

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

假设我声明了一个矩阵

mat X;

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

X.randn(500,500);

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

mat::fixed<n_rows, n_cols> 

其中，我引用：“通常比动态内存分配更快，但之后无法（直接或间接）更改矩阵的大小”。

不管以上——我应该使用这个：

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

尽管不使用new或，犰狳是否动态存储 X（即在堆上）delete？

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

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

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

之后，你应该使用

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

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

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

也可能没问题。

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

我最初错过了您对n. 正如 Mikhail 所说，处理 100000x100000 矩阵需要更多的关注，而不是简单地考虑实例化它们的方式。