int ludcmp(a, n, indx, d)
int n, *indx;
double **a, *d;
{
  int i, imax, j, k;
  double big, dum, sum, temp;
  double *vv, *vector();
  void free_vector();
  imax = 0;
  vv = vector(1, n);
  (*d) = 1.0;
  for (i = 1; i <= n; i++) {
    big = 0.0;
    for (j = 1; j <= n; j++)
      if ((temp = fabs(a[i][j])) > big)
        big = temp;
    if (!big)
      nrerror("Singular matrix in routine LUDCMP");
    vv[i] = 1.0 / big;
  }
  for (j = 1; j <= n; j++) {
    for (i = 1; i < j; i++) {
      sum = a[i][j];
      for (k = 1; k < i; k++)
        sum -= a[i][k] * a[k][j];
      a[i][j] = sum;
    }
    big = 0.0;
    for (i = j; i <= n; i++) {
      sum = a[i][j];
      for (k = 1; k < j; k++)
        sum -= a[i][k] * a[k][j];
      a[i][j] = sum;
      if ((dum = vv[i] * fabs(sum)) >= big) {
        big = dum;
        imax = i;
      }
    }
    if (j != imax) {
      for (k = 1; k <= n; k++) {
        dum = a[imax][k];
        a[imax][k] = a[j][k];
        a[j][k] = dum;
      }
      (*d) = -(*d);
      vv[imax] = vv[j];
    }
    indx[j] = imax;
    if (!a[j][j]) {
      a[j][j] = 1.0e-20;
    }
    if (j != n) {
      dum = 1.0 / (a[j][j]);
      for (i = j + 1; i <= n; i++)
        a[i][j] *= dum;
    }
  }
  free_vector(vv, 1, n);
  return (0);
}


x86 assembly (clang, -O0):
<ludcmp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0xc(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1050 <vector@plt>
mov    %rax,-0x58(%rbp)
mov    -0x20(%rbp),%rax
movsd  0xe93(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     123a <ludcmp+0x10a>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     11ec <ludcmp+0xbc>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe55(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
ucomisd -0x38(%rbp),%xmm0
jbe    11d9 <ludcmp+0xa9>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    11de <ludcmp+0xae>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1193 <ludcmp+0x63>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1212 <ludcmp+0xe2>
jp     1212 <ludcmp+0xe2>
lea    0xe15(%rip),%rdi
mov    $0x0,%al
call   1030 <nrerror@plt>
movsd  0xde6(%rip),%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1178 <ludcmp+0x48>
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     15d3 <ludcmp+0x4a3>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1316 <ludcmp+0x1e6>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    12ee <ludcmp+0x1be>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1281 <ludcmp+0x151>
movsd  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1254 <ludcmp+0x124>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     142c <ludcmp+0x2fc>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    13be <ludcmp+0x28e>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1351 <ludcmp+0x221>
movsd  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x48(%rbp),%xmm1
movaps 0xc1f(%rip),%xmm2
pand   %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
ucomisd -0x38(%rbp),%xmm0
jb     1419 <ludcmp+0x2e9>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
jmp    141e <ludcmp+0x2ee>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1324 <ludcmp+0x1f4>
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
je     14f8 <ludcmp+0x3c8>
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     14b7 <ludcmp+0x387>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    143f <ludcmp+0x30f>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x58(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%edx
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    154b <ludcmp+0x41b>
jp     154b <ludcmp+0x41b>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  0xac2(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
je     15c0 <ludcmp+0x490>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  0xa91(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     15bb <ludcmp+0x48b>
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1582 <ludcmp+0x452>
jmp    15c0 <ludcmp+0x490>
jmp    15c5 <ludcmp+0x495>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1241 <ludcmp+0x111>
mov    -0x58(%rbp),%rdi
mov    -0xc(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1040 <free_vector@plt>
xor    %eax,%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ludcmp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %rcx,%r14
mov    %rdx,0x70(%rsp)
mov    %esi,%ebp
mov    %rdi,%rbx
mov    $0x1,%edi
xor    %eax,%eax
call   1050 <vector@plt>
mov    %rax,%r13
movabs $0x3ff0000000000000,%rax
mov    %r14,0x68(%rsp)
mov    %rax,(%r14)
mov    %rbp,0x8(%rsp)
test   %ebp,%ebp
jle    1752 <ludcmp+0x622>
mov    0x8(%rsp),%ebp
lea    0x1(%rbp),%rax
mov    %rax,0x10(%rsp)
lea    -0x1(%rbp),%r12
mov    %ebp,%r14d
and    $0x3,%r14d
mov    %rbp,0x18(%rsp)
and    $0xfffffffc,%ebp
mov    $0x1,%r15d
xorpd  %xmm4,%xmm4
movsd  0xe75(%rip),%xmm7
movapd 0xe4c(%rip),%xmm8
movapd 0xe44(%rip),%xmm5
jmp    11de <ludcmp+0xae>
xchg   %ax,%ax
movapd %xmm7,%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x0(%r13,%r15,8)
add    $0x1,%r15
cmp    0x10(%rsp),%r15
je     12c7 <ludcmp+0x197>
mov    (%rbx,%r15,8),%rax
mov    $0x1,%ecx
xorpd  %xmm0,%xmm0
cmp    $0x3,%r12
jb     1245 <ludcmp+0x115>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movsd  0x8(%rax,%rcx,8),%xmm1
movsd  0x10(%rax,%rcx,8),%xmm2
movsd  0x18(%rax,%rcx,8),%xmm3
andpd  %xmm5,%xmm1
maxsd  %xmm0,%xmm1
movsd  0x20(%rax,%rcx,8),%xmm0
andpd  %xmm5,%xmm2
maxsd  %xmm1,%xmm2
andpd  %xmm5,%xmm3
maxsd  %xmm2,%xmm3
andpd  %xmm5,%xmm0
maxsd  %xmm3,%xmm0
add    $0x4,%rcx
cmp    %rcx,%rbp
jne    1200 <ludcmp+0xd0>
add    $0x1,%rcx
movapd %xmm0,%xmm1
test   %r14,%r14
je     127b <ludcmp+0x14b>
lea    (%rax,%rcx,8),%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rax,%rcx,8),%xmm1
andpd  %xmm8,%xmm1
maxsd  %xmm0,%xmm1
add    $0x1,%rcx
movapd %xmm1,%xmm0
cmp    %rcx,%r14
jne    1260 <ludcmp+0x130>
ucomisd %xmm4,%xmm1
jne    11c0 <ludcmp+0x90>
jp     11c0 <ludcmp+0x90>
lea    0xd9e(%rip),%rdi
xor    %eax,%eax
movapd %xmm1,0x20(%rsp)
call   1030 <nrerror@plt>
movapd 0x20(%rsp),%xmm1
movapd 0xd53(%rip),%xmm5
xorpd  %xmm4,%xmm4
movapd 0xd46(%rip),%xmm8
movsd  0xd5e(%rip),%xmm7
jmp    11c0 <ludcmp+0x90>
cmpl   $0x0,0x8(%rsp)
jle    1752 <ludcmp+0x622>
mov    0x8(%rsp),%rax
lea    0x1(%rax),%r12d
mov    0x10(%rsp),%rax
lea    -0x1(%rax),%rcx
add    $0xfffffffffffffffb,%rax
mov    %rax,0x58(%rsp)
mov    %rax,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %rcx,0x20(%rsp)
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x48(%rsp)
lea    0x1(%rcx),%rax
mov    %rax,0x40(%rsp)
mov    %rdx,%rax
mov    %rdx,0x50(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x38(%rsp)
mov    %r12,%rax
neg    %rax
mov    %rax,0x60(%rsp)
mov    0x18(%rsp),%rax
lea    -0x1(%rax),%r14
lea    0x10(%rbx),%r15
xor    %r11d,%r11d
mov    $0x1,%r9d
xorpd  %xmm0,%xmm0
movsd  0xcdb(%rip),%xmm1
movapd 0xcbb(%rip),%xmm2
xor    %r10d,%r10d
jmp    1379 <ludcmp+0x249>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r9
add    $0x1,%r11
add    $0xffffffffffffffff,%r14
add    $0x8,%r15
cmp    %r12,%r9
je     1752 <ludcmp+0x622>
mov    %r11,%r8
and    $0xfffffffffffffffe,%r8
cmp    $0x2,%r9
jb     1450 <ludcmp+0x320>
mov    $0x1,%ecx
xor    %edi,%edi
jmp    13b8 <ludcmp+0x288>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  %xmm3,0x0(%rbp,%r9,8)
add    $0x1,%rcx
add    $0x1,%rdi
cmp    %r9,%rcx
je     1450 <ludcmp+0x320>
mov    (%rbx,%rcx,8),%rbp
movsd  0x0(%rbp,%r9,8),%xmm3
cmp    $0x2,%rcx
jb     13a0 <ludcmp+0x270>
mov    $0x1,%edx
cmp    $0x1,%rdi
je     142b <ludcmp+0x2fb>
mov    %rdi,%rsi
and    $0xfffffffffffffffe,%rsi
neg    %rsi
mov    $0x1,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x0(%rbp,%rdx,8),%xmm4
movsd  0x8(%rbp,%rdx,8),%xmm5
mov    (%rbx,%rdx,8),%rax
mulsd  (%rax,%r9,8),%xmm4
mov    0x8(%rbx,%rdx,8),%rax
subsd  %xmm4,%xmm3
mulsd  (%rax,%r9,8),%xmm5
subsd  %xmm5,%xmm3
lea    (%rsi,%rdx,1),%rax
add    $0x2,%rax
add    $0x2,%rdx
cmp    $0x1,%rax
jne    13f0 <ludcmp+0x2c0>
test   $0x1,%dil
je     13a0 <ludcmp+0x270>
movsd  0x0(%rbp,%rdx,8),%xmm4
mov    (%rbx,%rdx,8),%rax
mulsd  (%rax,%r9,8),%xmm4
subsd  %xmm4,%xmm3
jmp    13a0 <ludcmp+0x270>
xchg   %ax,%ax
neg    %r8
xorpd  %xmm3,%xmm3
mov    %r9,%rcx
jmp    14a4 <ludcmp+0x374>
nopl   0x0(%rax)
movsd  %xmm4,(%rsi,%r9,8)
andpd  %xmm8,%xmm4
mulsd  0x0(%r13,%rcx,8),%xmm4
ucomisd %xmm3,%xmm4
movapd %xmm3,%xmm5
cmpnlesd %xmm4,%xmm5
movapd %xmm5,%xmm6
andnpd %xmm4,%xmm6
andpd  %xmm3,%xmm5
orpd   %xmm6,%xmm5
cmovae %ecx,%r10d
add    $0x1,%rcx
movapd %xmm5,%xmm3
cmp    %r12,%rcx
je     1530 <ludcmp+0x400>
mov    (%rbx,%rcx,8),%rsi
movsd  (%rsi,%r9,8),%xmm4
cmp    $0x2,%r9
jb     1460 <ludcmp+0x330>
mov    $0x1,%edx
cmp    $0x1,%r11
je     150a <ludcmp+0x3da>
mov    $0x1,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%rdx,8),%xmm5
movsd  0x8(%rsi,%rdx,8),%xmm6
mov    (%rbx,%rdx,8),%rax
mulsd  (%rax,%r9,8),%xmm5
mov    0x8(%rbx,%rdx,8),%rax
subsd  %xmm5,%xmm4
mulsd  (%rax,%r9,8),%xmm6
subsd  %xmm6,%xmm4
lea    (%r8,%rdx,1),%rax
add    $0x2,%rax
add    $0x2,%rdx
cmp    $0x1,%rax
jne    14d0 <ludcmp+0x3a0>
test   $0x1,%r11b
je     1460 <ludcmp+0x330>
movsd  (%rsi,%rdx,8),%xmm5
mov    (%rbx,%rdx,8),%rax
mulsd  (%rax,%r9,8),%xmm5
subsd  %xmm5,%xmm4
jmp    1460 <ludcmp+0x330>
nopl   0x0(%rax)
mov    %r10d,%eax
cmp    %rax,%r9
je     16d6 <ludcmp+0x5a6>
movslq %r10d,%r8
mov    (%rbx,%r8,8),%rbp
mov    (%rbx,%r9,8),%rcx
mov    $0x1,%edx
cmpq   $0x4,0x20(%rsp)
jb     1650 <ludcmp+0x520>
lea    0x8(%rbp),%rax
mov    0x10(%rsp),%rsi
lea    (%rcx,%rsi,8),%rsi
cmp    %rsi,%rax
jae    1587 <ludcmp+0x457>
mov    0x10(%rsp),%rax
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
lea    0x8(%rcx),%rsi
cmp    %rax,%rsi
jb     1650 <ludcmp+0x520>
cmpq   $0x4,0x58(%rsp)
jae    1593 <ludcmp+0x463>
xor    %esi,%esi
jmp    1602 <ludcmp+0x4d2>
mov    0x38(%rsp),%rdx
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movups 0x8(%rbp,%rsi,8),%xmm3
movups 0x18(%rbp,%rsi,8),%xmm4
movups 0x8(%rcx,%rsi,8),%xmm5
movups 0x18(%rcx,%rsi,8),%xmm6
movups %xmm5,0x8(%rbp,%rsi,8)
movups %xmm6,0x18(%rbp,%rsi,8)
movups %xmm3,0x8(%rcx,%rsi,8)
movups %xmm4,0x18(%rcx,%rsi,8)
movupd 0x28(%rbp,%rsi,8),%xmm3
movupd 0x38(%rbp,%rsi,8),%xmm4
movupd 0x28(%rcx,%rsi,8),%xmm5
movupd 0x38(%rcx,%rsi,8),%xmm6
movupd %xmm5,0x28(%rbp,%rsi,8)
movupd %xmm6,0x38(%rbp,%rsi,8)
movupd %xmm3,0x28(%rcx,%rsi,8)
movupd %xmm4,0x38(%rcx,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdx
jne    15a0 <ludcmp+0x470>
testb  $0x1,0x50(%rsp)
je     163f <ludcmp+0x50f>
shl    $0x3,%rsi
or     $0x8,%rsi
movupd 0x0(%rbp,%rsi,1),%xmm3
movupd 0x10(%rbp,%rsi,1),%xmm4
movupd (%rcx,%rsi,1),%xmm5
movupd 0x10(%rcx,%rsi,1),%xmm6
movupd %xmm5,0x0(%rbp,%rsi,1)
movupd %xmm6,0x10(%rbp,%rsi,1)
movupd %xmm3,(%rcx,%rsi,1)
movupd %xmm4,0x10(%rcx,%rsi,1)
mov    0x40(%rsp),%rdx
mov    0x48(%rsp),%rax
cmp    %rax,0x20(%rsp)
je     16b7 <ludcmp+0x587>
mov    %rdx,%rsi
test   $0x1,%r12b
jne    1673 <ludcmp+0x543>
movsd  0x0(%rbp,%rdx,8),%xmm3
movsd  (%rcx,%rdx,8),%xmm4
movsd  %xmm4,0x0(%rbp,%rdx,8)
movsd  %xmm3,(%rcx,%rdx,8)
lea    0x1(%rdx),%rsi
not    %rdx
cmp    0x60(%rsp),%rdx
je     16b7 <ludcmp+0x587>
nopl   (%rax)
movsd  0x0(%rbp,%rsi,8),%xmm3
movsd  (%rcx,%rsi,8),%xmm4
movsd  %xmm4,0x0(%rbp,%rsi,8)
movsd  %xmm3,(%rcx,%rsi,8)
movsd  0x8(%rbp,%rsi,8),%xmm3
movsd  0x8(%rcx,%rsi,8),%xmm4
movsd  %xmm4,0x8(%rbp,%rsi,8)
movsd  %xmm3,0x8(%rcx,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%r12
jne    1680 <ludcmp+0x550>
mov    0x68(%rsp),%rax
movsd  (%rax),%xmm3
xorpd  %xmm2,%xmm3
movlpd %xmm3,(%rax)
movsd  0x0(%r13,%r9,8),%xmm3
movsd  %xmm3,0x0(%r13,%r8,8)
mov    0x70(%rsp),%rax
mov    %r10d,(%rax,%r9,4)
mov    (%rbx,%r9,8),%rcx
movsd  (%rcx,%r9,8),%xmm4
ucomisd %xmm0,%xmm4
jne    16f1 <ludcmp+0x5c1>
jnp    1700 <ludcmp+0x5d0>
cmp    0x18(%rsp),%r9
je     1360 <ludcmp+0x230>
jmp    171d <ludcmp+0x5ed>
xchg   %ax,%ax
movabs $0x3bc79ca10c924223,%rax
9c c7 3b
mov    %rax,(%rcx,%r9,8)
movapd %xmm1,%xmm4
cmp    0x18(%rsp),%r9
je     1360 <ludcmp+0x230>
jae    1360 <ludcmp+0x230>
movapd %xmm7,%xmm3
divsd  %xmm4,%xmm3
xor    %ecx,%ecx
nopl   (%rax)
mov    (%r15,%rcx,8),%rax
movsd  (%rax,%r9,8),%xmm4
mulsd  %xmm3,%xmm4
movsd  %xmm4,(%rax,%r9,8)
add    $0x1,%rcx
cmp    %rcx,%r14
jne    1730 <ludcmp+0x600>
jmp    1360 <ludcmp+0x230>
mov    %r13,%rdi
mov    $0x1,%esi
mov    0x8(%rsp),%rdx
xor    %eax,%eax
call   1040 <free_vector@plt>
xor    %eax,%eax
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [34 x i8] c"Singular matrix in routine LUDCMP\00"
define i32 @ludcmp(double** %0, i32 %1, i32* %2, double* %3) {
%5 = alloca double**
%6 = alloca i32
%7 = alloca i32*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double*
store double** %0, double*** %5
store i32 %1, i32* %6
store i32* %2, i32** %7
store double* %3, double** %8
store i32 0, i32* %10
%18 = load i32, i32* %6
%19 = call double* (i32, i32, ...) bitcast (double* (...)* @vector to double* (i32, i32, ...)*)(i32 1, i32 %18)
store double* %19, double** %17
%20 = load double*, double** %8
store double 1.000000e+00, double* %20
store i32 1, i32* %9
br label %21
21: ; preds = %61, %4
%22 = load i32, i32* %9
%23 = load i32, i32* %6
%24 = icmp sle i32 %22, %23
br i1 %24, label %25, label %64
25: ; preds = %21
store double 0.000000e+00, double* %13
store i32 1, i32* %11
br label %26
26: ; preds = %46, %25
%27 = load i32, i32* %11
%28 = load i32, i32* %6
%29 = icmp sle i32 %27, %28
br i1 %29, label %30, label %49
30: ; preds = %26
%31 = load double**, double*** %5
%32 = load i32, i32* %9
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double*, double** %31, i64 %33
%35 = load double*, double** %34
%36 = load i32, i32* %11
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %35, i64 %37
%39 = load double, double* %38
%40 = call double @llvm.fabs.f64(double %39)
store double %40, double* %16
%41 = load double, double* %13
%42 = fcmp ogt double %40, %41
br i1 %42, label %43, label %45
43: ; preds = %30
%44 = load double, double* %16
store double %44, double* %13
br label %45
45: ; preds = %43, %30
br label %46
46: ; preds = %45
%47 = load i32, i32* %11
%48 = add nsw i32 %47, 1
store i32 %48, i32* %11
br label %26
49: ; preds = %26
%50 = load double, double* %13
%51 = fcmp une double %50, 0.000000e+00
br i1 %51, label %54, label %52
52: ; preds = %49
%53 = call i32 (i8*, ...) bitcast (i32 (...)* @nrerror to i32 (i8*, ...)*)(i8* getelementptr inbounds ([34 x i8], [34 x i8]* @.str, i64 0, i64 0))
br label %54
54: ; preds = %52, %49
%55 = load double, double* %13
%56 = fdiv double 1.000000e+00, %55
%57 = load double*, double** %17
%58 = load i32, i32* %9
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %57, i64 %59
store double %56, double* %60
br label %61
61: ; preds = %54
%62 = load i32, i32* %9
%63 = add nsw i32 %62, 1
store i32 %63, i32* %9
br label %21
64: ; preds = %21
store i32 1, i32* %11
br label %65
65: ; preds = %323, %64
%66 = load i32, i32* %11
%67 = load i32, i32* %6
%68 = icmp sle i32 %66, %67
br i1 %68, label %69, label %326
69: ; preds = %65
store i32 1, i32* %9
br label %70
70: ; preds = %123, %69
%71 = load i32, i32* %9
%72 = load i32, i32* %11
%73 = icmp slt i32 %71, %72
br i1 %73, label %74, label %126
74: ; preds = %70
%75 = load double**, double*** %5
%76 = load i32, i32* %9
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double*, double** %75, i64 %77
%79 = load double*, double** %78
%80 = load i32, i32* %11
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %79, i64 %81
%83 = load double, double* %82
store double %83, double* %15
store i32 1, i32* %12
br label %84
84: ; preds = %110, %74
%85 = load i32, i32* %12
%86 = load i32, i32* %9
%87 = icmp slt i32 %85, %86
br i1 %87, label %88, label %113
88: ; preds = %84
%89 = load double**, double*** %5
%90 = load i32, i32* %9
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double*, double** %89, i64 %91
%93 = load double*, double** %92
%94 = load i32, i32* %12
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double, double* %93, i64 %95
%97 = load double, double* %96
%98 = load double**, double*** %5
%99 = load i32, i32* %12
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double*, double** %98, i64 %100
%102 = load double*, double** %101
%103 = load i32, i32* %11
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds double, double* %102, i64 %104
%106 = load double, double* %105
%107 = load double, double* %15
%108 = fneg double %97
%109 = call double @llvm.fmuladd.f64(double %108, double %106, double %107)
store double %109, double* %15
br label %110
110: ; preds = %88
%111 = load i32, i32* %12
%112 = add nsw i32 %111, 1
store i32 %112, i32* %12
br label %84
113: ; preds = %84
%114 = load double, double* %15
%115 = load double**, double*** %5
%116 = load i32, i32* %9
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds double*, double** %115, i64 %117
%119 = load double*, double** %118
%120 = load i32, i32* %11
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %119, i64 %121
store double %114, double* %122
br label %123
123: ; preds = %113
%124 = load i32, i32* %9
%125 = add nsw i32 %124, 1
store i32 %125, i32* %9
br label %70
126: ; preds = %70
store double 0.000000e+00, double* %13
%127 = load i32, i32* %11
store i32 %127, i32* %9
br label %128
128: ; preds = %195, %126
%129 = load i32, i32* %9
%130 = load i32, i32* %6
%131 = icmp sle i32 %129, %130
br i1 %131, label %132, label %198
132: ; preds = %128
%133 = load double**, double*** %5
%134 = load i32, i32* %9
%135 = sext i32 %134 to i64
%136 = getelementptr inbounds double*, double** %133, i64 %135
%137 = load double*, double** %136
%138 = load i32, i32* %11
%139 = sext i32 %138 to i64
%140 = getelementptr inbounds double, double* %137, i64 %139
%141 = load double, double* %140
store double %141, double* %15
store i32 1, i32* %12
br label %142
142: ; preds = %168, %132
%143 = load i32, i32* %12
%144 = load i32, i32* %11
%145 = icmp slt i32 %143, %144
br i1 %145, label %146, label %171
146: ; preds = %142
%147 = load double**, double*** %5
%148 = load i32, i32* %9
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds double*, double** %147, i64 %149
%151 = load double*, double** %150
%152 = load i32, i32* %12
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds double, double* %151, i64 %153
%155 = load double, double* %154
%156 = load double**, double*** %5
%157 = load i32, i32* %12
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double*, double** %156, i64 %158
%160 = load double*, double** %159
%161 = load i32, i32* %11
%162 = sext i32 %161 to i64
%163 = getelementptr inbounds double, double* %160, i64 %162
%164 = load double, double* %163
%165 = load double, double* %15
%166 = fneg double %155
%167 = call double @llvm.fmuladd.f64(double %166, double %164, double %165)
store double %167, double* %15
br label %168
168: ; preds = %146
%169 = load i32, i32* %12
%170 = add nsw i32 %169, 1
store i32 %170, i32* %12
br label %142
171: ; preds = %142
%172 = load double, double* %15
%173 = load double**, double*** %5
%174 = load i32, i32* %9
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds double*, double** %173, i64 %175
%177 = load double*, double** %176
%178 = load i32, i32* %11
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds double, double* %177, i64 %179
store double %172, double* %180
%181 = load double*, double** %17
%182 = load i32, i32* %9
%183 = sext i32 %182 to i64
%184 = getelementptr inbounds double, double* %181, i64 %183
%185 = load double, double* %184
%186 = load double, double* %15
%187 = call double @llvm.fabs.f64(double %186)
%188 = fmul double %185, %187
store double %188, double* %14
%189 = load double, double* %13
%190 = fcmp oge double %188, %189
br i1 %190, label %191, label %194
191: ; preds = %171
%192 = load double, double* %14
store double %192, double* %13
%193 = load i32, i32* %9
store i32 %193, i32* %10
br label %194
194: ; preds = %191, %171
br label %195
195: ; preds = %194
%196 = load i32, i32* %9
%197 = add nsw i32 %196, 1
store i32 %197, i32* %9
br label %128
198: ; preds = %128
%199 = load i32, i32* %11
%200 = load i32, i32* %10
%201 = icmp ne i32 %199, %200
br i1 %201, label %202, label %260
202: ; preds = %198
store i32 1, i32* %12
br label %203
203: ; preds = %243, %202
%204 = load i32, i32* %12
%205 = load i32, i32* %6
%206 = icmp sle i32 %204, %205
br i1 %206, label %207, label %246
207: ; preds = %203
%208 = load double**, double*** %5
%209 = load i32, i32* %10
%210 = sext i32 %209 to i64
%211 = getelementptr inbounds double*, double** %208, i64 %210
%212 = load double*, double** %211
%213 = load i32, i32* %12
%214 = sext i32 %213 to i64
%215 = getelementptr inbounds double, double* %212, i64 %214
%216 = load double, double* %215
store double %216, double* %14
%217 = load double**, double*** %5
%218 = load i32, i32* %11
%219 = sext i32 %218 to i64
%220 = getelementptr inbounds double*, double** %217, i64 %219
%221 = load double*, double** %220
%222 = load i32, i32* %12
%223 = sext i32 %222 to i64
%224 = getelementptr inbounds double, double* %221, i64 %223
%225 = load double, double* %224
%226 = load double**, double*** %5
%227 = load i32, i32* %10
%228 = sext i32 %227 to i64
%229 = getelementptr inbounds double*, double** %226, i64 %228
%230 = load double*, double** %229
%231 = load i32, i32* %12
%232 = sext i32 %231 to i64
%233 = getelementptr inbounds double, double* %230, i64 %232
store double %225, double* %233
%234 = load double, double* %14
%235 = load double**, double*** %5
%236 = load i32, i32* %11
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds double*, double** %235, i64 %237
%239 = load double*, double** %238
%240 = load i32, i32* %12
%241 = sext i32 %240 to i64
%242 = getelementptr inbounds double, double* %239, i64 %241
store double %234, double* %242
br label %243
243: ; preds = %207
%244 = load i32, i32* %12
%245 = add nsw i32 %244, 1
store i32 %245, i32* %12
br label %203
246: ; preds = %203
%247 = load double*, double** %8
%248 = load double, double* %247
%249 = fneg double %248
%250 = load double*, double** %8
store double %249, double* %250
%251 = load double*, double** %17
%252 = load i32, i32* %11
%253 = sext i32 %252 to i64
%254 = getelementptr inbounds double, double* %251, i64 %253
%255 = load double, double* %254
%256 = load double*, double** %17
%257 = load i32, i32* %10
%258 = sext i32 %257 to i64
%259 = getelementptr inbounds double, double* %256, i64 %258
store double %255, double* %259
br label %260
260: ; preds = %246, %198
%261 = load i32, i32* %10
%262 = load i32*, i32** %7
%263 = load i32, i32* %11
%264 = sext i32 %263 to i64
%265 = getelementptr inbounds i32, i32* %262, i64 %264
store i32 %261, i32* %265
%266 = load double**, double*** %5
%267 = load i32, i32* %11
%268 = sext i32 %267 to i64
%269 = getelementptr inbounds double*, double** %266, i64 %268
%270 = load double*, double** %269
%271 = load i32, i32* %11
%272 = sext i32 %271 to i64
%273 = getelementptr inbounds double, double* %270, i64 %272
%274 = load double, double* %273
%275 = fcmp une double %274, 0.000000e+00
br i1 %275, label %285, label %276
276: ; preds = %260
%277 = load double**, double*** %5
%278 = load i32, i32* %11
%279 = sext i32 %278 to i64
%280 = getelementptr inbounds double*, double** %277, i64 %279
%281 = load double*, double** %280
%282 = load i32, i32* %11
%283 = sext i32 %282 to i64
%284 = getelementptr inbounds double, double* %281, i64 %283
store double 0x3BC79CA10C924223, double* %284
br label %285
285: ; preds = %276, %260
%286 = load i32, i32* %11
%287 = load i32, i32* %6
%288 = icmp ne i32 %286, %287
br i1 %288, label %289, label %322
289: ; preds = %285
%290 = load double**, double*** %5
%291 = load i32, i32* %11
%292 = sext i32 %291 to i64
%293 = getelementptr inbounds double*, double** %290, i64 %292
%294 = load double*, double** %293
%295 = load i32, i32* %11
%296 = sext i32 %295 to i64
%297 = getelementptr inbounds double, double* %294, i64 %296
%298 = load double, double* %297
%299 = fdiv double 1.000000e+00, %298
store double %299, double* %14
%300 = load i32, i32* %11
%301 = add nsw i32 %300, 1
store i32 %301, i32* %9
br label %302
302: ; preds = %318, %289
%303 = load i32, i32* %9
%304 = load i32, i32* %6
%305 = icmp sle i32 %303, %304
br i1 %305, label %306, label %321
306: ; preds = %302
%307 = load double, double* %14
%308 = load double**, double*** %5
%309 = load i32, i32* %9
%310 = sext i32 %309 to i64
%311 = getelementptr inbounds double*, double** %308, i64 %310
%312 = load double*, double** %311
%313 = load i32, i32* %11
%314 = sext i32 %313 to i64
%315 = getelementptr inbounds double, double* %312, i64 %314
%316 = load double, double* %315
%317 = fmul double %316, %307
store double %317, double* %315
br label %318
318: ; preds = %306
%319 = load i32, i32* %9
%320 = add nsw i32 %319, 1
store i32 %320, i32* %9
br label %302
321: ; preds = %302
br label %322
322: ; preds = %321, %285
br label %323
323: ; preds = %322
%324 = load i32, i32* %11
%325 = add nsw i32 %324, 1
store i32 %325, i32* %11
br label %65
326: ; preds = %65
%327 = load double*, double** %17
%328 = load i32, i32* %6
ret i32 0
}

LLVM IR (-O2):
@.str = [34 x i8] c"Singular matrix in routine LUDCMP\00"
define i32 @ludcmp(double** nocapture %0, i32 %1, i32* nocapture writeonly %2, double* nocapture %3) {
%5 = tail call double* (i32, i32, ...) bitcast (double* (...)* @vector to double* (i32, i32, ...)*)(i32 1, i32 %1)
store double 1.000000e+00, double* %3
%6 = icmp slt i32 %1, 1
br i1 %6, label %358, label %7
7: ; preds = %4
%8 = add nuw i32 %1, 1
%9 = zext i32 %8 to i64
%10 = add nsw i64 %9, -1
%11 = add nsw i64 %9, -2
%12 = and i64 %10, 3
%13 = icmp ult i64 %11, 3
%14 = and i64 %10, -4
%15 = icmp eq i64 %12, 0
br label %16
16: ; preds = %92, %7
%17 = phi i64 [ 1, %7 ], [ %95, %92 ]
%18 = getelementptr inbounds double*, double** %0, i64 %17
%19 = load double*, double** %18
br i1 %13, label %71, label %41
20: ; preds = %92
br i1 %6, label %358, label %21
21: ; preds = %20
%22 = add nuw i32 %1, 1
%23 = zext i32 %1 to i64
%24 = zext i32 %1 to i64
%25 = zext i32 %22 to i64
%26 = add nsw i64 %9, -1
%27 = add nsw i64 %9, -5
%28 = lshr i64 %27, 2
%29 = add nuw nsw i64 %28, 1
%30 = icmp ult i64 %26, 4
%31 = and i64 %26, -4
%32 = or i64 %31, 1
%33 = and i64 %29, 1
%34 = icmp ult i64 %27, 4
%35 = and i64 %29, 9223372036854775806
%36 = icmp eq i64 %33, 0
%37 = icmp eq i64 %26, %31
%38 = and i64 %25, 1
%39 = icmp eq i64 %38, 0
%40 = sub nsw i64 0, %25
br label %97
41: ; preds = %16, %41
%42 = phi i64 [ %68, %41 ], [ 1, %16 ]
%43 = phi double [ %67, %41 ], [ 0.000000e+00, %16 ]
%44 = phi i64 [ %69, %41 ], [ 0, %16 ]
%45 = getelementptr inbounds double, double* %19, i64 %42
%46 = load double, double* %45
%47 = tail call double @llvm.fabs.f64(double %46)
%48 = fcmp ogt double %47, %43
%49 = select i1 %48, double %47, double %43
%50 = add nuw nsw i64 %42, 1
%51 = getelementptr inbounds double, double* %19, i64 %50
%52 = load double, double* %51
%53 = tail call double @llvm.fabs.f64(double %52)
%54 = fcmp ogt double %53, %49
%55 = select i1 %54, double %53, double %49
%56 = add nuw nsw i64 %42, 2
%57 = getelementptr inbounds double, double* %19, i64 %56
%58 = load double, double* %57
%59 = tail call double @llvm.fabs.f64(double %58)
%60 = fcmp ogt double %59, %55
%61 = select i1 %60, double %59, double %55
%62 = add nuw nsw i64 %42, 3
%63 = getelementptr inbounds double, double* %19, i64 %62
%64 = load double, double* %63
%65 = tail call double @llvm.fabs.f64(double %64)
%66 = fcmp ogt double %65, %61
%67 = select i1 %66, double %65, double %61
%68 = add nuw nsw i64 %42, 4
%69 = add i64 %44, 4
%70 = icmp eq i64 %69, %14
br i1 %70, label %71, label %41
71: ; preds = %41, %16
%72 = phi double [ undef, %16 ], [ %67, %41 ]
%73 = phi i64 [ 1, %16 ], [ %68, %41 ]
%74 = phi double [ 0.000000e+00, %16 ], [ %67, %41 ]
br i1 %15, label %87, label %75
75: ; preds = %71, %75
%76 = phi i64 [ %84, %75 ], [ %73, %71 ]
%77 = phi double [ %83, %75 ], [ %74, %71 ]
%78 = phi i64 [ %85, %75 ], [ 0, %71 ]
%79 = getelementptr inbounds double, double* %19, i64 %76
%80 = load double, double* %79
%81 = tail call double @llvm.fabs.f64(double %80)
%82 = fcmp ogt double %81, %77
%83 = select i1 %82, double %81, double %77
%84 = add nuw nsw i64 %76, 1
%85 = add i64 %78, 1
%86 = icmp eq i64 %85, %12
br i1 %86, label %87, label %75
87: ; preds = %75, %71
%88 = phi double [ %72, %71 ], [ %83, %75 ]
%89 = fcmp une double %88, 0.000000e+00
br i1 %89, label %92, label %90
90: ; preds = %87
%91 = tail call i32 (i8*, ...) bitcast (i32 (...)* @nrerror to i32 (i8*, ...)*)(i8* getelementptr inbounds ([34 x i8], [34 x i8]* @.str, i64 0, i64 0))
br label %92
92: ; preds = %90, %87
%93 = fdiv double 1.000000e+00, %88
%94 = getelementptr inbounds double, double* %5, i64 %17
store double %93, double* %94
%95 = add nuw nsw i64 %17, 1
%96 = icmp eq i64 %95, %9
br i1 %96, label %20, label %16
97: ; preds = %21, %354
%98 = phi i64 [ 0, %21 ], [ %357, %354 ]
%99 = phi i64 [ 1, %21 ], [ %355, %354 ]
%100 = phi i32 [ 0, %21 ], [ %218, %354 ]
%101 = icmp ugt i64 %99, 1
br i1 %101, label %102, label %158
102: ; preds = %97, %153
%103 = phi i64 [ %157, %153 ], [ 0, %97 ]
%104 = phi i64 [ %155, %153 ], [ 1, %97 ]
%105 = getelementptr inbounds double*, double** %0, i64 %104
%106 = load double*, double** %105
%107 = getelementptr inbounds double, double* %106, i64 %99
%108 = load double, double* %107
%109 = icmp ugt i64 %104, 1
br i1 %109, label %110, label %153
110: ; preds = %102
%111 = and i64 %103, 1
%112 = icmp eq i64 %103, 1
br i1 %112, label %139, label %113
113: ; preds = %110
%114 = and i64 %103, -2
br label %115
115: ; preds = %115, %113
%116 = phi i64 [ 1, %113 ], [ %136, %115 ]
%117 = phi double [ %108, %113 ], [ %135, %115 ]
%118 = phi i64 [ 0, %113 ], [ %137, %115 ]
%119 = getelementptr inbounds double, double* %106, i64 %116
%120 = load double, double* %119
%121 = getelementptr inbounds double*, double** %0, i64 %116
%122 = load double*, double** %121
%123 = getelementptr inbounds double, double* %122, i64 %99
%124 = load double, double* %123
%125 = fneg double %120
%126 = tail call double @llvm.fmuladd.f64(double %125, double %124, double %117)
%127 = add nuw nsw i64 %116, 1
%128 = getelementptr inbounds double, double* %106, i64 %127
%129 = load double, double* %128
%130 = getelementptr inbounds double*, double** %0, i64 %127
%131 = load double*, double** %130
%132 = getelementptr inbounds double, double* %131, i64 %99
%133 = load double, double* %132
%134 = fneg double %129
%135 = tail call double @llvm.fmuladd.f64(double %134, double %133, double %126)
%136 = add nuw nsw i64 %116, 2
%137 = add i64 %118, 2
%138 = icmp eq i64 %137, %114
br i1 %138, label %139, label %115
139: ; preds = %115, %110
%140 = phi double [ undef, %110 ], [ %135, %115 ]
%141 = phi i64 [ 1, %110 ], [ %136, %115 ]
%142 = phi double [ %108, %110 ], [ %135, %115 ]
%143 = icmp eq i64 %111, 0
br i1 %143, label %153, label %144
144: ; preds = %139
%145 = getelementptr inbounds double, double* %106, i64 %141
%146 = load double, double* %145
%147 = getelementptr inbounds double*, double** %0, i64 %141
%148 = load double*, double** %147
%149 = getelementptr inbounds double, double* %148, i64 %99
%150 = load double, double* %149
%151 = fneg double %146
%152 = tail call double @llvm.fmuladd.f64(double %151, double %150, double %142)
br label %153
153: ; preds = %144, %139, %102
%154 = phi double [ %108, %102 ], [ %140, %139 ], [ %152, %144 ]
store double %154, double* %107
%155 = add nuw nsw i64 %104, 1
%156 = icmp eq i64 %155, %99
%157 = add i64 %103, 1
br i1 %156, label %158, label %102
158: ; preds = %153, %97
%159 = and i64 %98, 1
%160 = icmp eq i64 %98, 1
%161 = and i64 %98, -2
%162 = icmp eq i64 %159, 0
br label %163
163: ; preds = %158, %209
%164 = phi i64 [ %219, %209 ], [ %99, %158 ]
%165 = phi i32 [ %218, %209 ], [ %100, %158 ]
%166 = phi double [ %216, %209 ], [ 0.000000e+00, %158 ]
%167 = getelementptr inbounds double*, double** %0, i64 %164
%168 = load double*, double** %167
%169 = getelementptr inbounds double, double* %168, i64 %99
%170 = load double, double* %169
br i1 %101, label %171, label %209
171: ; preds = %163
br i1 %160, label %196, label %172
172: ; preds = %171, %172
%173 = phi i64 [ %193, %172 ], [ 1, %171 ]
%174 = phi double [ %192, %172 ], [ %170, %171 ]
%175 = phi i64 [ %194, %172 ], [ 0, %171 ]
%176 = getelementptr inbounds double, double* %168, i64 %173
%177 = load double, double* %176
%178 = getelementptr inbounds double*, double** %0, i64 %173
%179 = load double*, double** %178
%180 = getelementptr inbounds double, double* %179, i64 %99
%181 = load double, double* %180
%182 = fneg double %177
%183 = tail call double @llvm.fmuladd.f64(double %182, double %181, double %174)
%184 = add nuw nsw i64 %173, 1
%185 = getelementptr inbounds double, double* %168, i64 %184
%186 = load double, double* %185
%187 = getelementptr inbounds double*, double** %0, i64 %184
%188 = load double*, double** %187
%189 = getelementptr inbounds double, double* %188, i64 %99
%190 = load double, double* %189
%191 = fneg double %186
%192 = tail call double @llvm.fmuladd.f64(double %191, double %190, double %183)
%193 = add nuw nsw i64 %173, 2
%194 = add i64 %175, 2
%195 = icmp eq i64 %194, %161
br i1 %195, label %196, label %172
196: ; preds = %172, %171
%197 = phi double [ undef, %171 ], [ %192, %172 ]
%198 = phi i64 [ 1, %171 ], [ %193, %172 ]
%199 = phi double [ %170, %171 ], [ %192, %172 ]
br i1 %162, label %209, label %200
200: ; preds = %196
%201 = getelementptr inbounds double, double* %168, i64 %198
%202 = load double, double* %201
%203 = getelementptr inbounds double*, double** %0, i64 %198
%204 = load double*, double** %203
%205 = getelementptr inbounds double, double* %204, i64 %99
%206 = load double, double* %205
%207 = fneg double %202
%208 = tail call double @llvm.fmuladd.f64(double %207, double %206, double %199)
br label %209
209: ; preds = %200, %196, %163
%210 = phi double [ %170, %163 ], [ %197, %196 ], [ %208, %200 ]
store double %210, double* %169
%211 = getelementptr inbounds double, double* %5, i64 %164
%212 = load double, double* %211
%213 = tail call double @llvm.fabs.f64(double %210)
%214 = fmul double %213, %212
%215 = fcmp ult double %214, %166
%216 = select i1 %215, double %166, double %214
%217 = trunc i64 %164 to i32
%218 = select i1 %215, i32 %165, i32 %217
%219 = add nuw nsw i64 %164, 1
%220 = icmp eq i64 %219, %25
br i1 %220, label %221, label %163
221: ; preds = %209
%222 = zext i32 %218 to i64
%223 = icmp eq i64 %99, %222
br i1 %223, label %331, label %224
224: ; preds = %221
%225 = sext i32 %218 to i64
%226 = getelementptr inbounds double*, double** %0, i64 %225
%227 = load double*, double** %226
%228 = getelementptr inbounds double*, double** %0, i64 %99
%229 = load double*, double** %228
br i1 %30, label %300, label %230
230: ; preds = %224
%231 = getelementptr double, double* %227, i64 1
%232 = getelementptr double, double* %227, i64 %9
%233 = getelementptr double, double* %229, i64 1
%234 = getelementptr double, double* %229, i64 %9
%235 = icmp ult double* %231, %234
%236 = icmp ult double* %233, %232
%237 = and i1 %235, %236
br i1 %237, label %300, label %238
238: ; preds = %230
br i1 %34, label %279, label %239
239: ; preds = %238, %239
%240 = phi i64 [ %276, %239 ], [ 0, %238 ]
%241 = phi i64 [ %277, %239 ], [ 0, %238 ]
%242 = or i64 %240, 1
%243 = getelementptr inbounds double, double* %227, i64 %242
%244 = bitcast double* %243 to <2 x double>*
%245 = load <2 x double>, <2 x double>* %244
%246 = getelementptr inbounds double, double* %243, i64 2
%247 = bitcast double* %246 to <2 x double>*
%248 = load <2 x double>, <2 x double>* %247
%249 = getelementptr inbounds double, double* %229, i64 %242
%250 = bitcast double* %249 to <2 x double>*
%251 = load <2 x double>, <2 x double>* %250
%252 = getelementptr inbounds double, double* %249, i64 2
%253 = bitcast double* %252 to <2 x double>*
%254 = load <2 x double>, <2 x double>* %253
%255 = bitcast double* %243 to <2 x double>*
store <2 x double> %251, <2 x double>* %255
%256 = bitcast double* %246 to <2 x double>*
store <2 x double> %254, <2 x double>* %256
%257 = bitcast double* %249 to <2 x double>*
store <2 x double> %245, <2 x double>* %257
%258 = bitcast double* %252 to <2 x double>*
store <2 x double> %248, <2 x double>* %258
%259 = or i64 %240, 5
%260 = getelementptr inbounds double, double* %227, i64 %259
%261 = bitcast double* %260 to <2 x double>*
%262 = load <2 x double>, <2 x double>* %261
%263 = getelementptr inbounds double, double* %260, i64 2
%264 = bitcast double* %263 to <2 x double>*
%265 = load <2 x double>, <2 x double>* %264
%266 = getelementptr inbounds double, double* %229, i64 %259
%267 = bitcast double* %266 to <2 x double>*
%268 = load <2 x double>, <2 x double>* %267
%269 = getelementptr inbounds double, double* %266, i64 2
%270 = bitcast double* %269 to <2 x double>*
%271 = load <2 x double>, <2 x double>* %270
%272 = bitcast double* %260 to <2 x double>*
store <2 x double> %268, <2 x double>* %272
%273 = bitcast double* %263 to <2 x double>*
store <2 x double> %271, <2 x double>* %273
%274 = bitcast double* %266 to <2 x double>*
store <2 x double> %262, <2 x double>* %274
%275 = bitcast double* %269 to <2 x double>*
store <2 x double> %265, <2 x double>* %275
%276 = add nuw i64 %240, 8
%277 = add i64 %241, 2
%278 = icmp eq i64 %277, %35
br i1 %278, label %279, label %239
279: ; preds = %239, %238
%280 = phi i64 [ 0, %238 ], [ %276, %239 ]
br i1 %36, label %299, label %281
281: ; preds = %279
%282 = or i64 %280, 1
%283 = getelementptr inbounds double, double* %227, i64 %282
%284 = bitcast double* %283 to <2 x double>*
%285 = load <2 x double>, <2 x double>* %284
%286 = getelementptr inbounds double, double* %283, i64 2
%287 = bitcast double* %286 to <2 x double>*
%288 = load <2 x double>, <2 x double>* %287
%289 = getelementptr inbounds double, double* %229, i64 %282
%290 = bitcast double* %289 to <2 x double>*
%291 = load <2 x double>, <2 x double>* %290
%292 = getelementptr inbounds double, double* %289, i64 2
%293 = bitcast double* %292 to <2 x double>*
%294 = load <2 x double>, <2 x double>* %293
%295 = bitcast double* %283 to <2 x double>*
store <2 x double> %291, <2 x double>* %295
%296 = bitcast double* %286 to <2 x double>*
store <2 x double> %294, <2 x double>* %296
%297 = bitcast double* %289 to <2 x double>*
store <2 x double> %285, <2 x double>* %297
%298 = bitcast double* %292 to <2 x double>*
store <2 x double> %288, <2 x double>* %298
br label %299
299: ; preds = %279, %281
br i1 %37, label %325, label %300
300: ; preds = %230, %224, %299
%301 = phi i64 [ 1, %230 ], [ 1, %224 ], [ %32, %299 ]
%302 = xor i64 %301, -1
br i1 %39, label %303, label %309
303: ; preds = %300
%304 = getelementptr inbounds double, double* %227, i64 %301
%305 = load double, double* %304
%306 = getelementptr inbounds double, double* %229, i64 %301
%307 = load double, double* %306
store double %307, double* %304
store double %305, double* %306
%308 = add nuw nsw i64 %301, 1
br label %309
309: ; preds = %303, %300
%310 = phi i64 [ %301, %300 ], [ %308, %303 ]
%311 = icmp eq i64 %302, %40
br i1 %311, label %325, label %312
312: ; preds = %309, %312
%313 = phi i64 [ %323, %312 ], [ %310, %309 ]
%314 = getelementptr inbounds double, double* %227, i64 %313
%315 = load double, double* %314
%316 = getelementptr inbounds double, double* %229, i64 %313
%317 = load double, double* %316
store double %317, double* %314
store double %315, double* %316
%318 = add nuw nsw i64 %313, 1
%319 = getelementptr inbounds double, double* %227, i64 %318
%320 = load double, double* %319
%321 = getelementptr inbounds double, double* %229, i64 %318
%322 = load double, double* %321
store double %322, double* %319
store double %320, double* %321
%323 = add nuw nsw i64 %313, 2
%324 = icmp eq i64 %323, %25
br i1 %324, label %325, label %312
325: ; preds = %309, %312, %299
%326 = load double, double* %3
%327 = fneg double %326
store double %327, double* %3
%328 = getelementptr inbounds double, double* %5, i64 %99
%329 = load double, double* %328
%330 = getelementptr inbounds double, double* %5, i64 %225
store double %329, double* %330
br label %331
331: ; preds = %325, %221
%332 = getelementptr inbounds i32, i32* %2, i64 %99
store i32 %218, i32* %332
%333 = getelementptr inbounds double*, double** %0, i64 %99
%334 = load double*, double** %333
%335 = getelementptr inbounds double, double* %334, i64 %99
%336 = load double, double* %335
%337 = fcmp une double %336, 0.000000e+00
br i1 %337, label %339, label %338
338: ; preds = %331
store double 0x3BC79CA10C924223, double* %335
br label %339
339: ; preds = %338, %331
%340 = phi double [ 0x3BC79CA10C924223, %338 ], [ %336, %331 ]
%341 = icmp eq i64 %99, %24
br i1 %341, label %354, label %342
342: ; preds = %339
%343 = fdiv double 1.000000e+00, %340
%344 = icmp ult i64 %99, %23
br i1 %344, label %345, label %354
345: ; preds = %342, %345
%346 = phi i64 [ %347, %345 ], [ %99, %342 ]
%347 = add nuw nsw i64 %346, 1
%348 = getelementptr inbounds double*, double** %0, i64 %347
%349 = load double*, double** %348
%350 = getelementptr inbounds double, double* %349, i64 %99
%351 = load double, double* %350
%352 = fmul double %343, %351
store double %352, double* %350
%353 = icmp eq i64 %347, %24
br i1 %353, label %354, label %345
354: ; preds = %345, %342, %339
%355 = add nuw nsw i64 %99, 1
%356 = icmp eq i64 %355, %25
%357 = add i64 %98, 1
br i1 %356, label %358, label %97
358: ; preds = %354, %4, %20
tail call void (double*, i32, i32, ...) bitcast (void (...)* @free_vector to void (double*, i32, i32, ...)*)(double* %5, i32 1, i32 %1)
ret i32 0
}