|
199 | 199 | @muladd function perform_step!(integrator, cache::VerletLeapfrogConstantCache, |
200 | 200 | repeat_step = false) |
201 | 201 | @unpack t, dt, f, p = integrator |
202 | | - duprev, uprev, kduprev, kuprev = load_symp_state(integrator) |
| 202 | + duprev, uprev, kduprev, _ = load_symp_state(integrator) |
203 | 203 |
|
204 | | - # kick-drift-kick scheme of the Verlet Leapfrog method: |
| 204 | + # kick-drift-kick scheme of the Leapfrog method: |
205 | 205 | # update velocity |
206 | 206 | half = cache.half |
207 | 207 | du = duprev + dt * half * kduprev |
|
223 | 223 |
|
224 | 224 | @muladd function perform_step!(integrator, cache::VerletLeapfrogCache, repeat_step = false) |
225 | 225 | @unpack t, dt, f, p = integrator |
226 | | - duprev, uprev, kduprev, kuprev = load_symp_state(integrator) |
| 226 | + duprev, uprev, kduprev, _ = load_symp_state(integrator) |
227 | 227 | du, u, kdu, ku = alloc_symp_state(integrator) |
228 | 228 |
|
229 | 229 | # Kick-Drift-Kick scheme of the Verlet Leapfrog method: |
|
233 | 233 |
|
234 | 234 | # update position |
235 | 235 | tnew = t + half * dt |
236 | | - f.f2(ku, du, u, p, tnew) |
237 | | - @.. broadcast=false u=u + dt * ku |
| 236 | + f.f2(ku, du, uprev, p, tnew) |
| 237 | + @.. broadcast=false u=uprev + dt * ku |
238 | 238 |
|
239 | 239 | # update velocity |
240 | 240 | tnew = tnew + half * dt |
|
0 commit comments