#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_det2(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 * dim + 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 * dim + k]);
      if (temp > big) {
        big = temp;
        imax = i;
      }
    }
    if (k != imax) {
      for (j = 0; j < dim; j++) {
        temp = m[imax * dim + j];
        m[imax * dim + j] = m[k * dim + j];
        m[k * dim + j] = temp;
      }
      d = -d;
      vv[imax] = vv[k];
    }
    if (m[k * dim + k] == 0.0)
      m[k * dim + k] = 1.0e-25;
    for (i = k + 1; i < dim; i++) {
      temp = m[i * dim + k] /= m[k * dim + k];
      for (j = k + 1; j < dim; j++) {
        m[i * dim + j] -= temp * m[k * dim + j];
      }
    }
  }
  for (size_t i = 0; i < dim; i++) {
    d *= m[i * dim + i];
  }
  free(vv);
  return d;
}


x86 assembly (clang, -O0):
<LUDecompose_det2>:
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    1245 <LUDecompose_det2+0x115>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    11f6 <LUDecompose_det2+0xc6>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe4e(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
ucomisd -0x18(%rbp),%xmm0
jbe    11e0 <LUDecompose_det2+0xb0>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    11e5 <LUDecompose_det2+0xb5>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1191 <LUDecompose_det2+0x61>
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    121a <LUDecompose_det2+0xea>
jp     121a <LUDecompose_det2+0xea>
lea    0xe0b(%rip),%rdi
call   1050 <error@plt>
movsd  0xdde(%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_det2+0x40>
movq   $0x0,-0x60(%rbp)
mov    -0x60(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    153d <LUDecompose_det2+0x40d>
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    12fd <LUDecompose_det2+0x1cd>
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xd58(%rip),%xmm2
pand   %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    12e7 <LUDecompose_det2+0x1b7>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x48(%rbp),%rax
mov    %rax,-0x58(%rbp)
jmp    12ec <LUDecompose_det2+0x1bc>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    127a <LUDecompose_det2+0x14a>
mov    -0x60(%rbp),%rax
cmp    -0x58(%rbp),%rax
je     13de <LUDecompose_det2+0x2ae>
movq   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    13a3 <LUDecompose_det2+0x273>
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    1313 <LUDecompose_det2+0x1e3>
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
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    142b <LUDecompose_det2+0x2fb>
jp     142b <LUDecompose_det2+0x2fb>
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x60(%rbp),%rcx
movsd  0xbe2(%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    1527 <LUDecompose_det2+0x3f7>
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -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    1511 <LUDecompose_det2+0x3e1>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -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    1494 <LUDecompose_det2+0x364>
jmp    1516 <LUDecompose_det2+0x3e6>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    1437 <LUDecompose_det2+0x307>
jmp    152c <LUDecompose_det2+0x3fc>
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    124d <LUDecompose_det2+0x11d>
movq   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%rax
movslq -0xc(%rbp),%rcx
cmp    %rcx,%rax
jae    158a <LUDecompose_det2+0x45a>
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%rcx
movslq -0xc(%rbp),%rdx
imul   %rdx,%rcx
add    -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    1545 <LUDecompose_det2+0x415>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
movsd  -0x30(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<LUDecompose_det2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
mov    %esi,%ebx
mov    %rdi,%r12
mov    %esi,%edi
call   1040 <dvector@plt>
test   %ebx,%ebx
je     182f <LUDecompose_det2+0x6ff>
mov    %ebx,0xc(%rsp)
movslq %ebx,%rbp
cmp    $0x2,%rbp
mov    $0x1,%r8d
cmovae %rbp,%r8
lea    -0x1(%r8),%rcx
mov    %rcx,0x48(%rsp)
mov    %r8,%r15
and    $0xfffffffffffffffe,%r15
lea    0x8(%r12),%rbx
lea    0x0(,%rbp,8),%r13
xor    %r14d,%r14d
movsd  0xe6e(%rip),%xmm2
movapd 0xe85(%rip),%xmm10
movapd 0xe7d(%rip),%xmm3
mov    %rax,0x58(%rsp)
mov    %r8,0x50(%rsp)
jmp    11ce <LUDecompose_det2+0x9e>
nop
movapd %xmm2,%xmm0
divsd  %xmm4,%xmm0
movsd  %xmm0,(%rax,%r14,8)
add    $0x1,%r14
add    %r13,%rbx
cmp    %r8,%r14
je     128e <LUDecompose_det2+0x15e>
cmpq   $0x0,0x48(%rsp)
je     1280 <LUDecompose_det2+0x150>
xorpd  %xmm0,%xmm0
xor    %edx,%edx
movsd  -0x8(%rbx,%rdx,8),%xmm1
movsd  (%rbx,%rdx,8),%xmm4
andpd  %xmm3,%xmm1
maxsd  %xmm0,%xmm1
andpd  %xmm3,%xmm4
maxsd  %xmm1,%xmm4
add    $0x2,%rdx
movapd %xmm4,%xmm0
cmp    %rdx,%r15
jne    11e0 <LUDecompose_det2+0xb0>
test   $0x1,%r8b
je     122b <LUDecompose_det2+0xfb>
mov    %r14,%rcx
imul   %rbp,%rcx
add    %rcx,%rdx
movsd  (%r12,%rdx,8),%xmm0
andpd  %xmm10,%xmm0
maxsd  %xmm4,%xmm0
movapd %xmm0,%xmm4
ucomisd 0xddd(%rip),%xmm4
jne    11b0 <LUDecompose_det2+0x80>
jp     11b0 <LUDecompose_det2+0x80>
lea    0xdfa(%rip),%rdi
movapd %xmm4,0x10(%rsp)
call   1050 <error@plt>
movapd 0x10(%rsp),%xmm4
movapd 0xdc1(%rip),%xmm3
movapd 0xdb8(%rip),%xmm10
mov    0x50(%rsp),%r8
movsd  0xd8b(%rip),%xmm2
mov    0x58(%rsp),%rax
jmp    11b0 <LUDecompose_det2+0x80>
nop
xor    %edx,%edx
xorpd  %xmm4,%xmm4
test   $0x1,%r8b
jne    120e <LUDecompose_det2+0xde>
jmp    122b <LUDecompose_det2+0xfb>
cmpl   $0x0,0xc(%rsp)
je     18bd <LUDecompose_det2+0x78d>
cmp    $0x2,%rbp
mov    $0x1,%r11d
cmovae %rbp,%r11
lea    0x1(%rbp),%rcx
mov    %rcx,0x60(%rsp)
lea    0x8(,%rbp,8),%rcx
mov    %rcx,0x40(%rsp)
lea    (%r12,%r8,8),%rcx
mov    %rcx,0x98(%rsp)
lea    -0x4(%r11),%rcx
mov    %rcx,0x90(%rsp)
mov    %rcx,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %r8,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x80(%rsp)
mov    %rdx,%rcx
mov    %rdx,0x88(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x78(%rsp)
lea    (%r12,%rbp,8),%rcx
lea    0x30(%r12),%rdx
lea    0x8(%r12),%rsi
mov    %rsi,0xa0(%rsp)
lea    -0x1(%rbp),%rsi
mov    %rsi,0x30(%rsp)
lea    0x18(%r12),%rsi
mov    %rsi,0x28(%rsp)
lea    (%r12,%rbp,8),%rsi
add    $0x18,%rsi
mov    %rsi,0x20(%rsp)
movsd  0xcbf(%rip),%xmm2
xor    %ebx,%ebx
xorpd  %xmm9,%xmm9
movsd  0xcb7(%rip),%xmm8
movapd 0xcd7(%rip),%xmm0
mov    %r12,%r10
mov    %rdx,0x70(%rsp)
mov    %rdx,%r14
mov    %rcx,0x38(%rsp)
mov    %rcx,%rdi
mov    %r13,0x68(%rsp)
mov    %r11,0xa8(%rsp)
jmp    13dc <LUDecompose_det2+0x2ac>
nopl   0x0(%rax,%rax,1)
mov    0x40(%rsp),%rsi
mov    0xb0(%rsp),%rdi
add    %rsi,%rdi
mov    0x68(%rsp),%r13
mov    0xb8(%rsp),%r14
add    %r13,%r14
mov    0xc0(%rsp),%r10
add    %r13,%r10
addq   $0xffffffffffffffff,0x30(%rsp)
addq   $0x8,0x28(%rsp)
add    %rsi,0x20(%rsp)
add    %r13,0x38(%rsp)
mov    0xa8(%rsp),%r11
cmp    %r11,%rbx
mov    0x58(%rsp),%rax
mov    0x50(%rsp),%r8
je     1806 <LUDecompose_det2+0x6d6>
mov    %rbx,%rsi
add    $0x1,%rbx
mov    %rsi,%r15
cmp    %rbp,%rbx
jae    142a <LUDecompose_det2+0x2fa>
xorpd  %xmm3,%xmm3
mov    %rdi,%rdx
mov    %rsi,%r15
mov    %rbx,%rcx
nopl   0x0(%rax,%rax,1)
movsd  (%rdx),%xmm4
andpd  %xmm10,%xmm4
mulsd  (%rax,%rcx,8),%xmm4
ucomisd %xmm3,%xmm4
cmova  %rcx,%r15
maxsd  %xmm3,%xmm4
add    $0x1,%rcx
add    %r13,%rdx
movapd %xmm4,%xmm3
cmp    %rbp,%rcx
jb     1400 <LUDecompose_det2+0x2d0>
mov    %rsi,%r13
imul   %rbp,%r13
cmp    %r15,%rsi
mov    %r14,0xb8(%rsp)
mov    %rdi,0xb0(%rsp)
mov    %rsi,0x10(%rsp)
je     1622 <LUDecompose_det2+0x4f2>
mov    %rbx,%r9
mov    %r15,%rbx
imul   %rbp,%rbx
cmp    $0x4,%r8
jb     148d <LUDecompose_det2+0x35d>
lea    (%r8,%r13,1),%rcx
lea    (%r12,%rcx,8),%rcx
lea    (%r12,%rbx,8),%rdx
cmp    %rcx,%rdx
jae    151e <LUDecompose_det2+0x3ee>
lea    (%r12,%r13,8),%rcx
mov    0x98(%rsp),%rdx
lea    (%rdx,%rbx,8),%rdx
cmp    %rdx,%rcx
jae    151e <LUDecompose_det2+0x3ee>
xor    %ecx,%ecx
mov    %rcx,%rdx
or     $0x1,%rdx
test   $0x1,%r11b
je     14c4 <LUDecompose_det2+0x394>
mov    %rsi,%rdi
lea    (%rcx,%rbx,1),%rsi
movsd  (%r12,%rsi,8),%xmm3
add    %r13,%rcx
movsd  (%r12,%rcx,8),%xmm4
movsd  %xmm4,(%r12,%rsi,8)
mov    %rdi,%rsi
movsd  %xmm3,(%r12,%rcx,8)
mov    %rdx,%rcx
cmp    %rdx,%r11
je     1610 <LUDecompose_det2+0x4e0>
mov    0xa0(%rsp),%rdx
lea    (%rdx,%rbx,8),%rdx
nopl   0x0(%rax)
movsd  -0x8(%rdx,%rcx,8),%xmm3
movsd  (%r10,%rcx,8),%xmm4
movsd  %xmm4,-0x8(%rdx,%rcx,8)
movsd  %xmm3,(%r10,%rcx,8)
movsd  (%rdx,%rcx,8),%xmm3
movsd  0x8(%r10,%rcx,8),%xmm4
movsd  %xmm4,(%rdx,%rcx,8)
movsd  %xmm3,0x8(%r10,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r11
jne    14e0 <LUDecompose_det2+0x3b0>
jmp    1610 <LUDecompose_det2+0x4e0>
cmpq   $0x4,0x90(%rsp)
jae    152d <LUDecompose_det2+0x3fd>
xor    %ecx,%ecx
jmp    15a6 <LUDecompose_det2+0x476>
mov    0x70(%rsp),%rcx
lea    (%rcx,%rbx,8),%rsi
mov    0x78(%rsp),%rdi
xor    %ecx,%ecx
nopl   (%rax)
movups -0x30(%rsi,%rcx,8),%xmm3
movups -0x20(%rsi,%rcx,8),%xmm4
movups -0x30(%r14,%rcx,8),%xmm5
movups -0x20(%r14,%rcx,8),%xmm6
movups %xmm5,-0x30(%rsi,%rcx,8)
movups %xmm6,-0x20(%rsi,%rcx,8)
movups %xmm3,-0x30(%r14,%rcx,8)
movups %xmm4,-0x20(%r14,%rcx,8)
movupd -0x10(%rsi,%rcx,8),%xmm3
movupd (%rsi,%rcx,8),%xmm4
movupd -0x10(%r14,%rcx,8),%xmm5
movupd (%r14,%rcx,8),%xmm6
movupd %xmm5,-0x10(%rsi,%rcx,8)
movupd %xmm6,(%rsi,%rcx,8)
movupd %xmm3,-0x10(%r14,%rcx,8)
movupd %xmm4,(%r14,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdi
jne    1540 <LUDecompose_det2+0x410>
testb  $0x1,0x88(%rsp)
01
je     15eb <LUDecompose_det2+0x4bb>
lea    (%rcx,%rbx,1),%rdx
movupd (%r12,%rdx,8),%xmm3
movupd 0x10(%r12,%rdx,8),%xmm4
add    %r13,%rcx
movupd (%r12,%rcx,8),%xmm5
movupd 0x10(%r12,%rcx,8),%xmm6
movupd %xmm5,(%r12,%rdx,8)
movupd %xmm6,0x10(%r12,%rdx,8)
movupd %xmm3,(%r12,%rcx,8)
movupd %xmm4,0x10(%r12,%rcx,8)
mov    0x80(%rsp),%rdx
mov    %rdx,%rcx
cmp    %rdx,%r8
mov    0x10(%rsp),%rsi
jne    148f <LUDecompose_det2+0x35f>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xorpd  %xmm0,%xmm2
movsd  (%rax,%rsi,8),%xmm3
movsd  %xmm3,(%rax,%r15,8)
mov    %r9,%rbx
mov    %r10,0xc0(%rsp)
mov    %r13,%r8
add    %rsi,%r13
movsd  (%r12,%r13,8),%xmm4
ucomisd %xmm9,%xmm4
jne    1652 <LUDecompose_det2+0x522>
jp     1652 <LUDecompose_det2+0x522>
movabs $0x3abef2d0f5da7dd9,%rax
f2 be 3a
mov    %rax,(%r12,%r13,8)
movapd %xmm8,%xmm4
cmp    %rbp,%rbx
jae    1380 <LUDecompose_det2+0x250>
mov    0x10(%rsp),%rcx
mov    %rcx,%rdi
not    %rdi
add    %rbp,%rdi
lea    0x2(%rcx),%rdx
mov    0x30(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    0x40(%rsp),%rsi
imul   %rcx,%rsi
mov    %rsi,0xf0(%rsp)
imul   %rbp,%rdx
mov    %rdx,0xe0(%rsp)
lea    (%r12,%rdi,8),%rdx
lea    (%r8,%rbx,1),%rsi
lea    (%r12,%rsi,8),%rcx
mov    %rcx,0xd8(%rsp)
lea    (%rdx,%rsi,8),%rcx
mov    %rcx,0xe8(%rsp)
mov    %rdi,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0xd0(%rsp)
add    %rbx,%rcx
mov    %rcx,0xc8(%rsp)
mov    0x28(%rsp),%rdx
lea    (%rdx,%r8,8),%r9
lea    (%r12,%r8,8),%r10
mov    0x38(%rsp),%r15
mov    0x20(%rsp),%rsi
xor    %r14d,%r14d
mov    %rbx,%r11
nopw   0x0(%rax,%rax,1)
mov    %r11,%rdx
imul   %rbp,%rdx
add    0x10(%rsp),%rdx
movsd  (%r12,%rdx,8),%xmm3
divsd  %xmm4,%xmm3
movsd  %xmm3,(%r12,%rdx,8)
xorpd  %xmm0,%xmm3
mov    %rbx,%r8
cmp    $0x4,%rdi
jb     17c0 <LUDecompose_det2+0x690>
mov    %r14,%rdx
imul   %rbp,%rdx
mov    %rbx,%rcx
mov    0x60(%rsp),%rbx
add    %rdx,%rbx
lea    (%r12,%rbx,8),%rbx
add    0xf0(%rsp),%rbx
cmp    0xe8(%rsp),%rbx
mov    %rcx,%rbx
jae    1761 <LUDecompose_det2+0x631>
add    0xe0(%rsp),%rdx
lea    (%r12,%rdx,8),%rdx
mov    %rbx,%r8
cmp    %rdx,0xd8(%rsp)
jb     17c0 <LUDecompose_det2+0x690>
movapd %xmm3,%xmm4
unpcklpd %xmm3,%xmm4
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movupd -0x10(%r9,%rdx,8),%xmm5
movupd (%r9,%rdx,8),%xmm6
movupd -0x10(%rsi,%rdx,8),%xmm7
movupd (%rsi,%rdx,8),%xmm1
mulpd  %xmm4,%xmm5
addpd  %xmm7,%xmm5
mulpd  %xmm4,%xmm6
addpd  %xmm1,%xmm6
movupd %xmm5,-0x10(%rsi,%rdx,8)
movupd %xmm6,(%rsi,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rax
jne    1770 <LUDecompose_det2+0x640>
mov    0xc8(%rsp),%r8
cmp    0xd0(%rsp),%rdi
je     17df <LUDecompose_det2+0x6af>
xchg   %ax,%ax
movsd  (%r10,%r8,8),%xmm1
mulsd  %xmm3,%xmm1
addsd  (%r15,%r8,8),%xmm1
movsd  %xmm1,(%r15,%r8,8)
add    $0x1,%r8
cmp    %rbp,%r8
jb     17c0 <LUDecompose_det2+0x690>
add    $0x1,%r11
cmp    %rbp,%r11
jae    1380 <LUDecompose_det2+0x250>
movsd  (%r12,%r13,8),%xmm4
add    $0x1,%r14
mov    0x68(%rsp),%rcx
add    %rcx,%rsi
add    %rcx,%r15
jmp    16f0 <LUDecompose_det2+0x5c0>
cmpl   $0x0,0xc(%rsp)
je     18bd <LUDecompose_det2+0x78d>
mov    %r11d,%edi
and    $0x3,%edi
cmpq   $0x3,0x48(%rsp)
jae    183c <LUDecompose_det2+0x70c>
xor    %ecx,%ecx
test   %rdi,%rdi
jne    18a6 <LUDecompose_det2+0x776>
jmp    18bd <LUDecompose_det2+0x78d>
movsd  0x7c9(%rip),%xmm2
jmp    18bd <LUDecompose_det2+0x78d>
mov    %rbp,%r8
shl    $0x4,%r8
mov    %rbp,%r10
shl    $0x5,%r10
add    $0x20,%r10
lea    0x0(,%r13,2),%r9
add    %r13,%r9
and    $0xfffffffffffffffc,%r11
neg    %r11
mov    $0x1,%ecx
mov    %r12,%rbx
nopl   0x0(%rax,%rax,1)
mulsd  (%rbx),%xmm2
mulsd  0x8(%rbx,%rbp,8),%xmm2
mulsd  0x10(%rbx,%r8,1),%xmm2
mulsd  0x18(%rbx,%r9,1),%xmm2
add    %r10,%rbx
lea    (%r11,%rcx,1),%rdx
add    $0x4,%rdx
add    $0x4,%rcx
cmp    $0x1,%rdx
jne    1870 <LUDecompose_det2+0x740>
add    $0xffffffffffffffff,%rcx
test   %rdi,%rdi
je     18bd <LUDecompose_det2+0x78d>
imul   0x60(%rsp),%rcx
lea    (%r12,%rcx,8),%rcx
mulsd  (%rcx),%xmm2
add    %rsi,%rcx
add    $0xffffffffffffffff,%rdi
jne    18b0 <LUDecompose_det2+0x780>
movapd %xmm2,0x10(%rsp)
mov    %rax,%rdi
call   1030 <free@plt>
movaps 0x10(%rsp),%xmm0
add    $0xf8,%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_det2(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 = %59, %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 %62
23: ; preds = %18
store double 0.000000e+00, double* %5
store i64 0, i64* %10
br label %24
24: ; preds = %45, %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 %48
29: ; preds = %24
%30 = load double*, double** %3
%31 = load i64, i64* %9
%32 = load i32, i32* %4
%33 = sext i32 %32 to i64
%34 = mul i64 %31, %33
%35 = load i64, i64* %10
%36 = add i64 %34, %35
%37 = getelementptr inbounds double, double* %30, i64 %36
%38 = load double, double* %37
%39 = call double @llvm.fabs.f64(double %38)
store double %39, double* %6
%40 = load double, double* %5
%41 = fcmp ogt double %39, %40
br i1 %41, label %42, label %44
42: ; preds = %29
%43 = load double, double* %6
store double %43, double* %5
br label %44
44: ; preds = %42, %29
br label %45
45: ; preds = %44
%46 = load i64, i64* %10
%47 = add i64 %46, 1
store i64 %47, i64* %10
br label %24
48: ; preds = %24
%49 = load double, double* %5
%50 = fcmp oeq double %49, 0.000000e+00
br i1 %50, label %51, label %53
51: ; preds = %48
%52 = call i32 @error(i8* getelementptr inbounds ([49 x i8], [49 x i8]* @.str, i64 0, i64 0))
br label %53
53: ; preds = %51, %48
%54 = load double, double* %5
%55 = fdiv double 1.000000e+00, %54
%56 = load double*, double** %7
%57 = load i64, i64* %9
%58 = getelementptr inbounds double, double* %56, i64 %57
store double %55, double* %58
br label %59
59: ; preds = %53
%60 = load i64, i64* %9
%61 = add i64 %60, 1
store i64 %61, i64* %9
br label %18
62: ; preds = %18
store i64 0, i64* %14
br label %63
63: ; preds = %247, %62
%64 = load i64, i64* %14
%65 = load i32, i32* %4
%66 = sext i32 %65 to i64
%67 = icmp ult i64 %64, %66
br i1 %67, label %68, label %250
68: ; preds = %63
store double 0.000000e+00, double* %5
%69 = load i64, i64* %14
store i64 %69, i64* %13
%70 = load i64, i64* %14
%71 = add i64 %70, 1
store i64 %71, i64* %11
br label %72
72: ; preds = %100, %68
%73 = load i64, i64* %11
%74 = load i32, i32* %4
%75 = sext i32 %74 to i64
%76 = icmp ult i64 %73, %75
br i1 %76, label %77, label %103
77: ; preds = %72
%78 = load double*, double** %7
%79 = load i64, i64* %11
%80 = getelementptr inbounds double, double* %78, i64 %79
%81 = load double, double* %80
%82 = load double*, double** %3
%83 = load i64, i64* %11
%84 = load i32, i32* %4
%85 = sext i32 %84 to i64
%86 = mul i64 %83, %85
%87 = load i64, i64* %14
%88 = add i64 %86, %87
%89 = getelementptr inbounds double, double* %82, i64 %88
%90 = load double, double* %89
%91 = call double @llvm.fabs.f64(double %90)
%92 = fmul double %81, %91
store double %92, double* %6
%93 = load double, double* %6
%94 = load double, double* %5
%95 = fcmp ogt double %93, %94
br i1 %95, label %96, label %99
96: ; preds = %77
%97 = load double, double* %6
store double %97, double* %5
%98 = load i64, i64* %11
store i64 %98, i64* %13
br label %99
99: ; preds = %96, %77
br label %100
100: ; preds = %99
%101 = load i64, i64* %11
%102 = add i64 %101, 1
store i64 %102, i64* %11
br label %72
103: ; preds = %72
%104 = load i64, i64* %14
%105 = load i64, i64* %13
%106 = icmp ne i64 %104, %105
br i1 %106, label %107, label %162
107: ; preds = %103
store i64 0, i64* %12
br label %108
108: ; preds = %149, %107
%109 = load i64, i64* %12
%110 = load i32, i32* %4
%111 = sext i32 %110 to i64
%112 = icmp ult i64 %109, %111
br i1 %112, label %113, label %152
113: ; preds = %108
%114 = load double*, double** %3
%115 = load i64, i64* %13
%116 = load i32, i32* %4
%117 = sext i32 %116 to i64
%118 = mul i64 %115, %117
%119 = load i64, i64* %12
%120 = add i64 %118, %119
%121 = getelementptr inbounds double, double* %114, i64 %120
%122 = load double, double* %121
store double %122, double* %6
%123 = load double*, double** %3
%124 = load i64, i64* %14
%125 = load i32, i32* %4
%126 = sext i32 %125 to i64
%127 = mul i64 %124, %126
%128 = load i64, i64* %12
%129 = add i64 %127, %128
%130 = getelementptr inbounds double, double* %123, i64 %129
%131 = load double, double* %130
%132 = load double*, double** %3
%133 = load i64, i64* %13
%134 = load i32, i32* %4
%135 = sext i32 %134 to i64
%136 = mul i64 %133, %135
%137 = load i64, i64* %12
%138 = add i64 %136, %137
%139 = getelementptr inbounds double, double* %132, i64 %138
store double %131, double* %139
%140 = load double, double* %6
%141 = load double*, double** %3
%142 = load i64, i64* %14
%143 = load i32, i32* %4
%144 = sext i32 %143 to i64
%145 = mul i64 %142, %144
%146 = load i64, i64* %12
%147 = add i64 %145, %146
%148 = getelementptr inbounds double, double* %141, i64 %147
store double %140, double* %148
br label %149
149: ; preds = %113
%150 = load i64, i64* %12
%151 = add i64 %150, 1
store i64 %151, i64* %12
br label %108
152: ; preds = %108
%153 = load double, double* %8
%154 = fneg double %153
store double %154, double* %8
%155 = load double*, double** %7
%156 = load i64, i64* %14
%157 = getelementptr inbounds double, double* %155, i64 %156
%158 = load double, double* %157
%159 = load double*, double** %7
%160 = load i64, i64* %13
%161 = getelementptr inbounds double, double* %159, i64 %160
store double %158, double* %161
br label %162
162: ; preds = %152, %103
%163 = load double*, double** %3
%164 = load i64, i64* %14
%165 = load i32, i32* %4
%166 = sext i32 %165 to i64
%167 = mul i64 %164, %166
%168 = load i64, i64* %14
%169 = add i64 %167, %168
%170 = getelementptr inbounds double, double* %163, i64 %169
%171 = load double, double* %170
%172 = fcmp oeq double %171, 0.000000e+00
br i1 %172, label %173, label %182
173: ; preds = %162
%174 = load double*, double** %3
%175 = load i64, i64* %14
%176 = load i32, i32* %4
%177 = sext i32 %176 to i64
%178 = mul i64 %175, %177
%179 = load i64, i64* %14
%180 = add i64 %178, %179
%181 = getelementptr inbounds double, double* %174, i64 %180
store double 1.000000e-25, double* %181
br label %182
182: ; preds = %173, %162
%183 = load i64, i64* %14
%184 = add i64 %183, 1
store i64 %184, i64* %11
br label %185
185: ; preds = %243, %182
%186 = load i64, i64* %11
%187 = load i32, i32* %4
%188 = sext i32 %187 to i64
%189 = icmp ult i64 %186, %188
br i1 %189, label %190, label %246
190: ; preds = %185
%191 = load double*, double** %3
%192 = load i64, i64* %14
%193 = load i32, i32* %4
%194 = sext i32 %193 to i64
%195 = mul i64 %192, %194
%196 = load i64, i64* %14
%197 = add i64 %195, %196
%198 = getelementptr inbounds double, double* %191, i64 %197
%199 = load double, double* %198
%200 = load double*, double** %3
%201 = load i64, i64* %11
%202 = load i32, i32* %4
%203 = sext i32 %202 to i64
%204 = mul i64 %201, %203
%205 = load i64, i64* %14
%206 = add i64 %204, %205
%207 = getelementptr inbounds double, double* %200, i64 %206
%208 = load double, double* %207
%209 = fdiv double %208, %199
store double %209, double* %207
store double %209, double* %6
%210 = load i64, i64* %14
%211 = add i64 %210, 1
store i64 %211, i64* %12
br label %212
212: ; preds = %239, %190
%213 = load i64, i64* %12
%214 = load i32, i32* %4
%215 = sext i32 %214 to i64
%216 = icmp ult i64 %213, %215
br i1 %216, label %217, label %242
217: ; preds = %212
%218 = load double, double* %6
%219 = load double*, double** %3
%220 = load i64, i64* %14
%221 = load i32, i32* %4
%222 = sext i32 %221 to i64
%223 = mul i64 %220, %222
%224 = load i64, i64* %12
%225 = add i64 %223, %224
%226 = getelementptr inbounds double, double* %219, i64 %225
%227 = load double, double* %226
%228 = load double*, double** %3
%229 = load i64, i64* %11
%230 = load i32, i32* %4
%231 = sext i32 %230 to i64
%232 = mul i64 %229, %231
%233 = load i64, i64* %12
%234 = add i64 %232, %233
%235 = getelementptr inbounds double, double* %228, i64 %234
%236 = load double, double* %235
%237 = fneg double %218
%238 = call double @llvm.fmuladd.f64(double %237, double %227, double %236)
store double %238, double* %235
br label %239
239: ; preds = %217
%240 = load i64, i64* %12
%241 = add i64 %240, 1
store i64 %241, i64* %12
br label %212
242: ; preds = %212
br label %243
243: ; preds = %242
%244 = load i64, i64* %11
%245 = add i64 %244, 1
store i64 %245, i64* %11
br label %185
246: ; preds = %185
br label %247
247: ; preds = %246
%248 = load i64, i64* %14
%249 = add i64 %248, 1
store i64 %249, i64* %14
br label %63
250: ; preds = %63
store i64 0, i64* %15
br label %251
251: ; preds = %268, %250
%252 = load i64, i64* %15
%253 = load i32, i32* %4
%254 = sext i32 %253 to i64
%255 = icmp ult i64 %252, %254
br i1 %255, label %256, label %271
256: ; preds = %251
%257 = load double*, double** %3
%258 = load i64, i64* %15
%259 = load i32, i32* %4
%260 = sext i32 %259 to i64
%261 = mul i64 %258, %260
%262 = load i64, i64* %15
%263 = add i64 %261, %262
%264 = getelementptr inbounds double, double* %257, i64 %263
%265 = load double, double* %264
%266 = load double, double* %8
%267 = fmul double %266, %265
store double %267, double* %8
br label %268
268: ; preds = %256
%269 = load i64, i64* %15
%270 = add i64 %269, 1
store i64 %270, i64* %15
br label %251
271: ; preds = %251
%272 = load double*, double** %7
%273 = call i32 @free(double* %272)
%274 = load double, double* %8
ret double %274
}

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