在C++(犰狳)中存储大矩阵



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

假设我声明一个矩阵

mat X;

并将大小设置为(例如)500 行,500 列,带有随机条目:

X.randn(500,500);

犰狳是否动态存储X(即在堆上),尽管不使用newdelete?我问的原因是,似乎犰狳允许我将变量声明为:

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 之间,并且事先不知道。

犰狳是否动态存储 X(即在堆上),尽管没有 使用newdelete.?

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

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

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

之后,您应该使用

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

如果您已经知道此时的大小。在某些情况下,

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

可能也可以。

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

我最初错过了你对n大小的评论.正如 Mikhail 所说,处理 100000x100000 矩阵需要更加小心,而不仅仅是考虑实例化它们的方式。

最新更新