C++犰狳:从矢量初始化>破坏数据
问题描述:
我试图初始化矢量矢量到犰狳垫的数据。C++犰狳:从矢量初始化<vector <double>>破坏数据
发现这个网站解释如何做到这一点:
http://www.kaiyin.co.vu/2014/07/initialize-armadillo-matrix-with-std.html
打印向量的向量和所产生的垫我得出的结论是犰狳使我的数据的垃圾之后。
我的函数的主要部分:
vector<vector<double>> C_mat(num_of_functions,vector<double>(num_of_points));
for(int j=0; j< num_of_functions; ++j)
{
for(int i=0;i<num_of_points; ++i)
{
C_mat[j][i]=(functions[j](points[i].x()));// shouldn't bother you, this works.
}
}
mat C(&(C_mat.front()).front(),num_of_points, num_of_functions);
cout << C << endl << endl;
for(int i=0; i< num_of_functions;++i)
{
print_vec(C_mat[i]);
}
print_ve(矢量)是一个功能我写的打印的向量。
输出:
1.0000e+00 0e+00 1.9965e+01
1.0000e+00 4.7924e-322 1.2822e+00
1.0000e+00 1.1683e+01 0e+00
1.0000e+00 1.6936e+01 4.7924e-322
1.0000e+00 2.3361e-01 1.0237e+02
1.0000e+00 1.6445e+01 2.1512e+02
1.0000e+00 7.4271e+00 4.0931e-02
1.0000e+00 1.4162e+01 2.0284e+02
1.0000e+00 1.1670e+01 4.1371e+01
1.0000e+00 2.3633e+00 1.5042e+02
印刷矢量:
印刷矢量:
11.6828 16.9359 0.233613 16.4455 7.42708 14.1619 11.6701 2.36329 19.9653 1.28223
印刷矢量:
102.366 215.119 0.0409313 202.84 41.3711 150.421 102.143 4.18885 298.959 1.23308
UPDATE:
我试图从根本上让我自己的2D矢量出一维向量相对于改变代码的第一条评论。现在犰狳不破坏数据,但加扰它:
vector<double> C_mat(num_of_functions*num_of_points);
for(int i=0;i<num_of_points; ++i)
{
for(int j=0; j< num_of_functions; ++j)
{
C_mat[i*num_of_functions+j]= functions[j](points[i].x());
}
}
mat C(&C_mat.front(),num_of_points, num_of_functions);
输出:
1.0000e+00 1.1170e+01 1.6853e+01
4.5538e+00 9.3574e+01 1.0000e+00
1.5553e+01 1.0000e+00 1.6292e+01
1.0000e+00 8.7956e-01 1.9907e+02
6.0653e+00 5.8021e-01 1.0000e+00
2.7591e+01 1.0000e+00 1.0787e+01
1.0000e+00 1.8169e+01 8.7269e+01
1.3849e+01 2.4758e+02 1.0000e+00
1.4385e+02 1.0000e+00 1.7998e+01
1.0000e+00 4.7403e+00 2.4295e+02
1 4.5538 15.5528
1 6.06532 27.5911
1 13.8494 143.854
1 11.1698 93.574
1 0.879557 0.580215
1 18.1687 247.577
1 4.74031 16.8529
1 16.2919 199.069
1 10.787 87.2695
1 17.998 242.945
答
的std::vector
被equivalient数据指针和一个长度(和一些不相关的内容)。指针指向存储实际元素的一堆内存。第0个元素存储在地址pointer
,第N个元素存储在pointer+N
。这就是为什么犰狳垫构造函数可以将第一个元素和一个长度的地址和构造一个矩阵。
当你做一个向量时,你会得到num_of_functions
内部向量,每个向量都有自己的数据指针。这意味着数据不会存储在连续的内存中。
要初始化std::vector
的二维矩阵,您应该创建大小为num_of_functions*num_of_points
的矩阵,并按列主次序存储这些值。
犰狳不争抢任何东西。您的原始数据按行排列。为了与LAPACK兼容,Armadillo按列主要顺序存储数据。您传递了一个指向原始数据的指针,但未考虑到行和列已互换。这可能会更好:'mat C(&C_mat.front(),num_of_functions,num_of_points,false);'。要转置矩阵,请使用[.t()](http://arma.sourceforge.net/docs.html#t_st_members)或[trans()](http://arma.sourceforge.net/docs.html#反式)功能。 – hbrerkere
谢谢,解决了我的问题。 –