c--------------------------------------------------------------------- c--------------------------------------------------------------------- subroutine verify(no_time_steps, class, verified) c--------------------------------------------------------------------- c--------------------------------------------------------------------- c--------------------------------------------------------------------- c verification routine c--------------------------------------------------------------------- include 'header.h' include 'mpinpb.h' double precision xcrref(5),xceref(5),xcrdif(5),xcedif(5), > epsilon, xce(5), xcr(5), dtref integer m, no_time_steps character class logical verified c--------------------------------------------------------------------- c tolerance level c--------------------------------------------------------------------- epsilon = 1.0d-08 c--------------------------------------------------------------------- c compute the error norm and the residual norm, and exit if not printing c--------------------------------------------------------------------- call error_norm(xce) call copy_faces call rhs_norm(xcr) do m = 1, 5 xcr(m) = xcr(m) / dt enddo if (node .ne. 0) return class = 'U' verified = .true. do m = 1,5 xcrref(m) = 1.0 xceref(m) = 1.0 end do c--------------------------------------------------------------------- c reference data for 12X12X12 grids after 100 time steps, with DT = 1.50d-02 c--------------------------------------------------------------------- if ( (grid_points(1) .eq. 12 ) .and. > (grid_points(2) .eq. 12 ) .and. > (grid_points(3) .eq. 12 ) .and. > (no_time_steps .eq. 100 )) then class = 'S' dtref = 1.5d-2 c--------------------------------------------------------------------- c Reference values of RMS-norms of residual. c--------------------------------------------------------------------- xcrref(1) = 2.7470315451339479d-02 xcrref(2) = 1.0360746705285417d-02 xcrref(3) = 1.6235745065095532d-02 xcrref(4) = 1.5840557224455615d-02 xcrref(5) = 3.4849040609362460d-02 c--------------------------------------------------------------------- c Reference values of RMS-norms of solution error. c--------------------------------------------------------------------- xceref(1) = 2.7289258557377227d-05 xceref(2) = 1.0364446640837285d-05 xceref(3) = 1.6154798287166471d-05 xceref(4) = 1.5750704994480102d-05 xceref(5) = 3.4177666183390531d-05 c--------------------------------------------------------------------- c reference data for 36X36X36 grids after 400 time steps, with DT = 1.5d-03 c--------------------------------------------------------------------- elseif ( (grid_points(1) .eq. 36) .and. > (grid_points(2) .eq. 36) .and. > (grid_points(3) .eq. 36) .and. > (no_time_steps . eq. 400) ) then class = 'W' dtref = 1.5d-3 c--------------------------------------------------------------------- c Reference values of RMS-norms of residual. c--------------------------------------------------------------------- xcrref(1) = 0.1893253733584d-02 xcrref(2) = 0.1717075447775d-03 xcrref(3) = 0.2778153350936d-03 xcrref(4) = 0.2887475409984d-03 xcrref(5) = 0.3143611161242d-02 c--------------------------------------------------------------------- c Reference values of RMS-norms of solution error. c--------------------------------------------------------------------- xceref(1) = 0.7542088599534d-04 xceref(2) = 0.6512852253086d-05 xceref(3) = 0.1049092285688d-04 xceref(4) = 0.1128838671535d-04 xceref(5) = 0.1212845639773d-03 c--------------------------------------------------------------------- c reference data for 64X64X64 grids after 400 time steps, with DT = 1.5d-03 c--------------------------------------------------------------------- elseif ( (grid_points(1) .eq. 64) .and. > (grid_points(2) .eq. 64) .and. > (grid_points(3) .eq. 64) .and. > (no_time_steps . eq. 400) ) then class = 'A' dtref = 1.5d-3 c--------------------------------------------------------------------- c Reference values of RMS-norms of residual. c--------------------------------------------------------------------- xcrref(1) = 2.4799822399300195d0 xcrref(2) = 1.1276337964368832d0 xcrref(3) = 1.5028977888770491d0 xcrref(4) = 1.4217816211695179d0 xcrref(5) = 2.1292113035138280d0 c--------------------------------------------------------------------- c Reference values of RMS-norms of solution error. c--------------------------------------------------------------------- xceref(1) = 1.0900140297820550d-04 xceref(2) = 3.7343951769282091d-05 xceref(3) = 5.0092785406541633d-05 xceref(4) = 4.7671093939528255d-05 xceref(5) = 1.3621613399213001d-04 c--------------------------------------------------------------------- c reference data for 102X102X102 grids after 400 time steps, c with DT = 1.0d-03 c--------------------------------------------------------------------- elseif ( (grid_points(1) .eq. 102) .and. > (grid_points(2) .eq. 102) .and. > (grid_points(3) .eq. 102) .and. > (no_time_steps . eq. 400) ) then class = 'B' dtref = 1.0d-3 c--------------------------------------------------------------------- c Reference values of RMS-norms of residual. c--------------------------------------------------------------------- xcrref(1) = 0.6903293579998d+02 xcrref(2) = 0.3095134488084d+02 xcrref(3) = 0.4103336647017d+02 xcrref(4) = 0.3864769009604d+02 xcrref(5) = 0.5643482272596d+02 c--------------------------------------------------------------------- c Reference values of RMS-norms of solution error. c--------------------------------------------------------------------- xceref(1) = 0.9810006190188d-02 xceref(2) = 0.1022827905670d-02 xceref(3) = 0.1720597911692d-02 xceref(4) = 0.1694479428231d-02 xceref(5) = 0.1847456263981d-01 c--------------------------------------------------------------------- c reference data for 162X162X162 grids after 400 time steps, c with DT = 0.67d-03 c--------------------------------------------------------------------- elseif ( (grid_points(1) .eq. 162) .and. > (grid_points(2) .eq. 162) .and. > (grid_points(3) .eq. 162) .and. > (no_time_steps . eq. 400) ) then class = 'C' dtref = 0.67d-3 c--------------------------------------------------------------------- c Reference values of RMS-norms of residual. c--------------------------------------------------------------------- xcrref(1) = 0.5881691581829d+03 xcrref(2) = 0.2454417603569d+03 xcrref(3) = 0.3293829191851d+03 xcrref(4) = 0.3081924971891d+03 xcrref(5) = 0.4597223799176d+03 c--------------------------------------------------------------------- c Reference values of RMS-norms of solution error. c--------------------------------------------------------------------- xceref(1) = 0.2598120500183d+00 xceref(2) = 0.2590888922315d-01 xceref(3) = 0.5132886416320d-01 xceref(4) = 0.4806073419454d-01 xceref(5) = 0.5483377491301d+00 c--------------------------------------------------------------------- c reference data for 408X408X408 grids after 500 time steps, c with DT = 0.3d-03 c--------------------------------------------------------------------- elseif ( (grid_points(1) .eq. 408) .and. > (grid_points(2) .eq. 408) .and. > (grid_points(3) .eq. 408) .and. > (no_time_steps . eq. 500) ) then class = 'D' dtref = 0.30d-3 c--------------------------------------------------------------------- c Reference values of RMS-norms of residual. c--------------------------------------------------------------------- xcrref(1) = 0.1044696216887d+05 xcrref(2) = 0.3204427762578d+04 xcrref(3) = 0.4648680733032d+04 xcrref(4) = 0.4238923283697d+04 xcrref(5) = 0.7588412036136d+04 c--------------------------------------------------------------------- c Reference values of RMS-norms of solution error. c--------------------------------------------------------------------- xceref(1) = 0.5089471423669d+01 xceref(2) = 0.5323514855894d+00 xceref(3) = 0.1187051008971d+01 xceref(4) = 0.1083734951938d+01 xceref(5) = 0.1164108338568d+02 c--------------------------------------------------------------------- c reference data for 1020X1020X1020 grids after 500 time steps, c with DT = 0.1d-03 c--------------------------------------------------------------------- elseif ( (grid_points(1) .eq. 1020) .and. > (grid_points(2) .eq. 1020) .and. > (grid_points(3) .eq. 1020) .and. > (no_time_steps . eq. 500) ) then class = 'E' dtref = 0.10d-3 c--------------------------------------------------------------------- c Reference values of RMS-norms of residual. c--------------------------------------------------------------------- xcrref(1) = 0.6255387422609d+05 xcrref(2) = 0.1495317020012d+05 xcrref(3) = 0.2347595750586d+05 xcrref(4) = 0.2091099783534d+05 xcrref(5) = 0.4770412841218d+05 c--------------------------------------------------------------------- c Reference values of RMS-norms of solution error. c--------------------------------------------------------------------- xceref(1) = 0.6742735164909d+02 xceref(2) = 0.5390656036938d+01 xceref(3) = 0.1680647196477d+02 xceref(4) = 0.1536963126457d+02 xceref(5) = 0.1575330146156d+03 else verified = .false. endif c--------------------------------------------------------------------- c verification test for residuals if gridsize is one of c the defined grid sizes above (class .ne. 'U') c--------------------------------------------------------------------- c--------------------------------------------------------------------- c Compute the difference of solution values and the known reference values. c--------------------------------------------------------------------- do m = 1, 5 xcrdif(m) = dabs((xcr(m)-xcrref(m))/xcrref(m)) xcedif(m) = dabs((xce(m)-xceref(m))/xceref(m)) enddo c--------------------------------------------------------------------- c Output the comparison of computed results to known cases. c--------------------------------------------------------------------- if (class .ne. 'U') then write(*, 1990) class 1990 format(' Verification being performed for class ', a) write (*,2000) epsilon 2000 format(' accuracy setting for epsilon = ', E20.13) verified = (dabs(dt-dtref) .le. epsilon) if (.not.verified) then class = 'U' write (*,1000) dtref 1000 format(' DT does not match the reference value of ', > E15.8) endif else write(*, 1995) 1995 format(' Unknown class') endif if (class .ne. 'U') then write (*,2001) else write (*, 2005) endif 2001 format(' Comparison of RMS-norms of residual') 2005 format(' RMS-norms of residual') do m = 1, 5 if (class .eq. 'U') then write(*, 2015) m, xcr(m) else if (xcrdif(m) .le. epsilon) then write (*,2011) m,xcr(m),xcrref(m),xcrdif(m) else verified = .false. write (*,2010) m,xcr(m),xcrref(m),xcrdif(m) endif enddo if (class .ne. 'U') then write (*,2002) else write (*,2006) endif 2002 format(' Comparison of RMS-norms of solution error') 2006 format(' RMS-norms of solution error') do m = 1, 5 if (class .eq. 'U') then write(*, 2015) m, xce(m) else if (xcedif(m) .le. epsilon) then write (*,2011) m,xce(m),xceref(m),xcedif(m) else verified = .false. write (*,2010) m,xce(m),xceref(m),xcedif(m) endif enddo 2010 format(' FAILURE: ', i2, E20.13, E20.13, E20.13) 2011 format(' ', i2, E20.13, E20.13, E20.13) 2015 format(' ', i2, E20.13) if (class .eq. 'U') then write(*, 2022) write(*, 2023) 2022 format(' No reference values provided') 2023 format(' No verification performed') else if (verified) then write(*, 2020) 2020 format(' Verification Successful') else write(*, 2021) 2021 format(' Verification failed') endif return end