ptarray_dp_findsplit()

static void ptarray_dp_findsplit ( POINTARRAY *  pts, int  p1, int  p2, int *  split, double *  dist  )

static

Definition at line 1552 of file ptarray.c.

References distance2d_sqr_pt_seg(), getPoint2d_cp(), LWDEBUG, LWDEBUGF, POINT2D::x, and POINT2D::y.

Referenced by ptarray_simplify().

{
        int k;
        const POINT2D *pk, *pa, *pb;
        double tmp, d;

        LWDEBUG(4, "function called");

        *split = p1;
        d = -1;

        if (p1 + 1 < p2)
        {

                pa = getPoint2d_cp(pts, p1);
                pb = getPoint2d_cp(pts, p2);

                LWDEBUGF(4, "P%d(%f,%f) to P%d(%f,%f)",
                         p1, pa->x, pa->y, p2, pb->x, pb->y);

                for (k=p1+1; k<p2; k++)
                {
                        pk = getPoint2d_cp(pts, k);

                        LWDEBUGF(4, "P%d(%f,%f)", k, pk->x, pk->y);

                        /* distance computation */
                        tmp = distance2d_sqr_pt_seg(pk, pa, pb);

                        if (tmp > d)
                        {
                                d = tmp;        /* record the maximum */
                                *split = k;

                                LWDEBUGF(4, "P%d is farthest (%g)", k, d);
                        }
                }
                *dist = d;

        } /* length---should be redone if can == 0 */
        else
        {
                LWDEBUG(3, "segment too short, no split/no dist");
                *dist = -1;
        }

}

Here is the call graph for this function:

Here is the caller graph for this function: