void dgemm4(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 c0 = c[i * n + j];
      register double c1 = c[i * n + j + 1];
      register double c2 = c[i * n + j + 2];
      register double c3 = c[(i + 1) * n + j];
      register double c4 = c[(i + 1) * n + (j + 1)];
      register double c5 = c[(i + 1) * n + (j + 2)];
      register double c6 = c[(i + 2) * n + j];
      register double c7 = c[(i + 2) * n + (j + 1)];
      register double c8 = c[(i + 2) * n + (j + 2)];
      for (k = 0; k < n; k++) {
        register double a0 = a[i * n + k];
        register double a1 = a[(i + 1) * n + k];
        register double a2 = a[(i + 2) * n + k];
        register double b0 = b[k * n + j];
        register double b1 = b[k * n + (j + 1)];
        register double b2 = b[k * n + (j + 2)];
        c0 += a0 * b0;
        c1 += a0 * b1;
        c2 += a0 * b2;
        c3 += a1 * b0;
        c4 += a1 * b1;
        c5 += a1 * b2;
        c6 += a2 * b0;
        c7 += a2 * b1;
        c8 += a2 * b2;
      }
      c[i * n + j] = c0;
      c[i * n + j + 1] = c1;
      c[i * n + j + 2] = c2;
      c[(i + 1) * n + j] = c3;
      c[(i + 1) * n + (j + 1)] = c4;
      c[(i + 1) * n + (j + 2)] = c5;
      c[(i + 2) * n + j] = c6;
      c[(i + 2) * n + (j + 1)] = c7;
      c[(i + 2) * n + (j + 2)] = c8;
    }
  }
}


x86 assembly (clang, -O0):
<dgemm4>:
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    159b <dgemm4+0x49b>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1588 <dgemm4+0x488>
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
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%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
mov    -0x24(%rbp),%edx
add    $0x2,%edx
add    %edx,%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
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%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
mov    -0x24(%rbp),%edx
add    $0x2,%edx
add    %edx,%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    145b <dgemm4+0x35b>
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
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%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
mov    -0x24(%rbp),%edx
add    $0x2,%edx
add    %edx,%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    1263 <dgemm4+0x163>
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
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%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
mov    -0x24(%rbp),%edx
add    $0x2,%edx
add    %edx,%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
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%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
mov    -0x24(%rbp),%edx
add    $0x2,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x3,%eax
mov    %eax,-0x24(%rbp)
jmp    1131 <dgemm4+0x31>
jmp    158d <dgemm4+0x48d>
mov    -0x20(%rbp),%eax
add    $0x3,%eax
mov    %eax,-0x20(%rbp)
jmp    111e <dgemm4+0x1e>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dgemm4>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rsi,-0x30(%rsp)
test   %ecx,%ecx
jle    12de <dgemm4+0x1de>
mov    %ecx,%eax
addq   $0x10,-0x30(%rsp)
lea    0x0(,%rax,8),%r8
mov    %rax,%r10
shl    $0x4,%r10
add    %rdi,%r10
lea    (%r8,%r8,2),%rcx
mov    %rcx,-0x28(%rsp)
lea    (%rdi,%rax,8),%rcx
xor    %ebp,%ebp
mov    %rbp,%rsi
imul   %rax,%rsi
mov    %rsi,-0x8(%rsp)
lea    0x1(%rbp),%rsi
imul   %rax,%rsi
mov    %rsi,-0x10(%rsp)
mov    %rbp,-0x20(%rsp)
lea    0x2(%rbp),%rsi
imul   %rax,%rsi
mov    %rsi,-0x18(%rsp)
mov    -0x30(%rsp),%r15
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x8(%rsp),%rsi
lea    (%rbx,%rsi,1),%r9
mov    -0x10(%rsp),%rsi
lea    (%rbx,%rsi,1),%r12
lea    (%rbx,%rsi,1),%rbp
add    $0x2,%rbp
mov    -0x18(%rsp),%rsi
lea    (%rbx,%rsi,1),%r11
lea    0x2(%rbx,%rsi,1),%r14
mov    %rbp,%rsi
mov    %r12,%rbp
movsd  (%rdx,%r14,8),%xmm11
movupd (%rdx,%r11,8),%xmm0
movsd  (%rdx,%rsi,8),%xmm12
movupd (%rdx,%r12,8),%xmm3
movsd  0x10(%rdx,%r9,8),%xmm13
movupd (%rdx,%r9,8),%xmm7
mov    %r15,%r13
xor    %r12d,%r12d
nopl   0x0(%rax)
movapd %xmm7,%xmm8
movapd %xmm3,%xmm9
movapd %xmm0,%xmm10
movsd  (%rdi,%r12,8),%xmm0
movsd  (%rcx,%r12,8),%xmm5
movsd  (%r10,%r12,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,%r12
add    %r8,%r13
cmp    %r12,%rax
jne    11e0 <dgemm4+0xe0>
addpd  %xmm10,%xmm5
addpd  %xmm9,%xmm4
addpd  %xmm8,%xmm6
movupd %xmm6,(%rdx,%r9,8)
movsd  %xmm13,0x10(%rdx,%r9,8)
movupd %xmm4,(%rdx,%rbp,8)
movsd  %xmm12,(%rdx,%rsi,8)
movupd %xmm5,(%rdx,%r11,8)
movsd  %xmm11,(%rdx,%r14,8)
add    $0x3,%rbx
add    $0x18,%r15
cmp    %rax,%rbx
jb     1180 <dgemm4+0x80>
mov    -0x20(%rsp),%rbp
add    $0x3,%rbp
mov    -0x28(%rsp),%rsi
add    %rsi,%r10
add    %rsi,%rcx
add    %rsi,%rdi
cmp    %rax,%rbp
jb     1140 <dgemm4+0x40>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @dgemm4(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 %81, 1
%83 = add nsw i32 %80, %82
%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 %92, 2
%94 = add nsw i32 %91, %93
%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 %113, 1
%115 = add nsw i32 %112, %114
%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 %124, 2
%126 = add nsw i32 %123, %125
%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 %177, 1
%179 = add nsw i32 %176, %178
%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 %187, 2
%189 = add nsw i32 %186, %188
%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 %278, 1
%280 = add nsw i32 %277, %279
%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 %289, 2
%291 = add nsw i32 %288, %290
%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 %310, 1
%312 = add nsw i32 %309, %311
%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 %321, 2
%323 = add nsw i32 %320, %322
%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 void
}

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