 __global__ void k_sigmoid_f( floatType_t  *array,
                                      int    const size )
 {
-  int tid = blockDim.x * blockIdx.x + threadIdx.x;
+
+/* setup global threadID in X direction */
+
+  int tid = FIXME
+
   if( tid < size )
-    array[tid] = sigmoid_f( array[tid] );
+  {
+
+/* 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 */

                               floatType_t const *a3,
                               int         const  Xexamples )
{
  int tid = blockDim.x * blockIdx.x + threadIdx.x;
  if( tid < Xexamples )
  {
    delta3[INDX((int)Y[tid],tid,11)] = (floatType_t) 1.0;


    checkCUBLAS( cublasSgemm( cublasHandle, 
                              CUBLAS_OP_T, CUBLAS_OP_T,
                              FIXME, FIXME, FIXME,
			      &alpha, d_X, Xfeatures,
                              d_theta1, theta1Rows, &beta,
                              &d_z2[INDX(0,1,Xexamples)], Xexamples ) );                              
