 
     checkCUBLAS( cublasSgemm( cublasHandle, 
                               CUBLAS_OP_N, CUBLAS_OP_T,
-                              Xexamples, theta2Rows, theta2Cols,
+                              FIXME, FIXME, FIXME,
 			      &alpha, d_a2, Xexamples,
                               d_theta2, theta2Rows, &beta,
                               d_a3, Xexamples ) );                              


    checkCUBLAS( cublasSgemm( cublasHandle, 
                              CUBLAS_OP_T, CUBLAS_OP_N,
                              theta2Cols, Xexamples, theta2Rows,
			      &alpha, (float *)d_theta2, theta2Rows,
                              (float *)&d_delta3[1], 11, &beta,
                              (float *)d_delta2, theta1Rows+1 ) );   


__global__ void initOne( int size, floatType_t *array )
{

/* setup global threadID in X direction */

  int tid = FIXME;

  if( tid < size )
    array[tid] = (floatType_t) 1.0;
  return;
