searchsortedfirst of dual returning wrong result

[oscardssmith]

This is presumably caused by [email protected] and is causing test failures in SciML/DataInterpolations.jl#414

julia> searchsortedfirst([0.0, 1.0, 2.0, 3.0], 1.0)
2

julia> searchsortedfirst([0.0, 1.0, 2.0, 3.0], ForwardDiff.Dual(1.0, (1.0,)))
3

[devmotion]

This is both expected and (IMO) intended: a Dual with positive partial value is strictly greater than one with the same primal but zero partial value.

[devmotion]

The change is caused by #695 which fixed the broken isless etc definitions for Duals.

[oscardssmith]

Thanks for the quick response. I agree that DataInterpolations seems to be doing some fairly sketchy things here...

[ChrisRackauckas]

The DataInterpolations derivative implementation needs to take this into account based on the continuity argument. We can take it to that repo.