test_lw_dist2d_pt_ptarrayarc()

static void test_lw_dist2d_pt_ptarrayarc ( void  )

static

Definition at line 912 of file cu_measures.c.

References DIST_MIN, DISTPTS::distance, lw_dist2d_distpts_init(), lw_dist2d_pt_ptarrayarc(), LW_SUCCESS, lwgeom_as_lwline(), lwgeom_from_text(), lwline_free(), LWLINE::points, POINT2D::x, and POINT2D::y.

Referenced by measures_suite_setup().

{
        /* lw_dist2d_pt_ptarrayarc(const POINT2D *p, const POINTARRAY *pa, DISTPTS *dl) */
        DISTPTS dl;
        int rv;
        LWLINE *lwline;
        POINT2D P;

        /* Unit semi-circle above X axis */
        lwline = lwgeom_as_lwline(lwgeom_from_text("LINESTRING(-1 0, 0 1, 1 0)"));

        /* Point at origin */
        P.x = P.y = 0;
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1, 0.000001);

        /* Point above arc on Y axis */
        P.y = 2;
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1, 0.000001);

        /* Point 45 degrees off arc, 2 radii from center */
        P.y = P.x = 2 * cos(M_PI_4);
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1, 0.000001);

        /* Four unit semi-circles surrounding the 2x2 box around origin */
        lwline_free(lwline);
        lwline = lwgeom_as_lwline(lwgeom_from_text("LINESTRING(-1 -1, -2 0, -1 1, 0 2, 1 1, 2 0, 1 -1, 0 -2, -1 -1)"));

        /* Point at origin */
        P.x = P.y = 0;
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, sqrt(2.0), 0.000001);

        /* Point on box edge */
        P.x = -1; P.y = 0;
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, 1, 0.000001);

        /* Point within a semicircle lobe */
        P.x = -1.5; P.y = 0;
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.5, 0.000001);

        /* Point outside a semicircle lobe */
        P.x = -2.5; P.y = 0;
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.5, 0.000001);

        /* Point outside a semicircle lobe */
        P.y = -2.5; P.x = 0;
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, 0.5, 0.000001);

        /* Point outside a semicircle lobe */
        P.y = 2; P.x = 1;
        lw_dist2d_distpts_init(&dl, DIST_MIN);
        rv = lw_dist2d_pt_ptarrayarc(&P, lwline->points, &dl);
        CU_ASSERT_EQUAL( rv, LW_SUCCESS );
        CU_ASSERT_DOUBLE_EQUAL(dl.distance, sqrt(2.0)-1.0, 0.000001);

        /* Clean up */
        lwline_free(lwline);
}

Here is the call graph for this function:

Here is the caller graph for this function: