 
 __global__ void initOne( int size, floatType_t *array )
 {
-  int tid = blockDim.x * blockIdx.x + threadIdx.x;
+
+/* setup global threadID in X direction */
+
+  int tid = FIXME;
+
   if( tid < size )
     array[tid] = (floatType_t) 1.0;
   return;


   checkCUBLAS( cublasSgemm( cublasHandle, 
                              CUBLAS_OP_N, CUBLAS_OP_T,
                              theta1Rows, theta1Cols, Xexamples,
			      &recip, (float *)&d_delta2[1], theta1Rows+1,
                              (float *)d_X, Xfeatures, 
                              &beta, (float *)d_theta1Grad, theta1Rows ) );   

__global__ void k_sigmoid_f( floatType_t  *array,
                                     int    const size )
{

/* setup global threadID in X direction */

  int tid = FIXME

  if( tid < size )
  {

/* use the sigmoid_f function to complete this loop */

    array[tid] = FIXME

  } /* end if */
} /* end sigmoidGradient */

/* kernel for calculating the gradient of the sigmoid function */
