Added tests for generic operators to higher order - will need updating after the new derivative operator is finalized

Author: xtalax

test/generic_operator_check.jl

@@ -0,0 +1,45 @@
+using DiffEqOperators
+
+n = 100
+x=0.0:0.01:2π
+xprime = x[2:(end-1)]
+dx=diff(x)
+y = exp.(π*x)
+y_im = exp.(π*im*x)
+yim_ = y_im[2:(end-1)]
+y_ = y[2:(end-1)]
+
+
+
+for dor in 1:6, aor in 1:8
+
+    D1 = DerivativeOperator{Float64}(dor,aor,dx[1],length(x))
+    D2 = DiffEqOperators.FiniteDifference{Float64}(dor,aor,dx,length(x))
+    D = (D1,D2)
+    @test convert(Array, D1) ≈ convert(Array, D2)
+
+    #take derivatives
+    yprime1 = D1*y
+    yprime2 = D2*y
+
+    #test result
+    @test yprime1 ≈ (π^dor)*y_ # test operator with known derivative of exp(kx)
+    @test yprime2 ≈ (π^dor)*y_ # test operator with known derivative of exp(kx)
+
+    #test equivalance
+    @test yprime1 ≈ yprime2
+
+    #take derivatives
+    y_imprime1 = D1*y_im
+    y_imprime2 = D2*y_im
+
+    #test result
+    @test y_imprime1 ≈ ((pi*im)^dor)*yim_ # test operator with known derivative of exp(jkx)
+    @test y_imprime2 ≈ ((pi*im)^dor)*yim_ # test operator with known derivative of exp(jkx)
+
+    #test equivalance
+    @test y_imprime1 ≈ y_imprime2
+
+    #TODO: implement specific tests for the left and right boundary regions, waiting until after update
+    end
+end

test/generic_operator_validation.jl

@@ -5,7 +5,7 @@ import Base: isapprox
 @time @safetestset "Robin Boundary Condition Operators" begin include("robin.jl") end
 @time @safetestset "JacVec Operators Interface" begin include("jacvec_operators.jl") end
 @time @safetestset "Composite Operators Interface" begin include("composite_operators_interface.jl") end
-
+@time @safetestset "Validate and Compare Generic Operators" begin include("generic_operator_validation.jl") end
 @time @safetestset "Derivative Operators Interface" begin include("derivative_operators_interface.jl") end
 #@time @safetestset "2nd order check" begin include("2nd_order_check.jl") end
 #@time @safetestset "Finite Difference Operator" begin include("generic_operator_check.jl") end