diff --git a/src/ALE/Recon1d_PPM_CW.F90 b/src/ALE/Recon1d_PPM_CW.F90 index 6342c3ef3a..44894315f9 100644 --- a/src/ALE/Recon1d_PPM_CW.F90 +++ b/src/ALE/Recon1d_PPM_CW.F90 @@ -387,14 +387,14 @@ logical function unit_tests(this, verbose, stdout, stderr) ! f[i] = [ ( 3 i )^3 - ( 3 i - 3 )^3 ] i=1,2,3,4,5 ! means: 1, 7, 19, 37, 61 ! edges: 0, 3, 12, 27, 48, 75 -! call this%reconstruct( (/3.,3.,3.,3.,3./), (/1.,7.,19.,37.,61./) ) -! do k = 1, 5 -! ul(k) = this%f(k, 0.) -! um(k) = this%f(k, 0.5) -! ur(k) = this%f(k, 1.) -! enddo -! call test%real_arr(5, ul, (/1.,3.,12.,27.,61./), 'Return left edge', robits=2) -! call test%real_arr(5, ur, (/1.,12.,27.,48.,61./), 'Return right edge', robits=1) + call this%reconstruct( (/3.,3.,3.,3.,3./), (/1.,7.,19.,37.,61./) ) + do k = 1, 5 + ul(k) = this%f(k, 0.) + um(k) = this%f(k, 0.5) + ur(k) = this%f(k, 1.) + enddo + call test%real_arr(5, ul, (/1.,3.,12.,27.,61./), 'Return left edge', robits=2) + call test%real_arr(5, ur, (/1.,12.,27.,48.,61./), 'Return right edge', robits=1) call this%destroy() deallocate( um, ul, ur, ull, urr )