void dgemm(const int M, const int N, const int K, const double alpha,
           const double *A, const int lda, const double *B, const int ldb,
           const double beta, double *C, const int ldc) {
  int ib, jb, kb;
  int i, j, k;
  int ii, jj, kk;
  double a, b, c;
  double At[K * M];
  for (k = 0; k < K; ++k)
    for (i = 0; i < M; ++i)
      At[i * M + k] = A[k * lda + i];
  for (kb = 0; kb < K; kb += 64) {
    for (ib = 0; ib < M; ib += 64) {
      for (jb = 0; jb < N; jb += 64) {
        int Kb = (((64) < (K - kb)) ? (64) : (K - kb));
        int Nb = (((64) < (N - jb)) ? (64) : (N - jb));
        int Mb = (((64) < (M - ib)) ? (64) : (M - ib));
        for (i = 0; i < Mb; i = i + 1)
          for (j = 0; j < Nb; j = j + 1) {
            double cij = 0.;
            for (k = 0; k < K; ++k)
              cij += At[i * M + k] * B[k + j * ldb];
            C[j * ldc + i] += cij;
          }
      }
    }
  }
}


x86 assembly (clang, -O0):
<dgemm>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   %ecx,%eax
mov    %rsp,%rcx
mov    %rcx,-0x80(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0xa8(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x88(%rbp)
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11dd <dgemm+0xdd>
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11ca <dgemm+0xca>
mov    -0xa8(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x50(%rbp),%edx
imul   -0x24(%rbp),%edx
add    -0x48(%rbp),%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm0
mov    -0x48(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1181 <dgemm+0x81>
jmp    11cf <dgemm+0xcf>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    116e <dgemm+0x6e>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13e5 <dgemm+0x2e5>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13d2 <dgemm+0x2d2>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    13bf <dgemm+0x2bf>
mov    -0xc(%rbp),%ecx
sub    -0x44(%rbp),%ecx
mov    $0x40,%eax
cmp    %ecx,%eax
jge    1239 <dgemm+0x139>
mov    $0x40,%eax
mov    %eax,-0xac(%rbp)
jmp    1245 <dgemm+0x145>
mov    -0xc(%rbp),%eax
sub    -0x44(%rbp),%eax
mov    %eax,-0xac(%rbp)
mov    -0xac(%rbp),%eax
mov    %eax,-0x8c(%rbp)
mov    -0x8(%rbp),%ecx
sub    -0x40(%rbp),%ecx
mov    $0x40,%eax
cmp    %ecx,%eax
jge    1274 <dgemm+0x174>
mov    $0x40,%eax
mov    %eax,-0xb0(%rbp)
jmp    1280 <dgemm+0x180>
mov    -0x8(%rbp),%eax
sub    -0x40(%rbp),%eax
mov    %eax,-0xb0(%rbp)
mov    -0xb0(%rbp),%eax
mov    %eax,-0x90(%rbp)
mov    -0x4(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
mov    $0x40,%eax
cmp    %ecx,%eax
jge    12af <dgemm+0x1af>
mov    $0x40,%eax
mov    %eax,-0xb4(%rbp)
jmp    12bb <dgemm+0x1bb>
mov    -0x4(%rbp),%eax
sub    -0x3c(%rbp),%eax
mov    %eax,-0xb4(%rbp)
mov    -0xb4(%rbp),%eax
mov    %eax,-0x94(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x94(%rbp),%eax
jge    13ac <dgemm+0x2ac>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x90(%rbp),%eax
jge    1399 <dgemm+0x299>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1368 <dgemm+0x268>
mov    -0xa8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x50(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x4c(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1305 <dgemm+0x205>
movsd  -0xa0(%rbp),%xmm0
ff
mov    0x18(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   0x20(%rbp),%ecx
add    -0x48(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12e4 <dgemm+0x1e4>
jmp    139e <dgemm+0x29e>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    12ce <dgemm+0x1ce>
jmp    13b1 <dgemm+0x2b1>
mov    -0x40(%rbp),%eax
add    $0x40,%eax
mov    %eax,-0x40(%rbp)
jmp    120a <dgemm+0x10a>
jmp    13c4 <dgemm+0x2c4>
mov    -0x3c(%rbp),%eax
add    $0x40,%eax
mov    %eax,-0x3c(%rbp)
jmp    11f7 <dgemm+0xf7>
jmp    13d7 <dgemm+0x2d7>
mov    -0x44(%rbp),%eax
add    $0x40,%eax
mov    %eax,-0x44(%rbp)
jmp    11e4 <dgemm+0xe4>
mov    -0x80(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dgemm>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9,-0x68(%rbp)
mov    %rcx,%r10
mov    %esi,-0x34(%rbp)
mov    %edx,%eax
mov    %edi,-0x2c(%rbp)
imul   %edi,%eax
lea    0xf(,%rax,8),%rax
and    $0xfffffffffffffff0,%rax
mov    %rsp,%rcx
sub    %rax,%rcx
mov    %rcx,-0x40(%rbp)
mov    %rcx,%rsp
test   %edx,%edx
jle    140f <dgemm+0x30f>
mov    %edx,%r14d
mov    0x18(%rbp),%rcx
movslq -0x2c(%rbp),%rsi
movslq %r8d,%rdx
mov    %r14d,%eax
mov    %rax,-0x58(%rbp)
mov    %esi,%r15d
mov    %r15d,%r9d
and    $0xfffffffe,%r9d
mov    %r10,-0x50(%rbp)
lea    0x8(%r10),%rdi
mov    %rdx,%rbx
lea    0x0(,%rdx,8),%r12
mov    %rsi,%r8
shl    $0x4,%r8
xor    %r13d,%r13d
mov    -0x40(%rbp),%r10
mov    %rsi,-0x48(%rbp)
jmp    11a5 <dgemm+0xa5>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r13
add    %r12,%rdi
add    $0x8,%r10
cmp    -0x58(%rbp),%r13
je     1224 <dgemm+0x124>
cmpl   $0x0,-0x2c(%rbp)
jle    1190 <dgemm+0x90>
cmpl   $0x1,-0x2c(%rbp)
jne    11c0 <dgemm+0xc0>
xor    %eax,%eax
jmp    11f2 <dgemm+0xf2>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r10,%r11
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x8(%rdi,%rax,8),%xmm0
movsd  %xmm0,(%r11)
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%r11,%rsi,8)
add    $0x2,%rax
add    %r8,%r11
cmp    %rax,%r9
jne    11d0 <dgemm+0xd0>
test   $0x1,%r15b
je     1190 <dgemm+0x90>
mov    %r13,%rsi
imul   %rbx,%rsi
add    %rax,%rsi
mov    -0x50(%rbp),%rdx
movsd  (%rdx,%rsi,8),%xmm0
mov    -0x48(%rbp),%rsi
imul   %rsi,%rax
add    %r13,%rax
mov    -0x40(%rbp),%rdx
movsd  %xmm0,(%rdx,%rax,8)
jmp    1190 <dgemm+0x90>
test   %r14d,%r14d
jle    140f <dgemm+0x30f>
movslq 0x10(%rbp),%rsi
movslq 0x20(%rbp),%rdi
mov    -0x58(%rbp),%rax
and    $0xfffffffe,%eax
mov    -0x40(%rbp),%rdx
add    $0x8,%rdx
mov    %rdx,-0x78(%rbp)
mov    -0x48(%rbp),%rdx
lea    0x0(,%rdx,8),%rdx
mov    %rdx,-0x88(%rbp)
mov    -0x68(%rbp),%rdx
add    $0x8,%rdx
mov    %rdx,-0x80(%rbp)
mov    %rsi,-0x98(%rbp)
lea    0x0(,%rsi,8),%r15
xor    %edx,%edx
jmp    1290 <dgemm+0x190>
nopw   0x0(%rax,%rax,1)
mov    -0x70(%rbp),%rdx
add    $0x40,%edx
cmp    %r14d,%edx
jge    140f <dgemm+0x30f>
mov    %rdx,-0x70(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jle    1280 <dgemm+0x180>
movl   $0x0,-0x30(%rbp)
jmp    12be <dgemm+0x1be>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    -0x30(%rbp),%edx
add    $0x40,%edx
mov    %edx,-0x30(%rbp)
cmp    -0x2c(%rbp),%edx
jge    1280 <dgemm+0x180>
cmpl   $0x0,-0x34(%rbp)
jle    12b0 <dgemm+0x1b0>
mov    -0x2c(%rbp),%edx
sub    -0x30(%rbp),%edx
cmp    $0x40,%edx
mov    $0x40,%esi
mov    %edx,-0x5c(%rbp)
cmovl  %edx,%esi
movslq %esi,%rdx
mov    %rdx,-0x90(%rbp)
movl   $0x0,-0x38(%rbp)
jmp    12fe <dgemm+0x1fe>
nopl   0x0(%rax,%rax,1)
mov    -0x38(%rbp),%edx
add    $0x40,%edx
mov    %edx,-0x38(%rbp)
cmp    -0x34(%rbp),%edx
jge    12b0 <dgemm+0x1b0>
cmpl   $0x0,-0x5c(%rbp)
jle    12f0 <dgemm+0x1f0>
mov    -0x34(%rbp),%edx
sub    -0x38(%rbp),%edx
cmp    $0x40,%edx
mov    $0x40,%esi
mov    %edx,-0x60(%rbp)
cmovl  %edx,%esi
movslq %esi,%rsi
mov    -0x78(%rbp),%r13
xor    %r12d,%r12d
jmp    1344 <dgemm+0x244>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r12
add    -0x88(%rbp),%r13
cmp    -0x90(%rbp),%r12
jge    12f0 <dgemm+0x1f0>
cmpl   $0x0,-0x60(%rbp)
jle    1330 <dgemm+0x230>
mov    %r12,%rdx
imul   -0x48(%rbp),%rdx
mov    %rdx,-0x50(%rbp)
mov    -0x80(%rbp),%r11
xor    %r8d,%r8d
jmp    1380 <dgemm+0x280>
nop
mov    %r8,%rdx
imul   %rdi,%rdx
add    %r12,%rdx
addsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x1,%r8
add    %r15,%r11
cmp    %rsi,%r8
jge    1330 <dgemm+0x230>
cmp    $0x1,%r14d
jne    1390 <dgemm+0x290>
xorpd  %xmm0,%xmm0
xor    %r9d,%r9d
jmp    13d0 <dgemm+0x2d0>
nop
xorpd  %xmm1,%xmm1
xor    %r9d,%r9d
nopw   0x0(%rax,%rax,1)
movsd  -0x8(%r13,%r9,8),%xmm2
movsd  0x0(%r13,%r9,8),%xmm0
mulsd  -0x8(%r11,%r9,8),%xmm2
mulsd  (%r11,%r9,8),%xmm0
addsd  %xmm1,%xmm2
addsd  %xmm2,%xmm0
add    $0x2,%r9
movapd %xmm0,%xmm1
cmp    %r9,%rax
jne    13a0 <dgemm+0x2a0>
testb  $0x1,-0x58(%rbp)
je     1360 <dgemm+0x260>
mov    %r8,%r10
imul   -0x98(%rbp),%r10
ff
mov    %r14d,%ebx
mov    -0x50(%rbp),%rdx
lea    (%r9,%rdx,1),%r14
mov    -0x40(%rbp),%rdx
movsd  (%rdx,%r14,8),%xmm1
mov    %ebx,%r14d
add    %r9,%r10
mov    -0x68(%rbp),%rdx
mulsd  (%rdx,%r10,8),%xmm1
addsd  %xmm1,%xmm0
jmp    1360 <dgemm+0x260>
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @dgemm(i32 %0, i32 %1, i32 %2, double %3, double* %4, i32 %5, double* %6, i32 %7, double %8, double* %9, i32 %10) {
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca double
%16 = alloca double*
%17 = alloca i32
%18 = alloca double*
%19 = alloca i32
%20 = alloca double
%21 = alloca double*
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca i32
%31 = alloca i32
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca i8*
%36 = alloca i64
%37 = alloca i32
%38 = alloca i32
%39 = alloca i32
%40 = alloca double
store i32 %0, i32* %12
store i32 %1, i32* %13
store i32 %2, i32* %14
store double %3, double* %15
store double* %4, double** %16
store i32 %5, i32* %17
store double* %6, double** %18
store i32 %7, i32* %19
store double %8, double* %20
store double* %9, double** %21
store i32 %10, i32* %22
%41 = load i32, i32* %14
%42 = load i32, i32* %12
%43 = mul nsw i32 %41, %42
%44 = zext i32 %43 to i64
%45 = call i8* @llvm.stacksave()
store i8* %45, i8** %35
%46 = alloca double, i64 %44
store i64 %44, i64* %36
store i32 0, i32* %28
br label %47
47: ; preds = %77, %11
%48 = load i32, i32* %28
%49 = load i32, i32* %14
%50 = icmp slt i32 %48, %49
br i1 %50, label %51, label %80
51: ; preds = %47
store i32 0, i32* %26
br label %52
52: ; preds = %73, %51
%53 = load i32, i32* %26
%54 = load i32, i32* %12
%55 = icmp slt i32 %53, %54
br i1 %55, label %56, label %76
56: ; preds = %52
%57 = load double*, double** %16
%58 = load i32, i32* %28
%59 = load i32, i32* %17
%60 = mul nsw i32 %58, %59
%61 = load i32, i32* %26
%62 = add nsw i32 %60, %61
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %57, i64 %63
%65 = load double, double* %64
%66 = load i32, i32* %26
%67 = load i32, i32* %12
%68 = mul nsw i32 %66, %67
%69 = load i32, i32* %28
%70 = add nsw i32 %68, %69
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds double, double* %46, i64 %71
store double %65, double* %72
br label %73
73: ; preds = %56
%74 = load i32, i32* %26
%75 = add nsw i32 %74, 1
store i32 %75, i32* %26
br label %52
76: ; preds = %52
br label %77
77: ; preds = %76
%78 = load i32, i32* %28
%79 = add nsw i32 %78, 1
store i32 %79, i32* %28
br label %47
80: ; preds = %47
store i32 0, i32* %25
br label %81
81: ; preds = %194, %80
%82 = load i32, i32* %25
%83 = load i32, i32* %14
%84 = icmp slt i32 %82, %83
br i1 %84, label %85, label %197
85: ; preds = %81
store i32 0, i32* %23
br label %86
86: ; preds = %190, %85
%87 = load i32, i32* %23
%88 = load i32, i32* %12
%89 = icmp slt i32 %87, %88
br i1 %89, label %90, label %193
90: ; preds = %86
store i32 0, i32* %24
br label %91
91: ; preds = %186, %90
%92 = load i32, i32* %24
%93 = load i32, i32* %13
%94 = icmp slt i32 %92, %93
br i1 %94, label %95, label %189
95: ; preds = %91
%96 = load i32, i32* %14
%97 = load i32, i32* %25
%98 = sub nsw i32 %96, %97
%99 = icmp slt i32 64, %98
br i1 %99, label %100, label %101
100: ; preds = %95
br label %105
101: ; preds = %95
%102 = load i32, i32* %14
%103 = load i32, i32* %25
%104 = sub nsw i32 %102, %103
br label %105
105: ; preds = %101, %100
%106 = phi i32 [ 64, %100 ], [ %104, %101 ]
store i32 %106, i32* %37
%107 = load i32, i32* %13
%108 = load i32, i32* %24
%109 = sub nsw i32 %107, %108
%110 = icmp slt i32 64, %109
br i1 %110, label %111, label %112
111: ; preds = %105
br label %116
112: ; preds = %105
%113 = load i32, i32* %13
%114 = load i32, i32* %24
%115 = sub nsw i32 %113, %114
br label %116
116: ; preds = %112, %111
%117 = phi i32 [ 64, %111 ], [ %115, %112 ]
store i32 %117, i32* %38
%118 = load i32, i32* %12
%119 = load i32, i32* %23
%120 = sub nsw i32 %118, %119
%121 = icmp slt i32 64, %120
br i1 %121, label %122, label %123
122: ; preds = %116
br label %127
123: ; preds = %116
%124 = load i32, i32* %12
%125 = load i32, i32* %23
%126 = sub nsw i32 %124, %125
br label %127
127: ; preds = %123, %122
%128 = phi i32 [ 64, %122 ], [ %126, %123 ]
store i32 %128, i32* %39
store i32 0, i32* %26
br label %129
129: ; preds = %182, %127
%130 = load i32, i32* %26
%131 = load i32, i32* %39
%132 = icmp slt i32 %130, %131
br i1 %132, label %133, label %185
133: ; preds = %129
store i32 0, i32* %27
br label %134
134: ; preds = %178, %133
%135 = load i32, i32* %27
%136 = load i32, i32* %38
%137 = icmp slt i32 %135, %136
br i1 %137, label %138, label %181
138: ; preds = %134
store double 0.000000e+00, double* %40
store i32 0, i32* %28
br label %139
139: ; preds = %163, %138
%140 = load i32, i32* %28
%141 = load i32, i32* %14
%142 = icmp slt i32 %140, %141
br i1 %142, label %143, label %166
143: ; preds = %139
%144 = load i32, i32* %26
%145 = load i32, i32* %12
%146 = mul nsw i32 %144, %145
%147 = load i32, i32* %28
%148 = add nsw i32 %146, %147
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds double, double* %46, i64 %149
%151 = load double, double* %150
%152 = load double*, double** %18
%153 = load i32, i32* %28
%154 = load i32, i32* %27
%155 = load i32, i32* %19
%156 = mul nsw i32 %154, %155
%157 = add nsw i32 %153, %156
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double, double* %152, i64 %158
%160 = load double, double* %159
%161 = load double, double* %40
%162 = call double @llvm.fmuladd.f64(double %151, double %160, double %161)
store double %162, double* %40
br label %163
163: ; preds = %143
%164 = load i32, i32* %28
%165 = add nsw i32 %164, 1
store i32 %165, i32* %28
br label %139
166: ; preds = %139
%167 = load double, double* %40
%168 = load double*, double** %21
%169 = load i32, i32* %27
%170 = load i32, i32* %22
%171 = mul nsw i32 %169, %170
%172 = load i32, i32* %26
%173 = add nsw i32 %171, %172
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds double, double* %168, i64 %174
%176 = load double, double* %175
%177 = fadd double %176, %167
store double %177, double* %175
br label %178
178: ; preds = %166
%179 = load i32, i32* %27
%180 = add nsw i32 %179, 1
store i32 %180, i32* %27
br label %134
181: ; preds = %134
br label %182
182: ; preds = %181
%183 = load i32, i32* %26
%184 = add nsw i32 %183, 1
store i32 %184, i32* %26
br label %129
185: ; preds = %129
br label %186
186: ; preds = %185
%187 = load i32, i32* %24
%188 = add nsw i32 %187, 64
store i32 %188, i32* %24
br label %91
189: ; preds = %91
br label %190
190: ; preds = %189
%191 = load i32, i32* %23
%192 = add nsw i32 %191, 64
store i32 %192, i32* %23
br label %86
193: ; preds = %86
br label %194
194: ; preds = %193
%195 = load i32, i32* %25
%196 = add nsw i32 %195, 64
store i32 %196, i32* %25
br label %81
197: ; preds = %81
%198 = load i8*, i8** %35
ret void
}

LLVM IR (-O2):
define void @dgemm(i32 %0, i32 %1, i32 %2, double %3, double* nocapture %4, i32 %5, double* nocapture %6, i32 %7, double %8, double* nocapture %9, i32 %10) {
%12 = mul nsw i32 %2, %0
%13 = zext i32 %12 to i64
%14 = alloca double, i64 %13
%15 = icmp sgt i32 %2, 0
br i1 %15, label %16, label %152
16: ; preds = %11
%17 = icmp sgt i32 %0, 0
%18 = sext i32 %0 to i64
%19 = sext i32 %5 to i64
%20 = zext i32 %2 to i64
%21 = zext i32 %0 to i64
%22 = and i64 %21, 1
%23 = icmp eq i32 %0, 1
%24 = and i64 %21, 4294967294
%25 = icmp eq i64 %22, 0
br label %26
26: ; preds = %16, %69
%27 = phi i64 [ 0, %16 ], [ %70, %69 ]
br i1 %17, label %28, label %69
28: ; preds = %26
%29 = mul nsw i64 %27, %19
br i1 %23, label %60, label %41
30: ; preds = %69
br i1 %15, label %31, label %152
31: ; preds = %30
%32 = icmp sgt i32 %0, 0
%33 = icmp sgt i32 %1, 0
%34 = sext i32 %7 to i64
%35 = sext i32 %10 to i64
%36 = sext i32 %0 to i64
%37 = and i64 %20, 1
%38 = icmp eq i32 %2, 1
%39 = and i64 %20, 4294967294
%40 = icmp eq i64 %37, 0
br label %72
41: ; preds = %28, %41
%42 = phi i64 [ %57, %41 ], [ 0, %28 ]
%43 = phi i64 [ %58, %41 ], [ 0, %28 ]
%44 = add nsw i64 %42, %29
%45 = getelementptr inbounds double, double* %4, i64 %44
%46 = load double, double* %45
%47 = mul nsw i64 %42, %18
%48 = add nsw i64 %47, %27
%49 = getelementptr inbounds double, double* %14, i64 %48
store double %46, double* %49
%50 = or i64 %42, 1
%51 = add nsw i64 %50, %29
%52 = getelementptr inbounds double, double* %4, i64 %51
%53 = load double, double* %52
%54 = mul nsw i64 %50, %18
%55 = add nsw i64 %54, %27
%56 = getelementptr inbounds double, double* %14, i64 %55
store double %53, double* %56
%57 = add nuw nsw i64 %42, 2
%58 = add i64 %43, 2
%59 = icmp eq i64 %58, %24
br i1 %59, label %60, label %41
60: ; preds = %41, %28
%61 = phi i64 [ 0, %28 ], [ %57, %41 ]
br i1 %25, label %69, label %62
62: ; preds = %60
%63 = add nsw i64 %61, %29
%64 = getelementptr inbounds double, double* %4, i64 %63
%65 = load double, double* %64
%66 = mul nsw i64 %61, %18
%67 = add nsw i64 %66, %27
%68 = getelementptr inbounds double, double* %14, i64 %67
store double %65, double* %68
br label %69
69: ; preds = %62, %60, %26
%70 = add nuw nsw i64 %27, 1
%71 = icmp eq i64 %70, %20
br i1 %71, label %30, label %26
72: ; preds = %31, %149
%73 = phi i32 [ 0, %31 ], [ %150, %149 ]
br i1 %32, label %74, label %149
74: ; preds = %72, %146
%75 = phi i32 [ %147, %146 ], [ 0, %72 ]
br i1 %33, label %76, label %146
76: ; preds = %74
%77 = sub nsw i32 %0, %75
%78 = icmp slt i32 %77, 64
%79 = select i1 %78, i32 %77, i32 64
%80 = icmp sgt i32 %77, 0
%81 = sext i32 %79 to i64
br label %82
82: ; preds = %76, %143
%83 = phi i32 [ 0, %76 ], [ %144, %143 ]
%84 = sub nsw i32 %1, %83
br i1 %80, label %85, label %143
85: ; preds = %82
%86 = icmp slt i32 %84, 64
%87 = select i1 %86, i32 %84, i32 64
%88 = icmp sgt i32 %84, 0
%89 = sext i32 %87 to i64
br label %90
90: ; preds = %85, %140
%91 = phi i64 [ 0, %85 ], [ %141, %140 ]
br i1 %88, label %92, label %140
92: ; preds = %90
%93 = mul nsw i64 %91, %36
br label %94
94: ; preds = %131, %92
%95 = phi i64 [ 0, %92 ], [ %138, %131 ]
%96 = mul nsw i64 %95, %34
br i1 %38, label %119, label %97
97: ; preds = %94, %97
%98 = phi i64 [ %116, %97 ], [ 0, %94 ]
%99 = phi double [ %115, %97 ], [ 0.000000e+00, %94 ]
%100 = phi i64 [ %117, %97 ], [ 0, %94 ]
%101 = add nsw i64 %98, %93
%102 = getelementptr inbounds double, double* %14, i64 %101
%103 = load double, double* %102
%104 = add nsw i64 %98, %96
%105 = getelementptr inbounds double, double* %6, i64 %104
%106 = load double, double* %105
%107 = tail call double @llvm.fmuladd.f64(double %103, double %106, double %99)
%108 = or i64 %98, 1
%109 = add nsw i64 %108, %93
%110 = getelementptr inbounds double, double* %14, i64 %109
%111 = load double, double* %110
%112 = add nsw i64 %108, %96
%113 = getelementptr inbounds double, double* %6, i64 %112
%114 = load double, double* %113
%115 = tail call double @llvm.fmuladd.f64(double %111, double %114, double %107)
%116 = add nuw nsw i64 %98, 2
%117 = add i64 %100, 2
%118 = icmp eq i64 %117, %39
br i1 %118, label %119, label %97
119: ; preds = %97, %94
%120 = phi double [ undef, %94 ], [ %115, %97 ]
%121 = phi i64 [ 0, %94 ], [ %116, %97 ]
%122 = phi double [ 0.000000e+00, %94 ], [ %115, %97 ]
br i1 %40, label %131, label %123
123: ; preds = %119
%124 = add nsw i64 %121, %93
%125 = getelementptr inbounds double, double* %14, i64 %124
%126 = load double, double* %125
%127 = add nsw i64 %121, %96
%128 = getelementptr inbounds double, double* %6, i64 %127
%129 = load double, double* %128
%130 = tail call double @llvm.fmuladd.f64(double %126, double %129, double %122)
br label %131
131: ; preds = %119, %123
%132 = phi double [ %120, %119 ], [ %130, %123 ]
%133 = mul nsw i64 %95, %35
%134 = add nsw i64 %133, %91
%135 = getelementptr inbounds double, double* %9, i64 %134
%136 = load double, double* %135
%137 = fadd double %132, %136
store double %137, double* %135
%138 = add nuw nsw i64 %95, 1
%139 = icmp slt i64 %138, %89
br i1 %139, label %94, label %140
140: ; preds = %131, %90
%141 = add nuw nsw i64 %91, 1
%142 = icmp slt i64 %141, %81
br i1 %142, label %90, label %143
143: ; preds = %140, %82
%144 = add nuw nsw i32 %83, 64
%145 = icmp slt i32 %144, %1
br i1 %145, label %82, label %146
146: ; preds = %143, %74
%147 = add nuw nsw i32 %75, 64
%148 = icmp slt i32 %147, %0
br i1 %148, label %74, label %149
149: ; preds = %146, %72
%150 = add nuw nsw i32 %73, 64
%151 = icmp slt i32 %150, %2
br i1 %151, label %72, label %152
152: ; preds = %149, %11, %30
ret void
}