Finding quadrants is flawed #622
Description
As pointed out here #621 (comment), the function _quadrant used internally for trigonometric functions may yield wrong results.
See #621 (comment) for a proposal to fix the issue (at least partially).
Also, this issue can be related to #129 on the implementation of mod for intervals (see also PR #145 and #178).