Missing precision in the solution of combinations of ODE and jump processes using DifferentialEquations.jl and JumpProcesses.jl

Intuitive depends on who you’re talking to. It’s simply undefined because there’s two values of the solution at that t. Undefined can’t have a “better/more precise” answer since both are true. Given that it’s undefined, there’s an argument continuity for controlling which solution you get.

You can use sol(t, continuity = :right) if you want the right-continuous value. There are just as many cases (well, arguably more) where this choice is less intuitive which is why the default is set as it is.