examples/smpi/NAS/BT/add.f

   1 c---------------------------------------------------------------------
   2 c---------------------------------------------------------------------
   3
   4       subroutine  add
   5
   6 c---------------------------------------------------------------------
   7 c---------------------------------------------------------------------
   8
   9 c---------------------------------------------------------------------
  10 c     addition of update to the vector u
  11 c---------------------------------------------------------------------
  12
  13       include 'header.h'
  14
  15       integer  c, i, j, k, m
  16
  17       do     c = 1, ncells
  18          do     k = start(3,c), cell_size(3,c)-end(3,c)-1
  19             do     j = start(2,c), cell_size(2,c)-end(2,c)-1
  20                do     i = start(1,c), cell_size(1,c)-end(1,c)-1
  21                   do    m = 1, 5
  22                      u(m,i,j,k,c) = u(m,i,j,k,c) + rhs(m,i,j,k,c)
  23                   enddo
  24                enddo
  25             enddo
  26          enddo
  27       enddo
  28
  29       return
  30       end