@@ -59,12 +59,13 @@ generated code via:
59
59
generate_function (de)
60
60
61
61
# # Which returns:
62
- :((du, u, p, t)-> begin
63
- du .= let (x, y, z, σ, ρ, β) = (u[1 ], u[2 ], u[3 ], p[1 ], p[2 ], p[3 ])
64
- (σ * (y - x), x * (ρ - z) - y, x * y - β * z)
65
- end
62
+ :((# #363, u, p, t)->begin
63
+ let (x, y, z, σ, ρ, β) = (u[1 ], u[2 ], u[3 ], p[1 ], p[2 ], p[3 ])
64
+ # #363[1] = σ * (y - x)
65
+ # #363[2] = x * (ρ - z) - y
66
+ # #363[3] = x * y - β * z
67
+ end
66
68
end )
67
-
68
69
```
69
70
70
71
and get the generated function via:
@@ -94,10 +95,12 @@ nlsys_func = generate_function(ns)
94
95
which generates:
95
96
96
97
``` julia
97
- :((du, u, p)-> begin
98
- du .= let (y, x, z, σ, ρ, β) = (u[1 ], u[2 ], u[3 ], p[1 ], p[2 ], p[3 ])
99
- (σ * (y - x), x * (ρ - z) - y, x * y - β * z)
100
- end
98
+ :((# #366, u, p, t)->begin
99
+ let (y, z, x, ρ, σ, β) = (u[1 ], u[2 ], u[3 ], p[1 ], p[2 ], p[3 ])
100
+ # #366[1] = σ * (y - x)
101
+ # #366[2] = x * (ρ - z) - y
102
+ # #366[3] = x * y - β * z
103
+ end
101
104
end )
102
105
```
103
106
@@ -276,12 +279,13 @@ nlsys_func = generate_function(ns)
276
279
expands to:
277
280
278
281
``` julia
279
- :((du, u, p)-> begin
280
- du .= let (x, y, z, σ, ρ, β) = (u[1 ], u[2 ], u[3 ], p[1 ], p[2 ], p[3 ])
281
- (σ * (y - x), x * (ρ - z) - y, x * y - β * z)
282
- end
282
+ :((# #367, u, p, t)->begin
283
+ let (x, y, z, σ, ρ, β) = (u[1 ], u[2 ], u[3 ], p[1 ], p[2 ], p[3 ])
284
+ # #367[1] = σ * (y - x)
285
+ # #367[2] = x * (ρ - z) - y
286
+ # #367[3] = x * y - β * z
287
+ end
283
288
end )
284
-
285
289
```
286
290
287
291
In addition, the Jacobian calculations take into account intermediate variables
0 commit comments