Overlapping lines

[bjourne]

The algorithm does not handle the case in which lines are overlapping. Consider, f.e:

print(list(line_intersections([((0,0),(10,0)),((5,0),(15,0))])))

[adrianN]

There are probably many corner cases like this. The way I compute cross products is also not right because I use normal floats, so rounding errors can change the result.

It should work if the points are in General position. Typically you can add a small random perturbation to get to general position without changing the result significantly.

If you want proper battle tested code you should check out CGAL

I accept pull requests if you want handle some of the corner cases you encounter.