Problem: The symmetric difference of two geometry collections is
calculated as  difference(union(A, B), intersection(A, B)), and this
may cause some numerical robustness issues.

Fix: Rewrite the operation to instead calculate union(difference(A,
B), difference(B, A)) in order to make the algorithm more numerically
robust.