lw_dist3d_recursive()

int lw_dist3d_recursive	(	const LWGEOM *	lwg1,
		const LWGEOM *	lwg2,
		DISTPTS3D *	dl
	)

This is a recursive function delivering every possible combination of subgeometries.

Definition at line 399 of file measures3d.c.

References DIST_MIN, DISTPTS3D::distance, LWCOLLECTION::geoms, lw_dist3d_distribute_bruteforce(), lw_dist3d_recursive(), LW_FALSE, LW_TRUE, LWDEBUG, LWDEBUGF, lwgeom_as_lwcollection(), lwgeom_is_collection(), lwgeom_is_empty(), DISTPTS3D::mode, LWCOLLECTION::ngeoms, DISTPTS3D::tolerance, and LWGEOM::type.

Referenced by lw_dist3d_distanceline(), lw_dist3d_distancepoint(), lw_dist3d_recursive(), lwgeom_maxdistance3d_tolerance(), and lwgeom_mindistance3d_tolerance().

{
        int i, j;
        int n1=1;
        int n2=1;
        LWGEOM *g1 = NULL;
        LWGEOM *g2 = NULL;
        LWCOLLECTION *c1 = NULL;
        LWCOLLECTION *c2 = NULL;

        LWDEBUGF(2, "lw_dist3d_recursive is called with type1=%d, type2=%d", lwg1->type, lwg2->type);

        if (lwgeom_is_collection(lwg1))
        {
                LWDEBUG(3, "First geometry is collection");
                c1 = lwgeom_as_lwcollection(lwg1);
                n1 = c1->ngeoms;
        }
        if (lwgeom_is_collection(lwg2))
        {
                LWDEBUG(3, "Second geometry is collection");
                c2 = lwgeom_as_lwcollection(lwg2);
                n2 = c2->ngeoms;
        }

        for ( i = 0; i < n1; i++ )
        {

                if (lwgeom_is_collection(lwg1))
                {
                        g1 = c1->geoms[i];
                }
                else
                {
                        g1 = (LWGEOM*)lwg1;
                }

                if (lwgeom_is_empty(g1)) return LW_TRUE;

                if (lwgeom_is_collection(g1))
                {
                        LWDEBUG(3, "Found collection inside first geometry collection, recursing");
                        if (!lw_dist3d_recursive(g1, lwg2, dl)) return LW_FALSE;
                        continue;
                }
                for ( j = 0; j < n2; j++ )
                {
                        if (lwgeom_is_collection(lwg2))
                        {
                                g2 = c2->geoms[j];
                        }
                        else
                        {
                                g2 = (LWGEOM*)lwg2;
                        }
                        if (lwgeom_is_collection(g2))
                        {
                                LWDEBUG(3, "Found collection inside second geometry collection, recursing");
                                if (!lw_dist3d_recursive(g1, g2, dl)) return LW_FALSE;
                                continue;
                        }


                        /*If one of geometries is empty, return. True here only means continue searching. False would have stoped the process*/
                        if (lwgeom_is_empty(g1)||lwgeom_is_empty(g2)) return LW_TRUE;


                        if (!lw_dist3d_distribute_bruteforce(g1, g2, dl)) return LW_FALSE;
                        if (dl->distance<=dl->tolerance && dl->mode == DIST_MIN) return LW_TRUE; /*just a check if  the answer is already given*/
                }
        }
        return LW_TRUE;
}

Here is the call graph for this function:

Here is the caller graph for this function: