@@ -5,14 +5,14 @@ function timestep!(method::ExplicitRungeKuttaMethod, stepper, Δt; θ = nothing,
55 (; setup, psolver, u, temp, t, n) = stepper
66 (; closure_model, temperature) = setup
77 (; A, b, c) = method
8- (; uprev , ku, p, tempprev , ktemp, diff) = cache
8+ (; ustart , ku, p, tempstart , ktemp, diff) = cache
99 nstage = length (b)
1010 m = closure_model
1111
1212 # Update current solution
13- tprev = t
14- copyto! (uprev , u)
15- isnothing (temp) || copyto! (tempprev , temp)
13+ tstart = t
14+ copyto! (ustart , u)
15+ isnothing (temp) || copyto! (tempstart , temp)
1616
1717 for i = 1 : nstage
1818 # Compute force at current stage i
@@ -29,15 +29,15 @@ function timestep!(method::ExplicitRungeKuttaMethod, stepper, Δt; θ = nothing,
2929 isnothing (m) || (ku[i] .+ = m (u, θ))
3030
3131 # Intermediate time step
32- t = tprev + c[i] * Δt
32+ t = tstart + c[i] * Δt
3333
3434 # Apply stage forces
35- u .= uprev
35+ u .= ustart
3636 for j = 1 : i
3737 @. u += Δt * A[i, j] * ku[j]
3838 end
3939 if ! isnothing (temp)
40- temp .= tempprev
40+ temp .= tempstart
4141 for j = 1 : i
4242 @. temp += Δt * A[i, j] * ktemp[j]
4343 end
@@ -66,8 +66,8 @@ function timestep(method::ExplicitRungeKuttaMethod, stepper, Δt; θ = nothing)
6666 m = closure_model
6767
6868 # Update current solution (does not depend on previous step size)
69- tprev = t
70- uprev = u
69+ tstart = t
70+ ustart = u
7171 ku = ()
7272 ktemp = ()
7373
@@ -89,15 +89,15 @@ function timestep(method::ExplicitRungeKuttaMethod, stepper, Δt; θ = nothing)
8989 isnothing (temp) || (ktemp = (ktemp... , Ftemp))
9090
9191 # Intermediate time step
92- t = tprev + c[i] * Δt
92+ t = tstart + c[i] * Δt
9393
9494 # Apply stage forces
95- u = uprev
95+ u = ustart
9696 for j = 1 : i
9797 u = @. u + Δt * A[i, j] * ku[j]
9898 end
9999 if ! isnothing (temp)
100- temp = tempprev
100+ temp = tempstart
101101 for j = 1 : i
102102 temp = @. temp + Δt * A[i, j] * ktemp[j]
103103 end
0 commit comments