我在C++中使用犰狳库来存储/计算大型矩阵。我的理解是应该动态存储大型数组/矩阵(在堆上)。
假设我声明一个矩阵
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 之间,并且事先不知道。
犰狳是否动态存储 X(即在堆上),尽管没有 使用
new
或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
对象分配内存,但该对象肯定相当小。大部分可能隐藏在类mat
中的成员变量T* data
后面,您无法从外部影响其分配方式。
我最初错过了你对n
大小的评论.正如 Mikhail 所说,处理 100000x100000 矩阵需要更加小心,而不仅仅是考虑实例化它们的方式。