 void callStages(thrust::device_vector<float> &inout,
                 thrust::device_vector<float> &ybar,
                 thrust::device_vector<float> &vhat,
-                thrust::device_vector<float> &ysum, int compute_width,
-                int compute_height, int compute_m_size, int compute_n_size)
-{
-  const int nWm = (compute_width + MTS - 1) / MTS,
-            nHm = (compute_height + MTS - 1) / MTS;
+                thrust::device_vector<float> &ysum, int computeWidth,
+                int computeHeight, int compute_m_size, int compute_n_size)
+{
+  const int nWm = (computeWidth + MTS - 1) / MTS,
+            nHm = (computeHeight + MTS - 1) / MTS;
   const dim3 cg_img(compute_m_size, compute_n_size);
   const dim3 cg_ybar(nWm, 1);
   const dim3 cg_vhat(1, nHm);

             thrust::device_vector<float> &vhat,
             thrust::device_vector<float> &ysum)
{
  int compute_width = image_size;
  int compute_height = image_size;
  if (compute_width % 32 > 0)
    compute_width += (32 - (compute_width % 32));
  if (compute_height % 32 > 0)
    compute_height += (32 - (compute_height % 32));
  int compute_m_size = (compute_width + WS - 1) / WS;
  int compute_n_size = (compute_height + WS - 1) / WS;

  // set data structure sizes
  resizeIfNecessary(inout, compute_width * compute_height);
  resizeIfNecessary(ybar, compute_n_size * compute_width);
  resizeIfNecessary(vhat, compute_m_size * compute_height);
  resizeIfNecessary(ysum, compute_m_size * compute_n_size);



  cudaUnbindTexture(&textureReadDepth);
  inputImage->unmap();

  callStages(inout, ybar, vhat, ysum, computeWidth, computeHeight,
             compute_m_size, compute_n_size);
}

