1 c---------------------------------------------------------------------
2 c---------------------------------------------------------------------
6 c---------------------------------------------------------------------
7 c---------------------------------------------------------------------
9 c---------------------------------------------------------------------
11 c compute the solution error
13 c---------------------------------------------------------------------
20 c---------------------------------------------------------------------
22 c---------------------------------------------------------------------
26 double precision u000ijk(5), dummy(5)
41 call exact( iglob, jglob, k, u000ijk )
43 tmp = ( u000ijk(m) - u(m,i,j,k) )
44 dummy(m) = dummy(m) + tmp ** 2
50 c---------------------------------------------------------------------
51 c compute the global sum of individual contributions to dot product.
52 c---------------------------------------------------------------------
53 call MPI_ALLREDUCE( dummy,
62 errnm(m) = sqrt ( errnm(m) / ( (nx0-2)*(ny0-2)*(nz0-2) ) )
66 c write (*,1002) ( errnm(m), m = 1, 5 )
69 1002 format (1x/1x,'RMS-norm of error in soln. to ',
70 > 'first pde = ',1pe12.5/,
71 > 1x,'RMS-norm of error in soln. to ',
72 > 'second pde = ',1pe12.5/,
73 > 1x,'RMS-norm of error in soln. to ',
74 > 'third pde = ',1pe12.5/,
75 > 1x,'RMS-norm of error in soln. to ',
76 > 'fourth pde = ',1pe12.5/,
77 > 1x,'RMS-norm of error in soln. to ',
78 > 'fifth pde = ',1pe12.5)