296{
298 GBOX gbox, gbox_slow;
299 int i;
300
302 {
303#if 0
304
305 printf("\n\n------------\n");
307#endif
316#if 0
319 printf(
"line %d: diff %.9g\n", i, fabs(gbox.
xmin - gbox_slow.
xmin)+fabs(gbox.
ymin - gbox_slow.
ymin)+fabs(gbox.
zmin - gbox_slow.
zmin));
320 printf("------------\n");
321#endif
322 CU_ASSERT_DOUBLE_EQUAL(gbox.
xmin, gbox_slow.
xmin, 0.00000001);
323 CU_ASSERT_DOUBLE_EQUAL(gbox.
ymin, gbox_slow.
ymin, 0.00000001);
324 CU_ASSERT_DOUBLE_EQUAL(gbox.
zmin, gbox_slow.
zmin, 0.00000001);
325 CU_ASSERT_DOUBLE_EQUAL(gbox.
xmax, gbox_slow.
xmax, 0.00000001);
326 CU_ASSERT_DOUBLE_EQUAL(gbox.
ymax, gbox_slow.
ymax, 0.00000001);
327 CU_ASSERT_DOUBLE_EQUAL(gbox.
zmax, gbox_slow.
zmax, 0.00000001);
328 }
329
330}
char * gbox_to_string(const GBOX *gbox)
Allocate a string representation of the GBOX, based on dimensionality of flags.
void lwgeom_free(LWGEOM *geom)
#define LW_PARSER_CHECK_NONE
int lwgeom_calculate_gbox(const LWGEOM *lwgeom, GBOX *gbox)
Calculate bounding box of a geometry, automatically taking into account whether it is cartesian or ge...
#define FLAGS_SET_GEODETIC(flags, value)
#define LW_TRUE
Return types for functions with status returns.
LWGEOM * lwgeom_from_wkt(const char *wkt, const char check)
int gbox_geocentric_slow
For testing geodetic bounding box, we have a magic global variable.