cudaDeviceProp  prop;

     int count;

     ( cudaGetDeviceCount( &count ) );

     for (int i=; i< count; i++) {

         ( cudaGetDeviceProperties( &prop, i ) );

         printf( "   --- General Information for device %d ---\n", i );

         printf( "Name:  %s\n", prop.name );

         printf( "Compute capability:  %d.%d\n", prop.major, prop.minor );

         printf( "Clock rate:  %d\n", prop.clockRate );

         printf( "Device copy overlap:  " );

         if (prop.deviceOverlap)

             printf( "Enabled\n" );

         else

             printf( "Disabled\n");

         printf( "Kernel execution timeout :  " );

         if (prop.kernelExecTimeoutEnabled)

             printf( "Enabled\n" );

         else

             printf( "Disabled\n" );

         printf( "   --- Memory Information for device %d ---\n", i );

         printf( "Total global mem:  %ld\n", prop.totalGlobalMem );

         printf( "Total constant Mem:  %ld\n", prop.totalConstMem );

         printf( "Max mem pitch:  %ld\n", prop.memPitch );

         printf( "Texture Alignment:  %ld\n", prop.textureAlignment );

         printf( "   --- MP Information for device %d ---\n", i );

         printf( "Multiprocessor count:  %d\n",

             prop.multiProcessorCount );

         printf( "Shared mem per mp:  %ld\n", prop.sharedMemPerBlock );

         printf( "Registers per mp:  %d\n", prop.regsPerBlock );

         printf( "Threads in warp:  %d\n", prop.warpSize );

         printf( "Max threads per block:  %d\n",

             prop.maxThreadsPerBlock );

         printf( "Max thread dimensions:  (%d, %d, %d)\n",

             prop.maxThreadsDim[], prop.maxThreadsDim[],

             prop.maxThreadsDim[] );

         printf( "Max grid dimensions:  (%d, %d, %d)\n",

             prop.maxGridSize[], prop.maxGridSize[],

             prop.maxGridSize[] );

         printf( "\n" );

     }

     int n1 = rand() % base + base;

     int m1 = rand() % base + base;

     int n2 = m1;

     int m2 = rand() % base + base;

     int *g1 = new int[n1 * m1];

     int *g2 = new int[n2 * m2];

     printf("matrix A[%3d %3d]\n", n1, m1);

     for(int i = ; i < n1 * m1;i++)

     {

         g1[i] = rand() % large;

         //printf("%5d ", g1[i]);

         //if((i + 1) % m1 == 0)printf("\n");

     }

     printf("matrix B[%3d %3d]\n", n2, m2);

     for(int i = ; i < n2 * m2;i++)

     {

         g2[i] = rand() % large;

         //printf("%5d ", g2[i]);

         //if((i + 1) % m2 == 0)printf("\n");

     }

     int *g;

     g = matrixMultiply(g1,n1,m1,g2,n2,m2);

     printf("matrix C[%3d %3d]\n", n1, m2);

     for(int i = ; i< n1*m2;i++)

     {

         //printf("%5d ", g[i]);

         //if((i + 1) % m2 == 0) printf("\n");

     }