PDE example works

ChrisRackauckas · ChrisRackauckas · commit bd76e503dff2 · 2019-09-20T08:57:31.000-04:00
diff --git a/src/domains.jl b/src/domains.jl
@@ -1,3 +1,5 @@
+export IntervalDomain, ProductDomain, ⊗, CircleDomain
+
 abstract type AbstractDomain{T,N} end
 
 struct VarDomainPairing
diff --git a/src/equations.jl b/src/equations.jl
@@ -1,4 +1,4 @@
-export Equation
+export Equation, ConstrainedEquation
 
 
 """
diff --git a/src/systems/pde/pdesystem.jl b/src/systems/pde/pdesystem.jl
@@ -1,3 +1,5 @@
+export PDESystem
+
 struct PDESystem <: ModelingToolkit.AbstractSystem
   eq
   bcs
diff --git a/test/pde.jl b/test/pde.jl
@@ -14,15 +14,7 @@ domains = [t ∈ IntervalDomain(0.0,1.0),
 
 pdesys = PDESystem(eq,bcs,domains,[t,x],[u])
 discretization = MOLFiniteDifference(0.1)
-prob = discretize(pdesys,discretization)
+prob = discretize(pdesys,discretization) # This gives an ODEProblem since it's time-dependent
 
 using OrdinaryDiffEq
 sol = solve(prob,Tsit5(),saveat=0.1)
-
-#=
-using Plots
-plot(prob.space,Array(prob.extrapolation*sol[1]))
-plot!(prob.space,Array(prob.extrapolation*sol[2]))
-plot!(prob.space,Array(prob.extrapolation*sol[3]))
-plot!(prob.space,Array(prob.extrapolation*sol[4]))
-=#

Original file line number	Diff line number	Diff line change
`@@ -1,3 +1,5 @@`
	`1`	`+export IntervalDomain, ProductDomain, ⊗, CircleDomain`
	`2`	`+`
`1`	`3`	`abstract type AbstractDomain{T,N} end`
`2`	`4`
`3`	`5`	`struct VarDomainPairing`
Original file line number	Diff line number	Diff line change
`@@ -1,4 +1,4 @@`
`1`		`-export Equation`
	`1`	`+export Equation, ConstrainedEquation`
`2`	`2`
`3`	`3`
`4`	`4`	`"""`
Original file line number	Diff line number	Diff line change
`@@ -1,3 +1,5 @@`
	`1`	`+export PDESystem`
	`2`	`+`
`1`	`3`	`struct PDESystem <: ModelingToolkit.AbstractSystem`
`2`	`4`	`eq`
`3`	`5`	`bcs`