Support `is_interior_point` for mixed numeric types

[schillic]

I decided to also require that the tolerance is non-zero. That makes this a breaking change. I can undo that part if not desired.