Guest

Pi OMP Fortran

By: a guest on Apr 7th, 2010  |  syntax: Fortran  |  size: 0.41 KB  |  hits: 273  |  expires: Never
download  |  raw  |  embed  |  report abuse
Copied
  1. program pi_singlethread
  2.     implicit none    
  3.     integer limit
  4.     integer i
  5.     real :: pi=0
  6.    
  7.     limit = 20000000
  8.    
  9.     !$OMP PARALLEL      !inicio da paralelizacao do codigo
  10.    
  11.     !$OMP DO REDUCTION (+:pi)   !paralelizar o do
  12.     do i=0, limit
  13.         pi = pi + 4.0 / (4.0*i + 1.0);
  14.         pi = pi - 4.0 / (4.0*i + 3.0);
  15.     enddo
  16.     !$OMP END DO
  17.    
  18.     !$OMP END PARALLEL
  19.    
  20.     print*, pi
  21.    
  22. end