         }
     }
 
-
     const float xn      = n.x + d.x;
     const float yn      = n.y + d.y;
     const float sn      = n.z + d.z;

        // Either we return 1, and n has not been modified.
        // Or we quit the loop by exceeding the limit, and reject the point anyway.

        float3 t = make_float3( fabsf(d.x), fabsf(d.y), fabsf(d.z) );

        if( t.x < 0.5f && t.y < 0.5f && t.z < 0.5f ) {
            // return false, quit the loop, success


#undef PRINT_EXTREMA_DEBUG_INFO

#undef DEBUG_MEGA_LOOP

namespace popart{

/*************************************************************
