void dgemm2(double *a, double *b, double *c, int n) {
  int i, j, k;
  for (i = 0; i < n; i += 2) {
    for (j = 0; j < n; j += 2) {
      for (k = 0; k < n; k += 2) {
        c[i * n + j] = a[i * n + k] * b[k * n + j] +
                       a[i * n + (k + 1)] * b[(k + 1) * n + j] + c[i * n + j];
        c[(i + 1) * n + j] = a[(i + 1) * n + k] * b[k * n + j] +
                             a[(i + 1) * n + (k + 1)] * b[(k + 1) * n + j] +
                             c[(i + 1) * n + j];
        c[i * n + (j + 1)] = a[i * n + k] * b[k * n + (j + 1)] +
                             a[i * n + (k + 1)] * b[(k + 1) * n + (j + 1)] +
                             c[i * n + (j + 1)];
        c[(i + 1) * n + (j + 1)] =
            a[(i + 1) * n + k] * b[k * n + (j + 1)] +
            a[(i + 1) * n + (k + 1)] * b[(k + 1) * n + (j + 1)] +
            c[(i + 1) * n + (j + 1)];
      }
    }
  }
}


x86 assembly (clang, -O0):
<dgemm2>:
push   %rbp
mov    %rsp,%rbp
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    1410 <dgemm2+0x310>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13fd <dgemm2+0x2fd>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13ea <dgemm2+0x2ea>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x28(%rbp),%edx
add    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%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)
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
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x28(%rbp),%edx
add    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%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)
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
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),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x28(%rbp),%edx
add    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
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
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),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x28(%rbp),%edx
add    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x24(%rbp),%edx
add    $0x1,%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%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
addsd  (%rax,%rcx,8),%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)
mov    -0x28(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x28(%rbp)
jmp    1140 <dgemm2+0x40>
jmp    13ef <dgemm2+0x2ef>
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    112d <dgemm2+0x2d>
jmp    1402 <dgemm2+0x302>
mov    -0x20(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x20(%rbp)
jmp    111a <dgemm2+0x1a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dgemm2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rsi,-0x10(%rsp)
test   %ecx,%ecx
jle    1259 <dgemm2+0x159>
mov    %ecx,%eax
lea    (%rdi,%rax,8),%rcx
add    $0x8,%rcx
mov    %rax,%r8
shl    $0x4,%r8
addq   $0x8,-0x10(%rsp)
add    $0x8,%rdi
xor    %r11d,%r11d
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r11,%r10
imul   %rax,%r10
mov    %r11,-0x8(%rsp)
or     $0x1,%r11
imul   %rax,%r11
mov    -0x10(%rsp),%rsi
xor    %r15d,%r15d
nopl   0x0(%rax)
lea    (%r15,%r10,1),%r12
lea    (%r15,%r11,1),%r13
mov    %r15,%rbx
or     $0x1,%rbx
lea    (%rbx,%r10,1),%rbp
add    %r11,%rbx
mov    %rsi,%r14
xor    %r9d,%r9d
nopl   0x0(%rax)
movsd  -0x8(%rdi,%r9,8),%xmm0
movsd  (%rdi,%r9,8),%xmm1
mulsd  -0x8(%r14,%rax,8),%xmm1
mulsd  -0x8(%r14),%xmm0
addsd  %xmm1,%xmm0
addsd  (%rdx,%r12,8),%xmm0
movsd  %xmm0,(%rdx,%r12,8)
movsd  -0x8(%rcx,%r9,8),%xmm0
movsd  (%rcx,%r9,8),%xmm1
mulsd  -0x8(%r14,%rax,8),%xmm1
mulsd  -0x8(%r14),%xmm0
addsd  %xmm1,%xmm0
addsd  (%rdx,%r13,8),%xmm0
movsd  %xmm0,(%rdx,%r13,8)
movsd  -0x8(%rdi,%r9,8),%xmm0
movsd  (%rdi,%r9,8),%xmm1
mulsd  (%r14,%rax,8),%xmm1
mulsd  (%r14),%xmm0
addsd  %xmm1,%xmm0
addsd  (%rdx,%rbp,8),%xmm0
movsd  %xmm0,(%rdx,%rbp,8)
movsd  -0x8(%rcx,%r9,8),%xmm0
movsd  (%rcx,%r9,8),%xmm1
mulsd  (%r14,%rax,8),%xmm1
mulsd  (%r14),%xmm0
addsd  %xmm1,%xmm0
addsd  (%rdx,%rbx,8),%xmm0
movsd  %xmm0,(%rdx,%rbx,8)
add    $0x2,%r9
add    %r8,%r14
cmp    %rax,%r9
jb     1180 <dgemm2+0x80>
add    $0x2,%r15
add    $0x10,%rsi
cmp    %rax,%r15
jb     1160 <dgemm2+0x60>
mov    -0x8(%rsp),%r11
add    $0x2,%r11
add    %r8,%rcx
add    %r8,%rdi
cmp    %rax,%r11
jb     1140 <dgemm2+0x40>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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