WIP

Author: Keno

benchmark/thermalfluid.jl

@@ -103,6 +103,7 @@ function (this::CircularWall{Ne, Nn})(inner_heatport, outer_heatport) where {Ne,
     Qe = ntuple(_->continuous(), Ne)
 
     ntuple(i->initial!(Tn[i] - this.T0), Nn)
+    ntuple(i->initial!(Qe[i] - this.T0), Ne)
 
     dn = (this.d, (this.d .+ 2.0 .* cumsum(this.t_layer))...)
     Cni = Cn_circular_wall_inner.(dn[1:Ne], dn[2:Nn], this.cp, this.ρ)

docs/src/architecture.md

@@ -0,0 +1,43 @@
+```mermaid
+flowchart TD
+    US[Uncompressed State]
+    P[Parameters]
+    subgraph init
+        isicm([Init SICM])
+        ICS[Init Compressed State]
+        IB[Init Blob]
+        irhs([Init RHS])
+        IR[Init Residual]
+        nlsolve([nlsolve])
+        pr([Project])
+    end
+    subgraph dae
+        osicm([ODE/DAE SICM])
+        OCS[DAE/ODE Compressed State]
+        OB[DAE/ODE Blob]
+        orhs([ODE/DAE RHS])
+        OR[ODE/DAE Residual]
+        osolve([ODE solver])
+        cmp([Compress])
+    end
+    prep([Prepare])
+    US --> pr --> ICS
+    US --> cmp --> OCS
+    P --> isicm --> IB
+    P --> osicm --> OB
+    US --> osicm
+    ICS --> irhs
+    IB --> irhs
+    irhs --> IR
+    IR --> nlsolve
+    nlsolve --> ICS
+    OB --> orhs
+    OCS --> orhs
+    orhs --> OR
+    OR --> osolve
+    osolve --> OCS
+    ICS --> prep
+    IB --> prep
+    prep --> OCS
+    prep --> OB
+```

test/ipo.jl

@@ -11,13 +11,14 @@ using OrdinaryDiffEq
 #= Basic IPO: We need to read the incidence of the contained `-` =#
 @noinline function onecall!()
     x = continuous()
+    initial!(x - 1.0)
     always!(ddt(x) - x)
 end
 
 onecall!()
-sol = solve(DAECProblem(onecall!, (1,) .=> 1.), IDA())
+sol = solve(DAECProblem(onecall!), DFBDF(autodiff=false))
 @test all(map((x,y)->isapprox(x[], y, atol=1e-2), sol[1, :], exp.(sol.t)))
-sol = solve(ODECProblem(onecall!, (1,) .=> 1.), Rodas5(autodiff=false))
+sol = solve(ODECProblem(onecall!), Rodas5(autodiff=false))
 @test all(map((x,y)->isapprox(x[], y, atol=1e-2), sol[1, :], exp.(sol.t)))
 
 #==== Contained Equations ============#
@@ -27,23 +28,24 @@ function twocall!()
 end
 
 twocall!()
-dae_sol = solve(DAECProblem(twocall!, (1, 2) .=> 1.), IDA())
-ode_sol = solve(ODECProblem(twocall!, (1, 2) .=> 1.), Rodas5(autodiff=false))
+dae_sol = solve(DAECProblem(twocall!), DFBDF(autodiff=false))
+ode_sol = solve(ODECProblem(twocall!), Rodas5(autodiff=false))
 for (sol, i) in Iterators.product((dae_sol, ode_sol), 1:2)
     @test all(map((x,y)->isapprox(x[], y, atol=1e-2), sol[i, :], exp.(sol.t)))
 end
 
 #============== NonLinear ============#
 @noinline function sin!()
     x = continuous()
+    initial!(x - 1.0)
     always!(ddt(x) - sin(x))
 end
 function sin2!()
     sin!(); sin!();
     return nothing
 end
-dae_sol = solve(DAECProblem(sin2!, (1, 2) .=> 1.), IDA())
-ode_sol = solve(ODECProblem(sin2!, (1, 2) .=> 1.), Rodas5(autodiff=false))
+dae_sol = solve(DAECProblem(sin2!), DFBDF(autodiff=false))
+ode_sol = solve(ODECProblem(sin2!), Rodas5(autodiff=false))
 for (sol, i) in Iterators.product((dae_sol, ode_sol), 1:2)
     @test all(map((x,y)->isapprox(x[], y, atol=1e-2), sol[i, :], 2*acot.(exp.(-sol.t).*cot(1/2))))
 end
@@ -52,10 +54,11 @@ end
 @noinline sink!(x, v) = always!(x - v)
 function sinsink!()
     x = continuous()
+    initial!(x - 1.0)
     sink!(ddt(x), sin(x))
 end
-dae_sol = solve(DAECProblem(sinsink!, (1,) .=> 1.), IDA())
-ode_sol = solve(ODECProblem(sinsink!, (1,) .=> 1.), Rodas5(autodiff=false))
+dae_sol = solve(DAECProblem(sinsink!), DFBDF(autodiff=false))
+ode_sol = solve(ODECProblem(sinsink!), Rodas5(autodiff=false))
 for sol in (dae_sol, ode_sol)
     @test all(map((x,y)->isapprox(x[], y, atol=1e-2), sol[1, :], 2*acot.(exp.(-sol.t).*cot(1/2))))
 end
@@ -66,7 +69,7 @@ function sinsink2!()
     x = continuous()
     sink2!(x, sin(x))
 end
-dae_sol = solve(DAECProblem(sinsink2!, (1,) .=> 1.), IDA())
+dae_sol = solve(DAECProblem(sinsink2!, (1,) .=> 1.), DFBDF(autodiff=false))
 ode_sol = solve(ODECProblem(sinsink2!, (1,) .=> 1.), Rodas5(autodiff=false))
 for sol in (dae_sol, ode_sol)
     @test all(map((x,y)->isapprox(x[], y, atol=1e-2), sol[1, :], 2*acot.(exp.(-sol.t).*cot(1/2))))
@@ -91,7 +94,7 @@ function (this::sicm2!)()
     sicm!(this.a)(); sicm!(this.b)();
     return nothing
 end
-dae_sol = solve(DAECProblem(sicm2!(1., 1.), (1, 2) .=> 1.), IDA())
+dae_sol = solve(DAECProblem(sicm2!(1., 1.), (1, 2) .=> 1.), DFBDF(autodiff=false))
 ode_sol = solve(ODECProblem(sicm2!(1., 1.), (1, 2) .=> 1.), Rodas5(autodiff=false))
 for (sol, i) in Iterators.product((dae_sol, ode_sol), 1:2)
     @test all(map((x,y)->isapprox(x[], y, atol=1e-2), sol[i, :], 1. .+ sol.t))
@@ -104,6 +107,7 @@ end
 
 @noinline function (this::nlsicm!)()
     x = continuous()
+    initial!(x - 1.0)
     always!(ddt(x) - sin(this.arg))
 end
 
@@ -116,8 +120,8 @@ function (this::nlsicm2!)()
     nlsicm!(this.a)(); nlsicm!(this.b)();
     return nothing
 end
-dae_sol = solve(DAECProblem(nlsicm2!(1., 1.), (1, 2) .=> 1.), IDA())
-ode_sol = solve(ODECProblem(nlsicm2!(1., 1.), (1, 2) .=> 1.), Rodas5(autodiff=false))
+dae_sol = solve(DAECProblem(nlsicm2!(1., 1.)), DFBDF(autodiff=false))
+ode_sol = solve(ODECProblem(nlsicm2!(1., 1.)), Rodas5(autodiff=false))
 for (sol, i) in Iterators.product((dae_sol, ode_sol), 1:2)
     @test all(map((x,y)->isapprox(x[], y, atol=1e-2), sol[i, :], 1. .+ sin(1.)*sol.t))
 end