diff --git a/src/derivatives.jl b/src/derivatives.jl index 2cc5d33e..baf0d920 100644 --- a/src/derivatives.jl +++ b/src/derivatives.jl @@ -20,8 +20,6 @@ function _extrapolate_derivative_left(A, t, order) elseif extrapolation_left == ExtrapolationType.Constant zero(first(A.u) / one(A.t[1])) elseif extrapolation_left == ExtrapolationType.Linear - _derivative(A, first(A.t), 1) - zero(first(A.u) / one(A.t[1])) (order == 1) ? _derivative(A, first(A.t), 1) : zero(first(A.u) / one(A.t[1])) elseif extrapolation_left == ExtrapolationType.Extension (order == 1) ? _derivative(A, t, length(A.t)) : @@ -66,7 +64,7 @@ function _extrapolate_derivative_right(A, t, order) _derivative(A, t, A.iguesser) end, t_) else - # extrapolation_left == ExtrapolationType.Reflective + # extrapolation_right == ExtrapolationType.Reflective t_, n = transformation_reflective(A, t) sign = iseven(n) ? -1 : 1 (order == 1) ? sign * _derivative(A, t_, A.iguesser) : @@ -314,4 +312,4 @@ function _derivative( out += Δt₀^2 * (3c₁ + (3Δt₁ + Δt₀) * c₂ + (3Δt₁^2 + Δt₀ * 2Δt₁) * c₃) out -end \ No newline at end of file +end diff --git a/src/integral_inverses.jl b/src/integral_inverses.jl index 22fde725..075da29f 100644 --- a/src/integral_inverses.jl +++ b/src/integral_inverses.jl @@ -119,4 +119,4 @@ function _interpolate( idx_ = get_idx(A, t, idx; side = :first, lb = 1, ub_shift = 0) end A.u[idx] + (t - A.t[idx]) / A.itp.u[idx_] -end \ No newline at end of file +end diff --git a/src/interpolation_methods.jl b/src/interpolation_methods.jl index f8cdce89..73e240d1 100644 --- a/src/interpolation_methods.jl +++ b/src/interpolation_methods.jl @@ -335,4 +335,4 @@ function _interpolate( c₁, c₂, c₃ = get_parameters(A, idx) out += Δt₀^3 * (c₁ + Δt₁ * (c₂ + c₃ * Δt₁)) out -end \ No newline at end of file +end