void rec_matmul(float *A, float *B, float *C, int m, int n, int p, int ld,
                int add) {
  if ((m + n + p) <= 64) {
    int i, j, k;
    if (add) {
      for (i = 0; i < m; i++)
        for (k = 0; k < p; k++) {
          float c = 0.0;
          for (j = 0; j < n; j++)
            c += A[i * ld + j] * B[j * ld + k];
          C[i * ld + k] += c;
        }
    } else {
      for (i = 0; i < m; i++)
        for (k = 0; k < p; k++) {
          float c = 0.0;
          for (j = 0; j < n; j++)
            c += A[i * ld + j] * B[j * ld + k];
          C[i * ld + k] = c;
        }
    }
  } else if (m >= n && n >= p) {
    int m1 = m >> 1;
    rec_matmul(A, B, C, m1, n, p, ld, add);
    rec_matmul(A + m1 * ld, B, C + m1 * ld, m - m1, n, p, ld, add);
  } else if (n >= m && n >= p) {
    int n1 = n >> 1;
    rec_matmul(A, B, C, m, n1, p, ld, add);
    rec_matmul(A + n1, B + n1 * ld, C, m, n - n1, p, ld, 1);
  } else {
    int p1 = p >> 1;
    rec_matmul(A, B, C, m, n, p1, ld, add);
    rec_matmul(A, B + p1, C + p1, m, n, p - p1, ld, add);
  }
}


x86 assembly (clang, -O0):
<rec_matmul>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
mov    -0x1c(%rbp),%eax
add    -0x20(%rbp),%eax
add    -0x24(%rbp),%eax
cmp    $0x40,%eax
jg     12f7 <rec_matmul+0x1e7>
cmpl   $0x0,0x18(%rbp)
je     1224 <rec_matmul+0x114>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    121f <rec_matmul+0x10f>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    120c <rec_matmul+0xfc>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11de <rec_matmul+0xce>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   0x10(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x34(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1186 <rec_matmul+0x76>
movss  -0x34(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    116b <rec_matmul+0x5b>
jmp    1211 <rec_matmul+0x101>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1158 <rec_matmul+0x48>
jmp    12f2 <rec_matmul+0x1e2>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12ed <rec_matmul+0x1dd>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    12da <rec_matmul+0x1ca>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    12b1 <rec_matmul+0x1a1>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   0x10(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x38(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1259 <rec_matmul+0x149>
movss  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    123e <rec_matmul+0x12e>
jmp    12df <rec_matmul+0x1cf>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    122b <rec_matmul+0x11b>
jmp    12f2 <rec_matmul+0x1e2>
jmp    14b3 <rec_matmul+0x3a3>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jl     1396 <rec_matmul+0x286>
mov    -0x20(%rbp),%eax
cmp    -0x24(%rbp),%eax
jl     1396 <rec_matmul+0x286>
mov    -0x1c(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x3c(%rbp),%ecx
mov    -0x20(%rbp),%r8d
mov    -0x24(%rbp),%r9d
mov    0x10(%rbp),%r10d
mov    0x18(%rbp),%eax
mov    %r10d,(%rsp)
mov    %eax,0x8(%rsp)
call   1030 <rec_matmul@plt>
mov    -0x8(%rbp),%rdi
mov    -0x3c(%rbp),%eax
imul   0x10(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x3c(%rbp),%eax
imul   0x10(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rdx
mov    -0x1c(%rbp),%ecx
sub    -0x3c(%rbp),%ecx
mov    -0x20(%rbp),%r8d
mov    -0x24(%rbp),%r9d
mov    0x10(%rbp),%r10d
mov    0x18(%rbp),%eax
mov    %r10d,(%rsp)
mov    %eax,0x8(%rsp)
call   1030 <rec_matmul@plt>
jmp    14ae <rec_matmul+0x39e>
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     1430 <rec_matmul+0x320>
mov    -0x20(%rbp),%eax
cmp    -0x24(%rbp),%eax
jl     1430 <rec_matmul+0x320>
mov    -0x20(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
mov    -0x40(%rbp),%r8d
mov    -0x24(%rbp),%r9d
mov    0x10(%rbp),%r10d
mov    0x18(%rbp),%eax
mov    %r10d,(%rsp)
mov    %eax,0x8(%rsp)
call   1030 <rec_matmul@plt>
mov    -0x8(%rbp),%rdi
movslq -0x40(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x10(%rbp),%rsi
mov    -0x40(%rbp),%eax
imul   0x10(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%r8d
sub    -0x40(%rbp),%r8d
mov    -0x24(%rbp),%r9d
mov    0x10(%rbp),%eax
mov    %eax,(%rsp)
movl   $0x1,0x8(%rsp)
call   1030 <rec_matmul@plt>
jmp    14a9 <rec_matmul+0x399>
mov    -0x24(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%r8d
mov    -0x44(%rbp),%r9d
mov    0x10(%rbp),%r10d
mov    0x18(%rbp),%eax
mov    %r10d,(%rsp)
mov    %eax,0x8(%rsp)
call   1030 <rec_matmul@plt>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
movslq -0x44(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rsi
mov    -0x18(%rbp),%rdx
movslq -0x44(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%r8d
mov    -0x24(%rbp),%r9d
sub    -0x44(%rbp),%r9d
mov    0x10(%rbp),%r10d
mov    0x18(%rbp),%eax
mov    %r10d,(%rsp)
mov    %eax,0x8(%rsp)
call   1030 <rec_matmul@plt>
jmp    14ae <rec_matmul+0x39e>
jmp    14b3 <rec_matmul+0x3a3>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rec_matmul>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r9d,%r14d
mov    %r8d,%r15d
mov    %ecx,%r12d
mov    %rdx,%r13
mov    %rsi,0x18(%rsp)
mov    %rdi,0x10(%rsp)
mov    0x68(%rsp),%eax
mov    %rax,0x8(%rsp)
lea    (%r15,%r12,1),%eax
mov    %rax,(%rsp)
jmp    1198 <rec_matmul+0x88>
nopw   0x0(%rax,%rax,1)
mov    %r14d,%ebp
sar    %ebp
mov    0x10(%rsp),%rdi
mov    0x18(%rsp),%rbx
mov    %rbx,%rsi
mov    %r13,%rdx
mov    %r12d,%ecx
mov    %r15d,%r8d
mov    %ebp,%r9d
push   0x8(%rsp)
mov    0x68(%rsp),%eax
push   %rax
call   1030 <rec_matmul@plt>
add    $0x10,%rsp
movslq %ebp,%rax
lea    (%rbx,%rax,4),%rbx
mov    %rbx,0x18(%rsp)
lea    0x0(%r13,%rax,4),%r13
sub    %eax,%r14d
mov    (%rsp),%rax
add    %r14d,%eax
cmp    $0x41,%eax
jl     128f <rec_matmul+0x17f>
cmp    %r15d,%r12d
jl     11b2 <rec_matmul+0xa2>
cmp    %r14d,%r15d
jge    1240 <rec_matmul+0x130>
cmp    %r12d,%r15d
jl     1150 <rec_matmul+0x40>
cmp    %r14d,%r15d
jl     1150 <rec_matmul+0x40>
mov    %r15d,%r8d
sar    %r8d
mov    %r8d,(%rsp)
mov    0x10(%rsp),%rbx
mov    %rbx,%rdi
mov    0x18(%rsp),%rbp
mov    %rbp,%rsi
mov    %r13,%rdx
mov    %r12d,%ecx
mov    %r14d,%r9d
push   0x8(%rsp)
mov    0x68(%rsp),%eax
push   %rax
call   1030 <rec_matmul@plt>
add    $0x10,%rsp
movslq (%rsp),%rax
lea    (%rbx,%rax,4),%rbx
mov    %rbx,0x10(%rsp)
sub    %eax,%r15d
mov    0x60(%rsp),%ecx
imul   %ecx,%eax
cltq
lea    0x0(%rbp,%rax,4),%rbp
mov    %rbp,0x18(%rsp)
lea    (%r15,%r12,1),%eax
mov    %rax,(%rsp)
add    %r14d,%eax
mov    $0x1,%ecx
mov    %rcx,0x8(%rsp)
cmp    $0x41,%eax
jge    11a4 <rec_matmul+0x94>
jmp    14e8 <rec_matmul+0x3d8>
nopw   0x0(%rax,%rax,1)
mov    %r12d,%ebp
sar    %ebp
mov    0x10(%rsp),%rbx
mov    %rbx,%rdi
mov    0x18(%rsp),%rsi
mov    %r13,%rdx
mov    %ebp,%ecx
mov    %r15d,%r8d
mov    %r14d,%r9d
push   0x8(%rsp)
mov    0x68(%rsp),%eax
push   %rax
call   1030 <rec_matmul@plt>
add    $0x10,%rsp
sub    %ebp,%r12d
mov    0x60(%rsp),%eax
imul   %eax,%ebp
movslq %ebp,%rax
lea    (%rbx,%rax,4),%rbx
mov    %rbx,0x10(%rsp)
lea    0x0(%r13,%rax,4),%r13
jmp    113d <rec_matmul+0x2d>
cmpl   $0x0,0x8(%rsp)
mov    %r14,(%rsp)
je     13cf <rec_matmul+0x2bf>
test   %r12d,%r12d
jle    14f5 <rec_matmul+0x3e5>
mov    0x60(%rsp),%eax
movslq %eax,%rsi
mov    %r12d,%eax
mov    %rax,0x20(%rsp)
mov    (%rsp),%edx
mov    %r15d,%r11d
mov    %r11d,%edi
and    $0xfffffffe,%edi
mov    0x10(%rsp),%rax
lea    0x4(%rax),%rbp
lea    0x0(,%rsi,4),%rcx
mov    %rsi,0x8(%rsp)
lea    0x0(,%rsi,8),%rax
xor    %r10d,%r10d
mov    0x10(%rsp),%r9
jmp    1302 <rec_matmul+0x1f2>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r10
add    %rcx,%rbp
cmp    0x20(%rsp),%r10
je     14f5 <rec_matmul+0x3e5>
cmpl   $0x0,(%rsp)
jle    12f0 <rec_matmul+0x1e0>
mov    %r10,%rsi
imul   0x8(%rsp),%rsi
mov    0x18(%rsp),%r12
xor    %r8d,%r8d
jmp    1342 <rec_matmul+0x232>
nopl   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
lea    (%r8,%rsi,1),%rbx
addss  0x0(%r13,%rbx,4),%xmm0
movss  %xmm0,0x0(%r13,%rbx,4)
add    $0x1,%r8
add    $0x4,%r12
cmp    %rdx,%r8
je     12f0 <rec_matmul+0x1e0>
test   %r15d,%r15d
jle    1320 <rec_matmul+0x210>
cmp    $0x1,%r15d
jne    1360 <rec_matmul+0x250>
xorps  %xmm0,%xmm0
xor    %r14d,%r14d
jmp    139e <rec_matmul+0x28e>
cs nopw 0x0(%rax,%rax,1)
nop
xorps  %xmm1,%xmm1
mov    %r12,%rbx
xor    %r14d,%r14d
nopl   0x0(%rax)
movss  -0x4(%rbp,%r14,4),%xmm2
mulss  (%rbx),%xmm2
movss  0x0(%rbp,%r14,4),%xmm0
addss  %xmm1,%xmm2
mulss  (%rbx,%rcx,1),%xmm0
addss  %xmm2,%xmm0
add    $0x2,%r14
add    %rax,%rbx
movaps %xmm0,%xmm1
cmp    %r14,%rdi
jne    1370 <rec_matmul+0x260>
test   $0x1,%r11b
je     1323 <rec_matmul+0x213>
lea    (%r14,%rsi,1),%rbx
movss  (%r9,%rbx,4),%xmm1
imul   0x8(%rsp),%r14
add    %r8,%r14
mov    0x18(%rsp),%rbx
mulss  (%rbx,%r14,4),%xmm1
addss  %xmm1,%xmm0
jmp    1323 <rec_matmul+0x213>
test   %r12d,%r12d
jle    14f5 <rec_matmul+0x3e5>
mov    0x60(%rsp),%eax
movslq %eax,%rdx
mov    %r12d,%eax
mov    %rax,0x20(%rsp)
mov    (%rsp),%r12d
mov    %r15d,%r11d
mov    %r11d,%edi
and    $0xfffffffe,%edi
mov    0x10(%rsp),%rax
lea    0x4(%rax),%rbp
lea    0x0(,%rdx,4),%rcx
mov    %rdx,0x8(%rsp)
lea    0x0(,%rdx,8),%rax
xor    %r10d,%r10d
mov    0x10(%rsp),%r9
jmp    1432 <rec_matmul+0x322>
nopl   0x0(%rax)
add    $0x1,%r10
add    %rcx,%rbp
cmp    0x20(%rsp),%r10
je     14f5 <rec_matmul+0x3e5>
cmpl   $0x0,(%rsp)
jle    1420 <rec_matmul+0x310>
mov    %r10,%r14
imul   0x8(%rsp),%r14
mov    0x18(%rsp),%rsi
xor    %r8d,%r8d
jmp    146b <rec_matmul+0x35b>
nopl   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
lea    (%r8,%r14,1),%rdx
movss  %xmm0,0x0(%r13,%rdx,4)
add    $0x1,%r8
add    $0x4,%rsi
cmp    %r12,%r8
je     1420 <rec_matmul+0x310>
test   %r15d,%r15d
jle    1450 <rec_matmul+0x340>
cmp    $0x1,%r15d
jne    1480 <rec_matmul+0x370>
xorps  %xmm0,%xmm0
xor    %ebx,%ebx
jmp    14bc <rec_matmul+0x3ac>
nopl   (%rax)
xorps  %xmm1,%xmm1
mov    %rsi,%rdx
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movss  -0x4(%rbp,%rbx,4),%xmm2
mulss  (%rdx),%xmm2
movss  0x0(%rbp,%rbx,4),%xmm0
addss  %xmm1,%xmm2
mulss  (%rdx,%rcx,1),%xmm0
addss  %xmm2,%xmm0
add    $0x2,%rbx
add    %rax,%rdx
movaps %xmm0,%xmm1
cmp    %rbx,%rdi
jne    1490 <rec_matmul+0x380>
test   $0x1,%r11b
je     1453 <rec_matmul+0x343>
lea    (%rbx,%r14,1),%rdx
movss  (%r9,%rdx,4),%xmm1
imul   0x8(%rsp),%rbx
add    %r8,%rbx
mov    0x18(%rsp),%rdx
mulss  (%rdx,%rbx,4),%xmm1
addss  %xmm1,%xmm0
jmp    1453 <rec_matmul+0x343>
mov    %r14,(%rsp)
test   %r12d,%r12d
jg     12a7 <rec_matmul+0x197>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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