Frédéric BRIOL 480e9cf0f6 fix(formulas): geographic intersection calculation in sjoberg_intersection when a segment lies on the equator ...

Multiply the longitude offset term by `sign_lon_diff` when calculating
equator crossings in `sjoberg_geodesic::lon_of_equator_intersection()`.

Without this multiplication, when `sign_lon_diff` is negative, the offset
is applied in the wrong direction, leading to a wrong intersection
longitude (e.g. returning a coordinate outside of the segments' bounds).
This in turn causes overlay/intersection algorithms to miss crossings on
the boundary.

Including regression test for sjoberg_intersection

Fixes #1482

2026-07-09 20:04:41 +02:00

..

boost

fix(formulas): geographic intersection calculation in sjoberg_intersection when a segment lies on the equator

2026-07-09 20:04:41 +02:00