1 c---------------------------------------------------------------------
2 c---------------------------------------------------------------------
4 subroutine exact( i, j, k, u000ijk )
6 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
9 c---------------------------------------------------------------------
11 c compute the exact solution at (i,j,k)
13 c---------------------------------------------------------------------
19 c---------------------------------------------------------------------
21 c---------------------------------------------------------------------
23 double precision u000ijk(*)
25 c---------------------------------------------------------------------
27 c---------------------------------------------------------------------
29 double precision xi, eta, zeta
31 xi = ( dble ( i - 1 ) ) / ( nx0 - 1 )
32 eta = ( dble ( j - 1 ) ) / ( ny0 - 1 )
33 zeta = ( dble ( k - 1 ) ) / ( nz - 1 )
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