test_edge_distance_to_point()

static void test_edge_distance_to_point ( void  )

static

Definition at line 730 of file cu_geodetic.c.

{
        GEOGRAPHIC_EDGE e;
        GEOGRAPHIC_POINT g;
        GEOGRAPHIC_POINT closest;
        double d;
 
        /* closest point at origin, one degree away */
        edge_set(-50.0, 0.0, 50.0, 0.0, &e);
        point_set(0.0, 1.0, &g);
        d = edge_distance_to_point(&e, &g, 0);
        CU_ASSERT_DOUBLE_EQUAL(d, M_PI / 180.0, 0.00001);
 
        /* closest point at origin, one degree away */
        edge_set(-50.0, 0.0, 50.0, 0.0, &e);
        point_set(0.0, 2.0, &g);
        d = edge_distance_to_point(&e, &g, &closest);
#if 0
        printf("LINESTRING(%.8g %.8g, %.8g %.8g)\n", e.start.lon,  e.start.lat, e.end.lon,  e.end.lat);
        printf("POINT(%.9g %.9g)\n", g.lon, g.lat);
        printf("\nDISTANCE == %.8g\n", d);
#endif
        CU_ASSERT_DOUBLE_EQUAL(d, M_PI / 90.0, 0.00001);
        CU_ASSERT_DOUBLE_EQUAL(closest.lat, 0.0, 0.00001);
        CU_ASSERT_DOUBLE_EQUAL(closest.lon, 0.0, 0.00001);
 
        /* Ticket #2351 */
        edge_set(149.386990599235, -26.3567415843982, 149.386990599247, -26.3567415843965, &e);
        point_set(149.386990599235, -26.3567415843982, &g);
        d = edge_distance_to_point(&e, &g, &closest);
        CU_ASSERT_DOUBLE_EQUAL(d, 0.0, 0.00001);
        // printf("CLOSE POINT(%g %g)\n", closest.lon,  closest.lat);
        // printf(" ORIG POINT(%g %g)\n", g.lon, g.lat);
        CU_ASSERT_DOUBLE_EQUAL(g.lat, closest.lat, 0.00001);
        CU_ASSERT_DOUBLE_EQUAL(g.lon, closest.lon, 0.00001);
}

References edge_distance_to_point(), edge_set(), GEOGRAPHIC_EDGE::end, GEOGRAPHIC_POINT::lat, GEOGRAPHIC_POINT::lon, point_set(), and GEOGRAPHIC_EDGE::start.

Referenced by geodetic_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: