Hi everybody,
I’m using the GREAT StochasticDiffEq.jl package to define a multidimensional driftless SDE as follows:
using StochasticDiffEq
function f!(du, u, p, t)
du .= p[1]
return nothing
end
function g!(du, u, p, t)
du[1] = p[3] * exp(-p[2] * t)
du[2] = p[5] * exp(-p[4] * t)
return nothing
end
u0 = [1.0, 0.3]
tspan = (0.0, 1.0)
params = ([0.0, 0.0], 0.2, 0.2, 0.2, 0.1)
sde = SDEProblem(f!, g!, u0, tspan, params)
Is there a way of defining the SDE in order to avoid the unnecessary drift computations?
Thanks in advance for your help!