61 #define cbrt(x) \  argument
62     ((x) > 0.0 ? pow((double)(x), 1.0 / 3.0) : ((x) < 0.0 ? -pow((double)-(x), 1.0 / 3.0) : 0.0))
64 #define qbrt(x) \  argument
65     ((x) > 0.0 ? pow((double)(x), 1.0 / 4.0) : ((x) < 0.0 ? -pow((double)-(x), 1.0 / 4.0) : 0.0))
93     g->voxsize.x /= (flt)g->xsize;  in newgrid()
146     gmin->x = FHUGE;  in globalbound()
149     gmax->x = -FHUGE;  in globalbound()
155         min.x = -FHUGE;  in globalbound()
158         max.x = FHUGE;  in globalbound()
163             gmin->x = MYMIN(gmin->x, min.x);  in globalbound()
167             gmax->x = MYMAX(gmax->x, max.x);  in globalbound()
181     cur = g->cells[index->z * g->xsize * g->ysize + index->y * g->xsize + index->x];  in cellbound()
186     cellmin.x = voxel2x(g, index->x);  in cellbound()
190     cellmax.x = cellmin.x + g->voxsize.x;  in cellbound()
194     cmin->x = FHUGE;  in cellbound()
197     cmax->x = -FHUGE;  in cellbound()
202         min.x = -FHUGE;  in cellbound()
205         max.x = FHUGE;  in cellbound()
210             if ((min.x >= cellmin.x) && (max.x <= cellmax.x) && (min.y >= cellmin.y) &&  in cellbound()
212                 cmin->x = MYMIN(cmin->x, min.x);  in cellbound()
216                 cmax->x = MYMAX(cmax->x, max.x);  in cellbound()
229     if ((cmax->x - cmin->x) < EPSILON) {  in cellbound()
230         cmax->x += EPSILON;  in cellbound()
231         cmin->x -= EPSILON;  in cellbound()
301                 for (index.x = 0; index.x < g->xsize; index.x++) {  in engrid_scene()
340     list = &gold->cells[index->z * gold->xsize * gold->ysize + index->y * gold->xsize + index->x];  in engrid_cell()
348     len = 1.0 / (MYMAX(MYMAX(gsize.x, gsize.y), gsize.z));  in engrid_cell()
349     gsize.x *= len;  in engrid_cell()
356         xs = (int)((flt)numcbrt * gsize.x);  in engrid_cell()
412     int x, y, z, zindex, yindex, voxindex;  in engrid_object()  local
433             for (x = low.x; x <= high.x; x++) {  in engrid_object()
434                 voxindex = x + yindex + zindex;  in engrid_object()
447     index->x = (int)((pos->x - g->min.x) / g->voxsize.x);  in pos2grid()
451     if (index->x == g->xsize)  in pos2grid()
452         index->x--;  in pos2grid()
458     if (index->x < 0 || index->x > g->xsize || index->y < 0 || index->y > g->ysize ||  in pos2grid()
462     if (pos->x < g->min.x || pos->x > g->max.x || pos->y < g->min.y || pos->y > g->max.y ||  in pos2grid()
491     if (fabs(ry->d.x) < EPSILON) {  in grid_intersect()
492         tmax.x = FHUGE;  in grid_intersect()
493         tdelta.x = 0.0;  in grid_intersect()
494         step.x = 0;  in grid_intersect()
495         out.x = 0; /* never goes out of bounds on this axis */  in grid_intersect()
497     else if (ry->d.x < 0.0) {  in grid_intersect()
498         tmax.x = offset + ((voxel2x(g, curvox.x) - curpos.x) / ry->d.x);  in grid_intersect()
499         tdelta.x = g->voxsize.x / -ry->d.x;  in grid_intersect()
500         step.x = out.x = -1;  in grid_intersect()
503         tmax.x = offset + ((voxel2x(g, curvox.x + 1) - curpos.x) / ry->d.x);  in grid_intersect()
504         tdelta.x = g->voxsize.x / ry->d.x;  in grid_intersect()
505         step.x = 1;  in grid_intersect()
506         out.x = g->xsize;  in grid_intersect()
548     VScale(&pdeltaX, tdelta.x);  in grid_intersect()
554     nXp = Raypnt(ry, tmax.x);  in grid_intersect()
558     voxindex = curvox.z * g->xsize * g->ysize + curvox.y * g->xsize + curvox.x;  in grid_intersect()
560         if (tmax.x < tmax.y && tmax.x < tmax.z) {  in grid_intersect()
569             curvox.x += step.x;  in grid_intersect()
570             if (ry->maxdist < tmax.x || curvox.x == out.x)  in grid_intersect()
572             voxindex += step.x;  in grid_intersect()
573             tmax.x += tdelta.x;  in grid_intersect()
575             nXp.x += pdeltaX.x;  in grid_intersect()
594             nZp.x += pdeltaZ.x;  in grid_intersect()
613             nYp.x += pdeltaY.x;  in grid_intersect()
640     if (ry->d.x == 0.0) {  in grid_bounds_intersect()
641         if ((ry->o.x < g->min.x) || (ry->o.x > g->max.x))  in grid_bounds_intersect()
645         tx1 = (g->min.x - ry->o.x) / ry->d.x;  in grid_bounds_intersect()
646         tx2 = (g->max.x - ry->o.x) / ry->d.x;  in grid_bounds_intersect()