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 = 0; i < M * M; i ++ )
    {
        *matrixX ++ = ( float )( rand() % 2000 ) / 1000;
    }
}

void MulMatrix( float* matrixA, float* matrixB, float* matrixC )
{
    register int i, j, k;
    register float* p, *q, f;
    for ( j = 0; j < M; j ++ )
    {
        for ( i = 0; i < M; i ++ )
        {
            p = matrixA + j * M;
            q = matrixB + i;
            f = 0;
            for ( k = 0; 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( "
" );
    }*/

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


    return 0;
}

机器配置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 = 0; j < M; j ++ )
    {
        for ( i = 0; i < M; i ++ )
        {
            p = matrixA + j * M;
            q = matrixB + i;
            t = j * M + i;
            matrixC[t] = 0;
            for ( k = 0; k < M; k ++ )
            {
                matrixC[t] += *p * *q;
                p ++;
                q += M;
            }
        }
    }
}

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

难道也是CACHE。。