Skip to content

Commit 4957a10

Browse files
committed
use appropriate p, t in jacobian
1 parent 6f02500 commit 4957a10

File tree

2 files changed

+5
-3
lines changed

2 files changed

+5
-3
lines changed

Project.toml

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,7 +1,7 @@
11
name = "DynamicalSystemsBase"
22
uuid = "6e36e845-645a-534a-86f2-f5d4aa5a06b4"
33
repo = "https://github.com/JuliaDynamics/DynamicalSystemsBase.jl.git"
4-
version = "3.15.2"
4+
version = "3.15.3"
55

66
[deps]
77
ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"

src/core_systems/jacobian.jl

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -47,8 +47,10 @@ end
4747
function _jacobian(ds, ::Val{true})
4848
f = dynamic_rule(ds)
4949
u0 = current_state(ds)
50+
p = current_parameters(ds)
51+
t = initial_time(ds)
5052
cfg = ForwardDiff.JacobianConfig(
51-
(du, u) -> f(du, u, p, p), deepcopy(u0), deepcopy(u0)
53+
(du, u) -> f(du, u, p, t), deepcopy(u0), deepcopy(u0)
5254
)
5355
Jf! = (J0, u, p, t = 0) -> begin
5456
uv = @view u[:, 1]
@@ -63,7 +65,7 @@ end
6365

6466
function _jacobian(ds, ::Val{false})
6567
f = dynamic_rule(ds)
66-
Jf = (u, p, t = 0) -> ForwardDiff.jacobian((x) -> f(x, p, t), u)
68+
Jf = (u, p, t = 0) -> ForwardDiff.jacobian(x -> f(x, p, t), u)
6769
return Jf
6870
end
6971

0 commit comments

Comments
 (0)