   {
     sum = 0;
     for(j = 0; j <= i; ++j) {
-	  sum += *sigma++ * (data[j] - mean[j]);
+      sum += *sigma++ * (data[j] - mean[j]);
     }
     discrim += sum * sum;
   }

}

__device__ void copy_data(const PMatrix* data, float* sh_data,
						  int thidx, int thidy, int obs_num)
{
  if (obs_num >= data->rows)
	return;

  for (int chunk = 0; chunk < data->cols; chunk += blockDim.y)
  {

}

void get_tuned_layout(TuningInfo* info, PMatrix* data, PMatrix* params,
                      int max_block_params) {
  // query the device for smem / max # of threads
  int max_smem = smem_size();
  int max_threads = max_block_threads();
