Skip to content

Commit 82c1cc8

Browse files
author
Hadrien
committed
Actually OBABO need 2 noise processes
1 parent 2fe4848 commit 82c1cc8

File tree

2 files changed

+3
-2
lines changed

2 files changed

+3
-2
lines changed

src/alg_utils.jl

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -263,6 +263,7 @@ alg_needs_extra_process(alg::RS2) = true
263263
alg_needs_extra_process(alg::PL1WM) = true
264264
alg_needs_extra_process(alg::NON) = true
265265
alg_needs_extra_process(alg::NON2) = true
266+
alg_needs_extra_process(alg::OBABO) = true
266267

267268
OrdinaryDiffEq._alg_autodiff(alg::StochasticDiffEqNewtonAlgorithm{CS,AD,FDT,ST,CJ,Controller}) where {CS,AD,FDT,ST,CJ,Controller} = Val{AD}()
268269
OrdinaryDiffEq._alg_autodiff(alg::StochasticDiffEqNewtonAdaptiveAlgorithm{CS,AD,FDT,ST,CJ,Controller}) where {CS,AD,FDT,ST,CJ,Controller} = Val{AD}()

src/perform_step/dynamical.jl

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -200,7 +200,7 @@ end
200200

201201
# O
202202
cache.gt = integrator.g(u,p,t+dt)
203-
noise = cache.gt.*W.dW ./ sqdt # That should be a second noise
203+
noise = cache.gt.*W.dZ ./ sqdt # That should be a second noise
204204
if typeof(σ) <: AbstractMatrix || typeof(noise) <: Number
205205
du = c₂*du3 + σ*noise
206206
else
@@ -242,7 +242,7 @@ end
242242

243243
# O
244244
integrator.g(gtmp,u.x[2],p,t+dt)
245-
@.. noise = gtmp*W.dW / sqdt # That should be a second noise
245+
@.. noise = gtmp*W.dZ / sqdt # That should be a second noise
246246

247247
if typeof(σ) <: AbstractMatrix
248248
mul!(u.x[1],c₂,dutmp)

0 commit comments

Comments
 (0)