   float* d_Y;
   float* d_PC;
   float* d_Signals;
+
+  cudaMalloc(&d_PC, inCols * inCols * sizeof(float));
+  utilVectorSetByScalar(d_PC, 0, inCols * inCols);
 
   // Subtract mean for each dimension.
   pcaUpdateMean(d_inMatrix, inRows, inCols);


void pcaSVD(float* d_Y, unsigned int inRows, unsigned int inCols, float* d_PC)
{

}

void pcaCalculateSignals(float* d_PC, float* d_inMatrix, unsigned int inRows,

  }
  printf("&&&\n");

  // Calculate signals.
  pcaCalculateSignals(d_PC, d_inMatrix, inRows, inCols, d_Signals);

