OpenMP的效果 - 爱码网

今天简单地把#pragma omp parallel for private(j)语句加在一个练习时做的求解Laplace方程的有限差分程序上，4万个点的网格上，在P4超线程的CPU上，使得365秒减少到314秒，50秒的差距呀。

ps：然而在使用SOR方法时，4万个点的网格，单线程只用了14秒， OpenMP的效果，以后有时间再研究一下高斯赛德尔的SOR多线程版本，看有什么改进。

附带程序：

#include <omp.h>
2 OpenMP的效果

#include <iostream>
3 OpenMP的效果

#include <fstream>
4 OpenMP的效果

#include <cmath>
5 OpenMP的效果

#include <iomanip>
6 OpenMP的效果

using namespace std;
7 OpenMP的效果

//迭代终止误差
8 OpenMP的效果

const double eps=1e-7;
9 OpenMP的效果

//网格划分
10 OpenMP的效果

const int ndivx=200;
11 OpenMP的效果

const int ndivy=200;
12 OpenMP的效果

const int nNx=ndivx+1;
13 OpenMP的效果

const int nNy=ndivy+1;
14 OpenMP的效果

const double deltax=1.0/ndivx;
15 OpenMP的效果

const double deltay=1.0/ndivy;
16 OpenMP的效果

//常数
17

const double PI=3.14159265;
18 OpenMP的效果

const double omiga=1.902;//松弛因子
19 OpenMP的效果

int main(void)
24}