我正在使用boost lambda库找到给定的超过一百万点的最近点。当从文本中读取这些百万点时,它运行良好。但是为了减少编译时间,我已经将这些数据存储在数组中,并且我尝试使代码从数组中读取数据。尽管没有编译错误,但是对于存储在元组每个元素内的垃圾值,结果是错误的。代码是:
namespace
{
typedef double coord_t;
typedef boost::tuple<coord_t,coord_t,coord_t,coord_t,coord_t,coord_t> point_t;
使用Boost库计算距离:
coord_t distance_sq(const point_t& a, const point_t& b)
{
boost::geometry::distance
coord_t x = a.get<0>() - b.get<0>();
coord_t y = a.get<1>() - b.get<1>();
coord_t z = a.get<2>() - b.get<2>();
return x*x + y*y + z*z;
}
从一百万个点中可以找到最近点的输入点:
double px, py, pz;
px=4; py=4; pz=4;//-0.1958 0.8703 0.09787
point_t point(px, py, pz);
}
从文本文件读取数据:
ifstream f("MiniRotorCraft_02.txt");
while(f)
{
coord_t xa,ya,za,ua,va,wa;
f >> xa >> ya >> za >> ua >> va >> wa;
points.push_back(boost::make_tuple(xa,ya,za,ua,va,wa));//Million points entry into boost tuple
}
当我这样做时,它运行得很好。但是我希望我的代码不要从文本文件中读取,而是从代码本身的数组集中读取。所以我写了:
for (unsigned int i = 0; i < B[llp].size(); i++)
{
xa[i]=get<0>(B[llp][i]);
ya[i]=get<1>(B[llp][i]);
za[i]=get<2>(B[llp][i]);
ua[i]=get<3>(B[llp][i]);
va[i]=get<4>(B[llp][i]);
wa[i]=get<5>(B[llp][i]);
}
for (unsigned int i = 0; i < B[llp].size(); i++)
{
coord_t xa[i],ya[i],za[i],ua[i],va[i],wa[i];
//Instead of taking the data from text file I am trying to take from array that are stored in the code.
points.push_back(boost::make_tuple(xa[i],ya[i],za[i],ua[i],va[i],wa[i]));
}
但是它不起作用。它返回一些任意的垃圾值。 怎么做到呢。任何建议都将非常有益。
答案 0 :(得分:1)
以下代码
for (unsigned int i = 0; i < B[llp].size(); i++)
{
coord_t xa[i],ya[i],za[i],ua[i],va[i],wa[i];
//Instead of taking the data from text file I am trying to take from array that are stored in the code.
points.push_back(boost::make_tuple(xa[i],ya[i],za[i],ua[i],va[i],wa[i]));
}
不是C ++标准,因为使用
coord_t xa[i],ya[i],za[i],ua[i],va[i],wa[i];
您要声明大小为coord_t
的{{1}}(又名double
)的6个C样式数组,其中i
是运行时已知值(循环)。
但是还有更多:如果编译器接受使用运行时值初始化C样式数组(其中一些接受此作为语言扩展),则不会对其进行初始化(因此使用未定义的值),并且在下一行中
i
您正在传递大小为points.push_back(boost::make_tuple(xa[i],ya[i],za[i],ua[i],va[i],wa[i]));
的已初始化C样式数组的索引i
的值(第i + 1-n个值)。
结论:您的程序是未定义的行为,并且垃圾绝对是可以理解的。
可疑:您确定要排队吗
i
您是否不想推送在上一个 coord_t xa[i],ya[i],za[i],ua[i],va[i],wa[i];
循环中获得的值?