2

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

假设我声明了一个矩阵

mat X;

并将大小设置为(比如说)500 行,500 列随机条目:

X.randn(500,500);

犰狳是否X动态存储(即在堆上),尽管不使用newor 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 之间,并且事先不知道。

4

1 回答 1

5

犰狳是否动态存储 X(即在堆上),尽管不使用newor delete.?

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

我问的原因是,似乎犰狳允许我将变量声明为(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 矩阵比简单地考虑实例化它们的方式需要更多的关注。

于 2013-09-20T16:13:18.017 回答