OpenMP (3) - Fortran

c$omp paralleldo schedule(static)

c$omp&shared(omega,error,uold,u)

c$omp&private(i,j,resid)

c$omp&reduction(+:error)

do j = 2,m-1

do i = 2,n-1

resid = calcerror(uold,I,j)

u(i,j) = uold(i,j) - omega * resid

error = error + resid*resid

end do

enddo

c$omp end paralleldo

Previous slide

Next slide

Back to first slide

View graphic version