Nonlinear problem tests passing

Author: xtalax

test/components/DiscreteSpace.jl

@@ -49,7 +49,7 @@ using ModelingToolkit, MethodOfLines, DomainSets, Test, Symbolics, SymbolicUtils
     @test s.axies[x] == s.grid[x]
     @test s.axies[y] == s.grid[y]
 
-    @test all([all(I[i] .∈ (collect(s.Igrid),)) for I in values(s.Iedge), i in [1,2]])  
+    @test all([all(I[i] .∈ (collect(s.Igrid),)) for I in values(s.Iedges), i in [1,2]])  
 
     @test s.Iaxies == s.Igrid
 
@@ -101,7 +101,7 @@ end
     @test s.axies[x] != s.grid[x]
     @test s.axies[y] != s.grid[y]
 
-    @test all([all(I[i] .∈ (collect(s.Igrid),)) for I in values(s.Iedge), i in [1,2]])  
+    @test all([all(I[i] .∈ (collect(s.Igrid),)) for I in values(s.Iedges), i in [1,2]])  
 
     @test s.Iaxies != s.Igrid
 

test/pde_systems/MOL_1D_PDAE.jl

@@ -37,8 +37,9 @@ using ModelingToolkit: Differential
     # Method of lines discretization
     l = 20
     dx = range(0.0,1.0,length=l)
+    dx_ = dx[2] - dx[1]
     order = 2
-    discretization = MOLFiniteDifference([x=>dx],t)
+    discretization = MOLFiniteDifference([x=>dx_],t)
 
     # Convert the PDE problem into an ODE problem
     prob = discretize(pdesys,discretization)

test/pde_systems/MOL_NonlinearProblem.jl

@@ -102,8 +102,8 @@ end
     # test boundary
     @test all(abs.(u_sol[:,1]) .< eps(Float32))
     @test all(abs.(u_sol[1,:]) .< eps(Float32))
-    @test u_sol[:,end] ≈ 0:dy:1.0
-    @test u_sol[end,:] ≈ 0:dx:1.0
+    @test u_sol[2:end-1,end] ≈ (0.0:dy:1.0)[2:end-1]
+    @test u_sol[end,2:end-1] ≈ (0.0:dx:1.0)[2:end-1]
 
     # test interior with finite differences
     interior = CartesianIndices((axes(xs)[1], axes(ys)[1]))[2:end-1,2:end-1]