test_geography_substring()

static void test_geography_substring ( void  )

static

Definition at line 1597 of file cu_geodetic.c.

{
        LWGEOM *lwg;
        LWGEOM *result;
        double length;
        SPHEROID s;
 
        /* Init to Sphere */
        spheroid_init(&s, WGS84_RADIUS, WGS84_RADIUS);
        s.a = s.b = s.radius;
 
        /*
         * geography_substring(
         *   const LWLINE *lwline,
         *   const SPHEROID *s,
         *   double from, double to, double tolerance
         */
 
        lwg = lwgeom_from_wkt("LINESTRING(2 48, 2 49, 2 50)", LW_PARSER_CHECK_NONE);
        length = lwgeom_length(lwg);
        CU_ASSERT_DOUBLE_EQUAL(length, 2.0, 0.001);
        result = geography_substring((LWLINE*)lwg, &s, 0.0, 1.0, 0.01);
        length = lwgeom_length(result);
        CU_ASSERT_DOUBLE_EQUAL(length, 2.0, 0.001);
        lwgeom_free(lwg);
        lwgeom_free(result);
 
        lwg = lwgeom_from_wkt("LINESTRING(2 48, 2 49)", LW_PARSER_CHECK_NONE);
        length = lwgeom_length(lwg);
        CU_ASSERT_DOUBLE_EQUAL(length, 1.0, 0.001);
        result = geography_substring((LWLINE*)lwg, &s, 0.0, 1.0, 0.01);
        length = lwgeom_length(result);
        CU_ASSERT_DOUBLE_EQUAL(length, 1.0, 0.001);
        lwgeom_free(lwg);
        lwgeom_free(result);
 
        lwg = lwgeom_from_wkt("LINESTRING(2 48)", LW_PARSER_CHECK_NONE);
        length = lwgeom_length(lwg);
        CU_ASSERT_DOUBLE_EQUAL(length, 0.0, 0.001);
        result = geography_substring((LWLINE*)lwg, &s, 0.0, 1.0, 0.01);
        length = lwgeom_length(result);
        CU_ASSERT_DOUBLE_EQUAL(length, 0.0, 0.001);
        lwgeom_free(lwg);
        lwgeom_free(result);
}

References geography_substring(), LW_PARSER_CHECK_NONE, lwgeom_free(), lwgeom_from_wkt(), lwgeom_length(), result, s, spheroid_init(), and WGS84_RADIUS.

Referenced by geodetic_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: