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


__global__ void initOne( int size, floatType_t *array )
{
  int tid = blockDim.x * blockIdx.x + threadIdx.x;
  if( tid < size )
    array[tid] = (floatType_t) 1.0;
  return;

                                int         const Xexamples,
                                int         const size )
{

/* setup global threadID in X and Y directions */

  int tidx = blockDim.x * blockIdx.x + threadIdx.x;
  int tidy = blockDim.y * blockIdx.y + threadIdx.y;

  if( tidy < Xexamples && tidx < size )
  {

/* calculate the offset properly */

    delta2[INDX(tidx,tidy,size)] *= z2[INDX(tidy,tidx,Xexamples)];

  } /* end if */
} /* end k_updateDelta2 */

/* kernel for calculating the sigmoid of an array */
