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
New to julia: Convolution and juliadiff Numerics question	19	6344	August 21, 2017
How to find second and third derivative using ForwardDiff New to Julia	2	2213	February 19, 2019
ReverseDiff.jl with iterative solvers General Usage	14	1469	June 23, 2017
Simple numerical differentiation New to Julia question , forwarddiff , finitediff	5	4685	December 23, 2020
ForwardDiff failing on simple polynomial New to Julia	1	954	September 28, 2017