test_sphere_project()

static void test_sphere_project ( void  )

static

Definition at line 92 of file cu_geodetic.c.

References geographic_point_init(), GEOGRAPHIC_POINT::lat, ll2cart(), GEOGRAPHIC_POINT::lon, normalize(), robust_cross_product(), s, sphere_direction(), sphere_distance(), sphere_project(), unit_normal(), POINT2D::x, POINT3D::x, POINT2D::y, POINT3D::y, and POINT3D::z.

Referenced by geodetic_suite_setup().

{
        GEOGRAPHIC_POINT s, e;
        double dir1, dist1, dir2, dist2;

        dir1 = M_PI_2;
        dist1 = 0.1;

        geographic_point_init(0, 0, &s);
        sphere_project(&s, dist1, dir1, &e);

        CU_ASSERT_DOUBLE_EQUAL(e.lon, 0.1, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(e.lat, 0.0, 1e-14);

        /* Direct and inverse solutions agree */
        dist2 = sphere_distance(&s, &e);
        dir2 = sphere_direction(&s, &e, dist1);

        CU_ASSERT_DOUBLE_EQUAL(dist1, dist2, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(dir1, dir2, 1e-14);

        dist1 = sphere_distance(&e, &s);
        dir1 = sphere_direction(&e, &s, dist1);
        sphere_project(&e, dist1, dir1, &s);

        CU_ASSERT_DOUBLE_EQUAL(s.lon, 0.0, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(s.lat, 0.0, 1e-14);

        geographic_point_init(0, 0.2, &e);
        geographic_point_init(0, 0.4, &s);
        dist1 = sphere_distance(&s, &e);
        dir1 = sphere_direction(&e, &s, dist1);
        /* GeodSolve -i -E -p 16 -e 1 0 --input-string "0.2 0 0.4 0" */
        CU_ASSERT_DOUBLE_EQUAL(dir1, 0.0, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(dist1, 0.0034906585039887, 1e-14);

        geographic_point_init(0, 1, &s); /* same start point for remainder of tests */
        geographic_point_init(0, 2, &e);
        dist2 = sphere_distance(&s, &e);
        dir2 = sphere_direction(&s, &e, dist2);
        /* GeodSolve -i -E -p 16 -e 1 0 --input-string "1 0 2 0" */
        CU_ASSERT_DOUBLE_EQUAL(dir2, 0.0, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(dist2, 0.0174532925199433, 1e-14);

        geographic_point_init(1, 1, &e);
        dist2 = sphere_distance(&s, &e);
        dir2 = sphere_direction(&s, &e, dist2);
        /* GeodSolve -i -E -p 16 -e 1 0 --input-string "1 0 1 1" */
        CU_ASSERT_DOUBLE_EQUAL(dir2, 89.991273575329292895136 * M_PI / 180.0, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(dist2, 0.0174506342314906, 1e-14);

        geographic_point_init(0, 0, &e);
        dist2 = sphere_distance(&s, &e);
        dir2 = sphere_direction(&s, &e, dist2);
        /* GeodSolve -i -E -p 16 -e 1 0 --input-string "1 0 0 0" */
        CU_ASSERT_DOUBLE_EQUAL(dir2, M_PI, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(dist2, 0.0174532925199433, 1e-14);

        geographic_point_init(-1, 1, &e);
        dist2 = sphere_distance(&s, &e);
        dir2 = sphere_direction(&s, &e, dist2);
        /* GeodSolve -i -E -p 16 -e 1 0 --input-string "1 0 1 -1" */
        CU_ASSERT_DOUBLE_EQUAL(dir2, -89.991273575329292895136 * M_PI / 180.0, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(dist2, 0.0174506342314906, 1e-14);

        geographic_point_init(1, 2, &e);
        dist2 = sphere_distance(&s, &e);
        dir2 = sphere_direction(&s, &e, dist2);
        /* GeodSolve -i -E -p 16 -e 1 0 --input-string "1 0 2 1" */
        CU_ASSERT_DOUBLE_EQUAL(dir2, 44.978182941465044354783 * M_PI / 180.0, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(dist2, 0.0246782972905467, 1e-14);

        geographic_point_init(-1, 0, &e);
        dist2 = sphere_distance(&s, &e);
        dir2 = sphere_direction(&s, &e, dist2);
        /* GeodSolve -i -E -p 16 -e 1 0 --input-string "1 0 0 -1" */
        CU_ASSERT_DOUBLE_EQUAL(dir2, -134.995636455344851488216 * M_PI / 180.0, 1e-14);
        CU_ASSERT_DOUBLE_EQUAL(dist2, 0.0246820563917664, 1e-14);
}

Here is the call graph for this function:

Here is the caller graph for this function: