PostGIS  3.1.6dev-r@@SVN_REVISION@@

◆ do_median_test()

static void do_median_test ( char *  input,
char *  expected,
int  fail_if_not_converged,
int  iter_count 
)
static

Definition at line 1387 of file cu_algorithm.c.

1388 {
1391  LWPOINT* expected_result = NULL;
1392  POINT4D actual_pt;
1393  POINT4D expected_pt;
1394  const double tolerance = FP_TOLERANCE / 10.0;
1395 
1396  LWPOINT* result = lwgeom_median(g, tolerance, iter_count, fail_if_not_converged);
1397  int passed = LW_FALSE;
1398 
1399  if (expected != NULL)
1400  {
1401  expected_result = lwgeom_as_lwpoint(lwgeom_from_wkt(expected, LW_PARSER_CHECK_NONE));
1402  lwpoint_getPoint4d_p(expected_result, &expected_pt);
1403  }
1404  if (result != NULL)
1405  {
1406  lwpoint_getPoint4d_p(result, &actual_pt);
1407  }
1408 
1409  if (result != NULL && expected != NULL) /* got something, expecting something */
1410  {
1411  passed = LW_TRUE;
1412  passed = passed && lwgeom_is_empty((LWGEOM*) expected_result) == lwgeom_is_empty((LWGEOM*) result);
1413  passed = passed && (lwgeom_has_z((LWGEOM*) expected_result) == lwgeom_has_z((LWGEOM*) result));
1414 
1415  if (passed && !lwgeom_is_empty((LWGEOM*) result))
1416  {
1417  if (g->type == POINTTYPE)
1418  {
1419  passed &= fabs(actual_pt.x - expected_pt.x) < tolerance;
1420  passed &= fabs(actual_pt.y - expected_pt.y) < tolerance;
1421  passed &= (!lwgeom_has_z((LWGEOM*) expected_result) || fabs(actual_pt.z - expected_pt.z) < tolerance);
1422  passed &= (!lwgeom_has_m((LWGEOM*) expected_result) || fabs(actual_pt.m - expected_pt.m) < tolerance);
1423  }
1424  else
1425  {
1426  /* Check that the difference between the obtained geometric
1427  median and the expected point is within tolerance */
1428  uint32_t npoints = 1;
1429  int input_empty = LW_TRUE;
1430  POINT4D* points = lwmpoint_extract_points_4d(lwgeom_as_lwmpoint(g), &npoints, &input_empty);
1431  double distance_expected = test_weighted_distance(&expected_pt, points, npoints);
1432  double distance_result = test_weighted_distance(&actual_pt, points, npoints);
1433 
1434  passed = distance_result <= (1.0 + tolerance) * distance_expected;
1435  if (!passed)
1436  {
1437  printf("Diff: Got %.10f Expected %.10f\n", distance_result, distance_expected);
1438  }
1439  lwfree(points);
1440  }
1441  }
1442 
1443  if (!passed)
1444  {
1445  printf("median_test input %s (parsed %s) expected %s got %s\n",
1446  input, lwgeom_to_ewkt(g),
1447  lwgeom_to_ewkt((LWGEOM*) expected_result),
1449  }
1450 
1451  }
1452  else if (result == NULL && expected == NULL) /* got nothing, expecting nothing */
1453  {
1454  passed = LW_TRUE;
1455  }
1456  else if (result != NULL && expected == NULL) /* got something, expecting nothing */
1457  {
1458  passed = LW_FALSE;
1459  printf("median_test input %s (parsed %s) expected NULL got %s\n", input, lwgeom_to_ewkt(g), lwgeom_to_ewkt((LWGEOM*) result));
1460  }
1461  else if (result == NULL && expected != NULL) /* got nothing, expecting something */
1462  {
1463  passed = LW_FALSE;
1464  printf("%s", cu_error_msg);
1465  printf("median_test input %s (parsed %s) expected %s got NULL\n", input, lwgeom_to_ewkt(g), lwgeom_to_ewkt((LWGEOM*) expected_result));
1466  }
1467 
1468  CU_ASSERT_TRUE(passed);
1469 
1470  lwgeom_free(g);
1471  lwpoint_free(expected_result);
1473 }
static double test_weighted_distance(const POINT4D *curr, const POINT4D *points, uint32_t npoints)
char result[OUT_DOUBLE_BUFFER_SIZE]
Definition: cu_print.c:267
void cu_error_msg_reset()
char cu_error_msg[MAX_CUNIT_ERROR_LENGTH+1]
int lwpoint_getPoint4d_p(const LWPOINT *point, POINT4D *out)
Definition: lwpoint.c:57
#define LW_FALSE
Definition: liblwgeom.h:108
LWMPOINT * lwgeom_as_lwmpoint(const LWGEOM *lwgeom)
Definition: lwgeom.c:225
void lwpoint_free(LWPOINT *pt)
Definition: lwpoint.c:213
void lwgeom_free(LWGEOM *geom)
Definition: lwgeom.c:1138
#define LW_PARSER_CHECK_NONE
Definition: liblwgeom.h:2085
int lwgeom_has_z(const LWGEOM *geom)
Return LW_TRUE if geometry has Z ordinates.
Definition: lwgeom.c:917
#define POINTTYPE
LWTYPE numbers, used internally by PostGIS.
Definition: liblwgeom.h:116
char * lwgeom_to_ewkt(const LWGEOM *lwgeom)
Return an alloced string.
Definition: lwgeom.c:548
void lwfree(void *mem)
Definition: lwutil.c:242
LWGEOM * lwgeom_from_wkt(const char *wkt, const char check)
Definition: lwin_wkt.c:905
LWPOINT * lwgeom_median(const LWGEOM *g, double tol, uint32_t maxiter, char fail_if_not_converged)
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:107
int lwgeom_has_m(const LWGEOM *geom)
Return LW_TRUE if geometry has M ordinates.
Definition: lwgeom.c:924
POINT4D * lwmpoint_extract_points_4d(const LWMPOINT *g, uint32_t *npoints, int *input_empty)
#define FP_TOLERANCE
Floating point comparators.
static int lwgeom_is_empty(const LWGEOM *geom)
Return true or false depending on whether a geometry is an "empty" geometry (no vertices members)
Definition: lwinline.h:203
static LWPOINT * lwgeom_as_lwpoint(const LWGEOM *lwgeom)
Definition: lwinline.h:131
uint8_t type
Definition: liblwgeom.h:476
double m
Definition: liblwgeom.h:428
double x
Definition: liblwgeom.h:428
double z
Definition: liblwgeom.h:428
double y
Definition: liblwgeom.h:428

References cu_error_msg, cu_error_msg_reset(), FP_TOLERANCE, LW_FALSE, LW_PARSER_CHECK_NONE, LW_TRUE, lwfree(), lwgeom_as_lwmpoint(), lwgeom_as_lwpoint(), lwgeom_free(), lwgeom_from_wkt(), lwgeom_has_m(), lwgeom_has_z(), lwgeom_is_empty(), lwgeom_median(), lwgeom_to_ewkt(), lwmpoint_extract_points_4d(), lwpoint_free(), lwpoint_getPoint4d_p(), POINT4D::m, POINTTYPE, result, test_weighted_distance(), LWGEOM::type, POINT4D::x, POINT4D::y, and POINT4D::z.

Referenced by test_median_robustness().

Here is the call graph for this function:
Here is the caller graph for this function: