int dgemm3(double *a, double *b, double *c, int n) {
  int i, j, k;
  for (i = 0; i < n; i += 3)
    for (j = 0; j < n; j += 3) {
      register double c11 = c[i * n + j];
      register double c12 = c[i * n + j + 1];
      register double c13 = c[i * n + j + 2];
      register double c21 = c[(i + 1) * n + j];
      register double c22 = c[(i + 1) * n + j + 1];
      register double c23 = c[(i + 1) * n + j + 2];
      register double c31 = c[(i + 2) * n + j];
      register double c32 = c[(i + 2) * n + j + 1];
      register double c33 = c[(i + 2) * n + j + 2];
      for (k = 0; k < n; k++) {
        register double a1 = a[i * n + k];
        register double a2 = a[(i + 1) * n + k];
        register double a3 = a[(i + 2) * n + k];
        register double b1 = b[k * n + j];
        register double b2 = b[k * n + j + 1];
        register double b3 = b[k * n + j + 2];
        c11 += a1 * b1;
        c12 += a1 * b2;
        c13 += a1 * b3;
        c21 += a2 * b1;
        c22 += a2 * b2;
        c23 += a2 * b3;
        c31 += a3 * b1;
        c32 += a3 * b2;
        c33 += a3 * b3;
      }
      c[i * n + j] = c11;
      c[i * n + j + 1] = c12;
      c[i * n + j + 2] = c13;
      c[(i + 1) * n + j] = c21;
      c[(i + 1) * n + j + 1] = c22;
      c[(i + 1) * n + j + 2] = c23;
      c[(i + 2) * n + j] = c31;
      c[(i + 2) * n + j + 1] = c32;
      c[(i + 2) * n + j + 2] = c33;
    }
  return 0;
}


x86 assembly (clang, -O0):
<dgemm3>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1587 <dgemm3+0x487>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1574 <dgemm3+0x474>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    144f <dgemm3+0x34f>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x90(%rbp),%xmm2
ff
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0xa0(%rbp),%xmm2
ff
movsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    125b <dgemm3+0x15b>
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x3,%eax
mov    %eax,-0x24(%rbp)
jmp    1131 <dgemm3+0x31>
jmp    1579 <dgemm3+0x479>
mov    -0x20(%rbp),%eax
add    $0x3,%eax
mov    %eax,-0x20(%rbp)
jmp    111e <dgemm3+0x1e>
xor    %eax,%eax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dgemm3>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rsi,-0x20(%rsp)
test   %ecx,%ecx
jle    12bf <dgemm3+0x1bf>
mov    %ecx,%r10d
addq   $0x10,-0x20(%rsp)
lea    0x0(,%r10,8),%r8
mov    %r10,%rax
shl    $0x4,%rax
add    %rdi,%rax
lea    (%r8,%r8,2),%rcx
mov    %rcx,-0x18(%rsp)
lea    (%rdi,%r10,8),%rcx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbp,%rsi
imul   %r10,%rsi
mov    %rsi,-0x8(%rsp)
lea    0x1(%rbp),%r15
imul   %r10,%r15
mov    %rbp,-0x10(%rsp)
lea    0x2(%rbp),%r12
imul   %r10,%r12
mov    -0x20(%rsp),%r11
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
mov    -0x8(%rsp),%rsi
add    %rbx,%rsi
lea    (%rbx,%r15,1),%r9
lea    (%rbx,%r12,1),%rbp
movsd  0x10(%rdx,%rbp,8),%xmm11
movupd (%rdx,%rbp,8),%xmm0
movsd  0x10(%rdx,%r9,8),%xmm12
movupd (%rdx,%r9,8),%xmm3
movsd  0x10(%rdx,%rsi,8),%xmm13
movupd (%rdx,%rsi,8),%xmm7
mov    %r11,%r13
xor    %r14d,%r14d
nopl   0x0(%rax,%rax,1)
movapd %xmm7,%xmm8
movapd %xmm3,%xmm9
movapd %xmm0,%xmm10
movsd  (%rdi,%r14,8),%xmm0
movsd  (%rcx,%r14,8),%xmm5
movsd  (%rax,%r14,8),%xmm1
movupd -0x10(%r13),%xmm2
movsd  0x0(%r13),%xmm14
movapd %xmm0,%xmm6
unpcklpd %xmm0,%xmm6
mulpd  %xmm2,%xmm6
movapd %xmm6,%xmm7
addpd  %xmm8,%xmm7
mulsd  %xmm14,%xmm0
addsd  %xmm0,%xmm13
movapd %xmm5,%xmm4
unpcklpd %xmm5,%xmm4
mulpd  %xmm2,%xmm4
movapd %xmm4,%xmm3
addpd  %xmm9,%xmm3
mulsd  %xmm14,%xmm5
addsd  %xmm5,%xmm12
movapd %xmm1,%xmm5
unpcklpd %xmm1,%xmm5
mulpd  %xmm2,%xmm5
movapd %xmm5,%xmm0
addpd  %xmm10,%xmm0
mulsd  %xmm14,%xmm1
addsd  %xmm1,%xmm11
add    $0x1,%r14
add    %r8,%r13
cmp    %r14,%r10
jne    11c0 <dgemm3+0xc0>
addpd  %xmm10,%xmm5
addpd  %xmm9,%xmm4
addpd  %xmm8,%xmm6
movupd %xmm6,(%rdx,%rsi,8)
movsd  %xmm13,0x10(%rdx,%rsi,8)
movupd %xmm4,(%rdx,%r9,8)
movsd  %xmm12,0x10(%rdx,%r9,8)
movupd %xmm5,(%rdx,%rbp,8)
movsd  %xmm11,0x10(%rdx,%rbp,8)
add    $0x3,%rbx
add    $0x18,%r11
cmp    %r10,%rbx
jb     1180 <dgemm3+0x80>
mov    -0x10(%rsp),%rbp
add    $0x3,%rbp
mov    -0x18(%rsp),%rsi
add    %rsi,%rax
add    %rsi,%rcx
add    %rsi,%rdi
cmp    %r10,%rbp
jb     1150 <dgemm3+0x50>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @dgemm3(double* %0, double* %1, double* %2, i32 %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
store double* %0, double** %5
store double* %1, double** %6
store double* %2, double** %7
store i32 %3, i32* %8
store i32 0, i32* %9
br label %27
27: ; preds = %330, %4
%28 = load i32, i32* %9
%29 = load i32, i32* %8
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %333
31: ; preds = %27
store i32 0, i32* %10
br label %32
32: ; preds = %326, %31
%33 = load i32, i32* %10
%34 = load i32, i32* %8
%35 = icmp slt i32 %33, %34
br i1 %35, label %36, label %329
36: ; preds = %32
%37 = load double*, double** %7
%38 = load i32, i32* %9
%39 = load i32, i32* %8
%40 = mul nsw i32 %38, %39
%41 = load i32, i32* %10
%42 = add nsw i32 %40, %41
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %37, i64 %43
%45 = load double, double* %44
store double %45, double* %12
%46 = load double*, double** %7
%47 = load i32, i32* %9
%48 = load i32, i32* %8
%49 = mul nsw i32 %47, %48
%50 = load i32, i32* %10
%51 = add nsw i32 %49, %50
%52 = add nsw i32 %51, 1
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %46, i64 %53
%55 = load double, double* %54
store double %55, double* %13
%56 = load double*, double** %7
%57 = load i32, i32* %9
%58 = load i32, i32* %8
%59 = mul nsw i32 %57, %58
%60 = load i32, i32* %10
%61 = add nsw i32 %59, %60
%62 = add nsw i32 %61, 2
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %56, i64 %63
%65 = load double, double* %64
store double %65, double* %14
%66 = load double*, double** %7
%67 = load i32, i32* %9
%68 = add nsw i32 %67, 1
%69 = load i32, i32* %8
%70 = mul nsw i32 %68, %69
%71 = load i32, i32* %10
%72 = add nsw i32 %70, %71
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %66, i64 %73
%75 = load double, double* %74
store double %75, double* %15
%76 = load double*, double** %7
%77 = load i32, i32* %9
%78 = add nsw i32 %77, 1
%79 = load i32, i32* %8
%80 = mul nsw i32 %78, %79
%81 = load i32, i32* %10
%82 = add nsw i32 %80, %81
%83 = add nsw i32 %82, 1
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %76, i64 %84
%86 = load double, double* %85
store double %86, double* %16
%87 = load double*, double** %7
%88 = load i32, i32* %9
%89 = add nsw i32 %88, 1
%90 = load i32, i32* %8
%91 = mul nsw i32 %89, %90
%92 = load i32, i32* %10
%93 = add nsw i32 %91, %92
%94 = add nsw i32 %93, 2
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double, double* %87, i64 %95
%97 = load double, double* %96
store double %97, double* %17
%98 = load double*, double** %7
%99 = load i32, i32* %9
%100 = add nsw i32 %99, 2
%101 = load i32, i32* %8
%102 = mul nsw i32 %100, %101
%103 = load i32, i32* %10
%104 = add nsw i32 %102, %103
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds double, double* %98, i64 %105
%107 = load double, double* %106
store double %107, double* %18
%108 = load double*, double** %7
%109 = load i32, i32* %9
%110 = add nsw i32 %109, 2
%111 = load i32, i32* %8
%112 = mul nsw i32 %110, %111
%113 = load i32, i32* %10
%114 = add nsw i32 %112, %113
%115 = add nsw i32 %114, 1
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %108, i64 %116
%118 = load double, double* %117
store double %118, double* %19
%119 = load double*, double** %7
%120 = load i32, i32* %9
%121 = add nsw i32 %120, 2
%122 = load i32, i32* %8
%123 = mul nsw i32 %121, %122
%124 = load i32, i32* %10
%125 = add nsw i32 %123, %124
%126 = add nsw i32 %125, 2
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds double, double* %119, i64 %127
%129 = load double, double* %128
store double %129, double* %20
store i32 0, i32* %11
br label %130
130: ; preds = %229, %36
%131 = load i32, i32* %11
%132 = load i32, i32* %8
%133 = icmp slt i32 %131, %132
br i1 %133, label %134, label %232
134: ; preds = %130
%135 = load double*, double** %5
%136 = load i32, i32* %9
%137 = load i32, i32* %8
%138 = mul nsw i32 %136, %137
%139 = load i32, i32* %11
%140 = add nsw i32 %138, %139
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds double, double* %135, i64 %141
%143 = load double, double* %142
store double %143, double* %21
%144 = load double*, double** %5
%145 = load i32, i32* %9
%146 = add nsw i32 %145, 1
%147 = load i32, i32* %8
%148 = mul nsw i32 %146, %147
%149 = load i32, i32* %11
%150 = add nsw i32 %148, %149
%151 = sext i32 %150 to i64
%152 = getelementptr inbounds double, double* %144, i64 %151
%153 = load double, double* %152
store double %153, double* %22
%154 = load double*, double** %5
%155 = load i32, i32* %9
%156 = add nsw i32 %155, 2
%157 = load i32, i32* %8
%158 = mul nsw i32 %156, %157
%159 = load i32, i32* %11
%160 = add nsw i32 %158, %159
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds double, double* %154, i64 %161
%163 = load double, double* %162
store double %163, double* %23
%164 = load double*, double** %6
%165 = load i32, i32* %11
%166 = load i32, i32* %8
%167 = mul nsw i32 %165, %166
%168 = load i32, i32* %10
%169 = add nsw i32 %167, %168
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds double, double* %164, i64 %170
%172 = load double, double* %171
store double %172, double* %24
%173 = load double*, double** %6
%174 = load i32, i32* %11
%175 = load i32, i32* %8
%176 = mul nsw i32 %174, %175
%177 = load i32, i32* %10
%178 = add nsw i32 %176, %177
%179 = add nsw i32 %178, 1
%180 = sext i32 %179 to i64
%181 = getelementptr inbounds double, double* %173, i64 %180
%182 = load double, double* %181
store double %182, double* %25
%183 = load double*, double** %6
%184 = load i32, i32* %11
%185 = load i32, i32* %8
%186 = mul nsw i32 %184, %185
%187 = load i32, i32* %10
%188 = add nsw i32 %186, %187
%189 = add nsw i32 %188, 2
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %183, i64 %190
%192 = load double, double* %191
store double %192, double* %26
%193 = load double, double* %21
%194 = load double, double* %24
%195 = load double, double* %12
%196 = call double @llvm.fmuladd.f64(double %193, double %194, double %195)
store double %196, double* %12
%197 = load double, double* %21
%198 = load double, double* %25
%199 = load double, double* %13
%200 = call double @llvm.fmuladd.f64(double %197, double %198, double %199)
store double %200, double* %13
%201 = load double, double* %21
%202 = load double, double* %26
%203 = load double, double* %14
%204 = call double @llvm.fmuladd.f64(double %201, double %202, double %203)
store double %204, double* %14
%205 = load double, double* %22
%206 = load double, double* %24
%207 = load double, double* %15
%208 = call double @llvm.fmuladd.f64(double %205, double %206, double %207)
store double %208, double* %15
%209 = load double, double* %22
%210 = load double, double* %25
%211 = load double, double* %16
%212 = call double @llvm.fmuladd.f64(double %209, double %210, double %211)
store double %212, double* %16
%213 = load double, double* %22
%214 = load double, double* %26
%215 = load double, double* %17
%216 = call double @llvm.fmuladd.f64(double %213, double %214, double %215)
store double %216, double* %17
%217 = load double, double* %23
%218 = load double, double* %24
%219 = load double, double* %18
%220 = call double @llvm.fmuladd.f64(double %217, double %218, double %219)
store double %220, double* %18
%221 = load double, double* %23
%222 = load double, double* %25
%223 = load double, double* %19
%224 = call double @llvm.fmuladd.f64(double %221, double %222, double %223)
store double %224, double* %19
%225 = load double, double* %23
%226 = load double, double* %26
%227 = load double, double* %20
%228 = call double @llvm.fmuladd.f64(double %225, double %226, double %227)
store double %228, double* %20
br label %229
229: ; preds = %134
%230 = load i32, i32* %11
%231 = add nsw i32 %230, 1
store i32 %231, i32* %11
br label %130
232: ; preds = %130
%233 = load double, double* %12
%234 = load double*, double** %7
%235 = load i32, i32* %9
%236 = load i32, i32* %8
%237 = mul nsw i32 %235, %236
%238 = load i32, i32* %10
%239 = add nsw i32 %237, %238
%240 = sext i32 %239 to i64
%241 = getelementptr inbounds double, double* %234, i64 %240
store double %233, double* %241
%242 = load double, double* %13
%243 = load double*, double** %7
%244 = load i32, i32* %9
%245 = load i32, i32* %8
%246 = mul nsw i32 %244, %245
%247 = load i32, i32* %10
%248 = add nsw i32 %246, %247
%249 = add nsw i32 %248, 1
%250 = sext i32 %249 to i64
%251 = getelementptr inbounds double, double* %243, i64 %250
store double %242, double* %251
%252 = load double, double* %14
%253 = load double*, double** %7
%254 = load i32, i32* %9
%255 = load i32, i32* %8
%256 = mul nsw i32 %254, %255
%257 = load i32, i32* %10
%258 = add nsw i32 %256, %257
%259 = add nsw i32 %258, 2
%260 = sext i32 %259 to i64
%261 = getelementptr inbounds double, double* %253, i64 %260
store double %252, double* %261
%262 = load double, double* %15
%263 = load double*, double** %7
%264 = load i32, i32* %9
%265 = add nsw i32 %264, 1
%266 = load i32, i32* %8
%267 = mul nsw i32 %265, %266
%268 = load i32, i32* %10
%269 = add nsw i32 %267, %268
%270 = sext i32 %269 to i64
%271 = getelementptr inbounds double, double* %263, i64 %270
store double %262, double* %271
%272 = load double, double* %16
%273 = load double*, double** %7
%274 = load i32, i32* %9
%275 = add nsw i32 %274, 1
%276 = load i32, i32* %8
%277 = mul nsw i32 %275, %276
%278 = load i32, i32* %10
%279 = add nsw i32 %277, %278
%280 = add nsw i32 %279, 1
%281 = sext i32 %280 to i64
%282 = getelementptr inbounds double, double* %273, i64 %281
store double %272, double* %282
%283 = load double, double* %17
%284 = load double*, double** %7
%285 = load i32, i32* %9
%286 = add nsw i32 %285, 1
%287 = load i32, i32* %8
%288 = mul nsw i32 %286, %287
%289 = load i32, i32* %10
%290 = add nsw i32 %288, %289
%291 = add nsw i32 %290, 2
%292 = sext i32 %291 to i64
%293 = getelementptr inbounds double, double* %284, i64 %292
store double %283, double* %293
%294 = load double, double* %18
%295 = load double*, double** %7
%296 = load i32, i32* %9
%297 = add nsw i32 %296, 2
%298 = load i32, i32* %8
%299 = mul nsw i32 %297, %298
%300 = load i32, i32* %10
%301 = add nsw i32 %299, %300
%302 = sext i32 %301 to i64
%303 = getelementptr inbounds double, double* %295, i64 %302
store double %294, double* %303
%304 = load double, double* %19
%305 = load double*, double** %7
%306 = load i32, i32* %9
%307 = add nsw i32 %306, 2
%308 = load i32, i32* %8
%309 = mul nsw i32 %307, %308
%310 = load i32, i32* %10
%311 = add nsw i32 %309, %310
%312 = add nsw i32 %311, 1
%313 = sext i32 %312 to i64
%314 = getelementptr inbounds double, double* %305, i64 %313
store double %304, double* %314
%315 = load double, double* %20
%316 = load double*, double** %7
%317 = load i32, i32* %9
%318 = add nsw i32 %317, 2
%319 = load i32, i32* %8
%320 = mul nsw i32 %318, %319
%321 = load i32, i32* %10
%322 = add nsw i32 %320, %321
%323 = add nsw i32 %322, 2
%324 = sext i32 %323 to i64
%325 = getelementptr inbounds double, double* %316, i64 %324
store double %315, double* %325
br label %326
326: ; preds = %232
%327 = load i32, i32* %10
%328 = add nsw i32 %327, 3
store i32 %328, i32* %10
br label %32
329: ; preds = %32
br label %330
330: ; preds = %329
%331 = load i32, i32* %9
%332 = add nsw i32 %331, 3
store i32 %332, i32* %9
br label %27
333: ; preds = %27
ret i32 0
}

LLVM IR (-O2):
define i32 @dgemm3(double* nocapture %0, double* nocapture %1, double* nocapture %2, i32 %3) {
%5 = icmp sgt i32 %3, 0
br i1 %5, label %6, label %90
6: ; preds = %4
%7 = zext i32 %3 to i64
%8 = zext i32 %3 to i64
br label %9
9: ; preds = %87, %6
%10 = phi i64 [ 0, %6 ], [ %88, %87 ]
%11 = mul nsw i64 %10, %7
%12 = add nuw nsw i64 %10, 1
%13 = mul nsw i64 %12, %7
%14 = add nuw nsw i64 %10, 2
%15 = mul nsw i64 %14, %7
br label %16
16: ; preds = %78, %9
%17 = phi i64 [ 0, %9 ], [ %85, %78 ]
%18 = add nuw nsw i64 %17, %11
%19 = getelementptr inbounds double, double* %2, i64 %18
%20 = add nuw nsw i64 %18, 2
%21 = getelementptr inbounds double, double* %2, i64 %20
%22 = add nuw nsw i64 %17, %13
%23 = getelementptr inbounds double, double* %2, i64 %22
%24 = add nuw nsw i64 %22, 2
%25 = getelementptr inbounds double, double* %2, i64 %24
%26 = add nuw nsw i64 %17, %15
%27 = getelementptr inbounds double, double* %2, i64 %26
%28 = add nuw nsw i64 %26, 2
%29 = getelementptr inbounds double, double* %2, i64 %28
%30 = load double, double* %29
%31 = bitcast double* %27 to <2 x double>*
%32 = load <2 x double>, <2 x double>* %31
%33 = load double, double* %25
%34 = bitcast double* %23 to <2 x double>*
%35 = load <2 x double>, <2 x double>* %34
%36 = load double, double* %21
%37 = bitcast double* %19 to <2 x double>*
%38 = load <2 x double>, <2 x double>* %37
br label %39
39: ; preds = %16, %39
%40 = phi i64 [ 0, %16 ], [ %76, %39 ]
%41 = phi double [ %36, %16 ], [ %67, %39 ]
%42 = phi double [ %30, %16 ], [ %75, %39 ]
%43 = phi double [ %33, %16 ], [ %71, %39 ]
%44 = phi <2 x double> [ %32, %16 ], [ %74, %39 ]
%45 = phi <2 x double> [ %35, %16 ], [ %70, %39 ]
%46 = phi <2 x double> [ %38, %16 ], [ %66, %39 ]
%47 = add nuw nsw i64 %40, %11
%48 = getelementptr inbounds double, double* %0, i64 %47
%49 = load double, double* %48
%50 = add nuw nsw i64 %40, %13
%51 = getelementptr inbounds double, double* %0, i64 %50
%52 = load double, double* %51
%53 = add nuw nsw i64 %40, %15
%54 = getelementptr inbounds double, double* %0, i64 %53
%55 = load double, double* %54
%56 = mul nsw i64 %40, %7
%57 = add nuw nsw i64 %56, %17
%58 = getelementptr inbounds double, double* %1, i64 %57
%59 = bitcast double* %58 to <2 x double>*
%60 = load <2 x double>, <2 x double>* %59
%61 = add nuw nsw i64 %57, 2
%62 = getelementptr inbounds double, double* %1, i64 %61
%63 = load double, double* %62
%64 = insertelement <2 x double> poison, double %49, i64 0
%65 = shufflevector <2 x double> %64, <2 x double> poison, <2 x i32> zeroinitializer
%66 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %65, <2 x double> %60, <2 x double> %46)
%67 = tail call double @llvm.fmuladd.f64(double %49, double %63, double %41)
%68 = insertelement <2 x double> poison, double %52, i64 0
%69 = shufflevector <2 x double> %68, <2 x double> poison, <2 x i32> zeroinitializer
%70 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %69, <2 x double> %60, <2 x double> %45)
%71 = tail call double @llvm.fmuladd.f64(double %52, double %63, double %43)
%72 = insertelement <2 x double> poison, double %55, i64 0
%73 = shufflevector <2 x double> %72, <2 x double> poison, <2 x i32> zeroinitializer
%74 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %73, <2 x double> %60, <2 x double> %44)
%75 = tail call double @llvm.fmuladd.f64(double %55, double %63, double %42)
%76 = add nuw nsw i64 %40, 1
%77 = icmp eq i64 %76, %8
br i1 %77, label %78, label %39
78: ; preds = %39
%79 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %73, <2 x double> %60, <2 x double> %44)
%80 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %69, <2 x double> %60, <2 x double> %45)
%81 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %65, <2 x double> %60, <2 x double> %46)
%82 = bitcast double* %19 to <2 x double>*
store <2 x double> %81, <2 x double>* %82
store double %67, double* %21
%83 = bitcast double* %23 to <2 x double>*
store <2 x double> %80, <2 x double>* %83
store double %71, double* %25
%84 = bitcast double* %27 to <2 x double>*
store <2 x double> %79, <2 x double>* %84
store double %75, double* %29
%85 = add nuw nsw i64 %17, 3
%86 = icmp ult i64 %85, %7
br i1 %86, label %16, label %87
87: ; preds = %78
%88 = add nuw nsw i64 %10, 3
%89 = icmp ult i64 %88, %7
br i1 %89, label %9, label %90
90: ; preds = %87, %4
ret i32 0
}