void dlevmar_trans_mat_mat_mult(double *a, double *b, int n, int m) {
  register int i, j, k, jj, kk;
  register double sum, *bim, *akm;
  const int bsize = 32;
  for (jj = 0; jj < m; jj += bsize) {
    for (i = 0; i < m; ++i) {
      bim = b + i * m;
      for (j = (((jj) >= (i)) ? (jj) : (i));
           j < (((jj + bsize) <= (m)) ? (jj + bsize) : (m)); ++j)
        bim[j] = 0.0;
    }
    for (kk = 0; kk < n; kk += bsize) {
      for (i = 0; i < m; ++i) {
        bim = b + i * m;
        for (j = (((jj) >= (i)) ? (jj) : (i));
             j < (((jj + bsize) <= (m)) ? (jj + bsize) : (m)); ++j) {
          sum = 0.0;
          for (k = kk; k < (((kk + bsize) <= (n)) ? (kk + bsize) : (n)); ++k) {
            akm = a + k * m;
            sum += akm[i] * akm[j];
          }
          bim[j] += sum;
        }
      }
    }
  }
  for (i = 0; i < m; ++i)
    for (j = 0; j < i; ++j)
      b[i * m + j] = b[j * m + i];
}


x86 assembly (clang, -O0):
<dlevmar_trans_mat_mat_mult>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x20,-0x4c(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1372 <dlevmar_trans_mat_mat_mult+0x272>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11e3 <dlevmar_trans_mat_mat_mult+0xe3>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     116f <dlevmar_trans_mat_mat_mult+0x6f>
mov    -0x28(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    1175 <dlevmar_trans_mat_mat_mult+0x75>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x28(%rbp),%eax
add    $0x20,%eax
cmp    -0x18(%rbp),%eax
jg     119e <dlevmar_trans_mat_mat_mult+0x9e>
mov    -0x28(%rbp),%eax
add    $0x20,%eax
mov    %eax,-0x58(%rbp)
jmp    11a4 <dlevmar_trans_mat_mat_mult+0xa4>
mov    -0x18(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x58(%rbp),%ecx
cmp    %ecx,%eax
jge    11d0 <dlevmar_trans_mat_mat_mult+0xd0>
mov    -0x40(%rbp),%rax
movslq -0x20(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    117b <dlevmar_trans_mat_mat_mult+0x7b>
jmp    11d5 <dlevmar_trans_mat_mat_mult+0xd5>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1133 <dlevmar_trans_mat_mat_mult+0x33>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    135f <dlevmar_trans_mat_mat_mult+0x25f>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    134c <dlevmar_trans_mat_mat_mult+0x24c>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     1239 <dlevmar_trans_mat_mat_mult+0x139>
mov    -0x28(%rbp),%eax
mov    %eax,-0x5c(%rbp)
jmp    123f <dlevmar_trans_mat_mat_mult+0x13f>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x28(%rbp),%eax
add    $0x20,%eax
cmp    -0x18(%rbp),%eax
jg     1268 <dlevmar_trans_mat_mat_mult+0x168>
mov    -0x28(%rbp),%eax
add    $0x20,%eax
mov    %eax,-0x64(%rbp)
jmp    126e <dlevmar_trans_mat_mat_mult+0x16e>
mov    -0x18(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x60(%rbp),%eax
mov    -0x64(%rbp),%ecx
cmp    %ecx,%eax
jge    1339 <dlevmar_trans_mat_mat_mult+0x239>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x20,%eax
cmp    -0x14(%rbp),%eax
jg     12ad <dlevmar_trans_mat_mat_mult+0x1ad>
mov    -0x2c(%rbp),%eax
add    $0x20,%eax
mov    %eax,-0x6c(%rbp)
jmp    12b3 <dlevmar_trans_mat_mat_mult+0x1b3>
mov    -0x14(%rbp),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x68(%rbp),%eax
mov    -0x6c(%rbp),%ecx
cmp    %ecx,%eax
jge    1314 <dlevmar_trans_mat_mat_mult+0x214>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    128a <dlevmar_trans_mat_mat_mult+0x18a>
movsd  -0x38(%rbp),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x20(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1245 <dlevmar_trans_mat_mat_mult+0x145>
jmp    133e <dlevmar_trans_mat_mat_mult+0x23e>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    11fd <dlevmar_trans_mat_mat_mult+0xfd>
jmp    1351 <dlevmar_trans_mat_mat_mult+0x251>
mov    -0x2c(%rbp),%eax
add    $0x20,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ea <dlevmar_trans_mat_mat_mult+0xea>
jmp    1364 <dlevmar_trans_mat_mat_mult+0x264>
mov    -0x28(%rbp),%eax
add    $0x20,%eax
mov    %eax,-0x28(%rbp)
jmp    1120 <dlevmar_trans_mat_mat_mult+0x20>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    13e5 <dlevmar_trans_mat_mat_mult+0x2e5>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13d2 <dlevmar_trans_mat_mat_mult+0x2d2>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    138c <dlevmar_trans_mat_mat_mult+0x28c>
jmp    13d7 <dlevmar_trans_mat_mat_mult+0x2d7>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1379 <dlevmar_trans_mat_mat_mult+0x279>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dlevmar_trans_mat_mat_mult>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xb8,%rsp
mov    %rsi,0x8(%rsp)
mov    %rdi,0x60(%rsp)
test   %ecx,%ecx
jle    1587 <dlevmar_trans_mat_mat_mult+0x477>
mov    %ecx,%r12d
mov    %ecx,%ebx
mov    %rbx,%rax
shl    $0x8,%rax
mov    %rax,0x70(%rsp)
lea    0x0(,%rbx,8),%rax
mov    %rax,0x58(%rsp)
mov    %rbx,%rax
shl    $0x5,%rax
mov    %rax,0x88(%rsp)
xor    %eax,%eax
mov    %rbx,0x48(%rsp)
mov    %ecx,0x28(%rsp)
mov    %edx,0x4(%rsp)
jmp    1193 <dlevmar_trans_mat_mat_mult+0x83>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x38(%rsp),%rax
cmp    %rbx,%rax
mov    0x28(%rsp),%r12d
jae    14e1 <dlevmar_trans_mat_mat_mult+0x3d1>
mov    %rax,0x18(%rsp)
add    $0x20,%rax
cmp    %eax,%r12d
mov    %eax,%r14d
cmovl  %r12d,%r14d
cmp    %rbx,%rax
mov    %rax,0x38(%rsp)
mov    %eax,%r15d
cmova  %r12d,%r15d
xor    %r13d,%r13d
xor    %ebp,%ebp
jmp    11cc <dlevmar_trans_mat_mat_mult+0xbc>
nopl   0x0(%rax)
add    $0x1,%rbp
add    %r12d,%r13d
cmp    %rbp,%rbx
je     1210 <dlevmar_trans_mat_mat_mult+0x100>
mov    0x18(%rsp),%rcx
cmp    %ebp,%ecx
mov    %ebp,%eax
cmova  %ecx,%eax
cmp    %rbp,%rcx
cmovbe %ebp,%ecx
cmp    %ecx,%r15d
jle    11c0 <dlevmar_trans_mat_mat_mult+0xb0>
mov    %eax,%ecx
not    %eax
add    %r14d,%eax
lea    0x8(,%rax,8),%rdx
movslq %r13d,%rax
add    %rax,%rcx
mov    0x8(%rsp),%rax
lea    (%rax,%rcx,8),%rdi
xor    %esi,%esi
call   1030 <memset@plt>
jmp    11c0 <dlevmar_trans_mat_mat_mult+0xb0>
nopw   0x0(%rax,%rax,1)
mov    0x4(%rsp),%edx
test   %edx,%edx
jle    1180 <dlevmar_trans_mat_mat_mult+0x70>
mov    0x38(%rsp),%rax
cmp    %rbx,%rax
cmova  %ebx,%eax
mov    %eax,0x30(%rsp)
cltq
mov    %rax,0xa8(%rsp)
mov    $0x20,%eax
xor    %r10d,%r10d
xor    %ecx,%ecx
xor    %esi,%esi
jmp    127a <dlevmar_trans_mat_mat_mult+0x16a>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x20,%rsi
mov    0x78(%rsp),%rcx
add    $0x1,%rcx
mov    0x40(%rsp),%r10
add    0x70(%rsp),%r10
mov    0x2c(%rsp),%eax
add    $0x20,%eax
mov    0x4(%rsp),%edx
cmp    %edx,%esi
jge    1180 <dlevmar_trans_mat_mat_mult+0x70>
cmp    %eax,%edx
mov    %eax,0x2c(%rsp)
cmovl  %edx,%eax
and    $0x3,%eax
neg    %rax
mov    %rax,0x80(%rsp)
mov    %rcx,0x78(%rsp)
shl    $0x5,%rcx
lea    0x20(%rsi),%edi
cmp    %edi,%edx
mov    %edi,%eax
cmovl  %edx,%eax
movslq %eax,%rbp
not    %rcx
mov    %rbp,%rax
mov    %rbp,0x98(%rsp)
add    %rbp,%rcx
mov    %rcx,0xa0(%rsp)
cmp    %edx,%edi
cmovg  %edx,%edi
mov    %edi,0x34(%rsp)
movslq %edi,%rax
mov    %rax,0x90(%rsp)
xor    %eax,%eax
mov    %rax,0x10(%rsp)
mov    %r10,0x40(%rsp)
xor    %ecx,%ecx
mov    %rsi,0x20(%rsp)
jmp    1311 <dlevmar_trans_mat_mat_mult+0x201>
nopl   0x0(%rax)
mov    0x50(%rsp),%rcx
add    $0x1,%rcx
add    $0x8,%r10
addq   $0x8,0x10(%rsp)
mov    0x48(%rsp),%rbx
cmp    %rbx,%rcx
je     1250 <dlevmar_trans_mat_mat_mult+0x140>
mov    0x18(%rsp),%rax
cmp    %ecx,%eax
mov    %ecx,%edx
cmova  %eax,%edx
cmp    %rcx,%rax
mov    %rcx,0x50(%rsp)
cmovbe %ecx,%eax
cmp    %eax,0x30(%rsp)
jle    12f0 <dlevmar_trans_mat_mat_mult+0x1e0>
mov    0x40(%rsp),%rax
lea    (%rax,%rdx,8),%r15
shl    $0x3,%rdx
mov    %rdx,0x68(%rsp)
mov    0x50(%rsp),%r11
mov    %r11,%rax
imul   0x48(%rsp),%rax
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
mov    %rax,0xb0(%rsp)
mov    0x18(%rsp),%rax
cmp    %r11d,%eax
cmova  %eax,%r11d
jmp    13a5 <dlevmar_trans_mat_mat_mult+0x295>
nopl   (%rax)
mov    0xb0(%rsp),%rax
addsd  (%rax,%r11,8),%xmm0
movsd  %xmm0,(%rax,%r11,8)
add    $0x1,%r11
add    $0x8,%r15
addq   $0x8,0x68(%rsp)
cmp    0xa8(%rsp),%r11
mov    0x20(%rsp),%rsi
jge    12f0 <dlevmar_trans_mat_mat_mult+0x1e0>
xorpd  %xmm0,%xmm0
cmp    %esi,0x34(%rsp)
jle    1370 <dlevmar_trans_mat_mat_mult+0x260>
testb  $0x3,0x98(%rsp)
03
je     1406 <dlevmar_trans_mat_mat_mult+0x2f6>
xorpd  %xmm0,%xmm0
mov    0x60(%rsp),%rcx
xor    %eax,%eax
mov    0x58(%rsp),%rsi
mov    0x80(%rsp),%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%r10,1),%xmm1
mulsd  (%rcx,%r15,1),%xmm1
addsd  %xmm1,%xmm0
add    $0xffffffffffffffff,%rax
add    %rsi,%rcx
cmp    %rax,%rdx
jne    13e0 <dlevmar_trans_mat_mat_mult+0x2d0>
mov    0x20(%rsp),%rdi
sub    %rax,%rdi
jmp    140f <dlevmar_trans_mat_mat_mult+0x2ff>
xorpd  %xmm0,%xmm0
mov    0x20(%rsp),%rdi
cmpq   $0x3,0xa0(%rsp)
mov    0x88(%rsp),%rdx
mov    0x90(%rsp),%r14
jb     1370 <dlevmar_trans_mat_mat_mult+0x260>
lea    0x3(%rdi),%rcx
mov    0x58(%rsp),%r9
imul   %r9,%rcx
mov    0x68(%rsp),%r12
lea    (%r12,%rcx,1),%rsi
mov    0x10(%rsp),%rbx
mov    %rbx,0x10(%rsp)
add    %rbx,%rcx
lea    0x2(%rdi),%rax
imul   %r9,%rax
lea    (%r12,%rax,1),%r8
add    %rbx,%rax
mov    %r9,%r13
imul   %rdi,%r13
lea    (%rbx,%r13,1),%rbp
add    %r12,%r13
lea    0x1(%rdi),%rbx
imul   %r9,%rbx
add    %rbx,%r12
add    0x10(%rsp),%rbx
mov    0x60(%rsp),%r9
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%r9,%rbp,1),%xmm1
mulsd  (%r9,%r13,1),%xmm1
addsd  %xmm0,%xmm1
movsd  (%r9,%rbx,1),%xmm0
mulsd  (%r9,%r12,1),%xmm0
addsd  %xmm1,%xmm0
movsd  (%r9,%rax,1),%xmm1
mulsd  (%r9,%r8,1),%xmm1
addsd  %xmm0,%xmm1
movsd  (%r9,%rcx,1),%xmm0
mulsd  (%r9,%rsi,1),%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rdi
add    %rdx,%r9
cmp    %r14,%rdi
jl     1490 <dlevmar_trans_mat_mat_mult+0x380>
jmp    1370 <dlevmar_trans_mat_mat_mult+0x260>
test   %r12d,%r12d
jle    1587 <dlevmar_trans_mat_mat_mult+0x477>
lea    0x0(,%rbx,8),%rax
mov    %rbx,%r10
shl    $0x4,%r10
mov    0x8(%rsp),%r8
lea    0x8(%r8),%rdx
xor    %ecx,%ecx
mov    0x8(%rsp),%r9
jmp    1520 <dlevmar_trans_mat_mat_mult+0x410>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rcx
add    $0x8,%r8
add    %rax,%rdx
cmp    %rbx,%rcx
je     1587 <dlevmar_trans_mat_mat_mult+0x477>
test   %rcx,%rcx
je     1510 <dlevmar_trans_mat_mat_mult+0x400>
cmp    $0x1,%rcx
jne    1530 <dlevmar_trans_mat_mat_mult+0x420>
xor    %ebp,%ebp
jmp    1560 <dlevmar_trans_mat_mat_mult+0x450>
nop
mov    %rcx,%rsi
and    $0xfffffffffffffffe,%rsi
mov    %r8,%rdi
xor    %ebp,%ebp
nopl   0x0(%rax)
movsd  (%rdi),%xmm0
movsd  %xmm0,-0x8(%rdx,%rbp,8)
movsd  (%rdi,%rax,1),%xmm0
movsd  %xmm0,(%rdx,%rbp,8)
add    $0x2,%rbp
add    %r10,%rdi
cmp    %rbp,%rsi
jne    1540 <dlevmar_trans_mat_mat_mult+0x430>
test   $0x1,%cl
je     1510 <dlevmar_trans_mat_mat_mult+0x400>
mov    %rcx,%rsi
imul   %rbx,%rsi
mov    %rbp,%rdi
imul   %rbx,%rdi
add    %rcx,%rdi
movsd  (%r9,%rdi,8),%xmm0
add    %rbp,%rsi
movsd  %xmm0,(%r9,%rsi,8)
jmp    1510 <dlevmar_trans_mat_mat_mult+0x400>
add    $0xb8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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