home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Professional / OS2PRO194.ISO / os2 / prgramer / adaptor / examples / overlap / jacobi.f < prev   
Encoding:
Text File  |  1993-06-28  |  1.7 KB  |  68 lines
  1.       PROGRAM LAPLACE
  2.       REAL F(:,:), DF(:,:), HF(:,:[1:1])
  3.       INTEGER MAXX, MAXY
  4.       REAL FMAX
  5.       INTEGER ITER
  6. c 
  7. c     read in sizes 
  8. c
  9.       PRINT *,'MAXX = (z.B. 64) '
  10.       READ *,MAXX
  11.       PRINT *,'MAXY = (z.B. 64) '
  12.       READ *,MAXY
  13.       ALLOCATE (F(MAXX,MAXY), DF(MAXX,MAXY))
  14.       ALLOCATE (HF(MAXX,MAXY))
  15. C
  16. !HPF$ INDEPENDENT, LOCAL_ACCESS
  17.       do j = 1, maxy
  18. !HPF$    INDEPENDENT, LOCAL_ACCESS
  19.          do i = 1, maxx 
  20.            if (j == maxy) then
  21.               F(i,j) = 1.0
  22.            else if (j == 1) then
  23.               F(i,j) = 2.0
  24.            else if ((i==1) .or. (i==maxx)) then
  25.               F(i,j) = 2.0
  26.            else
  27.               F(i,j) = 0.0
  28.            end if
  29.          end do
  30.       end do
  31.  
  32.       ITER = 0
  33.       FMAX = 1
  34.       DF = 0
  35.       DO WHILE (FMAX .gt. 0.001)
  36.          ITER = ITER + 1
  37. c        print *, 'iter = ', iter
  38. c        do i = 1, maxx
  39. c           do j = 1, maxy
  40. c              print *, 'f(',i,',',j,') = ',f(i,j)
  41. c           end do
  42. c        end do
  43.          hf = f          ! communicate with overlap
  44. c        do i = 1, maxx
  45. c           do j = 1, maxy
  46. c              print *, 'hf(',i,',',j,') = ', hf(i,j)
  47. c           end do
  48. c        end do
  49.          forall (j=2:maxy-1,i=2:maxx-1)
  50.              df(i,j) = (hf(i,j-1) + hf(i,j+1) + 
  51.      $                  hf(i-1,j) + hf(i+1,j) ) * 0.25  - hf(i,j)
  52.          end forall
  53. !HPF$    INDEPENDENT, LOCAL_ACCESS
  54.          do j = 2, maxy-1
  55. !HPF$      INDEPENDENT, LOCAL_ACCESS
  56.            do i = 2, maxx-1 
  57.              f(i,j) = f(i,j) + df(i,j)
  58.            end do
  59.          end do
  60.          DF = ABS(DF)
  61.          FMAX = MAXVAL (DF)
  62.          PRINT *,'Iteration ',ITER,'  Max = ',FMAX
  63.       END DO
  64.       PRINT *, ITER, ' Iterationen benoetigt'
  65.       deallocate (hf, df, f)
  66.       END 
  67.  
  68.