Skip to content

Commit

Permalink
Update to match the OrdinaryDiffEq interpolation change
Browse files Browse the repository at this point in the history
  • Loading branch information
ChrisRackauckas committed Jan 7, 2024
1 parent 917c348 commit 065a1a9
Showing 1 changed file with 7 additions and 14 deletions.
21 changes: 7 additions & 14 deletions src/solve.jl
Original file line number Diff line number Diff line change
Expand Up @@ -168,20 +168,13 @@ function DiffEqBase.__init(prob::DiffEqBase.AbstractDDEProblem,
stats = DDEStats(0)

# create solution
if iscomposite(alg)
id = OrdinaryDiffEq.CompositeInterpolationData(f_with_history, timeseries, ts, ks,
Int[], dense, cache, differential_vars)
sol = DiffEqBase.build_solution(prob, alg.alg, ts, timeseries;
dense = dense, k = ks, interp = id,
alg_choice = id.alg_choice, calculate_error = false,
stats = stats)
else
id = OrdinaryDiffEq.InterpolationData(f_with_history, timeseries, ts, ks, dense,
cache, differential_vars)
sol = DiffEqBase.build_solution(prob, alg.alg, ts, timeseries;
dense = dense, k = ks, interp = id,
calculate_error = false, stats = stats)
end
alg_choice = iscomposite(alg) ? Int[] : nothing
id = OrdinaryDiffEq.CompositeInterpolationData(f_with_history, timeseries, ts, ks,
alg_choice, dense, cache, differential_vars, false)
sol = DiffEqBase.build_solution(prob, alg.alg, ts, timeseries;
dense = dense, k = ks, interp = id,
alg_choice = id.alg_choice, calculate_error = false,
stats = stats)

# retrieve time stops, time points at which solutions is saved, and discontinuities
tstops_internal = OrdinaryDiffEq.initialize_tstops(tType, tstops, d_discontinuities,
Expand Down

0 comments on commit 065a1a9

Please sign in to comment.