examples/smpi/NAS/LU/exact.f

   1 c---------------------------------------------------------------------
   2 c---------------------------------------------------------------------
   3
   4       subroutine exact( i, j, k, u000ijk )
   5
   6 c---------------------------------------------------------------------
   7 c---------------------------------------------------------------------
   8
   9 c---------------------------------------------------------------------
  10 c
  11 c   compute the exact solution at (i,j,k)
  12 c
  13 c---------------------------------------------------------------------
  14
  15       implicit none
  16
  17       include 'applu.incl'
  18
  19 c---------------------------------------------------------------------
  20 c  input parameters
  21 c---------------------------------------------------------------------
  22       integer i, j, k
  23       double precision u000ijk(*)
  24
  25 c---------------------------------------------------------------------
  26 c  local variables
  27 c---------------------------------------------------------------------
  28       integer m
  29       double precision xi, eta, zeta
  30
  31       xi  = ( dble ( i - 1 ) ) / ( nx0 - 1 )
  32       eta  = ( dble ( j - 1 ) ) / ( ny0 - 1 )
  33       zeta = ( dble ( k - 1 ) ) / ( nz - 1 )
  34
  35
  36       do m = 1, 5
  37          u000ijk(m) =  ce(m,1)
  38      >        + ce(m,2) * xi
  39      >        + ce(m,3) * eta
  40      >        + ce(m,4) * zeta
  41      >        + ce(m,5) * xi * xi
  42      >        + ce(m,6) * eta * eta
  43      >        + ce(m,7) * zeta * zeta
  44      >        + ce(m,8) * xi * xi * xi
  45      >        + ce(m,9) * eta * eta * eta
  46      >        + ce(m,10) * zeta * zeta * zeta
  47      >        + ce(m,11) * xi * xi * xi * xi
  48      >        + ce(m,12) * eta * eta * eta * eta
  49      >        + ce(m,13) * zeta * zeta * zeta * zeta
  50       end do
  51
  52       return
  53       end