diff --git a/test/test_FiniteDifferences2D.jl b/test/test_FiniteDifferences2D.jl index 5553953f..92287849 100644 --- a/test/test_FiniteDifferences2D.jl +++ b/test/test_FiniteDifferences2D.jl @@ -36,12 +36,16 @@ end @parallel d_ya!(R, Ay) = (@all(R) = @d_ya(Ay); return) @parallel d_xi!(R, Axyy) = (@all(R) = @d_xi(Axyy); return) @parallel d_yi!(R, Axxy) = (@all(R) = @d_yi(Axxy); return) + @parallel d2_xa!(R, Axx) = (@all(R) = @d2_xa(Axx); return) + @parallel d2_ya!(R, Ayy) = (@all(R) = @d2_ya(Ayy); return) @parallel d2_xi!(R, Axxyy) = (@all(R) = @d2_xi(Axxyy); return) @parallel d2_yi!(R, Axxyy) = (@all(R) = @d2_yi(Axxyy); return) R.=0; @parallel d_xa!(R, Ax); @test all(Array(R .== Ax[2:end, :].-Ax[1:end-1, :])) R.=0; @parallel d_ya!(R, Ay); @test all(Array(R .== Ay[ :,2:end].-Ay[ :,1:end-1])) R.=0; @parallel d_xi!(R, Axyy); @test all(Array(R .== Axyy[2:end ,2:end-1].-Axyy[1:end-1,2:end-1])) R.=0; @parallel d_yi!(R, Axxy); @test all(Array(R .== Axxy[2:end-1,2:end ].-Axxy[2:end-1,1:end-1])) + R.=0; @parallel d2_xa!(R, Axx); @test all(Array(R .== (Axx[3:end, :].-Axx[2:end-1, :]).-(Axx[2:end-1, :].-Axx[1:end-2, :]))) + R.=0; @parallel d2_ya!(R, Ayy); @test all(Array(R .== (Ayy[ :,3:end].-Ayy[ :,2:end-1]).-(Ayy[ :,2:end-1].-Ayy[ :,1:end-2]))) R.=0; @parallel d2_xi!(R, Axxyy); @test all(Array(R .== (Axxyy[3:end,2:end-1].-Axxyy[2:end-1,2:end-1]).-(Axxyy[2:end-1,2:end-1].-Axxyy[1:end-2,2:end-1]))) R.=0; @parallel d2_yi!(R, Axxyy); @test all(Array(R .== (Axxyy[2:end-1,3:end].-Axxyy[2:end-1,2:end-1]).-(Axxyy[2:end-1,2:end-1].-Axxyy[2:end-1,1:end-2]))) end;