ptarray_dp_findsplit_in_place()

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

static

Definition at line 1519 of file ptarray.c.

{
        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;
        }
        else
        {
                LWDEBUG(3, "segment too short, no split/no dist");
                *dist = -1;
        }
}

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

Referenced by ptarray_simplify_in_place().

Here is the call graph for this function:

Here is the caller graph for this function: