Update solve.jl

Author: ChrisRackauckas

src/solve.jl

@@ -2,6 +2,8 @@
 # when the solver (i.e. StochasticDiffEq.jl) allows for jumps in __init
 struct ForceJumpDispatch end
 
+
+
 function DiffEqBase.__solve(jump_prob::DiffEqBase.AbstractJumpProblem{P},
         alg::DiffEqBase.DEAlgorithm;
         kwargs...) where {P}
@@ -10,6 +12,15 @@ function DiffEqBase.__solve(jump_prob::DiffEqBase.AbstractJumpProblem{P},
     integrator.sol
 end
 
+#Ambiguity Fix
+function DiffEqBase.__solve(jump_prob::DiffEqBase.AbstractJumpProblem{P},
+        alg::Union{SciMLBase.AbstractRODEAlgorithm, SciMLBase.AbstractSDEAlgorithm};
+        kwargs...) where {P}
+    integrator = DiffEqBase.__init(jump_prob, alg, ForceJumpDispatch(); kwargs...)
+    solve!(integrator)
+    integrator.sol
+end
+
 # if passed a JumpProblem over a DiscreteProblem, and no aggregator is selected use
 # SSAStepper
 function DiffEqBase.__solve(jump_prob::DiffEqBase.AbstractJumpProblem{P};