#352 more 2D tests

Author: valentinsulzer

test/MOL/MOL_1D_Linear_Diffusion.jl test/MOL/MOL_1D_Diffusion.jltest/MOL/MOL_1D_Linear_Diffusion.jl renamed to test/MOL/MOL_1D_Diffusion.jl

@@ -37,12 +37,15 @@ using ModelingToolkit: Differential
 
     # Method of lines discretization
     dx = 0.1; dy = 0.2
-    discretization = MOLFiniteDifference([x=>dx,y=>dy],t)
-    prob = ModelingToolkit.discretize(pdesys,discretization)
-    sol = solve(prob,Tsit5())
-
-    # Test against exact solution
-    # TODO: do this properly when sol[u] with reshape etc works
-    @test sol.u[1][1] ≈ analytic_sol_func(sol.t[1],0.1,0.2)
-    @test sol.u[1][2] ≈ analytic_sol_func(sol.t[1],0.2,0.2)
+    for order in [2,4,6]
+        discretization = MOLFiniteDifference([x=>dx,y=>dy],t;centered_order=order)
+        prob = ModelingToolkit.discretize(pdesys,discretization)
+        sol = solve(prob,Tsit5())
+
+        # Test against exact solution
+        # TODO: do this properly when sol[u] with reshape etc works
+        @test sol.u[1][1] ≈ analytic_sol_func(sol.t[1],0.1,0.2)
+        @test sol.u[1][2] ≈ analytic_sol_func(sol.t[1],0.2,0.2)
+
+    end
 end

test/MOL/MOL_2D_Diffusion.jl

@@ -38,7 +38,8 @@ end
 if GROUP == "All" || GROUP == "MOLFiniteDifference"
     @time @safetestset "MOLFiniteDifference Interface" begin include("MOL/MOLtest.jl") end
     #@time @safetestset "MOLFiniteDifference Interface: Linear Convection" begin include("MOL/MOL_1D_Linear_Convection.jl") end
-    @time @safetestset "MOLFiniteDifference Interface: Linear Diffusion" begin include("MOL/MOL_1D_Linear_Diffusion.jl") end
+    @time @safetestset "MOLFiniteDifference Interface: 1D Diffusion" begin include("MOL/MOL_1D_Diffusion.jl") end
+    @time @safetestset "MOLFiniteDifference Interface: 2D Diffusion" begin include("MOL/MOL_2D_Diffusion.jl") end
 end
 
 if GROUP == "All" || GROUP == "Misc"