intel compiler的表现

好久没弄这个东西，今天突然想试下，代码没写完，以后补。

 #include <stdio.h>

 #include <stdlib.h>

 #include <time.h>

 #include <math.h>

 #include <Windows.h>

 #define M      1024

 float matA[M][M];

 float matB[M][M];

 float matC[M][M];

 void InitMatrix( float* matrixX )

 {

     register int i;

     for ( i = ; i < M * M; i ++ )

     {

         *matrixX ++ = ( float )( rand() %  ) / ;

     }

 }

 void MulMatrix( float* matrixA, float* matrixB, float* matrixC )

 {

     register int i, j, k;

     register float* p, *q, f;

     for ( j = ; j < M; j ++ )

     {

         for ( i = ; i < M; i ++ )

         {

             p = matrixA + j * M;

             q = matrixB + i;

             f = ;

             for ( k = ; k < M; k ++ )

             {

                 f += *p * *q;

                 p ++;

                 q += M;

             }

             matrixC[j * M + i] = f;

         }

     }

 }

 int main()

 {

     DWORD t;

     //register int i, j;

     srand( ( unsigned int )time( NULL ) );

     InitMatrix( ( float* )matA );

     InitMatrix( ( float* )matB );

     t = ::GetTickCount();

     MulMatrix( ( float* )matA, ( float* )matA, ( float* )matC );

     t = ::GetTickCount() - t;

     /*for ( j = 0; j < M; j ++ )

     {

         for ( i = 0; i < M; i ++ )

         {

             printf( "%.2f ", matC[j][i] );

         }

         printf( "\n" );

     }*/

     printf( "TIME:%d\n", t );

     return ;

 }

机器配置E3 1231V3 MEM:16G VS2010SP1 ICC 2015XE GTX660 将来把CUDA带来一起测下

1. CPU单线程仅一个O2

4750ms 大体如此

多线程原来测过，这次代码没加上。按物理核计算 4核应该6秒左右超线程估计会好些。应该能到5秒左右。

2. 单文件转为ICC编译额外添加优化项/Qipo /Qparallel

2600ms左右

多线程依然没测，以后补

3. CUDA也没测

4.MKL没测。。有点对不住这个CPU了。。呵呵，心血来潮，以后一定补上。

5. 比较搞笑的是，我在收工的时候突发奇想，要把MatrixC相关代码改为本地，试下有没有效果，这个还真有，平均少100ms

看来，高手们教的CACHE命中还是很有道理的。

上面代码是改后的，改之前为

void MulMatrix( float* matrixA, float* matrixB, float* matrixC )

{

    register int i, j, k, t;

    register float* p, *q;

    for ( j = ; j < M; j ++ )

    {

        for ( i = ; i < M; i ++ )

        {

            p = matrixA + j * M;

            q = matrixB + i;

            t = j * M + i;

            matrixC[t] = ;

            for ( k = ; k < M; k ++ )

            {

                matrixC[t] += *p * *q;

                p ++;

                q += M;

            }

        }

    }

}

6. 更搞笑的是，把q += M; 中的M改为100.。。。。竟然变为了原来的1/10

难道也是CACHE。。

巴特西

intel compiler的表现

最新文章

热门文章