Support for GEOS 3.13 and RelateNG.

Most functionality remains the same, but new GEOS predicate implementation has a few small changes.

- Boundary Node Rule relate matrices might be different when
  using the "multi-valent end point" rule.
- Relate matrices for situations with invalid MultiPolygons
  with shared boundaries might be different. Run ST_MakeValid
  to get valid inputs to feed to the calculation.
- Zero length LineStrings are treated as if they are the
  equivalent Point object.

References #5770