Handle edge-cases of inverted root pairs in the condition rootfinding

[DaniGlez]

Closes #1290 (see the issue discussion for context)

Checklist

• Appropriate tests were added (requires Make ModAB iteration always shrink the bracket even if there's numerical noise NonlinearSolve.jl#861 to be merged)
• Any code changes were done in a way that does not break public API
• All documentation related to code changes were updated
• The new code follows the
  contributor guidelines, in particular the SciML Style Guide and
  COLPRAC.
• Any new documentation only uses public API

Additional context

Add any other context about the problem here.

[DaniGlez]

Also please squash if merging, I was a bit sloppy with the WIP commits 😅

[ChrisRackauckas]

Needs a test

[DaniGlez]

Test added; I would still like SciML/NonlinearSolve.jl#861 to be merged first, so that the ModAB behavior is consistent (if it hits the buggy path, a slightly different original root might be found)

[DaniGlez]

Fuzzer script for posterity, in case the test needs to be recreated in the future: https://gist.github.com/DaniGlez/6100f2530fc76d4f9673a8f9f7422564

[DaniGlez]

I think the JumpProcesses test failure should be solved by SciML/NonlinearSolve.jl#861 as well