lw_dist3d_ptarray_ptarray()

int lw_dist3d_ptarray_ptarray	(	POINTARRAY *	l1,
		POINTARRAY *	l2,
		DISTPTS3D *	dl
	)

Finds all combinationes of segments between two pointarrays.

Definition at line 867 of file measures3d.c.

References DIST_MAX, DIST_MIN, DISTPTS3D::distance, getPoint3dz_p(), lw_dist3d_pt_pt(), lw_dist3d_seg_seg(), LW_TRUE, LWDEBUGF, DISTPTS3D::mode, POINTARRAY::npoints, DISTPTS3D::tolerance, and DISTPTS3D::twisted.

Referenced by lw_dist3d_line_line(), lw_dist3d_line_poly(), lw_dist3d_poly_poly(), and lw_dist3d_ptarray_poly().

{
        int t,u;
        POINT3DZ        start, end;
        POINT3DZ        start2, end2;
        int twist = dl->twisted;
        LWDEBUGF(2, "lw_dist3d_ptarray_ptarray called (points: %d-%d)",l1->npoints, l2->npoints);



        if (dl->mode == DIST_MAX)/*If we are searching for maxdistance we go straight to point-point calculation since the maxdistance have to be between two vertexes*/
        {
                for (t=0; t<l1->npoints; t++) /*for each segment in L1 */
                {
                        getPoint3dz_p(l1, t, &start);
                        for (u=0; u<l2->npoints; u++) /*for each segment in L2 */
                        {
                                getPoint3dz_p(l2, u, &start2);
                                lw_dist3d_pt_pt(&start,&start2,dl);
                                LWDEBUGF(4, "maxdist_ptarray_ptarray; seg %i * seg %i, dist = %g\n",t,u,dl->distance);
                                LWDEBUGF(3, " seg%d-seg%d dist: %f, mindist: %f",
                                         t, u, dl->distance, dl->tolerance);
                        }
                }
        }
        else
        {
                getPoint3dz_p(l1, 0, &start);
                for (t=1; t<l1->npoints; t++) /*for each segment in L1 */
                {
                        getPoint3dz_p(l1, t, &end);
                        getPoint3dz_p(l2, 0, &start2);
                        for (u=1; u<l2->npoints; u++) /*for each segment in L2 */
                        {
                                getPoint3dz_p(l2, u, &end2);
                                dl->twisted=twist;
                                lw_dist3d_seg_seg(&start, &end, &start2, &end2,dl);
                                LWDEBUGF(4, "mindist_ptarray_ptarray; seg %i * seg %i, dist = %g\n",t,u,dl->distance);
                                LWDEBUGF(3, " seg%d-seg%d dist: %f, mindist: %f",
                                         t, u, dl->distance, dl->tolerance);
                                if (dl->distance<=dl->tolerance && dl->mode == DIST_MIN) return LW_TRUE; /*just a check if  the answer is already given*/
                                start2 = end2;
                        }
                        start = end;
                }
        }
        return LW_TRUE;
}

Here is the call graph for this function:

Here is the caller graph for this function: