#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double *dvector(int);
double fabs(double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
int /*<<< orphan*/ free_dvector(double *);
int /*<<< orphan*/ stderr;
int ludcmp(double *a[], int n, int indx[], double *d) {
  int i, imax = 0, j, k;
  double big, dum, sum, temp;
  double *vv;
  vv = dvector(n);
  *d = 1.0;
  for (i = 0; i < n; i++) {
    big = 0;
    for (j = 0; j < n; j++)
      if ((temp = fabs(a[i][j])) > big)
        big = temp;
    if (big == 0.0) {
      fprintf(stderr, "ludcmp: singular matrix\n");
      return 1;
    }
    vv[i] = 1.0 / big;
  }
  for (j = 0; j < n; j++) {
    for (i = 0; i < j; i++) {
      sum = a[i][j];
      for (k = 0; 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 = 0; 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 = 0; 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] == 0.0)
      a[j][j] = 1.0e-32;
    if (j != n) {
      dum = 1.0 / a[j][j];
      for (i = j + 1; i < n; i++)
        a[i][j] *= dum;
    }
  }
  free_dvector(vv);
  return 0;
}


x86 assembly (clang, -O0):
<ludcmp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x14(%rbp),%edi
call   1040 <dvector@plt>
mov    %rax,-0x60(%rbp)
mov    -0x28(%rbp),%rax
movsd  0xe9a(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1246 <ludcmp+0x116>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11e5 <ludcmp+0xb5>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe5c(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
ucomisd -0x40(%rbp),%xmm0
jbe    11d2 <ludcmp+0xa2>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    11d7 <ludcmp+0xa7>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    118c <ludcmp+0x5c>
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    121e <ludcmp+0xee>
jp     121e <ludcmp+0xee>
mov    0x2ddc(%rip),%rax
mov    (%rax),%edi
lea    0xe13(%rip),%rsi
call   1030 <fprintf@plt>
movl   $0x1,-0x4(%rbp)
jmp    15ef <ludcmp+0x4bf>
movsd  0xdda(%rip),%xmm0
divsd  -0x40(%rbp),%xmm0
mov    -0x60(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1171 <ludcmp+0x41>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    15df <ludcmp+0x4af>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1322 <ludcmp+0x1f2>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    12fa <ludcmp+0x1ca>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    128d <ludcmp+0x15d>
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1260 <ludcmp+0x130>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1438 <ludcmp+0x308>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    13ca <ludcmp+0x29a>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    135d <ludcmp+0x22d>
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm1
movaps 0xc13(%rip),%xmm2
pand   %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
ucomisd -0x40(%rbp),%xmm0
jb     1425 <ludcmp+0x2f5>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    142a <ludcmp+0x2fa>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1330 <ludcmp+0x200>
mov    -0x34(%rbp),%eax
cmp    -0x30(%rbp),%eax
je     1504 <ludcmp+0x3d4>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    14c3 <ludcmp+0x393>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    144b <ludcmp+0x31b>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x60(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x60(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%edx
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1557 <ludcmp+0x427>
jp     1557 <ludcmp+0x427>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  0xab6(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
je     15cc <ludcmp+0x49c>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  0xa85(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    15c7 <ludcmp+0x497>
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    158e <ludcmp+0x45e>
jmp    15cc <ludcmp+0x49c>
jmp    15d1 <ludcmp+0x4a1>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    124d <ludcmp+0x11d>
mov    -0x60(%rbp),%rdi
call   1050 <free_dvector@plt>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%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    $0x48,%rsp
mov    %rcx,%rbx
mov    %rdx,0x40(%rsp)
mov    %esi,%ebp
mov    %rdi,%r15
mov    %esi,%edi
call   1040 <dvector@plt>
movabs $0x3ff0000000000000,%rcx
mov    %rbx,0x38(%rsp)
mov    %rcx,(%rbx)
mov    %ebp,0xc(%rsp)
test   %ebp,%ebp
jle    16a2 <ludcmp+0x572>
mov    0xc(%rsp),%r13d
lea    -0x1(%r13),%r8
mov    %r13d,%esi
and    $0x3,%esi
mov    %r13d,%edi
and    $0xfffffffc,%edi
xor    %ebx,%ebx
xorpd  %xmm2,%xmm2
movsd  0xe8c(%rip),%xmm8
movapd 0xe64(%rip),%xmm1
movapd 0xe5c(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    (%r15,%rbx,8),%rdx
cmp    $0x3,%r8
jae    11d0 <ludcmp+0xa0>
xor    %ecx,%ecx
xorpd  %xmm4,%xmm4
jmp    1224 <ludcmp+0xf4>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorpd  %xmm5,%xmm5
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rdx,%rcx,8),%xmm6
movsd  0x8(%rdx,%rcx,8),%xmm7
movsd  0x10(%rdx,%rcx,8),%xmm0
movsd  0x18(%rdx,%rcx,8),%xmm4
andpd  %xmm3,%xmm6
maxsd  %xmm5,%xmm6
andpd  %xmm3,%xmm7
maxsd  %xmm6,%xmm7
andpd  %xmm3,%xmm0
maxsd  %xmm7,%xmm0
andpd  %xmm3,%xmm4
maxsd  %xmm0,%xmm4
add    $0x4,%rcx
movapd %xmm4,%xmm5
cmp    %rcx,%rdi
jne    11e0 <ludcmp+0xb0>
movapd %xmm4,%xmm5
test   %rsi,%rsi
je     125a <ludcmp+0x12a>
lea    (%rdx,%rcx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rcx,%rdx,8),%xmm5
andpd  %xmm1,%xmm5
maxsd  %xmm4,%xmm5
add    $0x1,%rdx
movapd %xmm5,%xmm4
cmp    %rdx,%rsi
jne    1240 <ludcmp+0x110>
ucomisd %xmm2,%xmm5
jne    1266 <ludcmp+0x136>
jnp    16ae <ludcmp+0x57e>
movapd %xmm8,%xmm0
divsd  %xmm5,%xmm0
movsd  %xmm0,(%rax,%rbx,8)
add    $0x1,%rbx
cmp    %r13,%rbx
jne    11b0 <ludcmp+0x80>
cmpl   $0x0,0xc(%rsp)
jle    16a2 <ludcmp+0x572>
mov    %r13d,%ecx
and    $0xfffffffc,%ecx
add    $0xfffffffffffffffc,%rcx
mov    %rcx,0x28(%rsp)
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %r13d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,0x18(%rsp)
mov    %rdx,%rcx
mov    %rdx,0x20(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x10(%rsp)
mov    %r13,%rcx
neg    %rcx
mov    %rcx,0x30(%rsp)
xor    %edx,%edx
xorpd  %xmm2,%xmm2
movsd  0xd4d(%rip),%xmm3
movapd 0xd2d(%rip),%xmm4
mov    %r13,%r14
mov    %r15,%r12
xor    %r8d,%r8d
jmp    1304 <ludcmp+0x1d4>
xchg   %ax,%ax
add    $0x8,%r12
add    $0xffffffffffffffff,%r14
mov    %rcx,%rdx
cmp    %r13,%rcx
je     16a2 <ludcmp+0x572>
test   %rdx,%rdx
je     13b0 <ludcmp+0x280>
xor    %r9d,%r9d
jmp    1332 <ludcmp+0x202>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  %xmm5,(%rbx,%rdx,8)
add    $0x1,%r9
cmp    %rdx,%r9
je     13b0 <ludcmp+0x280>
mov    (%r15,%r9,8),%rbx
movsd  (%rbx,%rdx,8),%xmm5
test   %r9,%r9
je     1320 <ludcmp+0x1f0>
cmp    $0x1,%r9
jne    1350 <ludcmp+0x220>
xor    %edi,%edi
jmp    1390 <ludcmp+0x260>
nopw   0x0(%rax,%rax,1)
mov    %r9,%rcx
and    $0xfffffffffffffffe,%rcx
xor    %edi,%edi
nopl   0x0(%rax)
movsd  (%rbx,%rdi,8),%xmm0
movsd  0x8(%rbx,%rdi,8),%xmm6
mov    (%r15,%rdi,8),%rbp
mov    0x8(%r15,%rdi,8),%rsi
mulsd  0x0(%rbp,%rdx,8),%xmm0
mulsd  (%rsi,%rdx,8),%xmm6
subsd  %xmm0,%xmm5
subsd  %xmm6,%xmm5
add    $0x2,%rdi
cmp    %rdi,%rcx
jne    1360 <ludcmp+0x230>
test   $0x1,%r9b
je     1320 <ludcmp+0x1f0>
movsd  (%rbx,%rdi,8),%xmm0
mov    (%r15,%rdi,8),%rcx
mulsd  (%rcx,%rdx,8),%xmm0
subsd  %xmm0,%xmm5
jmp    1320 <ludcmp+0x1f0>
nopl   (%rax)
mov    %rdx,%r9
and    $0xfffffffffffffffe,%r9
xorpd  %xmm5,%xmm5
mov    %rdx,%rsi
jmp    13fc <ludcmp+0x2cc>
movsd  %xmm6,(%rdi,%rdx,8)
andpd  %xmm1,%xmm6
mulsd  (%rax,%rsi,8),%xmm6
ucomisd %xmm5,%xmm6
cmovae %esi,%r8d
movapd %xmm5,%xmm0
cmpnlesd %xmm6,%xmm0
movapd %xmm0,%xmm7
andnpd %xmm6,%xmm7
andpd  %xmm5,%xmm0
orpd   %xmm7,%xmm0
add    $0x1,%rsi
movapd %xmm0,%xmm5
cmp    %r13,%rsi
je     1470 <ludcmp+0x340>
mov    (%r15,%rsi,8),%rdi
movsd  (%rdi,%rdx,8),%xmm6
test   %rdx,%rdx
je     13c0 <ludcmp+0x290>
xor    %ecx,%ecx
cmp    $0x1,%rdx
je     1450 <ludcmp+0x320>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rdi,%rcx,8),%xmm0
movsd  0x8(%rdi,%rcx,8),%xmm7
mov    (%r15,%rcx,8),%rbp
mov    0x8(%r15,%rcx,8),%rbx
mulsd  0x0(%rbp,%rdx,8),%xmm0
mulsd  (%rbx,%rdx,8),%xmm7
subsd  %xmm0,%xmm6
subsd  %xmm7,%xmm6
add    $0x2,%rcx
cmp    %rcx,%r9
jne    1420 <ludcmp+0x2f0>
test   $0x1,%dl
je     13c0 <ludcmp+0x290>
movsd  (%rdi,%rcx,8),%xmm0
mov    (%r15,%rcx,8),%rcx
mulsd  (%rcx,%rdx,8),%xmm0
subsd  %xmm0,%xmm6
jmp    13c0 <ludcmp+0x290>
mov    %r8d,%ecx
cmp    %rcx,%rdx
je     161c <ludcmp+0x4ec>
movslq %r8d,%r9
mov    (%r15,%r9,8),%r10
mov    (%r15,%rdx,8),%r11
cmpl   $0x4,0xc(%rsp)
jb     14a8 <ludcmp+0x378>
lea    (%r11,%r13,8),%rcx
cmp    %rcx,%r10
jae    1522 <ludcmp+0x3f2>
lea    (%r10,%r13,8),%rcx
cmp    %rcx,%r11
jae    1522 <ludcmp+0x3f2>
xor    %esi,%esi
mov    %rsi,%rcx
test   $0x1,%r13b
je     14d2 <ludcmp+0x3a2>
movsd  (%r10,%rsi,8),%xmm0
movsd  (%r11,%rsi,8),%xmm5
movsd  %xmm5,(%r10,%rsi,8)
movsd  %xmm0,(%r11,%rsi,8)
mov    %rsi,%rcx
or     $0x1,%rcx
not    %rsi
cmp    0x30(%rsp),%rsi
je     1600 <ludcmp+0x4d0>
movsd  (%r10,%rcx,8),%xmm0
movsd  (%r11,%rcx,8),%xmm5
movsd  %xmm5,(%r10,%rcx,8)
movsd  %xmm0,(%r11,%rcx,8)
movsd  0x8(%r10,%rcx,8),%xmm0
movsd  0x8(%r11,%rcx,8),%xmm5
movsd  %xmm5,0x8(%r10,%rcx,8)
movsd  %xmm0,0x8(%r11,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r13
jne    14e0 <ludcmp+0x3b0>
jmp    1600 <ludcmp+0x4d0>
cmpq   $0x0,0x28(%rsp)
je     1690 <ludcmp+0x560>
mov    0x10(%rsp),%rcx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movups (%r10,%rsi,8),%xmm0
movups 0x10(%r10,%rsi,8),%xmm5
movups (%r11,%rsi,8),%xmm6
movups 0x10(%r11,%rsi,8),%xmm7
movups %xmm6,(%r10,%rsi,8)
movups %xmm7,0x10(%r10,%rsi,8)
movups %xmm0,(%r11,%rsi,8)
movups %xmm5,0x10(%r11,%rsi,8)
movupd 0x20(%r10,%rsi,8),%xmm0
movupd 0x30(%r10,%rsi,8),%xmm5
movupd 0x20(%r11,%rsi,8),%xmm6
movupd 0x30(%r11,%rsi,8),%xmm7
movupd %xmm6,0x20(%r10,%rsi,8)
movupd %xmm7,0x30(%r10,%rsi,8)
movupd %xmm0,0x20(%r11,%rsi,8)
movupd %xmm5,0x30(%r11,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rcx
jne    1540 <ludcmp+0x410>
testb  $0x1,0x20(%rsp)
je     15e9 <ludcmp+0x4b9>
movupd (%r10,%rsi,8),%xmm0
movupd 0x10(%r10,%rsi,8),%xmm5
movupd (%r11,%rsi,8),%xmm6
movupd 0x10(%r11,%rsi,8),%xmm7
movupd %xmm6,(%r10,%rsi,8)
movupd %xmm7,0x10(%r10,%rsi,8)
movupd %xmm0,(%r11,%rsi,8)
movupd %xmm5,0x10(%r11,%rsi,8)
mov    0x18(%rsp),%rcx
mov    %rcx,%rsi
cmp    %r13,%rcx
jne    14aa <ludcmp+0x37a>
nopw   0x0(%rax,%rax,1)
mov    0x38(%rsp),%rcx
movsd  (%rcx),%xmm0
xorpd  %xmm4,%xmm0
movlpd %xmm0,(%rcx)
movsd  (%rax,%rdx,8),%xmm0
movsd  %xmm0,(%rax,%r9,8)
mov    0x40(%rsp),%rcx
mov    %r8d,(%rcx,%rdx,4)
mov    (%r15,%rdx,8),%rcx
movsd  (%rcx,%rdx,8),%xmm6
ucomisd %xmm2,%xmm6
jne    1648 <ludcmp+0x518>
jp     1648 <ludcmp+0x518>
movabs $0x3949f623d5a8a733,%rsi
f6 49 39
mov    %rsi,(%rcx,%rdx,8)
movapd %xmm3,%xmm6
lea    0x1(%rdx),%rcx
cmp    %r13,%rcx
jae    12f0 <ludcmp+0x1c0>
movapd %xmm8,%xmm5
divsd  %xmm6,%xmm5
mov    $0x1,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    (%r12,%rsi,8),%rdi
movsd  (%rdi,%rdx,8),%xmm0
mulsd  %xmm5,%xmm0
movsd  %xmm0,(%rdi,%rdx,8)
add    $0x1,%rsi
cmp    %rsi,%r14
jne    1670 <ludcmp+0x540>
jmp    12f0 <ludcmp+0x1c0>
xor    %esi,%esi
testb  $0x1,0x20(%rsp)
jne    15b5 <ludcmp+0x485>
jmp    15e9 <ludcmp+0x4b9>
mov    %rax,%rdi
call   1050 <free_dvector@plt>
xor    %eax,%eax
jmp    16c8 <ludcmp+0x598>
mov    0x292b(%rip),%rax
mov    (%rax),%edi
lea    0x972(%rip),%rsi
call   1030 <fprintf@plt>
mov    $0x1,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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