Merge pull request #30 from xtalax/master

Fix tests

Author: ChrisRackauckas

Project.toml

@@ -5,6 +5,7 @@ version = "0.1.0"
 
 [deps]
 DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"
+DiffEqOperators = "9fdde737-9c7f-55bf-ade8-46b3f136cc48"
 DomainSets = "5b8099bc-c8ec-5219-889f-1d9e522a28bf"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
@@ -17,11 +18,11 @@ Symbolics = "0c5d862f-8b57-4792-8d23-62f2024744c7"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
 [compat]
-OrdinaryDiffEq = "6"
-ModelingToolkit = "8"
-SymbolicUtils = "0.19"
-SafeTestsets = "0.0.1"
 DomainSets = "0.5"
+ModelingToolkit = "8"
 NonlinearSolve = "0.3"
+OrdinaryDiffEq = "6"
+SafeTestsets = "0.0.1"
 SciMLBase = "1"
+SymbolicUtils = "0.19"
 Symbolics = "4"

test/runtests.jl

@@ -17,6 +17,7 @@ if GROUP == "All" || GROUP == "MOLFiniteDifference"
     @time @safetestset "MOLFiniteDifference Interface: 1D HigherOrder" begin include("pde_systems/MOL_1D_HigherOrder.jl") end
     @time @safetestset "MOLFiniteDifference Interface: 1D Partial DAE" begin include("pde_systems/MOL_1D_PDAE.jl") end
     @time @safetestset "MOLFiniteDifference Interface: Stationary Nonlinear Problems" begin include("pde_systems/MOL_NonlinearProblem.jl") end
+    @time @safetestset "MOLFiniteDifference Utils" begin include("utils_test.jl") end
 end
 
 end

test/utils_test.jl

@@ -0,0 +1,74 @@
+using MethodOfLines, Test, ModelingToolkit
+
+@testset "count differentials 1D" begin
+    @parameters t x
+    @variables u(..)
+    Dt = Differential(t)
+
+    Dx = Differential(x)
+    eq  = Dt(u(t,x)) ~ -Dx(u(t,x))
+    @test first(MethodOfLines.differential_order(eq.rhs, x.val)) == 1
+    @test isempty(MethodOfLines.differential_order(eq.rhs, t.val))
+    @test first(MethodOfLines.differential_order(eq.lhs, t.val)) == 1
+    @test isempty(MethodOfLines.differential_order(eq.lhs, x.val))
+
+    Dxx = Differential(x)^2
+    eq  = Dt(u(t,x)) ~ Dxx(u(t,x))
+    @test first(MethodOfLines.differential_order(eq.rhs, x.val)) == 2
+    @test isempty(MethodOfLines.differential_order(eq.rhs, t.val))
+    @test first(MethodOfLines.differential_order(eq.lhs, t.val)) == 1
+    @test isempty(MethodOfLines.differential_order(eq.lhs, x.val))
+
+    Dxxxx = Differential(x)^4
+    eq  = Dt(u(t,x)) ~ -Dxxxx(u(t,x))
+    @test first(MethodOfLines.differential_order(eq.rhs, x.val)) == 4
+    @test isempty(MethodOfLines.differential_order(eq.rhs, t.val))
+    @test first(MethodOfLines.differential_order(eq.lhs, t.val)) == 1
+    @test isempty(MethodOfLines.differential_order(eq.lhs, x.val))
+end
+
+@testset "count differentials 2D" begin
+    @parameters t x y
+    @variables u(..)
+    Dxx = Differential(x)^2
+    Dyy = Differential(y)^2
+    Dt = Differential(t)
+
+    eq  = Dt(u(t,x,y)) ~ Dxx(u(t,x,y)) + Dyy(u(t,x,y))
+    @test first(MethodOfLines.differential_order(eq.rhs, x.val)) == 2
+    @test first(MethodOfLines.differential_order(eq.rhs, y.val)) == 2
+    @test isempty(MethodOfLines.differential_order(eq.rhs, t.val))
+    @test first(MethodOfLines.differential_order(eq.lhs, t.val)) == 1
+    @test isempty(MethodOfLines.differential_order(eq.lhs, x.val))
+    @test isempty(MethodOfLines.differential_order(eq.lhs, y.val))
+end
+
+@testset "count with mixed terms" begin
+    @parameters t x y
+    @variables u(..)
+    Dxx = Differential(x)^2
+    Dyy = Differential(y)^2
+    Dx = Differential(x)
+    Dy = Differential(y)
+    Dt = Differential(t)
+
+    eq  = Dt(u(t,x,y)) ~ Dxx(u(t,x,y)) + Dyy(u(t,x,y)) + Dx(Dy(u(t,x,y)))
+    @test MethodOfLines.differential_order(eq.rhs, x.val) == Set([2, 1])
+    @test MethodOfLines.differential_order(eq.rhs, y.val) == Set([2, 1])
+end
+
+@testset "Kuramoto–Sivashinsky equation" begin
+    @parameters x, t
+    @variables u(..)
+    Dt = Differential(t)
+    Dx = Differential(x)
+    Dx2 = Differential(x)^2
+    Dx3 = Differential(x)^3
+    Dx4 = Differential(x)^4
+
+    α = 1
+    β = 4
+    γ = 1
+    eq = Dt(u(x,t)) + u(x,t)*Dx(u(x,t)) + α*Dx2(u(x,t)) + β*Dx3(u(x,t)) + γ*Dx4(u(x,t)) ~ 0
+    @test MethodOfLines.differential_order(eq.lhs, x.val) == Set([4, 3, 2, 1])
+end