#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 const);
int /*<<< orphan*/ error(char *);
double fabs(double);
int /*<<< orphan*/ free(double *);
double LUDecompose_det(double **m, const int dim) {
  double big = 0.0;
  double temp = 0.0;
  double *vv = dvector(dim);
  double d = 1.0;
  for (size_t i = 0; i < dim; i++) {
    big = 0.0;
    for (size_t j = 0; j < dim; j++) {
      if ((temp = fabs(m[i][j])) > big)
        big = temp;
    }
    if (big == 0.0) {
      error("LU decomposition does not like singular matrices");
    }
    vv[i] = 1.0 / big;
  }
  size_t i, j;
  size_t imax;
  for (size_t k = 0; k < dim; k++) {
    big = 0.0;
    imax = k;
    for (i = k + 1; i < dim; i++) {
      temp = vv[i] * fabs(m[i][k]);
      if (temp > big) {
        big = temp;
        imax = i;
      }
    }
    if (k != imax) {
      for (j = 0; j < dim; j++) {
        temp = m[imax][j];
        m[imax][j] = m[k][j];
        m[k][j] = temp;
      }
      d = -d;
      vv[imax] = vv[k];
    }
    if (m[k][k] == 0.0)
      m[k][k] = 1.0e-25;
    for (i = k + 1; i < dim; i++) {
      temp = m[i][k] /= m[k][k];
      for (j = k + 1; j < dim; j++) {
        m[i][j] -= temp * m[k][j];
      }
    }
  }
  for (size_t i = 0; i < dim; i++) {
    d *= m[i][i];
  }
  free(vv);
  return d;
}


x86 assembly (clang, -O0):
<LUDecompose_det>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0xc(%rbp),%edi
call   1040 <dvector@plt>
mov    %rax,-0x28(%rbp)
movsd  0xe9d(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    1241 <LUDecompose_det+0x111>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    11f2 <LUDecompose_det+0xc2>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe52(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
ucomisd -0x18(%rbp),%xmm0
jbe    11dc <LUDecompose_det+0xac>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    11e1 <LUDecompose_det+0xb1>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1191 <LUDecompose_det+0x61>
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1216 <LUDecompose_det+0xe6>
jp     1216 <LUDecompose_det+0xe6>
lea    0xe0f(%rip),%rdi
call   1050 <error@plt>
movsd  0xde2(%rip),%xmm0
divsd  -0x18(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1170 <LUDecompose_det+0x40>
movq   $0x0,-0x60(%rbp)
mov    -0x60(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    150d <LUDecompose_det+0x3dd>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x60(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    12f5 <LUDecompose_det+0x1c5>
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xd60(%rip),%xmm2
pand   %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    12df <LUDecompose_det+0x1af>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x48(%rbp),%rax
mov    %rax,-0x58(%rbp)
jmp    12e4 <LUDecompose_det+0x1b4>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    1276 <LUDecompose_det+0x146>
mov    -0x60(%rbp),%rax
cmp    -0x58(%rbp),%rax
je     13c6 <LUDecompose_det+0x296>
movq   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    138b <LUDecompose_det+0x25b>
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    130b <LUDecompose_det+0x1db>
movsd  -0x30(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x58(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    140b <LUDecompose_det+0x2db>
jp     140b <LUDecompose_det+0x2db>
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x60(%rbp),%rcx
movsd  0xc02(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    14f7 <LUDecompose_det+0x3c7>
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  %xmm0,-0x20(%rbp)
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    14e1 <LUDecompose_det+0x3b1>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    146c <LUDecompose_det+0x33c>
jmp    14e6 <LUDecompose_det+0x3b6>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    1417 <LUDecompose_det+0x2e7>
jmp    14fc <LUDecompose_det+0x3cc>
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    1249 <LUDecompose_det+0x119>
movq   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    1556 <LUDecompose_det+0x426>
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x68(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x68(%rbp)
jmp    1515 <LUDecompose_det+0x3e5>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
movsd  -0x30(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<LUDecompose_det>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %esi,%r15d
mov    %rdi,%r12
mov    %esi,%edi
call   1040 <dvector@plt>
test   %r15d,%r15d
je     174f <LUDecompose_det+0x61f>
movslq %r15d,%rbp
cmp    $0x2,%rbp
mov    $0x1,%esi
cmovae %rbp,%rsi
lea    -0x1(%rsi),%rcx
mov    %rcx,0x40(%rsp)
mov    %esi,%ebx
and    $0x3,%ebx
mov    %rsi,%rdi
and    $0xfffffffffffffffc,%rdi
xor    %r14d,%r14d
lea    0xeba(%rip),%r13
movsd  0xe71(%rip),%xmm10
movapd 0xe88(%rip),%xmm11
movapd 0xe80(%rip),%xmm5
mov    %rax,0x48(%rsp)
mov    %rsi,0x8(%rsp)
mov    %rdi,0x30(%rsp)
jmp    11dc <LUDecompose_det+0xac>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd %xmm10,%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%r14,8)
add    $0x1,%r14
cmp    %rsi,%r14
je     12d1 <LUDecompose_det+0x1a1>
mov    (%r12,%r14,8),%rdx
cmpq   $0x3,0x40(%rsp)
jae    11f0 <LUDecompose_det+0xc0>
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
jmp    1244 <LUDecompose_det+0x114>
xorpd  %xmm1,%xmm1
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rdx,%rcx,8),%xmm2
movsd  0x8(%rdx,%rcx,8),%xmm3
movsd  0x10(%rdx,%rcx,8),%xmm4
movsd  0x18(%rdx,%rcx,8),%xmm0
andpd  %xmm5,%xmm2
maxsd  %xmm1,%xmm2
andpd  %xmm5,%xmm3
maxsd  %xmm2,%xmm3
andpd  %xmm5,%xmm4
maxsd  %xmm3,%xmm4
andpd  %xmm5,%xmm0
maxsd  %xmm4,%xmm0
add    $0x4,%rcx
movapd %xmm0,%xmm1
cmp    %rcx,%rdi
jne    1200 <LUDecompose_det+0xd0>
movapd %xmm0,%xmm1
test   %rbx,%rbx
je     127b <LUDecompose_det+0x14b>
lea    (%rdx,%rcx,8),%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rdx,%rcx,8),%xmm1
andpd  %xmm11,%xmm1
maxsd  %xmm0,%xmm1
add    $0x1,%rcx
movapd %xmm1,%xmm0
cmp    %rcx,%rbx
jne    1260 <LUDecompose_det+0x130>
ucomisd 0xd8d(%rip),%xmm1
jne    11c0 <LUDecompose_det+0x90>
jp     11c0 <LUDecompose_det+0x90>
mov    %r13,%rdi
movapd %xmm1,0x20(%rsp)
call   1050 <error@plt>
movapd 0x20(%rsp),%xmm1
movapd 0xd75(%rip),%xmm5
movapd 0xd6c(%rip),%xmm11
mov    0x30(%rsp),%rdi
mov    0x8(%rsp),%rsi
movsd  0xd39(%rip),%xmm10
mov    0x48(%rsp),%rax
jmp    11c0 <LUDecompose_det+0x90>
test   %r15d,%r15d
je     17b4 <LUDecompose_det+0x684>
mov    %r15d,0x1c(%rsp)
cmp    $0x2,%rbp
mov    $0x1,%r8d
cmovae %rbp,%r8
lea    -0x4(%r8),%rcx
mov    %rcx,0x60(%rsp)
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %rdx,%rcx
mov    %rdx,0x58(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x50(%rsp)
lea    -0x5(%rbp),%rcx
mov    %rcx,0x38(%rsp)
movsd  0xcdc(%rip),%xmm10
mov    $0x7,%ecx
mov    %rcx,0x10(%rsp)
xor    %r15d,%r15d
xorpd  %xmm9,%xmm9
movsd  0xcc9(%rip),%xmm8
movapd 0xce9(%rip),%xmm2
mov    %r8,0x68(%rsp)
jmp    136f <LUDecompose_det+0x23f>
xchg   %ax,%ax
addq   $0x1,0x10(%rsp)
addq   $0xffffffffffffffff,0x38(%rsp)
mov    0x68(%rsp),%r8
cmp    %r8,%r15
mov    0x48(%rsp),%rax
je     1736 <LUDecompose_det+0x606>
mov    %r15,%r11
add    $0x1,%r15
mov    %r11,%rdx
cmp    %rbp,%r15
jae    13bd <LUDecompose_det+0x28d>
xorpd  %xmm3,%xmm3
mov    %r11,%rdx
mov    %r15,%rdi
nopl   0x0(%rax,%rax,1)
mov    (%r12,%rdi,8),%rsi
movsd  (%rsi,%r11,8),%xmm4
andpd  %xmm11,%xmm4
mulsd  (%rax,%rdi,8),%xmm4
ucomisd %xmm3,%xmm4
cmova  %rdi,%rdx
maxsd  %xmm3,%xmm4
add    $0x1,%rdi
movapd %xmm4,%xmm3
cmp    %rbp,%rdi
jb     1390 <LUDecompose_det+0x260>
cmp    %rdx,%r11
jne    13d0 <LUDecompose_det+0x2a0>
mov    (%r12,%r11,8),%rdi
jmp    1540 <LUDecompose_det+0x410>
nopl   0x0(%rax,%rax,1)
mov    (%r12,%rdx,8),%rcx
mov    (%r12,%r11,8),%rdi
mov    0x8(%rsp),%rbx
cmp    $0x4,%rbx
jb     13f5 <LUDecompose_det+0x2c5>
lea    (%rdi,%rbx,8),%rsi
cmp    %rsi,%rcx
jae    146a <LUDecompose_det+0x33a>
lea    (%rcx,%rbx,8),%rsi
cmp    %rsi,%rdi
jae    146a <LUDecompose_det+0x33a>
xor    %esi,%esi
mov    %rsi,%rbx
or     $0x1,%rbx
test   $0x1,%r8b
je     141b <LUDecompose_det+0x2eb>
movsd  (%rcx,%rsi,8),%xmm3
movsd  (%rdi,%rsi,8),%xmm4
movsd  %xmm4,(%rcx,%rsi,8)
movsd  %xmm3,(%rdi,%rsi,8)
mov    %rbx,%rsi
cmp    %rbx,%r8
je     1530 <LUDecompose_det+0x400>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rcx,%rsi,8),%xmm3
movsd  (%rdi,%rsi,8),%xmm4
movsd  %xmm4,(%rcx,%rsi,8)
movsd  %xmm3,(%rdi,%rsi,8)
movsd  0x8(%rcx,%rsi,8),%xmm3
movsd  0x8(%rdi,%rsi,8),%xmm4
movsd  %xmm4,0x8(%rcx,%rsi,8)
movsd  %xmm3,0x8(%rdi,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%r8
jne    1430 <LUDecompose_det+0x300>
jmp    1530 <LUDecompose_det+0x400>
cmpq   $0x4,0x60(%rsp)
jae    1476 <LUDecompose_det+0x346>
xor    %esi,%esi
jmp    14de <LUDecompose_det+0x3ae>
mov    0x50(%rsp),%rbx
xor    %esi,%esi
nopl   (%rax)
movups (%rcx,%rsi,8),%xmm3
movups 0x10(%rcx,%rsi,8),%xmm4
movups (%rdi,%rsi,8),%xmm5
movups 0x10(%rdi,%rsi,8),%xmm6
movups %xmm5,(%rcx,%rsi,8)
movups %xmm6,0x10(%rcx,%rsi,8)
movups %xmm3,(%rdi,%rsi,8)
movups %xmm4,0x10(%rdi,%rsi,8)
movupd 0x20(%rcx,%rsi,8),%xmm3
movupd 0x30(%rcx,%rsi,8),%xmm4
movupd 0x20(%rdi,%rsi,8),%xmm5
movupd 0x30(%rdi,%rsi,8),%xmm6
movupd %xmm5,0x20(%rcx,%rsi,8)
movupd %xmm6,0x30(%rcx,%rsi,8)
movupd %xmm3,0x20(%rdi,%rsi,8)
movupd %xmm4,0x30(%rdi,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rbx
jne    1480 <LUDecompose_det+0x350>
testb  $0x1,0x58(%rsp)
je     1511 <LUDecompose_det+0x3e1>
movupd (%rcx,%rsi,8),%xmm3
movupd 0x10(%rcx,%rsi,8),%xmm4
movupd (%rdi,%rsi,8),%xmm5
movupd 0x10(%rdi,%rsi,8),%xmm6
movupd %xmm5,(%rcx,%rsi,8)
movupd %xmm6,0x10(%rcx,%rsi,8)
movupd %xmm3,(%rdi,%rsi,8)
movupd %xmm4,0x10(%rdi,%rsi,8)
mov    0x30(%rsp),%rbx
mov    %rbx,%rsi
cmp    %rbx,0x8(%rsp)
jne    13f7 <LUDecompose_det+0x2c7>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xorpd  %xmm2,%xmm10
movsd  (%rax,%r11,8),%xmm3
movsd  %xmm3,(%rax,%rdx,8)
movsd  (%rdi,%r11,8),%xmm4
ucomisd %xmm9,%xmm4
jne    1562 <LUDecompose_det+0x432>
jp     1562 <LUDecompose_det+0x432>
movabs $0x3abef2d0f5da7dd9,%rax
f2 be 3a
mov    %rax,(%rdi,%r11,8)
movapd %xmm8,%xmm4
cmp    %rbp,%r15
jae    1350 <LUDecompose_det+0x220>
mov    0x38(%rsp),%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %r11,%r9
not    %r9
lea    (%r9,%rbp,1),%rax
add    $0xfffffffffffffffc,%rax
mov    %rax,0x20(%rsp)
shr    $0x2,%rax
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x70(%rsp)
add    %rbp,%r9
add    $0x1,%rax
mov    %rax,0x80(%rsp)
lea    (%rdi,%r15,8),%rbx
lea    (%rdi,%rbp,8),%r13
mov    %r9,%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,%rcx
add    %r15,%rax
mov    %rax,0x78(%rsp)
mov    0x10(%rsp),%rax
lea    (%rdi,%rax,8),%rsi
mov    %r15,%r8
nopl   (%rax)
mov    (%r12,%r8,8),%r14
movsd  (%r14,%r11,8),%xmm3
divsd  %xmm4,%xmm3
movsd  %xmm3,(%r14,%r11,8)
xorpd  %xmm2,%xmm3
mov    %r15,%rax
cmp    $0x4,%r9
jb     1700 <LUDecompose_det+0x5d0>
lea    (%r14,%r15,8),%rax
cmp    %r13,%rax
jae    160e <LUDecompose_det+0x4de>
lea    (%r14,%rbp,8),%rdx
mov    %r15,%rax
cmp    %rdx,%rbx
jb     1700 <LUDecompose_det+0x5d0>
movapd %xmm3,%xmm4
unpcklpd %xmm3,%xmm4
cmpq   $0x4,0x20(%rsp)
jae    1625 <LUDecompose_det+0x4f5>
xor    %edx,%edx
jmp    16af <LUDecompose_det+0x57f>
mov    0x10(%rsp),%rax
lea    (%r14,%rax,8),%rax
mov    0x70(%rsp),%r10
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movupd -0x30(%rsi,%rdx,8),%xmm5
movupd -0x20(%rsi,%rdx,8),%xmm6
mulpd  %xmm4,%xmm5
movupd -0x30(%rax,%rdx,8),%xmm7
addpd  %xmm5,%xmm7
movupd -0x20(%rax,%rdx,8),%xmm5
movupd -0x10(%rax,%rdx,8),%xmm1
movupd (%rax,%rdx,8),%xmm0
mulpd  %xmm4,%xmm6
addpd  %xmm5,%xmm6
movupd %xmm7,-0x30(%rax,%rdx,8)
movupd %xmm6,-0x20(%rax,%rdx,8)
movupd -0x10(%rsi,%rdx,8),%xmm5
movupd (%rsi,%rdx,8),%xmm6
mulpd  %xmm4,%xmm5
addpd  %xmm1,%xmm5
mulpd  %xmm4,%xmm6
addpd  %xmm0,%xmm6
movupd %xmm5,-0x10(%rax,%rdx,8)
movupd %xmm6,(%rax,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%r10
jne    1640 <LUDecompose_det+0x510>
testb  $0x1,0x80(%rsp)
01
je     16f1 <LUDecompose_det+0x5c1>
add    %r15,%rdx
movupd (%rdi,%rdx,8),%xmm0
movupd 0x10(%rdi,%rdx,8),%xmm1
movupd (%r14,%rdx,8),%xmm5
movupd 0x10(%r14,%rdx,8),%xmm6
mulpd  %xmm4,%xmm0
addpd  %xmm5,%xmm0
mulpd  %xmm4,%xmm1
addpd  %xmm6,%xmm1
movupd %xmm0,(%r14,%rdx,8)
movupd %xmm1,0x10(%r14,%rdx,8)
mov    0x78(%rsp),%rax
cmp    %rcx,%r9
je     171e <LUDecompose_det+0x5ee>
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm0
mulsd  %xmm3,%xmm0
addsd  (%r14,%rax,8),%xmm0
movsd  %xmm0,(%r14,%rax,8)
add    $0x1,%rax
cmp    %rbp,%rax
jb     1700 <LUDecompose_det+0x5d0>
add    $0x1,%r8
cmp    %rbp,%r8
jae    1350 <LUDecompose_det+0x220>
movsd  (%rdi,%r11,8),%xmm4
jmp    15d0 <LUDecompose_det+0x4a0>
cmpl   $0x0,0x1c(%rsp)
je     17b4 <LUDecompose_det+0x684>
mov    %r8d,%edi
and    $0x3,%edi
cmpq   $0x3,0x40(%rsp)
jae    175a <LUDecompose_det+0x62a>
xor    %ecx,%ecx
jmp    1797 <LUDecompose_det+0x667>
movsd  0x8a8(%rip),%xmm10
jmp    17b4 <LUDecompose_det+0x684>
and    $0xfffffffffffffffc,%r8
xor    %ecx,%ecx
mov    (%r12,%rcx,8),%rdx
mov    0x8(%r12,%rcx,8),%rsi
mulsd  (%rdx,%rcx,8),%xmm10
mulsd  0x8(%rsi,%rcx,8),%xmm10
mov    0x10(%r12,%rcx,8),%rdx
mulsd  0x10(%rdx,%rcx,8),%xmm10
mov    0x18(%r12,%rcx,8),%rdx
mulsd  0x18(%rdx,%rcx,8),%xmm10
add    $0x4,%rcx
cmp    %rcx,%r8
jne    1760 <LUDecompose_det+0x630>
test   %rdi,%rdi
je     17b4 <LUDecompose_det+0x684>
nopl   0x0(%rax)
mov    (%r12,%rcx,8),%rdx
mulsd  (%rdx,%rcx,8),%xmm10
add    $0x1,%rcx
add    $0xffffffffffffffff,%rdi
jne    17a0 <LUDecompose_det+0x670>
movapd %xmm10,0x20(%rsp)
mov    %rax,%rdi
call   1030 <free@plt>
movaps 0x20(%rsp),%xmm0
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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