examples/smpi/NAS/BT/exact_solution.f

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