OpenMP 2: FORTRAN Specific

parameter (n=5000,l=100)

real*8,dimension(n,n):: a,b,c
real*8,dimension(n):: d

a=0.5d0
b=1.d0

do i=1,n
a(i,i)=real(i)
end do

do i=1,n
b(i,i)=sin(i*0.5)
end do

c=matmul(a,b)

d(1:l)=d(1:l)+c(1:l,1)*b(1,1:l)

write(*,'(5e13.6)') (d(i),i=5,25,5)
stop
end