diff --git a/src/integrals.jl b/src/integrals.jl index 5e2e5907..de1b4633 100644 --- a/src/integrals.jl +++ b/src/integrals.jl @@ -74,7 +74,7 @@ function _integral(A::QuadraticSpline{<:AbstractVector{<:Number}}, idx::Number, α, β = get_parameters(A, idx) uᵢ = A.u[idx] Δt = t - A.t[idx] - Δt_full = (A.t[idx + 1] - A.t[idx]) + Δt_full = A.t[idx + 1] - A.t[idx] Δt * (α * Δt^2 / (3Δt_full^2) + β * Δt / (2Δt_full) + uᵢ) end diff --git a/src/interpolation_utils.jl b/src/interpolation_utils.jl index 748cf662..2d1c6432 100644 --- a/src/interpolation_utils.jl +++ b/src/interpolation_utils.jl @@ -81,9 +81,9 @@ function quadratic_spline_params(t::AbstractVector, sc::AbstractVector) k[(end - 2):end] .= t[end] k[4:(end - 3)] .= t[2:(end - 2)] - # Create and solve linear system Ac = u, where: + # Create linear system Ac = u, where: # - A consists of basis function evaulations in t - # - c are the 1D control points + # - c are 1D control points n = length(t) dtype_sc = typeof(t[1] / t[1])