Skip to content

Commit 0ca92b3

Browse files
format
1 parent dc4ec92 commit 0ca92b3

File tree

49 files changed

+221
-178
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

49 files changed

+221
-178
lines changed

lib/OrdinaryDiffEqAdamsBashforthMoulton/src/adams_bashforth_moulton_caches.jl

+4-2
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,9 @@
11
abstract type ABMMutableCache <: OrdinaryDiffEqMutableCache end
22
abstract type ABMVariableCoefficientMutableCache <: OrdinaryDiffEqMutableCache end
3-
get_fsalfirstlast(cache::ABMMutableCache,u) = (cache.fsalfirst, cache.k)
4-
get_fsalfirstlast(cache::ABMVariableCoefficientMutableCache,u) = (cache.fsalfirst, cache.k4)
3+
get_fsalfirstlast(cache::ABMMutableCache, u) = (cache.fsalfirst, cache.k)
4+
function get_fsalfirstlast(cache::ABMVariableCoefficientMutableCache, u)
5+
(cache.fsalfirst, cache.k4)
6+
end
57
@cache mutable struct AB3Cache{uType, rateType} <: ABMMutableCache
68
u::uType
79
uprev::uType

lib/OrdinaryDiffEqBDF/src/algorithms.jl

+3-2
Original file line numberDiff line numberDiff line change
@@ -150,7 +150,7 @@ end
150150
function QNDF1(; chunk_size = Val{0}(), autodiff = Val{true}(), standardtag = Val{true}(),
151151
concrete_jac = nothing, diff_type = Val{:forward},
152152
linsolve = nothing, precs = DEFAULT_PRECS, nlsolve = NLNewton(),
153-
extrapolant = :linear, kappa = -37//200,
153+
extrapolant = :linear, kappa = -37 // 200,
154154
controller = :Standard, step_limiter! = trivial_limiter!)
155155
QNDF1{
156156
_unwrap_val(chunk_size), _unwrap_val(autodiff), typeof(linsolve), typeof(nlsolve),
@@ -233,7 +233,8 @@ function QNDF(; max_order::Val{MO} = Val{5}(), chunk_size = Val{0}(),
233233
diff_type = Val{:forward},
234234
linsolve = nothing, precs = DEFAULT_PRECS, nlsolve = NLNewton(), κ = nothing,
235235
tol = nothing,
236-
extrapolant = :linear, kappa = (-37//200, -1//9, -823//10000, -83//2000, 0//1),
236+
extrapolant = :linear, kappa = (
237+
-37 // 200, -1 // 9, -823 // 10000, -83 // 2000, 0 // 1),
237238
controller = :Standard, step_limiter! = trivial_limiter!) where {MO}
238239
QNDF{MO, _unwrap_val(chunk_size), _unwrap_val(autodiff), typeof(linsolve),
239240
typeof(nlsolve), typeof(precs), diff_type, _unwrap_val(standardtag),

lib/OrdinaryDiffEqBDF/src/bdf_caches.jl

+3-1
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,7 @@
11
abstract type BDFMutableCache <: OrdinaryDiffEqMutableCache end
2-
get_fsalfirstlast(cache::BDFMutableCache,u) = (cache.fsalfirst, du_alias_or_new(cache.nlsolver, cache.fsalfirst))
2+
function get_fsalfirstlast(cache::BDFMutableCache, u)
3+
(cache.fsalfirst, du_alias_or_new(cache.nlsolver, cache.fsalfirst))
4+
end
35

46
@cache mutable struct ABDF2ConstantCache{N, dtType, rate_prototype} <:
57
OrdinaryDiffEqConstantCache

lib/OrdinaryDiffEqBDF/src/controllers.jl

+2-2
Original file line numberDiff line numberDiff line change
@@ -206,7 +206,7 @@ function choose_order!(alg::FBDF, integrator,
206206
Val(max_order))
207207
local terk_tmp
208208
if u isa Number
209-
terk_tmp = fd_weights[k - 2, 1]*u
209+
terk_tmp = fd_weights[k - 2, 1] * u
210210
for i in 2:(k - 2)
211211
terk_tmp += fd_weights[i, k - 2] * u_history[i - 1]
212212
end
@@ -215,7 +215,7 @@ function choose_order!(alg::FBDF, integrator,
215215
# we need terk_tmp to be mutable.
216216
# so it can be updated
217217
terk_tmp = similar(u)
218-
@.. terk_tmp = fd_weights[k - 2, 1]*_vec(u)
218+
@.. terk_tmp = fd_weights[k - 2, 1] * _vec(u)
219219
for i in 2:(k - 2)
220220
@.. @views terk_tmp += fd_weights[i, k - 2] * u_history[:, i - 1]
221221
end

lib/OrdinaryDiffEqBDF/src/dae_caches.jl

+4-2
Original file line numberDiff line numberDiff line change
@@ -1,5 +1,7 @@
11
abstract type DAEBDFMutableCache <: OrdinaryDiffEqMutableCache end
2-
get_fsalfirstlast(cache::DAEBDFMutableCache,u) = (cache.fsalfirst, du_alias_or_new(cache.nlsolver, cache.fsalfirst))
2+
function get_fsalfirstlast(cache::DAEBDFMutableCache, u)
3+
(cache.fsalfirst, du_alias_or_new(cache.nlsolver, cache.fsalfirst))
4+
end
35

46
@cache mutable struct DImplicitEulerCache{uType, rateType, uNoUnitsType, N} <:
57
DAEBDFMutableCache
@@ -13,7 +15,7 @@ get_fsalfirstlast(cache::DAEBDFMutableCache,u) = (cache.fsalfirst, du_alias_or_n
1315
end
1416

1517
# Not FSAL
16-
get_fsalfirstlast(cache::DImplicitEulerCache,u) = (u,u)
18+
get_fsalfirstlast(cache::DImplicitEulerCache, u) = (u, u)
1719

1820
mutable struct DImplicitEulerConstantCache{N} <: OrdinaryDiffEqConstantCache
1921
nlsolver::N

lib/OrdinaryDiffEqBDF/src/dae_perform_step.jl

-2
Original file line numberDiff line numberDiff line change
@@ -160,8 +160,6 @@ end
160160
end
161161

162162
function initialize!(integrator, cache::DABDF2Cache)
163-
164-
165163
integrator.kshortsize = 2
166164
@unpack k₁, k₂ = cache.eulercache
167165
resize!(integrator.k, integrator.kshortsize)

lib/OrdinaryDiffEqCore/ext/OrdinaryDiffEqCoreEnzymeCoreExt.jl

+16-4
Original file line numberDiff line numberDiff line change
@@ -2,10 +2,22 @@ module OrdinaryDiffEqCoreEnzymeCoreExt
22
import OrdinaryDiffEqCore, EnzymeCore
33

44
EnzymeCore.EnzymeRules.inactive(::typeof(OrdinaryDiffEqCore.increment_nf!), args...) = true
5-
EnzymeCore.EnzymeRules.inactive(::typeof(OrdinaryDiffEqCore.fixed_t_for_floatingpoint_error!), args...) = true
6-
EnzymeCore.EnzymeRules.inactive(::typeof(OrdinaryDiffEqCore.increment_accept!), args...) = true
7-
EnzymeCore.EnzymeRules.inactive(::typeof(OrdinaryDiffEqCore.increment_reject!), args...) = true
8-
EnzymeCore.EnzymeRules.inactive(::typeof(OrdinaryDiffEqCore.increment_nf_perform_step!), args...) = true
5+
function EnzymeCore.EnzymeRules.inactive(
6+
::typeof(OrdinaryDiffEqCore.fixed_t_for_floatingpoint_error!), args...)
7+
true
8+
end
9+
function EnzymeCore.EnzymeRules.inactive(
10+
::typeof(OrdinaryDiffEqCore.increment_accept!), args...)
11+
true
12+
end
13+
function EnzymeCore.EnzymeRules.inactive(
14+
::typeof(OrdinaryDiffEqCore.increment_reject!), args...)
15+
true
16+
end
17+
function EnzymeCore.EnzymeRules.inactive(
18+
::typeof(OrdinaryDiffEqCore.increment_nf_perform_step!), args...)
19+
true
20+
end
921
EnzymeCore.EnzymeRules.inactive(::typeof(OrdinaryDiffEqCore.check_error!), args...) = true
1022
EnzymeCore.EnzymeRules.inactive(::typeof(OrdinaryDiffEqCore.log_step!), args...) = true
1123

lib/OrdinaryDiffEqCore/src/OrdinaryDiffEqCore.jl

+1-2
Original file line numberDiff line numberDiff line change
@@ -41,7 +41,7 @@ using SimpleUnPack, RecursiveArrayTools, DataStructures, ArrayInterface
4141
import TruncatedStacktraces
4242

4343
import StaticArraysCore: SArray, MVector, SVector, StaticArray, MMatrix,
44-
StaticMatrix
44+
StaticMatrix
4545

4646
# Integrator Interface
4747
import DiffEqBase: resize!, deleteat!, addat!, full_cache, user_cache, u_cache, du_cache,
@@ -54,7 +54,6 @@ import DiffEqBase: resize!, deleteat!, addat!, full_cache, user_cache, u_cache,
5454
isautodifferentiable,
5555
get_tstops, get_tstops_array, get_tstops_max
5656

57-
5857
using DiffEqBase: check_error!, @def, _vec, _reshape
5958

6059
using FastBroadcast: @.., True, False

lib/OrdinaryDiffEqCore/src/caches/basic_caches.jl

+6-5
Original file line numberDiff line numberDiff line change
@@ -5,15 +5,15 @@ struct ODEEmptyCache <: OrdinaryDiffEqConstantCache end
55
struct ODEChunkCache{CS} <: OrdinaryDiffEqConstantCache end
66

77
# Don't worry about the potential alloc on a constant cache
8-
get_fsalfirstlast(cache::OrdinaryDiffEqConstantCache,u) = (zero(u), zero(u))
8+
get_fsalfirstlast(cache::OrdinaryDiffEqConstantCache, u) = (zero(u), zero(u))
99

1010
mutable struct CompositeCache{T, F} <: OrdinaryDiffEqCache
1111
caches::T
1212
choice_function::F
1313
current::Int
1414
end
1515

16-
get_fsalfirstlast(cache::CompositeCache,u) = get_fsalfirstlast(cache.caches[1],u)
16+
get_fsalfirstlast(cache::CompositeCache, u) = get_fsalfirstlast(cache.caches[1], u)
1717

1818
mutable struct DefaultCache{T1, T2, T3, T4, T5, T6, A, F, uType} <: OrdinaryDiffEqCache
1919
args::A
@@ -28,12 +28,13 @@ mutable struct DefaultCache{T1, T2, T3, T4, T5, T6, A, F, uType} <: OrdinaryDiff
2828
cache6::T6
2929
function DefaultCache{T1, T2, T3, T4, T5, T6, F, uType}(
3030
args, choice_function, current, u) where {T1, T2, T3, T4, T5, T6, F, uType}
31-
new{T1, T2, T3, T4, T5, T6, typeof(args), F, uType}(args, choice_function, current, u)
31+
new{T1, T2, T3, T4, T5, T6, typeof(args), F, uType}(
32+
args, choice_function, current, u)
3233
end
3334
end
3435

35-
function get_fsalfirstlast(cache::DefaultCache,u)
36-
(cache.u,cache.u) # will be overwritten by the cache choice
36+
function get_fsalfirstlast(cache::DefaultCache, u)
37+
(cache.u, cache.u) # will be overwritten by the cache choice
3738
end
3839

3940
function alg_cache(alg::CompositeAlgorithm, u, rate_prototype, ::Type{uEltypeNoUnits},

lib/OrdinaryDiffEqCore/src/integrators/integrator_utils.jl

+8-8
Original file line numberDiff line numberDiff line change
@@ -255,8 +255,8 @@ function _loopfooter!(integrator)
255255
end
256256
if integrator.opts.progress && integrator.iter % integrator.opts.progress_steps == 0
257257
log_step!(integrator.opts.progress_name, integrator.opts.progress_id,
258-
integrator.opts.progress_message, integrator.dt, integrator.u,
259-
integrator.p, integrator.t, integrator.sol.prob.tspan)
258+
integrator.opts.progress_message, integrator.dt, integrator.u,
259+
integrator.p, integrator.t, integrator.sol.prob.tspan)
260260
end
261261

262262
# Take value because if t is dual then maxeig can be dual
@@ -280,18 +280,18 @@ end
280280

281281
function log_step!(progress_name, progress_id, progress_message, dt, u, p, t, tspan)
282282
t1, t2 = tspan
283-
@logmsg(LogLevel(-1),progress_name,
284-
_id=progress_id,
285-
message=progress_message(dt, u, p, t),
286-
progress=(t - t1) / (t2 - t1))
283+
@logmsg(LogLevel(-1), progress_name,
284+
_id=progress_id,
285+
message=progress_message(dt, u, p, t),
286+
progress=(t - t1) / (t2 - t1))
287287
end
288288

289289
function fixed_t_for_floatingpoint_error!(integrator, ttmp)
290290
if has_tstop(integrator)
291291
tstop = integrator.tdir * first_tstop(integrator)
292292
if abs(ttmp - tstop) <
293-
100eps(float(max(integrator.t, tstop) / oneunit(integrator.t))) *
294-
oneunit(integrator.t)
293+
100eps(float(max(integrator.t, tstop) / oneunit(integrator.t))) *
294+
oneunit(integrator.t)
295295
tstop
296296
else
297297
ttmp

lib/OrdinaryDiffEqCore/src/interp_func.jl

+1-2
Original file line numberDiff line numberDiff line change
@@ -62,7 +62,6 @@ end
6262

6363
# strip interpolation of function information
6464
function SciMLBase.strip_interpolation(id::InterpolationData)
65-
6665
cache = strip_cache(id.cache)
6766

6867
InterpolationData(nothing, id.timeseries,
@@ -78,7 +77,7 @@ end
7877
function strip_cache(cache)
7978
if hasfield(typeof(cache), :jac_config) || hasfield(typeof(cache), :grad_config)
8079
fieldnums = length(fieldnames(typeof(cache)))
81-
noth_list = fill(nothing,fieldnums)
80+
noth_list = fill(nothing, fieldnums)
8281
cache_type_name = Base.typename(typeof(cache)).wrapper
8382
cache_type_name(noth_list...)
8483
else

lib/OrdinaryDiffEqCore/src/perform_step/composite_perform_step.jl

+19-19
Original file line numberDiff line numberDiff line change
@@ -32,40 +32,40 @@ function initialize!(integrator, cache::DefaultCache)
3232
init_ith_default_cache(cache, algs, cache.current)
3333
u = integrator.u
3434
if cache.current == 1
35-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache1,u)
35+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache1, u)
3636
integrator.fsalfirst = fsalfirst
3737
integrator.fsallast = fsallast
3838
initialize!(integrator, cache.cache1)
3939
elseif cache.current == 2
40-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache2,u)
40+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache2, u)
4141
integrator.fsalfirst = fsalfirst
4242
integrator.fsallast = fsallast
4343
initialize!(integrator, cache.cache2)
4444
# the controller was initialized by default for algs[1]
4545
reset_alg_dependent_opts!(integrator.opts.controller, algs[1], algs[2])
4646
elseif cache.current == 3
47-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache3,u)
47+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache3, u)
4848
integrator.fsalfirst = fsalfirst
4949
integrator.fsallast = fsallast
5050
initialize!(integrator, cache.cache3)
5151
# the controller was initialized by default for algs[1]
5252
reset_alg_dependent_opts!(integrator.opts.controller, algs[1], algs[3])
5353
elseif cache.current == 4
54-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache4,u)
54+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache4, u)
5555
integrator.fsalfirst = fsalfirst
5656
integrator.fsallast = fsallast
5757
initialize!(integrator, cache.cache4)
5858
# the controller was initialized by default for algs[1]
5959
reset_alg_dependent_opts!(integrator.opts.controller, algs[1], algs[4])
6060
elseif cache.current == 5
61-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache5,u)
61+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache5, u)
6262
integrator.fsalfirst = fsalfirst
6363
integrator.fsallast = fsallast
6464
initialize!(integrator, cache.cache5)
6565
# the controller was initialized by default for algs[1]
6666
reset_alg_dependent_opts!(integrator.opts.controller, algs[1], algs[5])
6767
elseif cache.current == 6
68-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache6,u)
68+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache6, u)
6969
integrator.fsalfirst = fsalfirst
7070
integrator.fsallast = fsallast
7171
initialize!(integrator, cache.cache6)
@@ -79,20 +79,20 @@ function initialize!(integrator, cache::CompositeCache)
7979
cache.current = cache.choice_function(integrator)
8080
u = integrator.u
8181
if cache.current == 1
82-
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[1],u)
82+
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[1], u)
8383
integrator.fsalfirst = fsalfirst
8484
integrator.fsallast = fsallast
8585
initialize!(integrator, @inbounds(cache.caches[1]))
8686
elseif cache.current == 2
87-
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[2],u)
87+
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[2], u)
8888
integrator.fsalfirst = fsalfirst
8989
integrator.fsallast = fsallast
9090
initialize!(integrator, @inbounds(cache.caches[2]))
9191
# the controller was initialized by default for integrator.alg.algs[1]
9292
reset_alg_dependent_opts!(integrator.opts.controller, integrator.alg.algs[1],
9393
integrator.alg.algs[2])
9494
else
95-
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[cache.current],u)
95+
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[cache.current], u)
9696
integrator.fsalfirst = fsalfirst
9797
integrator.fsallast = fsallast
9898
initialize!(integrator, @inbounds(cache.caches[cache.current]))
@@ -106,12 +106,12 @@ function initialize!(integrator, cache::CompositeCache{Tuple{T1, T2}, F}) where
106106
cache.current = cache.choice_function(integrator)
107107
u = integrator.u
108108
if cache.current == 1
109-
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[1],u)
109+
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[1], u)
110110
integrator.fsalfirst = fsalfirst
111111
integrator.fsallast = fsallast
112112
initialize!(integrator, @inbounds(cache.caches[1]))
113113
elseif cache.current == 2
114-
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[2],u)
114+
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[2], u)
115115
integrator.fsalfirst = fsalfirst
116116
integrator.fsallast = fsallast
117117
initialize!(integrator, @inbounds(cache.caches[2]))
@@ -172,12 +172,12 @@ function choose_algorithm!(integrator,
172172
@inbounds if new_current != old_current
173173
cache.current = new_current
174174
if new_current == 1
175-
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[1],u)
175+
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[1], u)
176176
integrator.fsalfirst = fsalfirst
177177
integrator.fsallast = fsallast
178178
initialize!(integrator, @inbounds(cache.caches[1]))
179179
elseif new_current == 2
180-
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[2],u)
180+
fsalfirst, fsallast = get_fsalfirstlast(cache.caches[2], u)
181181
integrator.fsalfirst = fsalfirst
182182
integrator.fsallast = fsallast
183183
initialize!(integrator, @inbounds(cache.caches[2]))
@@ -205,37 +205,37 @@ function choose_algorithm!(integrator, cache::DefaultCache)
205205
cache.current = new_current
206206
init_ith_default_cache(cache, algs, new_current)
207207
if new_current == 1
208-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache1,u)
208+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache1, u)
209209
integrator.fsalfirst = fsalfirst
210210
integrator.fsallast = fsallast
211211
initialize!(integrator, @inbounds(cache.cache1))
212212
new_cache = cache.cache1
213213
elseif new_current == 2
214-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache2,u)
214+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache2, u)
215215
integrator.fsalfirst = fsalfirst
216216
integrator.fsallast = fsallast
217217
initialize!(integrator, @inbounds(cache.cache2))
218218
new_cache = cache.cache2
219219
elseif new_current == 3
220-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache3,u)
220+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache3, u)
221221
integrator.fsalfirst = fsalfirst
222222
integrator.fsallast = fsallast
223223
initialize!(integrator, @inbounds(cache.cache3))
224224
new_cache = cache.cache3
225225
elseif new_current == 4
226-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache4,u)
226+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache4, u)
227227
integrator.fsalfirst = fsalfirst
228228
integrator.fsallast = fsallast
229229
initialize!(integrator, @inbounds(cache.cache4))
230230
new_cache = cache.cache4
231231
elseif new_current == 5
232-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache5,u)
232+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache5, u)
233233
integrator.fsalfirst = fsalfirst
234234
integrator.fsallast = fsallast
235235
initialize!(integrator, @inbounds(cache.cache5))
236236
new_cache = cache.cache5
237237
elseif new_current == 6
238-
fsalfirst, fsallast = get_fsalfirstlast(cache.cache6,u)
238+
fsalfirst, fsallast = get_fsalfirstlast(cache.cache6, u)
239239
integrator.fsalfirst = fsalfirst
240240
integrator.fsallast = fsallast
241241
initialize!(integrator, @inbounds(cache.cache6))

lib/OrdinaryDiffEqCore/src/precompilation_setup.jl

+5-5
Original file line numberDiff line numberDiff line change
@@ -12,16 +12,16 @@ PrecompileTools.@compile_workload begin
1212
ODEProblem(lorenz, [1.0; 0.0; 0.0], (0.0, 1.0))
1313
ODEProblem(lorenz, [1.0; 0.0; 0.0], (0.0, 1.0), Float64[])
1414
ODEProblem{true, SciMLBase.AutoSpecialize}(lorenz, [1.0; 0.0; 0.0],
15-
(0.0, 1.0))
15+
(0.0, 1.0))
1616
ODEProblem{true, SciMLBase.AutoSpecialize}(lorenz, [1.0; 0.0; 0.0],
17-
(0.0, 1.0), Float64[])
17+
(0.0, 1.0), Float64[])
1818
ODEProblem{true, SciMLBase.FunctionWrapperSpecialize}(lorenz, [1.0; 0.0; 0.0],
19-
(0.0, 1.0))
19+
(0.0, 1.0))
2020
ODEProblem{true, SciMLBase.FunctionWrapperSpecialize}(lorenz, [1.0; 0.0; 0.0],
21-
(0.0, 1.0), Float64[])
21+
(0.0, 1.0), Float64[])
2222
ODEProblem{true, SciMLBase.NoSpecialize}(lorenz, [1.0; 0.0; 0.0], (0.0, 1.0))
2323
ODEProblem{true, SciMLBase.NoSpecialize}(lorenz, [1.0; 0.0; 0.0], (0.0, 1.0),
24-
Float64[])
24+
Float64[])
2525

2626
lorenz([1.0; 0.0; 0.0], [1.0; 0.0; 0.0], DiffEqBase.NullParameters(), 0.0)
2727
lorenz([1.0; 0.0; 0.0], [1.0; 0.0; 0.0], Float64[], 0.0)

lib/OrdinaryDiffEqCore/src/solve.jl

+1-1
Original file line numberDiff line numberDiff line change
@@ -469,7 +469,7 @@ function DiffEqBase.__init(
469469
reinitiailize = true
470470
saveiter = 0 # Starts at 0 so first save is at 1
471471
saveiter_dense = 0
472-
faslfirst, fsallast = get_fsalfirstlast(cache,rate_prototype)
472+
faslfirst, fsallast = get_fsalfirstlast(cache, rate_prototype)
473473

474474
integrator = ODEIntegrator{typeof(_alg), isinplace(prob), uType, typeof(du),
475475
tType, typeof(p),

0 commit comments

Comments
 (0)