待编辑

原始代码

取自深入理解计算机系统p182页，这里使用了在线编译器Compiler Explorer(强烈推荐，可以实时看汇编代码，网址Compiler Explorer (godbolt.org)
)

#define N 16
typedef int fix_martrix[N][N];
int  fix_prod_ele(fix_martrix A,fix_martrix B,long i,long k)
{
    long j;
    int result =0;

    for (j=0;j<N;++j){
        result+=A[i][j]*B[j][k];
    }
    return result;
}

未优化时，即使用-Og的汇编代码

fix_prod_ele(int (*) [16], int (*) [16], long, long):
        movq    %rdi, %r9
        movq    %rsi, %r10
        movq    %rdx, %r8
        movq    %rcx, %rdi
        movl    $0, %esi
        movl    $0, %eax
        jmp     .L3
.L4:
        movq    %r8, %rcx
        salq    $6, %rcx
        addq    %r9, %rcx
        movq    %rax, %rdx
        salq    $6, %rdx
        addq    %r10, %rdx
        movl    (%rdx,%rdi,4), %edx     ;注意这里和下面使用了乘法指令伸缩，在一些处理器中会招致严重的性能处罚
        imull   (%rcx,%rax,4), %edx
        addl    %edx, %esi
        addq    $1, %rax
.L3:
        cmpq    $15, %rax
        jle     .L4
        movl    %esi, %eax
        ret

使用-O1优化后代码

fix_prod_ele(int (*) [16], int (*) [16], long, long):
        salq    $6, %rdx
        addq    %rdx, %rdi
        leaq    (%rsi,%rcx,4), %rax
        leaq    1024(%rax), %rsi
        movl    $0, %ecx
.L2:
        movl    (%rdi), %edx
        imull   (%rax), %edx
        addl    %edx, %ecx
        addq    $4, %rdi          ;优化后的寻址使用了加法和立即数来寻址，更快
        addq    $64, %rax
        cmpq    %rsi, %rax
        jne     .L2
        movl    %ecx, %eax
        ret