使用高级构造函数将阵列转换为Armadillo Vec

Using advanced constructors to convert array to armadillo vec

本文关键字：转换 Armadillo Vec 阵列 高级 构造函数      更新时间：2023-10-16

我正在尝试将类型ODE_VECTOR的对象转换为Armadillo VEC，以便新对象将继承Armadillo库的线性代数功能。作为问题的一个简单示例，我正在使用C 数组的演示：

#include <iostream>
#include <random>
#include <armadillo>
using namespace std;
using namespace arma;
int main() {
double b1 = 0.2;
double b2 = 0.1;
const double state[2] = {b1, b2};
rowvec B(&state[0], 2);
cout << B << endl;
mat A(2,2); A.fill(0.2); A.diag().ones();
cout << A << endl;
cout << B(0) * A(0,0) << endl;    // x
return 0;
}

B1和B2中包含的数据保存在B的元素中，以使运行此脚本返回以下内容：

   0.2000
   0.1000
   1.0000   0.2000
   0.2000   1.0000
0.2

但是，如果我用

替换x 
cout << B * A << endl;

我遇到的错误表明b尚未继承标准的腋窝线性代数功能：

[ 50%] Building CXX object CMakeFiles/untitled.dir/main.cpp.o
[100%] Linking CXX executable untitled
CMakeFiles/untitled.dir/main.cpp.o: In function `void arma::blas::gemv<double>(char const*, int const*, int const*, double const*, double const*, int const*, double const*, int const*, double const*, double*, int const*)':
/usr/include/armadillo_bits/wrapper_blas.hpp:36: undefined reference to `wrapper_dgemv_'
CMakeFiles/untitled.dir/main.cpp.o: In function `void arma::blas::gemm<double>(char const*, char const*, int const*, int const*, int const*, double const*, double const*, int const*, double const*, int const*, double const*, double*, int const*)':
/usr/include/armadillo_bits/wrapper_blas.hpp:71: undefined reference to `wrapper_dgemm_'
collect2: error: ld returned 1 exit status
CMakeFiles/untitled.dir/build.make:94: recipe for target 'untitled' failed
make[3]: *** [untitled] Error 1
CMakeFiles/Makefile2:67: recipe for target 'CMakeFiles/untitled.dir/all' failed
make[2]: *** [CMakeFiles/untitled.dir/all] Error 2
CMakeFiles/Makefile2:79: recipe for target 'CMakeFiles/untitled.dir/rule' failed
make[1]: *** [CMakeFiles/untitled.dir/rule] Error 2
Makefile:118: recipe for target 'untitled' failed
make: *** [untitled] Error 2

任何人都可以解释如何正确使用高级构造函数(ARMA文档(来使此操作？

谢谢

您需要在乘法中具有正确的格式。B型是2x1，A为2x2。更改为

mat B(&state[0], 1,2);

或

rowvec B(&state[0], 2);
相关文章：