From d1d59c43b5fdcccb64eec64605282bb843ae6334 Mon Sep 17 00:00:00 2001 From: Pierre Martinon Date: Thu, 21 Nov 2024 16:27:02 +0100 Subject: [PATCH] todo: check gauss legendre s=2 --- src/irk.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/irk.jl b/src/irk.jl index 0479a15..4b00f06 100644 --- a/src/irk.jl +++ b/src/irk.jl @@ -227,7 +227,7 @@ function setWorkArray(docp::DOCP{ <: GenericIRK}, xu, time_grid, v) if docp.dim_OCP_x == 1 xij = work[end] else - @views xij = work[end-docp.dim_OCP_x+1:end] + @views xij = work[end-docp.dim_OCP_x+1:end] #seems to displace the allocs somehow... end # dynamics for stage equation k_i^j = f(t_i^j, x_i^j, u_i, v) docp.ocp.dynamics((@view work[offset+1:offset+docp.dim_OCP_x]), tij, xij, ui, v)