Dear all, I am successfully playing with sol = solve(ensembleprob,SRIW1(),EnsembleThreads(),trajectories=ϕ,saveat=δt); which generates a nice set of values. Among those values sol[1,:,:]; is a position, sol[2,:,:]; a velocity. unfortunately, I need an acceleration… how can I nicely proceed ? t…

sadly ended with a five-point stencil scheme : for i in 3:(lmnt-2) kt_a[i] = (-kt[i+2]+8*kt[i+1]-8*kt[i-1]+kt[i-2])/(12*δt) end; :roll_eyes:

Derivative for a set of ... points

rafael.guerra March 16, 2021, 2:36pm 3

@Frederic_Gilt, you may want to take a look at this discourse thread.

Topic		Replies	Views
Simple numerical differentiation New to Julia question , forwarddiff , finitediff	5	4561	December 23, 2020
Problem with ForardDiff New to Julia package , forwarddiff	1	251	June 8, 2022
Derivative of user-defined function with variables as arguments in a PDE system General Usage question , package , diffeq , modelingtoolkit , methodoflines	1	276	March 21, 2024
How to get sol.du in DifferentialEquations.jl General Usage diffeq	3	1346	February 23, 2019
Vectorization of ForwardDiff.gradient function New to Julia	7	1839	October 29, 2017