Merge pull request #96 from ericagol/fix_transit

Fixing starting guess for transit time
ericagol · Mar 16, 2023 · 3402d0b · 3402d0b
2 parents ca63ded + 08cc7e4
commit 3402d0b
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/transits/timing.jl b/src/transits/timing.jl
@@ -17,7 +17,7 @@ function detect_transits!(s::State{T},d::Derivatives{T},tt::TransitOutput{T},int
             # A transit has occurred between the time steps - integrate ahl21!
             tt.count[i] += 1
             if tt.count[i] <= tt.ntt
-                dt0 = -tt.gsave[i]*intr.h/(gi-tt.gsave[i]) # Starting estimate
+                dt0 = -gi*intr.h/(gi-tt.gsave[i]) # Starting estimate
                 set_state!(s,tt.s_prior)
                 findtransit!(tt.ti,i,dt0,s,d,tt,intr;grad=grad) # Search for transit time (integrating 'backward')
             end