test_spheroid_distance()

static void test_spheroid_distance ( void  )

static

Definition at line 1275 of file cu_geodetic.c.

References point_set(), s, spheroid_distance(), and spheroid_init().

Referenced by geodetic_suite_setup().

{
        GEOGRAPHIC_POINT g1, g2;
        double d;
#if ! PROJ_GEODESIC
        double epsilon; /* irregular */
#else
        const double epsilon = 1e-8; /* at least 10 nm precision */
#endif
        SPHEROID s;

        /* Init to WGS84 */
        spheroid_init(&s, 6378137.0, 6356752.314245179498);

        /* One vertical degree
        $ GeodSolve -E -i -p 16 --input-string "0 0 1 0" */
        point_set(0.0, 0.0, &g1);
        point_set(0.0, 1.0, &g2);
        d = spheroid_distance(&g1, &g2, &s);
#if ! PROJ_GEODESIC
        epsilon = 1e-6;
#endif
        CU_ASSERT_DOUBLE_EQUAL(d, 110574.3885577987957342, epsilon);

        /* Ten horizontal degrees
        $ GeodSolve -E -i -p 16 --input-string "0 -10 0 0" */
        point_set(-10.0, 0.0, &g1);
        point_set(0.0, 0.0, &g2);
        d = spheroid_distance(&g1, &g2, &s);
#if ! PROJ_GEODESIC
        epsilon = 1e-3;
#endif
        CU_ASSERT_DOUBLE_EQUAL(d, 1113194.9079327357264771, epsilon);

        /* One horizonal degree
        $ GeodSolve -E -i -p 16 --input-string "0 -1 0 0" */
        point_set(-1.0, 0.0, &g1);
        point_set(0.0, 0.0, &g2);
        d = spheroid_distance(&g1, &g2, &s);
#if ! PROJ_GEODESIC
        epsilon = 1e-4;
#endif
        CU_ASSERT_DOUBLE_EQUAL(d, 111319.4907932735726477, epsilon);

        /* Around world w/ slight bend
        $ GeodSolve -E -i -p 16 --input-string "0 -180 1 0" */
        point_set(-180.0, 0.0, &g1);
        point_set(0.0, 1.0, &g2);
        d = spheroid_distance(&g1, &g2, &s);
#if ! PROJ_GEODESIC
        epsilon = 1e-5;
#endif
        CU_ASSERT_DOUBLE_EQUAL(d, 19893357.0700676468277450, epsilon);

        /* Up to pole
        $ GeodSolve -E -i -p 16 --input-string "0 -180 90 0" */
        point_set(-180.0, 0.0, &g1);
        point_set(0.0, 90.0, &g2);
        d = spheroid_distance(&g1, &g2, &s);
#if ! PROJ_GEODESIC
        epsilon = 1e-6;
#endif
        CU_ASSERT_DOUBLE_EQUAL(d, 10001965.7293127228117396, epsilon);

}

Here is the call graph for this function:

Here is the caller graph for this function: