examples/smpi/NAS/SP/exact_solution.f

   1
   2 c---------------------------------------------------------------------
   3 c---------------------------------------------------------------------
   4
   5        subroutine exact_solution(xi,eta,zeta,dtemp)
   6
   7 c---------------------------------------------------------------------
   8 c---------------------------------------------------------------------
   9
  10 c---------------------------------------------------------------------
  11 c this function returns the exact solution at point xi, eta, zeta
  12 c---------------------------------------------------------------------
  13
  14        include 'header.h'
  15
  16        double precision  xi, eta, zeta, dtemp(5)
  17        integer m
  18
  19        do  m = 1, 5
  20           dtemp(m) =  ce(m,1) +
  21      >    xi*(ce(m,2) + xi*(ce(m,5) + xi*(ce(m,8) + xi*ce(m,11)))) +
  22      >    eta*(ce(m,3) + eta*(ce(m,6) + eta*(ce(m,9) + eta*ce(m,12))))+
  23      >    zeta*(ce(m,4) + zeta*(ce(m,7) + zeta*(ce(m,10) +
  24      >    zeta*ce(m,13))))
  25        end do
  26
  27        return
  28        end
  29
  30