home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Professional / OS2PRO194.ISO / os2 / prgramer / adaptor / examples / laplace / graphic.f < prev    next >
Encoding:
Text File  |  1993-09-28  |  1.5 KB  |  62 lines
  1.       PROGRAM LAPLACE
  2.       REAL F(:,:), DF(:,:)
  3.       REAL HF(:,:)
  4.       INTEGER BILD (:,:), HOSTBILD(:,:)
  5. cmf$  layout hostbild (:host)
  6.       LOGICAL CMASK (:,:)
  7.       INTEGER MAXX, MAXY
  8.       REAL FMAX
  9.       INTEGER ITER, K
  10. c 
  11. c     read in sizes 
  12. c
  13.       PRINT *,'MAXX = (z.B. 64) '
  14.       READ *,MAXX
  15.       PRINT *,'MAXY = (z.B. 64) '
  16.       READ *,MAXY
  17.       ALLOCATE (F(MAXX,MAXY), DF(1:MAXX,MAXY))
  18.       ALLOCATE (HF(MAXX,1:MAXY))
  19.       ALLOCATE (CMASK(MAXX,1:MAXY))
  20.       allocate (bild(maxx,maxy), hostbild(maxx,maxy))
  21.       call x_display_init (maxx, maxy)
  22. C
  23.       CMASK = .FALSE.
  24.       CMASK (2:MAXX-1,2:MAXY-1) = .TRUE.
  25.       F = 2.
  26.       F(:,MAXY) = 1.
  27.       WHERE (CMASK)
  28.           F = 0.0
  29.       ENDWHERE
  30.       ITER = 0
  31.       FMAX = 1
  32.       DO WHILE (FMAX .gt. 0.001)
  33.         do k = 1, 10
  34.          ITER = ITER + 1
  35.          HF = CSHIFT (F,1,1)
  36.          DF = CSHIFT (F,1,-1)
  37.          HF = HF + DF
  38.          DF = CSHIFT (F,2,1)
  39.          HF = HF + DF
  40.          DF = CSHIFT (F,2,-1)
  41.          HF = (HF + DF) * 0.25
  42.          DF = 0
  43.          WHERE (CMASK)
  44.             DF = HF - F
  45.             F  = F + DF
  46.          ENDWHERE
  47.         end do
  48.         bild = f * 128.0
  49.         hostbild = bild
  50.         call x_show_bild (hostbild)
  51.         DF = ABS(DF)
  52.         FMAX = MAXVAL (DF)
  53.         PRINT *,'Iteration ',ITER,'  Max = ',FMAX
  54.       END DO
  55.       PRINT *, ITER, ' Iterationen benoetigt'
  56.       read *, iter
  57.       call x_display_exit ()
  58.       deallocate (hostbild, bild)
  59.       deallocate (cmask, hf, df, f)
  60.       END 
  61.  
  62.