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 回答