home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / adaptor.zip / adapt.zip / adaptor / examples / dalib / replic / rep4.f < prev   
Text File  |  1993-06-25  |  563b  |  31 lines
  1.       program replicate
  2.  
  3.       real a(10,10,10,10)
  4.       real ra(5,5,5,10), x
  5.       integer i, j, k, l, n
  6.  
  7. cmf$  layout ra(:replicated)
  8.  
  9.       forall (i=1:10,j=1:10,k=1:10,l=1:10)
  10.         a(i,j,k,l) = i + 2*j - k + 3*l
  11.       end forall
  12.  
  13.       ra = 0
  14.       ra = a(4:8,2:6,5:9,1:10)
  15.  
  16.       n = 0
  17.       do i = 4, 8
  18.         do j = 2, 6
  19.           do k = 5,9
  20.             do l = 1, 10
  21.                n = n + i + 2*j - k + 3*l
  22.             end do
  23.           end do
  24.         end do
  25.       end do
  26.  
  27.       x = sum (ra)
  28.       print *, 'Result is ', x, ' should be ', n
  29.  
  30.       end
  31.