之前的一篇文章给出了用OpenMP编程求PI的四种方法,算是OpenMP编程的初级实例,本文将继续探讨用OpenMP并行编程解决更多的问题。
关于求PI的实现方法可见我的前一篇博客:http://blog.csdn.net/u014768676/article/details/45065341
1.使用临界区,并行规约方式来实现两个向量点积的运算
首先是临界区的方法:
#include <stdio.h>
#include <omp.h>
int main(){
int i,n;
float a[100],b[100],result;
//some initializations
n=100;
result=0.0;
for(i=0;i<n;i++){
a[i]=i*1.0;
b[i]=i*2.0;
}
#pragma omp parallel for
for(i=0;i<n;i++){
#pragma omp critical
result=result+(a[i]*b[i]);
}
printf("%f\n",result);
}