我有一个简单的程序,它生成(使用Boost
)一些初始速度和位置,并计算传播一定距离所需的时间。根据横向距离(x, y)
,将最终的轴向 ( z
) 速度添加到向量中。这是一个简单的程序:
#include <iostream>
#include <boost/random.hpp>
#include <boost/random/normal_distribution.hpp>
using namespace std;
int main()
{
boost::mt19937 engine(static_cast<unsigned int>(time(0)));
boost::normal_distribution<double> nd(0.0, 1.0);
boost::variate_generator< boost::mt19937, boost::normal_distribution<double> > normal_std_one(engine, nd);
double coordX, coordY, coordZ, time;
double velX, velY, velZ;
const double factor = 0.01;
const double distance = 15.0;
vector<double> cont;
int i;
for(i=0; i<1000000000; i++)
{
coordX = factor*normal_std_one();
coordY = factor*normal_std_one();
coordZ = 0.0;
velX = normal_std_one();
velY = normal_std_one();
velZ = 20.0*normal_std_one()+300;
time = distance/velZ;
coordX += velX*time;
coordY += velY*time;
if(sqrt(coordX*coordX + coordY*coordY) < 0.02)
{
cont.push_back(velZ);
}
}
cout << cont.size() << endl;
return 0;
}
我认为一个很好的补充是for
使用OpenMP
. 我通过在循环启动之前添加以下行来做到这一点:
#pragma omp parallel for
此外,我已将-fopenmp
编译器选项和 `-fopenmp* 添加到链接器设置中。我的程序编译和链接没有错误,但是当我执行文件时,我收到消息:
Process terminated with status -1073741819 (0 minutes, 2 seconds)
我不清楚我在这里做错了什么。我正在使用 Windows 和 g++(通过 Code::Blocks IDE)。