extern void matmul(const char *tr, int n, int k, int m, double alpha,
                   const double *A, const double *B, double beta, double *C) {
  double d;
  int i, j, x,
      f = tr[0] == 'N' ? (tr[1] == 'N' ? 1 : 2) : (tr[1] == 'N' ? 3 : 4);
  for (i = 0; i < n; i++)
    for (j = 0; j < k; j++) {
      d = 0.0;
      switch (f) {
      case 1:
        for (x = 0; x < m; x++)
          d += A[i + x * n] * B[x + j * m];
        break;
      case 2:
        for (x = 0; x < m; x++)
          d += A[i + x * n] * B[j + x * k];
        break;
      case 3:
        for (x = 0; x < m; x++)
          d += A[x + i * m] * B[x + j * m];
        break;
      case 4:
        for (x = 0; x < m; x++)
          d += A[x + i * m] * B[j + x * k];
        break;
      }
      if (beta == 0.0)
        C[i + j * n] = alpha * d;
      else
        C[i + j * n] = alpha * d + beta * C[i + j * n];
    }
}


x86 assembly (clang, -O0):
<matmul>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movsbl (%rax),%eax
cmp    $0x4e,%eax
jne    1157 <matmul+0x57>
mov    -0x8(%rbp),%rax
movsbl 0x1(%rax),%edx
mov    $0x2,%eax
mov    $0x1,%ecx
cmp    $0x4e,%edx
cmove  %ecx,%eax
mov    %eax,-0x54(%rbp)
jmp    1172 <matmul+0x72>
mov    -0x8(%rbp),%rax
movsbl 0x1(%rax),%edx
mov    $0x4,%eax
mov    $0x3,%ecx
cmp    $0x4e,%edx
cmove  %ecx,%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x50(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    141b <matmul+0x31b>
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1408 <matmul+0x308>
movq   $0x0,-0x40(%rbp)
mov    -0x50(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,%ecx
mov    %rcx,-0x60(%rbp)
sub    $0x3,%eax
ja     136f <matmul+0x26f>
mov    -0x60(%rbp),%rax
lea    0xe3a(%rip),%rcx
movslq (%rcx,%rax,4),%rax
add    %rcx,%rax
jmp    *%rax
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1232 <matmul+0x132>
mov    -0x28(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4c(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x4c(%rbp),%ecx
mov    -0x48(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11d6 <matmul+0xd6>
jmp    136f <matmul+0x26f>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    129a <matmul+0x19a>
mov    -0x28(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x4c(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x48(%rbp),%ecx
mov    -0x4c(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    123e <matmul+0x13e>
jmp    136f <matmul+0x26f>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1302 <matmul+0x202>
mov    -0x28(%rbp),%rax
mov    -0x4c(%rbp),%ecx
mov    -0x44(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x4c(%rbp),%ecx
mov    -0x48(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12a6 <matmul+0x1a6>
jmp    136f <matmul+0x26f>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    136a <matmul+0x26a>
mov    -0x28(%rbp),%rax
mov    -0x4c(%rbp),%ecx
mov    -0x44(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x48(%rbp),%ecx
mov    -0x4c(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    130e <matmul+0x20e>
jmp    136f <matmul+0x26f>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    13ae <matmul+0x2ae>
jp     13ae <matmul+0x2ae>
movsd  -0x20(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
mov    0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x48(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    13f5 <matmul+0x2f5>
movsd  -0x20(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mov    0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x48(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x48(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    13fa <matmul+0x2fa>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1192 <matmul+0x92>
jmp    140d <matmul+0x30d>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    117f <matmul+0x7f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<matmul>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r9,-0x68(%rsp)
mov    %r8,-0x60(%rsp)
mov    %edx,-0x4c(%rsp)
xor    %r9d,%r9d
cmpb   $0x4e,(%rdi)
setne  %bl
xor    %eax,%eax
cmpb   $0x4e,0x1(%rdi)
sete   %dl
test   %esi,%esi
jle    1528 <matmul+0x428>
mov    %dl,%al
mov    %bl,%r9b
add    %r9d,%r9d
sub    %eax,%r9d
add    $0x2,%r9d
movslq -0x4c(%rsp),%r8
mov    %esi,%esi
movslq %ecx,%rdx
mov    %r8d,%eax
mov    %rax,-0x8(%rsp)
mov    %edx,%ebx
lea    -0x1(%rbx),%rax
mov    %rax,-0x48(%rsp)
mov    %rbx,-0x40(%rsp)
and    $0xfffffffe,%ebx
mov    -0x60(%rsp),%rdi
lea    0x8(%rdi),%rax
mov    %rdx,-0x38(%rsp)
lea    0x0(,%rdx,8),%rdx
mov    %rdx,-0x18(%rsp)
mov    %r8,%r11
shl    $0x4,%r11
mov    -0x68(%rsp),%rdx
add    $0x8,%rdx
mov    %rdx,-0x10(%rsp)
mov    %rsi,%r15
shl    $0x4,%r15
xor    %edx,%edx
mov    %rdx,-0x70(%rsp)
add    $0xffffffff,%r9d
xorpd  %xmm2,%xmm2
mov    %rdi,-0x30(%rsp)
mov    %r9,-0x58(%rsp)
mov    %ecx,-0x24(%rsp)
jmp    11e5 <matmul+0xe5>
nopl   0x0(%rax)
mov    -0x70(%rsp),%rdx
add    $0x1,%rdx
add    -0x18(%rsp),%rax
addq   $0x8,-0x30(%rsp)
mov    %rdx,%rdi
mov    %rdx,-0x70(%rsp)
cmp    %rsi,%rdx
je     1528 <matmul+0x428>
cmpl   $0x0,-0x4c(%rsp)
jle    11c0 <matmul+0xc0>
mov    -0x70(%rsp),%rdx
imul   -0x38(%rsp),%rdx
mov    %rdx,-0x20(%rsp)
mov    -0x10(%rsp),%r13
mov    -0x68(%rsp),%r14
xor    %r10d,%r10d
jmp    124f <matmul+0x14f>
nopl   0x0(%rax,%rax,1)
mov    %r10,%rdx
imul   %rsi,%rdx
add    -0x70(%rsp),%rdx
mov    0x38(%rsp),%rdi
movsd  (%rdi,%rdx,8),%xmm4
mulsd  %xmm1,%xmm4
mulsd  %xmm0,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,(%rdi,%rdx,8)
add    $0x1,%r10
add    $0x8,%r14
add    -0x18(%rsp),%r13
cmp    -0x8(%rsp),%r10
je     11c0 <matmul+0xc0>
xorpd  %xmm3,%xmm3
cmp    $0x3,%r9d
ja     14fe <matmul+0x3fe>
lea    0xd9c(%rip),%rdi
movslq (%rdi,%r9,4),%rdx
add    %rdi,%rdx
jmp    *%rdx
test   %ecx,%ecx
jle    14fe <matmul+0x3fe>
cmpq   $0x0,-0x48(%rsp)
je     13f4 <matmul+0x2f4>
xorpd  %xmm3,%xmm3
mov    -0x30(%rsp),%rbp
xor    %edx,%edx
nopl   0x0(%rax)
movsd  0x0(%rbp),%xmm4
mulsd  -0x8(%r13,%rdx,8),%xmm4
addsd  %xmm3,%xmm4
movsd  0x0(%rbp,%rsi,8),%xmm3
mulsd  0x0(%r13,%rdx,8),%xmm3
addsd  %xmm4,%xmm3
add    $0x2,%rdx
add    %r15,%rbp
cmp    %rdx,%rbx
jne    1290 <matmul+0x190>
jmp    13fa <matmul+0x2fa>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
test   %ecx,%ecx
jle    14fe <matmul+0x3fe>
cmpq   $0x0,-0x48(%rsp)
je     1442 <matmul+0x342>
xorpd  %xmm3,%xmm3
mov    %r14,%rcx
mov    -0x30(%rsp),%r12
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%r12),%xmm4
mulsd  (%rcx),%xmm4
addsd  %xmm3,%xmm4
movsd  (%r12,%rsi,8),%xmm3
mulsd  (%rcx,%r8,8),%xmm3
addsd  %xmm4,%xmm3
add    $0x2,%rdi
add    %r15,%r12
add    %r11,%rcx
cmp    %rdi,%rbx
jne    1300 <matmul+0x200>
jmp    1448 <matmul+0x348>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
test   %ecx,%ecx
jle    14fe <matmul+0x3fe>
cmpq   $0x0,-0x48(%rsp)
je     1486 <matmul+0x386>
xorpd  %xmm4,%xmm4
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movsd  -0x8(%rax,%rdx,8),%xmm5
movsd  (%rax,%rdx,8),%xmm3
mulsd  -0x8(%r13,%rdx,8),%xmm5
mulsd  0x0(%r13,%rdx,8),%xmm3
addsd  %xmm4,%xmm5
addsd  %xmm5,%xmm3
add    $0x2,%rdx
movapd %xmm3,%xmm4
cmp    %rdx,%rbx
jne    1360 <matmul+0x260>
testb  $0x1,-0x40(%rsp)
jne    1493 <matmul+0x393>
jmp    14fe <matmul+0x3fe>
xchg   %ax,%ax
test   %ecx,%ecx
jle    14fe <matmul+0x3fe>
cmpq   $0x0,-0x48(%rsp)
je     14c4 <matmul+0x3c4>
xorpd  %xmm4,%xmm4
mov    %r14,%rbp
xor    %edx,%edx
nopl   (%rax)
movsd  -0x8(%rax,%rdx,8),%xmm5
mulsd  0x0(%rbp),%xmm5
movsd  (%rax,%rdx,8),%xmm3
addsd  %xmm4,%xmm5
mulsd  0x0(%rbp,%r8,8),%xmm3
addsd  %xmm5,%xmm3
add    $0x2,%rdx
add    %r11,%rbp
movapd %xmm3,%xmm4
cmp    %rdx,%rbx
jne    13c0 <matmul+0x2c0>
jmp    14ca <matmul+0x3ca>
xorpd  %xmm3,%xmm3
xor    %edx,%edx
testb  $0x1,-0x40(%rsp)
mov    -0x58(%rsp),%r9
je     14fe <matmul+0x3fe>
mov    %r10,%r12
imul   -0x38(%rsp),%r12
mov    %rdx,%rdi
imul   %rsi,%rdi
add    -0x70(%rsp),%rdi
mov    -0x60(%rsp),%rbp
movsd  0x0(%rbp,%rdi,8),%xmm4
add    %rdx,%r12
mov    -0x68(%rsp),%rdx
mulsd  (%rdx,%r12,8),%xmm4
mov    -0x58(%rsp),%r9
jmp    14fa <matmul+0x3fa>
xorpd  %xmm3,%xmm3
xor    %edi,%edi
testb  $0x1,-0x40(%rsp)
mov    -0x24(%rsp),%ecx
mov    -0x58(%rsp),%r9
je     14fe <matmul+0x3fe>
mov    %rdi,%rdx
imul   %rsi,%rdx
add    -0x70(%rsp),%rdx
mov    -0x60(%rsp),%rbp
movsd  0x0(%rbp,%rdx,8),%xmm4
imul   %r8,%rdi
add    %r10,%rdi
mov    -0x68(%rsp),%rdx
mulsd  (%rdx,%rdi,8),%xmm4
jmp    14fa <matmul+0x3fa>
xorpd  %xmm3,%xmm3
xor    %edx,%edx
testb  $0x1,-0x40(%rsp)
je     14fe <matmul+0x3fe>
mov    %r10,%r12
imul   -0x38(%rsp),%r12
mov    -0x20(%rsp),%rdi
lea    (%rdx,%rdi,1),%rbp
mov    -0x60(%rsp),%rdi
movsd  (%rdi,%rbp,8),%xmm4
mov    -0x58(%rsp),%r9
add    %rdx,%r12
mov    -0x68(%rsp),%rdx
mulsd  (%rdx,%r12,8),%xmm4
jmp    14fa <matmul+0x3fa>
xorpd  %xmm3,%xmm3
xor    %edx,%edx
testb  $0x1,-0x40(%rsp)
mov    -0x58(%rsp),%r9
je     14fe <matmul+0x3fe>
mov    -0x20(%rsp),%rdi
add    %rdx,%rdi
mov    -0x60(%rsp),%r12
movsd  (%r12,%rdi,8),%xmm4
mov    -0x68(%rsp),%rdi
imul   %r8,%rdx
add    %r10,%rdx
mulsd  (%rdi,%rdx,8),%xmm4
addsd  %xmm4,%xmm3
ucomisd %xmm2,%xmm1
jne    1210 <matmul+0x110>
jp     1210 <matmul+0x110>
mulsd  %xmm0,%xmm3
mov    %r10,%rdx
imul   %rsi,%rdx
add    -0x70(%rsp),%rdx
mov    0x38(%rsp),%rdi
jmp    1232 <matmul+0x132>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @matmul(i8* nocapture %0, i32 %1, i32 %2, i32 %3, double %4, double* nocapture %5, double* nocapture %6, double %7, double* nocapture %8) {
%10 = load i8, i8* %0
%11 = icmp eq i8 %10, 78
%12 = getelementptr inbounds i8, i8* %0, i64 1
%13 = load i8, i8* %12
%14 = icmp eq i8 %13, 78
%15 = select i1 %14, i32 1, i32 2
%16 = select i1 %14, i32 3, i32 4
%17 = select i1 %11, i32 %15, i32 %16
%18 = icmp sgt i32 %1, 0
br i1 %18, label %19, label %230
19: ; preds = %9
%20 = icmp sgt i32 %2, 0
%21 = icmp sgt i32 %3, 0
%22 = fcmp oeq double %7, 0.000000e+00
%23 = sext i32 %2 to i64
%24 = zext i32 %1 to i64
%25 = sext i32 %3 to i64
%26 = zext i32 %1 to i64
%27 = zext i32 %2 to i64
%28 = zext i32 %3 to i64
%29 = add nsw i64 %28, -1
%30 = and i64 %28, 1
%31 = icmp eq i64 %29, 0
%32 = and i64 %28, 4294967294
%33 = icmp eq i64 %30, 0
%34 = and i64 %28, 1
%35 = icmp eq i64 %29, 0
%36 = and i64 %28, 4294967294
%37 = icmp eq i64 %34, 0
%38 = and i64 %28, 1
%39 = icmp eq i64 %29, 0
%40 = and i64 %28, 4294967294
%41 = icmp eq i64 %38, 0
%42 = and i64 %28, 1
%43 = icmp eq i64 %29, 0
%44 = and i64 %28, 4294967294
%45 = icmp eq i64 %42, 0
br label %46
46: ; preds = %19, %227
%47 = phi i64 [ 0, %19 ], [ %228, %227 ]
br i1 %20, label %48, label %227
48: ; preds = %46
%49 = mul nsw i64 %47, %25
br label %50
50: ; preds = %48, %224
%51 = phi i64 [ 0, %48 ], [ %225, %224 ]
switch i32 %17, label %210 [
i32 1, label %59
i32 2, label %57
i32 3, label %54
i32 4, label %52
]
52: ; preds = %50
br i1 %21, label %53, label %210
53: ; preds = %52
br i1 %31, label %197, label %134
54: ; preds = %50
br i1 %21, label %55, label %210
55: ; preds = %54
%56 = mul nsw i64 %51, %25
br i1 %35, label %185, label %112
57: ; preds = %50
br i1 %21, label %58, label %210
58: ; preds = %57
br i1 %39, label %171, label %86
59: ; preds = %50
br i1 %21, label %60, label %210
60: ; preds = %59
%61 = mul nsw i64 %51, %25
br i1 %43, label %158, label %62
62: ; preds = %60, %62
%63 = phi i64 [ %83, %62 ], [ 0, %60 ]
%64 = phi double [ %82, %62 ], [ 0.000000e+00, %60 ]
%65 = phi i64 [ %84, %62 ], [ 0, %60 ]
%66 = mul nsw i64 %63, %24
%67 = add nuw nsw i64 %66, %47
%68 = getelementptr inbounds double, double* %5, i64 %67
%69 = load double, double* %68
%70 = add nsw i64 %63, %61
%71 = getelementptr inbounds double, double* %6, i64 %70
%72 = load double, double* %71
%73 = tail call double @llvm.fmuladd.f64(double %69, double %72, double %64)
%74 = or i64 %63, 1
%75 = mul nsw i64 %74, %24
%76 = add nuw nsw i64 %75, %47
%77 = getelementptr inbounds double, double* %5, i64 %76
%78 = load double, double* %77
%79 = add nsw i64 %74, %61
%80 = getelementptr inbounds double, double* %6, i64 %79
%81 = load double, double* %80
%82 = tail call double @llvm.fmuladd.f64(double %78, double %81, double %73)
%83 = add nuw nsw i64 %63, 2
%84 = add i64 %65, 2
%85 = icmp eq i64 %84, %44
br i1 %85, label %158, label %62
86: ; preds = %58, %86
%87 = phi i64 [ %109, %86 ], [ 0, %58 ]
%88 = phi double [ %108, %86 ], [ 0.000000e+00, %58 ]
%89 = phi i64 [ %110, %86 ], [ 0, %58 ]
%90 = mul nsw i64 %87, %24
%91 = add nuw nsw i64 %90, %47
%92 = getelementptr inbounds double, double* %5, i64 %91
%93 = load double, double* %92
%94 = mul nsw i64 %87, %23
%95 = add nsw i64 %94, %51
%96 = getelementptr inbounds double, double* %6, i64 %95
%97 = load double, double* %96
%98 = tail call double @llvm.fmuladd.f64(double %93, double %97, double %88)
%99 = or i64 %87, 1
%100 = mul nsw i64 %99, %24
%101 = add nuw nsw i64 %100, %47
%102 = getelementptr inbounds double, double* %5, i64 %101
%103 = load double, double* %102
%104 = mul nsw i64 %99, %23
%105 = add nsw i64 %104, %51
%106 = getelementptr inbounds double, double* %6, i64 %105
%107 = load double, double* %106
%108 = tail call double @llvm.fmuladd.f64(double %103, double %107, double %98)
%109 = add nuw nsw i64 %87, 2
%110 = add i64 %89, 2
%111 = icmp eq i64 %110, %40
br i1 %111, label %171, label %86
112: ; preds = %55, %112
%113 = phi i64 [ %131, %112 ], [ 0, %55 ]
%114 = phi double [ %130, %112 ], [ 0.000000e+00, %55 ]
%115 = phi i64 [ %132, %112 ], [ 0, %55 ]
%116 = add nsw i64 %113, %49
%117 = getelementptr inbounds double, double* %5, i64 %116
%118 = load double, double* %117
%119 = add nsw i64 %113, %56
%120 = getelementptr inbounds double, double* %6, i64 %119
%121 = load double, double* %120
%122 = tail call double @llvm.fmuladd.f64(double %118, double %121, double %114)
%123 = or i64 %113, 1
%124 = add nsw i64 %123, %49
%125 = getelementptr inbounds double, double* %5, i64 %124
%126 = load double, double* %125
%127 = add nsw i64 %123, %56
%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)
%131 = add nuw nsw i64 %113, 2
%132 = add i64 %115, 2
%133 = icmp eq i64 %132, %36
br i1 %133, label %185, label %112
134: ; preds = %53, %134
%135 = phi i64 [ %155, %134 ], [ 0, %53 ]
%136 = phi double [ %154, %134 ], [ 0.000000e+00, %53 ]
%137 = phi i64 [ %156, %134 ], [ 0, %53 ]
%138 = add nsw i64 %135, %49
%139 = getelementptr inbounds double, double* %5, i64 %138
%140 = load double, double* %139
%141 = mul nsw i64 %135, %23
%142 = add nsw i64 %141, %51
%143 = getelementptr inbounds double, double* %6, i64 %142
%144 = load double, double* %143
%145 = tail call double @llvm.fmuladd.f64(double %140, double %144, double %136)
%146 = or i64 %135, 1
%147 = add nsw i64 %146, %49
%148 = getelementptr inbounds double, double* %5, i64 %147
%149 = load double, double* %148
%150 = mul nsw i64 %146, %23
%151 = add nsw i64 %150, %51
%152 = getelementptr inbounds double, double* %6, i64 %151
%153 = load double, double* %152
%154 = tail call double @llvm.fmuladd.f64(double %149, double %153, double %145)
%155 = add nuw nsw i64 %135, 2
%156 = add i64 %137, 2
%157 = icmp eq i64 %156, %32
br i1 %157, label %197, label %134
158: ; preds = %62, %60
%159 = phi double [ undef, %60 ], [ %82, %62 ]
%160 = phi i64 [ 0, %60 ], [ %83, %62 ]
%161 = phi double [ 0.000000e+00, %60 ], [ %82, %62 ]
br i1 %45, label %210, label %162
162: ; preds = %158
%163 = mul nsw i64 %160, %24
%164 = add nuw nsw i64 %163, %47
%165 = getelementptr inbounds double, double* %5, i64 %164
%166 = load double, double* %165
%167 = add nsw i64 %160, %61
%168 = getelementptr inbounds double, double* %6, i64 %167
%169 = load double, double* %168
%170 = tail call double @llvm.fmuladd.f64(double %166, double %169, double %161)
br label %210
171: ; preds = %86, %58
%172 = phi double [ undef, %58 ], [ %108, %86 ]
%173 = phi i64 [ 0, %58 ], [ %109, %86 ]
%174 = phi double [ 0.000000e+00, %58 ], [ %108, %86 ]
br i1 %41, label %210, label %175
175: ; preds = %171
%176 = mul nsw i64 %173, %24
%177 = add nuw nsw i64 %176, %47
%178 = getelementptr inbounds double, double* %5, i64 %177
%179 = load double, double* %178
%180 = mul nsw i64 %173, %23
%181 = add nsw i64 %180, %51
%182 = getelementptr inbounds double, double* %6, i64 %181
%183 = load double, double* %182
%184 = tail call double @llvm.fmuladd.f64(double %179, double %183, double %174)
br label %210
185: ; preds = %112, %55
%186 = phi double [ undef, %55 ], [ %130, %112 ]
%187 = phi i64 [ 0, %55 ], [ %131, %112 ]
%188 = phi double [ 0.000000e+00, %55 ], [ %130, %112 ]
br i1 %37, label %210, label %189
189: ; preds = %185
%190 = add nsw i64 %187, %49
%191 = getelementptr inbounds double, double* %5, i64 %190
%192 = load double, double* %191
%193 = add nsw i64 %187, %56
%194 = getelementptr inbounds double, double* %6, i64 %193
%195 = load double, double* %194
%196 = tail call double @llvm.fmuladd.f64(double %192, double %195, double %188)
br label %210
197: ; preds = %134, %53
%198 = phi double [ undef, %53 ], [ %154, %134 ]
%199 = phi i64 [ 0, %53 ], [ %155, %134 ]
%200 = phi double [ 0.000000e+00, %53 ], [ %154, %134 ]
br i1 %33, label %210, label %201
201: ; preds = %197
%202 = add nsw i64 %199, %49
%203 = getelementptr inbounds double, double* %5, i64 %202
%204 = load double, double* %203
%205 = mul nsw i64 %199, %23
%206 = add nsw i64 %205, %51
%207 = getelementptr inbounds double, double* %6, i64 %206
%208 = load double, double* %207
%209 = tail call double @llvm.fmuladd.f64(double %204, double %208, double %200)
br label %210
210: ; preds = %201, %197, %189, %185, %175, %171, %162, %158, %52, %54, %57, %59, %50
%211 = phi double [ 0.000000e+00, %50 ], [ 0.000000e+00, %59 ], [ 0.000000e+00, %57 ], [ 0.000000e+00, %54 ], [ 0.000000e+00, %52 ], [ %159, %158 ], [ %170, %162 ], [ %172, %171 ], [ %184, %175 ], [ %186, %185 ], [ %196, %189 ], [ %198, %197 ], [ %209, %201 ]
br i1 %22, label %212, label %217
212: ; preds = %210
%213 = fmul double %211, %4
%214 = mul nsw i64 %51, %24
%215 = add nuw nsw i64 %214, %47
%216 = getelementptr inbounds double, double* %8, i64 %215
store double %213, double* %216
br label %224
217: ; preds = %210
%218 = mul nsw i64 %51, %24
%219 = add nuw nsw i64 %218, %47
%220 = getelementptr inbounds double, double* %8, i64 %219
%221 = load double, double* %220
%222 = fmul double %221, %7
%223 = tail call double @llvm.fmuladd.f64(double %4, double %211, double %222)
store double %223, double* %220
br label %224
224: ; preds = %212, %217
%225 = add nuw nsw i64 %51, 1
%226 = icmp eq i64 %225, %27
br i1 %226, label %227, label %50
227: ; preds = %224, %46
%228 = add nuw nsw i64 %47, 1
%229 = icmp eq i64 %228, %26
br i1 %229, label %230, label %46
230: ; preds = %227, %9
ret void
}