         {
 /* set the stream for the cublas call */
 
-          CUBLAS_CALL( cublasSetStream( handle, stream[i] ) );
+          FIXME
 
 /* call the cublas dgemm function */
 
-          CUBLAS_CALL( cublasDgemm( handle, CUBLAS_OP_N, CUBLAS_OP_N,
-                 tileSize, tileSize, tileSize,
-                 &alpha,
-                 d_a[i], tileSize,
-                 d_b[i], tileSize,
-                 &beta,
-                 d_c[i], tileSize ) );
+          FIXME
 
         } /* end for */
 


        int aOffset = INDX( rowTile[i], k, linearTiles ) *
                            tileSize * tileSize;
        CUDA_CALL( cudaMemcpyAsync( d_a[i], &p_a[aOffset], tileBytes,
                     cudaMemcpyHostToDevice, stream[i] ) );
        } /* end for */

/* stream B into device */

    for( int i = 0; i < nstreams; i++ )
    {
/* destroy the streams here */
      FIXME
    } /* end for */

    CUDA_CALL( cudaFreeHost( p_a ) );
