do_median_test()

static void do_median_test ( char *  input, char *  expected  )

static

Definition at line 1022 of file cu_algorithm.c.

References FP_EQUALS, FP_TOLERANCE, LW_PARSER_CHECK_NONE, LW_TRUE, lwgeom_as_lwpoint(), lwgeom_free(), lwgeom_from_wkt(), lwgeom_has_z(), lwgeom_is_empty(), lwgeom_median(), lwgeom_to_ewkt(), lwpoint_free(), lwpoint_getPoint3dz_p(), POINT3DZ::x, POINT3DZ::y, and POINT3DZ::z.

Referenced by test_median_robustness().

{
        LWGEOM* g = lwgeom_from_wkt(input, LW_PARSER_CHECK_NONE);
        LWPOINT* expected_result = lwgeom_as_lwpoint(lwgeom_from_wkt(expected, LW_PARSER_CHECK_NONE));
        POINT3DZ actual_pt;
        POINT3DZ expected_pt;

        LWPOINT* result = lwgeom_median(g, FP_TOLERANCE / 10.0, 1000, LW_TRUE);
        int passed = LW_TRUE;

        lwpoint_getPoint3dz_p(result, &actual_pt);
        lwpoint_getPoint3dz_p(expected_result, &expected_pt);

        passed = passed && lwgeom_is_empty((LWGEOM*) expected_result) == lwgeom_is_empty((LWGEOM*) result);
        passed = passed && (lwgeom_has_z((LWGEOM*) expected_result) == lwgeom_has_z((LWGEOM*) result));

        if (!lwgeom_is_empty((LWGEOM*) result))
        {
                passed = passed && FP_EQUALS(actual_pt.x, expected_pt.x);
                passed = passed && FP_EQUALS(actual_pt.y, expected_pt.y);
                passed = passed && (!lwgeom_has_z((LWGEOM*) expected_result) || FP_EQUALS(actual_pt.z, expected_pt.z));
        }

        if (!passed)
                printf("median_test expected %s got %s\n", lwgeom_to_ewkt((LWGEOM*) expected_result), lwgeom_to_ewkt((LWGEOM*) result));

        CU_ASSERT_TRUE(passed);

        lwgeom_free(g);
        lwpoint_free(expected_result);
        lwpoint_free(result);
}

Here is the call graph for this function:

Here is the caller graph for this function: