We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
step
1 parent e485198 commit c369c9eCopy full SHA for c369c9e
src/integrator.jl
@@ -218,10 +218,8 @@ function step(
218
ϵ = fwd ? step_size(lf) : -step_size(lf)
219
ϵ = ϵ'
220
221
- res = if FullTraj
222
- Vector{P}(undef, n_steps)
223
- else
224
- Vector{P}(undef, 1)
+ if FullTraj
+ res = Vector{P}(undef, n_steps)
225
end
226
227
(; θ, r) = z
@@ -244,20 +242,18 @@ function step(
244
242
# Update result
245
243
if FullTraj
246
res[i] = z
247
248
- res[1] = z
249
250
if !isfinite(z)
251
# Remove undef
252
253
- res = res[isassigned.(Ref(res), 1:n_steps)]
+ resize!(res, i)
254
255
break
256
257
258
return if FullTraj
259
res
260
else
261
- first(res)
+ z
262
263
0 commit comments