         {
         int boffset = INDX( k, colTile[i], linearTiles ) *
                             tileSize * tileSize;
-#if 0
-        CUBLAS_CALL( cublasSetMatrixAsync( tileSize, tileSize, 8,
-          &p_b[boffset],
-          tileSize,
-          d_b[i],
-          tileSize,
-          stream[i] ));
-#endif
         CUDA_CALL( cudaMemcpyAsync( d_b[i], &p_b[boffset], tileBytes,
                      cudaMemcpyHostToDevice, stream[i] ) );
 
-//        printMat<<<1,1,0,stream[i]>>>(d_b[i], tileSize*tileSize);
         } /* end for */
 
 /* do the dgemm */

      {
        int coffset = INDX( rowTile[i], colTile[i], linearTiles ) *
                            tileSize * tileSize;
#if 0
        CUBLAS_CALL( cublasSetMatrixAsync( tileSize, tileSize, 8,
          &p_c[coffset],
          tileSize,
          d_c[i],
          tileSize,
          stream[i] ));
#endif
        CUDA_CALL( cudaMemcpyAsync( d_c[i], &p_c[coffset], tileBytes,
                     cudaMemcpyHostToDevice, stream[i] ) );


      {
        int coffset = INDX( rowTile[i], colTile[i], linearTiles ) *
                            tileSize * tileSize;
        CUDA_CALL( cudaMemcpyAsync( &p_c[coffset], d_c[i], tileBytes,
                        cudaMemcpyDeviceToHost, stream[i] ) );
      } /* end for */
