                                            vd->vox[1] * vd->vox[1] +
                                            vd->vox[2] * vd->vox[2]));
   int numSlices = max(1, static_cast<int>(_renderState._quality * diagonalVoxels));
-  //std::cerr << diagonalVoxels << " " << numSlices << std::endl;
 
   // Inverse modelview-projection matrix.
   vvMatrix mvp, pr;

    justClippedPlane = false;
    justClippedSphere = false;

    // pre-multiply alpha
    src.x *= src.w;
    src.y *= src.w;

         true, // Front to back.
         0, // Mip mode.
         true, // Local illumination.
         true, // Opacity correction.
         false, // Jittering.
         false, // Clip sphere.
         false // Clip plane.
        ><<<gridSize, blockSize>>>(d_output, width, height,
                                   diagonalVoxels / (float)numSlices,
                                   volSize * 0.5f,
                                   L, H,
                                   center, radius * radius,
