change variable name

Author: Hadrien

src/caches/dynamical_caches.jl

@@ -125,7 +125,7 @@ end
 
 mutable struct OBABOConstantCache{uType,rateNoiseType, uEltypeNoUnits, uCoeffType, uCoeffMType} <: StochasticDynamicalEqConstantCache
   k::uType
-  sig::rateNoiseType
+  gt::rateNoiseType
   half::uEltypeNoUnits
   c₂::uCoeffType
   σ::uCoeffMType
@@ -147,7 +147,7 @@ end
 
 function alg_cache(alg::OBABO,prob,u,ΔW,ΔZ,p,rate_prototype,noise_rate_prototype,jump_rate_prototype,::Type{uEltypeNoUnits},::Type{uBottomEltypeNoUnits},::Type{tTypeNoUnits},uprev,f,t,dt,::Type{Val{false}}) where {uEltypeNoUnits,uBottomEltypeNoUnits,tTypeNoUnits}
   k = zero(rate_prototype.x[1])
-  sig = zero(noise_rate_prototype)
+  gt = zero(noise_rate_prototype)
   half=uEltypeNoUnits(1//2)
 
   if typeof(alg.gamma) <: AbstractMatrix
@@ -158,7 +158,7 @@ function alg_cache(alg::OBABO,prob,u,ΔW,ΔZ,p,rate_prototype,noise_rate_prototy
       σ = sqrt.(1 .- alg.scale_noise*c₂.^2)
   end
    # if scale_noise == false, c2 = 1
-  OBABOConstantCache(k, sig, half, c₂, σ)
+  OBABOConstantCache(k, gt, half, c₂, σ)
 end
 
 function alg_cache(alg::OBABO,prob,u,ΔW,ΔZ,p,rate_prototype, noise_rate_prototype,jump_rate_prototype,::Type{uEltypeNoUnits},::Type{uBottomEltypeNoUnits},::Type{tTypeNoUnits},uprev,f,t,dt,::Type{Val{true}}) where {uEltypeNoUnits,uBottomEltypeNoUnits,tTypeNoUnits}

src/perform_step/dynamical.jl

@@ -160,7 +160,7 @@ function initialize!(integrator, cache::OBABOConstantCache)
 
   verify_f2(integrator.f.f2, du1, u1, p, t, integrator, cache)
   cache.k = integrator.f.f1(du1,u1,p,t)
-  cache.sig = integrator.g(u1,p,t)
+  cache.gt = integrator.g(u1,p,t)
 end
 
 function initialize!(integrator, cache::OBABOCache)
@@ -181,7 +181,7 @@ end
   u1 = uprev.x[2]
 
   # O
-  noise = cache.sig.*W.dW ./ sqdt
+  noise = cache.gt.*W.dW ./ sqdt
   if typeof(σ) <: AbstractMatrix || typeof(noise) <: Number
     du2 = c₂*du1 + σ*noise
   else
@@ -199,8 +199,8 @@ end
   du3 = dumid + half*dt*cache.k
 
   # O
-  cache.sig = integrator.g(u,p,t+dt)
-  noise = cache.sig.*W.dW ./ sqdt # That should be a second noise
+  cache.gt = integrator.g(u,p,t+dt)
+  noise = cache.gt.*W.dW ./ sqdt # That should be a second noise
   if typeof(σ) <: AbstractMatrix || typeof(noise) <: Number
     du = c₂*du3 + σ*noise
   else