dis-err.c

并行计算的演示程序

#include <omp.h>
#include <stdio.h>

/*

  for(i=4;i<100;i++){               
       a[i] = b[i-2] + 1;
       c[i] = b[i-1] + f[i];  
       b[i] = a[i-1] + 2;    
       d[i] = d[i+1] + b[i-1];     
  }
  
*/

#define Iter 100000

main()
{
  int i;

  int a[Iter],b[Iter],c[Iter],d[Iter],f[Iter];
  int a1[Iter],b1[Iter],c1[Iter],d1[Iter],f1[Iter];


  for(i=0;i<Iter;i++) a[i]=b[i]=c[i]=d[i]=f[i]=i;
  for(i=0;i<Iter;i++) a1[i]=b1[i]=c1[i]=d1[i]=f1[i]=i;
  

  for(i=4;i<Iter;i++){
     a1[i] = b1[i-2] + 1;
     c1[i] = b1[i-1] + f1[i];
     b1[i] = a1[i-1] + 2;
     d1[i] = d1[i+1] + b1[i-1];
  }

#pragma omp parallel for shared(a,b,c,d,f) private(i)   
  for(i=4;i<Iter;i++){ 
     a[i] = b[i-2] + 1;
     c[i] = b[i-1] + f[i];
     b[i] = a[i-1] + 2;
     d[i] = d[i+1] + b[i-1];
  }
  

for(i=4;i<Iter;i++) {

  if ( a[i]!=a1[i]) printf("a[%d] = %d , a1[%d] = %d\n",i,a[i],i,a1[i]);
  if ( b[i]!=b1[i]) printf("b[%d] = %d , b1[%d] = %d\n",i,b[i],i,b1[i]); 
  if ( c[i]!=c1[i]) printf("c[%d] = %d , c1[%d] = %d\n",i,c[i],i,c1[i]);
  if ( d[i]!=d1[i]) printf("d[%d] = %d , d1[%d] = %d\n\n\n",i,d[i],i,d1[i]);
  

}  
  
  
  
 return 0;  
}