home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Professional / OS2PRO194.ISO / os2 / prgramer / adaptor / examples / bsp2 / flow.f9 < prev    next >
Encoding:
Text File  |  1993-06-15  |  1.6 KB  |  55 lines
  1.       program flow
  2.       implicit none
  3.  
  4.       integer x_dim, y_dim, points
  5. c     dimension of image: x_dim x y_dim
  6.       parameter (x_dim = 128, y_dim = 128)
  7.       integer x_h, y_h
  8.       parameter (x_h = x_dim / 2, y_h = y_dim / 2)
  9.       parameter (points = 4096)
  10.       integer i, j
  11.  
  12.       real x_coord(points), y_coord(points), dx (points), dy (points)
  13.       integer image_buffer (x_dim, y_dim )
  14. cmf$  layout image_buffer (:host)
  15.       integer image (x_dim, y_dim)
  16.       integer x (points), y (points)
  17.       integer color (points)
  18.       integer dummy
  19.  
  20.       i = x_dim
  21.       j = y_dim
  22.       call x_display_init (x_dim,y_dim)
  23.  
  24.       x_coord = 0
  25.       y_coord = 0
  26.       do i = 1, points !parallel
  27.          x_coord(i)=sin(i*2*3.1415926*8/real(points))
  28.      $              *(real(i)/real(points))
  29.       end do
  30.       do i = 1,points
  31.          y_coord(i)=cos(i*2*3.1415926*8/real(points))
  32.      $              *(real(i)/real(points))
  33.       end do
  34.       x_coord = 0.5*x_coord
  35.       y_coord = 0.5*y_coord
  36.       do i = 1, 250
  37.         dx = x_coord + (i/25.0) * y_coord * x_coord
  38.         dy = - y_coord - (i/50.0) * x_coord * (y_coord-x_coord)
  39.         x_coord = x_coord + 0.01 * dx
  40.         y_coord = y_coord + 0.01 * dy
  41. c       display all points
  42.         image = 50
  43.         do j = 1, points   !parallel
  44.            x(j) = x_h + x_coord(j) * x_h
  45.            y(j) = y_h + y_coord(j) * y_h
  46.            color (j) = 255
  47.         end do
  48. c       scatter for a two dimensional structure
  49.         call global_send (image, x, y, color)
  50.         image_buffer = image
  51.         call x_show_bild (image_buffer)
  52.       end do
  53.       read *, dummy
  54.       end 
  55.