141 bool use_spheroid =
true;
145 g1 = PG_GETARG_GSERIALIZED_P(0);
146 g2 = PG_GETARG_GSERIALIZED_P(1);
149 if ( PG_NARGS() > 2 && ! PG_ARGISNULL(2) )
150 tolerance = PG_GETARG_FLOAT8(2);
153 if ( PG_NARGS() > 3 && ! PG_ARGISNULL(3) )
154 use_spheroid = PG_GETARG_BOOL(3);
162 if ( ! use_spheroid )
163 s.a =
s.b =
s.radius;
171 PG_FREE_IF_COPY(g1, 0);
172 PG_FREE_IF_COPY(g2, 1);
182 POSTGIS_DEBUGF(2,
"[GIST] '%s' got distance %g", __func__,
distance);
187 PG_FREE_IF_COPY(g1, 0);
188 PG_FREE_IF_COPY(g2, 1);
int32_t gserialized_get_srid(const GSERIALIZED *s)
Extract the SRID from the serialized form (it is packed into three bytes so this is a handy function)...
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
double lwgeom_distance_spheroid(const LWGEOM *lwgeom1, const LWGEOM *lwgeom2, const SPHEROID *spheroid, double tolerance)
Calculate the geodetic distance from lwgeom1 to lwgeom2 on the spheroid.
void lwgeom_free(LWGEOM *geom)
int lwgeom_is_empty(const LWGEOM *geom)
Return true or false depending on whether a geometry is an "empty" geometry (no vertices members)
void lwgeom_add_bbox_deep(LWGEOM *lwgeom, GBOX *gbox)
Compute a box for geom and all sub-geometries, if not already computed.
void error_if_srid_mismatch(int srid1, int srid2)
#define FP_TOLERANCE
Floating point comparators.
Datum distance(PG_FUNCTION_ARGS)