   int *indicesPtr = thrust::raw_pointer_cast(seedIndices.data());
   float4 *seedBufferPtr = thrust::raw_pointer_cast(seedBuffer.data());
 
-  seed << <dimGrid, dimBlock >> >(outputSurface, imageSize, labelCount,
-    seedBufferPtr, computePtr, idPtr, indicesPtr);
+  seed<<<dimGrid, dimBlock>>>(outputSurface, imageSize, labelCount,
+                              seedBufferPtr, computePtr, idPtr, indicesPtr);
   HANDLE_ERROR(cudaThreadSynchronize());
 }
 
 void Apollonius::runStepsKernels()
 {
-  apolloniusStep << <dimGrid, dimBlock >> >
-    (distancesTexture, thrust::raw_pointer_cast(computeVector.data()),
-    1, imageSize, imageSize);
+  apolloniusStep<<<dimGrid, dimBlock>>>(
+      distancesTexture, thrust::raw_pointer_cast(computeVector.data()), 1,
+      imageSize, imageSize);
 
   for (int k = (imageSize / 2); k > 0; k /= 2)
   {
-    apolloniusStep << <dimGrid, dimBlock >> >
-      (distancesTexture, thrust::raw_pointer_cast(computeVector.data()),
-      k, imageSize, imageSize);
+    apolloniusStep<<<dimGrid, dimBlock>>>(
+        distancesTexture, thrust::raw_pointer_cast(computeVector.data()), k,
+        imageSize, imageSize);
   }
   HANDLE_ERROR(cudaThreadSynchronize());
 }

 */
__device__ float4 rgb(int maxValue, int value)
{
    float ratio = 2.0f * value / maxValue;
    float b = max(0.0f, 1 - ratio);  // NOLINT
    float r = max(0.0f, ratio - 1);  // NOLINT
    float g = max(0.0f, 1.0f - b - r);  // NOLINT

    return make_float4(r, g, b, 1.0);
}

__global__ void gather(cudaSurfaceObject_t output, int imageSize,

  int *computePtr = thrust::raw_pointer_cast(computeVector.data());
  int *seedIdsPtr = thrust::raw_pointer_cast(seedIds.data());
  int *seedIndicesPtr = thrust::raw_pointer_cast(seedIndices.data());
  gather<<<dimGrid, dimBlock>>>(outputSurface, imageSize, labelCount,
                                computePtr, seedIdsPtr, seedIndicesPtr);
  HANDLE_ERROR(cudaThreadSynchronize());
}

