Skip to content
Snippets Groups Projects
Commit 842183de authored by Brian Groenke's avatar Brian Groenke
Browse files

Add missing import for adstrip

parent cad2f38c
No related branches found
No related tags found
1 merge request!36Fix CFL step limiter callback choking on autodiff types
......@@ -50,7 +50,8 @@ p.soil.n .= 2.0
p.soil.Tₘ .= 273.15 # K
# CryoGrid front-end for ODEProblem
prob = CryoGridProblem(model,tspan,p)
# OPTIONAL: solve for short time period to
# OPTIONAL: solve for short time period to "warm up" (i.e. let Julia compile everything)
solve(remake(prob, tspan=(prob.tspan[1],prob.tspan[1]+60.0)), Euler(), dt=10.0);
# solve with forward Euler (w/ CFL) and construct CryoGridOutput from solution
out = @time solve(prob, Euler(), dt=2*60.0, callback=CFLStepLimiter(model), saveat=24*3600.0, progress=true) |> CryoGridOutput;
# Plot it!
......
......@@ -4,7 +4,7 @@ import CryoGrid
using CryoGrid.Numerics
using CryoGrid.Setup: CryoGridSetup, HeatOnlySetup, getvar
using CryoGrid.Utils
using CryoGrid.Utils: adstrip
using DiffEqCallbacks
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment