examples/smpi/NAS/SP/add.f

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