test_gbox_from_spherical_coordinates()

static void test_gbox_from_spherical_coordinates ( void  )

static

Definition at line 224 of file cu_geodetic.c.

{
#if RANDOM_TEST
        const double gtolerance = 0.000001;
        const int loops = RANDOM_TEST;
        int i;
        double ll[64];
        GBOX gbox;
        GBOX gbox_slow;
        int rndlat;
        int rndlon;
 
        POINTARRAY *pa;
        LWGEOM *lwline;
 
        ll[0] = -3.083333333333333333333333333333333;
        ll[1] = 9.83333333333333333333333333333333;
        ll[2] = 15.5;
        ll[3] = -5.25;
 
        pa = ptarray_construct_reference_data(0, 0, 2, (uint8_t*)ll);
 
        lwline = lwline_as_lwgeom(lwline_construct(SRID_UNKNOWN, 0, pa));
        FLAGS_SET_GEODETIC(lwline->flags, 1);
 
        srandomdev();
 
        for ( i = 0; i < loops; i++ )
        {
                rndlat = (int)(90.0 - 180.0 * (double)random() / pow(2.0, 31.0));
                rndlon = (int)(180.0 - 360.0 * (double)random() / pow(2.0, 31.0));
                ll[0] = (double)rndlon;
                ll[1] = (double)rndlat;
 
                rndlat = (int)(90.0 - 180.0 * (double)random() / pow(2.0, 31.0));
                rndlon = (int)(180.0 - 360.0 * (double)random() / pow(2.0, 31.0));
                ll[2] = (double)rndlon;
                ll[3] = (double)rndlat;
 
                gbox_geocentric_slow = LW_FALSE;
                lwgeom_calculate_gbox_geodetic(lwline, &gbox);
                gbox_geocentric_slow = LW_TRUE;
                lwgeom_calculate_gbox_geodetic(lwline, &gbox_slow);
                gbox_geocentric_slow = LW_FALSE;
 
                if (
                        ( fabs( gbox.xmin - gbox_slow.xmin ) > gtolerance ) ||
                        ( fabs( gbox.xmax - gbox_slow.xmax ) > gtolerance ) ||
                        ( fabs( gbox.ymin - gbox_slow.ymin ) > gtolerance ) ||
                        ( fabs( gbox.ymax - gbox_slow.ymax ) > gtolerance ) ||
                        ( fabs( gbox.zmin - gbox_slow.zmin ) > gtolerance ) ||
                        ( fabs( gbox.zmax - gbox_slow.zmax ) > gtolerance ) )
                {
                        printf("\n-------\n");
                        printf("If you are seeing this, cut and paste it, it is a randomly generated test case!\n");
                        printf("LOOP: %d\n", i);
                        printf("SEGMENT (Lon Lat): (%.9g %.9g) (%.9g %.9g)\n", ll[0], ll[1], ll[2], ll[3]);
                        printf("CALC: %s\n", gbox_to_string(&gbox));
                        printf("SLOW: %s\n", gbox_to_string(&gbox_slow));
                        printf("-------\n\n");
                        CU_FAIL_FATAL(Slow (GOOD) and fast (CALC) box calculations returned different values!!);
                }
 
        }
 
        lwgeom_free(lwline);
#endif /* RANDOM_TEST */
}

References LWGEOM::flags, FLAGS_SET_GEODETIC, gbox_geocentric_slow, gbox_to_string(), LW_FALSE, LW_TRUE, lwgeom_calculate_gbox_geodetic(), lwgeom_free(), lwline_as_lwgeom(), lwline_construct(), ptarray_construct_reference_data(), RANDOM_TEST, SRID_UNKNOWN, GBOX::xmax, GBOX::xmin, GBOX::ymax, GBOX::ymin, GBOX::zmax, and GBOX::zmin.

Referenced by geodetic_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: