16

我是使用犰狳的新手,尽管尝试/搜索了很多,但无法获得以下信息。

我需要对两个巨大的(动态)数组(不是向量)进行关联。我决定为此使用犰狳。我了解如何使用向量初始化 arma::mat,但我可以使用数组来这样做吗?我不明白,因为我在文档中没有看到任何提及。出于内部设计原因,我试图避免使用向量。我尝试使用示例数组手动初始化每个元素(作为一个愚蠢但起点)。像下面的代码这样的东西是行不通的。

using namespace std;
using namespace arma;   

mat A(SIZE, 1), B(SIZE, 1);

for(int i = 0; i < SIZE; i++)
{
    A << v[i] << endr;
    B << c[i] << endr;
}

cout << "A: " << endl;
A.print();
cout << "B: " << endl;
B.print();

对于输入数组 v = {1, 2, 0, -1, .9} 和 c = {0, .5, 1, -2, -5}。输出将是:

A:
        0
B:
  -5.0000

这是可以理解的。使用数组初始化 arma::mat 或 arma::colvector 的任何解决方法?提前致谢!

4

1 回答 1

10

假设您的数组 v 和 c 是双数组,您可以只使用辅助内存构造函数:

来自犰狳文档

  • mat(aux_mem*, n_rows, n_cols, copy_aux_mem = true, strict = true)

    使用可写辅助存储器中的数据创建矩阵。默认情况下,矩阵分配自己的内存并从辅助内存中复制数据(为了安全)。但是,如果 copy_aux_mem 设置为 false,矩阵将改为直接使用辅助内存(即不复制)。这更快,但除非您知道自己在做什么,否则可能很危险!

仅当 copy_aux_mem 设置为 false(即矩阵直接使用辅助内存)时,严格变量才会生效。如果 strict 设置为 true,则矩阵将在其生命周期内绑定到辅助内存;矩阵中元素的数量不能(直接或间接)改变。如果 strict 设置为 false,则矩阵在其生命周期内不会绑定到辅助存储器,即可以更改矩阵的大小。如果请求的元素数量与辅助内存的大小不同,则会分配新的内存并且不再使用辅助内存。

  • 垫子(常量 aux_mem*,n_rows,n_cols)

通过从只读辅助存储器复制数据来创建矩阵。

这意味着您可以通过复制源数据来创建矩阵,如下所示:

mat A_Copy(v, SIZE, 1);
mat B_Copy(c, SIZE, 1);

或者,您实际上可以重用已为数组分配的内存来创建只读矩阵,如下所示:

mat A_InPlace(v, SIZE, 1, /*copy_aux_mem*/false, /*strict*/true);
mat B_InPlace(c, SIZE, 1, /*copy_aux_mem*/false, /*strict*/true);

如果您使用向量,这将变得更加简单

vec A_Vec_Copy(v, SIZE);
vec B_Vec_Copy(c, SIZE);

或者:

vec A_Vec_InPlace(v, SIZE, false, true);
vec B_Vec_InPlace(c, SIZE, false, true);
于 2013-10-15T03:45:47.293 回答