Small fix to function edge_distance_to_point

[ezimanyi]

The third parameter closest can be NULL. This is tested at the end of the function but not at the beginning. The code below is extracted from PostGIS 3.1.2.

double edge_distance_to_point(const GEOGRAPHIC_EDGE *e, const GEOGRAPHIC_POINT *gp, GEOGRAPHIC_POINT *closest)
{
  double d1 = 1000000000.0, d2, d3, d_nearest;
  POINT3D n, p, k;
  GEOGRAPHIC_POINT gk, g_nearest;

  /* Zero length edge, */
  if ( geographic_point_equals(&(e->start), &(e->end)) )
  {
    if (closest) <----------------- Missing test
      *closest = e->start;
    return sphere_distance(&(e->start), gp);
  }

  [...]

  if (closest)
    *closest = g_nearest;

  return d_nearest;
}

Regards

Esteban

[Paul Ramsey <pramsey@…>]

In d3f1c2b5/git:

Add guard to return pointer, closes #4953