Better implementation of `xy_ordered_union` by hashing (x,y) coordinates

[beddalumia]

Refer to this discussion on fortran-lang.

• This would greatly increase efficiency, so it is something we surely want to do at some point...

• For now the priority is to get a minimal implementation working asap, so I just defer this to a future patch.