#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 */
float fabs(double);
int matLUdcmp(double *aa, int *indx, double *dd, double *r) {
  double a[4][4];
  int i, imax, j, k;
  double big, dum, sum, temp;
  float vv[4];
  (*dd) = 1.0;
  for (i = 0; i < 4; i++)
    for (j = 0; j < 4; j++)
      a[i][j] = aa[i * 4 + j];
  for (i = 0; i < 4; i++) {
    big = 0.0;
    for (j = 0; j < 4; j++)
      if ((temp = fabs(a[i][j])) > big)
        big = temp;
    if (big == 0.0)
      return -1;
    vv[i] = 1.0 / big;
  }
  for (j = 0; j < 4; 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 < 4; 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 < 4; k++) {
        dum = a[imax][k];
        a[imax][k] = a[j][k];
        a[j][k] = dum;
      }
      (*dd) = -(*dd);
      vv[imax] = vv[j];
    }
    indx[j] = imax;
    if (a[j][j] == 0.0)
      a[j][j] = 1e-20;
    if (j < 3) {
      dum = 1.0 / a[j][j];
      for (i = j + 1; i < 4; i++)
        a[i][j] *= dum;
    }
  }
  for (i = 0; i < 4; i++)
    for (j = 0; j < 4; j++)
      r[i * 4 + j] = a[i][j];
  return 0;
}


x86 assembly (clang, -O0):
<matLUdcmp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x100,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    -0x20(%rbp),%rax
movsd  0xec9(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0xb4(%rbp)
cmpl   $0x4,-0xb4(%rbp)
jge    11d2 <matLUdcmp+0xc2>
movl   $0x0,-0xbc(%rbp)
cmpl   $0x4,-0xbc(%rbp)
jge    11b9 <matLUdcmp+0xa9>
mov    -0x10(%rbp),%rax
mov    -0xb4(%rbp),%ecx
shl    $0x2,%ecx
add    -0xbc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xbc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xbc(%rbp)
jmp    115c <matLUdcmp+0x4c>
jmp    11be <matLUdcmp+0xae>
mov    -0xb4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
jmp    1145 <matLUdcmp+0x35>
movl   $0x0,-0xb4(%rbp)
cmpl   $0x4,-0xb4(%rbp)
jge    12d3 <matLUdcmp+0x1c3>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movl   $0x0,-0xbc(%rbp)
cmpl   $0x4,-0xbc(%rbp)
jge    1274 <matLUdcmp+0x164>
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
ucomisd -0xc8(%rbp),%xmm0
ff
jbe    125b <matLUdcmp+0x14b>
movsd  -0xe0(%rbp),%xmm0
ff
movsd  %xmm0,-0xc8(%rbp)
ff
jmp    1260 <matLUdcmp+0x150>
mov    -0xbc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xbc(%rbp)
jmp    11fe <matLUdcmp+0xee>
movsd  -0xc8(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    129b <matLUdcmp+0x18b>
jp     129b <matLUdcmp+0x18b>
movl   $0xffffffff,-0x4(%rbp)
jmp    1888 <matLUdcmp+0x778>
movsd  0xd5d(%rip),%xmm0
divsd  -0xc8(%rbp),%xmm0
ff
cvtsd2ss %xmm0,%xmm0
movslq -0xb4(%rbp),%rax
movss  %xmm0,-0xf0(%rbp,%rax,4)
ff ff
mov    -0xb4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
jmp    11dc <matLUdcmp+0xcc>
movl   $0x0,-0xbc(%rbp)
cmpl   $0x4,-0xbc(%rbp)
jge    17ea <matLUdcmp+0x6da>
movl   $0x0,-0xb4(%rbp)
mov    -0xb4(%rbp),%eax
cmp    -0xbc(%rbp),%eax
jge    140d <matLUdcmp+0x2fd>
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movl   $0x0,-0xc0(%rbp)
mov    -0xc0(%rbp),%eax
cmp    -0xb4(%rbp),%eax
jge    13d0 <matLUdcmp+0x2c0>
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xc0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0xc0(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xd8(%rbp),%xmm1
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
mov    -0xc0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc0(%rbp)
jmp    1339 <matLUdcmp+0x229>
movsd  -0xd8(%rbp),%xmm0
ff
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xb4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
jmp    12f4 <matLUdcmp+0x1e4>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
mov    -0xbc(%rbp),%eax
mov    %eax,-0xb4(%rbp)
cmpl   $0x4,-0xb4(%rbp)
jge    15a7 <matLUdcmp+0x497>
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movl   $0x0,-0xc0(%rbp)
mov    -0xc0(%rbp),%eax
cmp    -0xbc(%rbp),%eax
jge    14fb <matLUdcmp+0x3eb>
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xc0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0xc0(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xd8(%rbp),%xmm1
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
mov    -0xc0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc0(%rbp)
jmp    1464 <matLUdcmp+0x354>
movsd  -0xd8(%rbp),%xmm0
ff
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movslq -0xb4(%rbp),%rax
movss  -0xf0(%rbp,%rax,4),%xmm0
ff ff
movss  %xmm0,-0xf4(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1030 <fabs@plt>
movaps %xmm0,%xmm1
movss  -0xf4(%rbp),%xmm0
ff
mulss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
ucomisd -0xc8(%rbp),%xmm0
ff
jb     158e <matLUdcmp+0x47e>
movsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0xc8(%rbp)
ff
mov    -0xb4(%rbp),%eax
mov    %eax,-0xb8(%rbp)
jmp    1593 <matLUdcmp+0x483>
mov    -0xb4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
jmp    1424 <matLUdcmp+0x314>
mov    -0xbc(%rbp),%eax
cmp    -0xb8(%rbp),%eax
je     16bf <matLUdcmp+0x5af>
movl   $0x0,-0xc0(%rbp)
cmpl   $0x4,-0xc0(%rbp)
jge    1678 <matLUdcmp+0x568>
movslq -0xb8(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xc0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movslq -0xbc(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xc0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0xb8(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xc0(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xd0(%rbp),%xmm0
ff
movslq -0xbc(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xc0(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc0(%rbp)
jmp    15c3 <matLUdcmp+0x4b3>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movslq -0xbc(%rbp),%rax
movss  -0xf0(%rbp,%rax,4),%xmm0
ff ff
movslq -0xb8(%rbp),%rax
movss  %xmm0,-0xf0(%rbp,%rax,4)
ff ff
mov    -0xb8(%rbp),%edx
mov    -0x18(%rbp),%rax
movslq -0xbc(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
movslq -0xbc(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1730 <matLUdcmp+0x620>
jp     1730 <matLUdcmp+0x620>
movslq -0xbc(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  0x8dd(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
cmpl   $0x3,-0xbc(%rbp)
jge    17d1 <matLUdcmp+0x6c1>
movslq -0xbc(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  0x89f(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
mov    -0xbc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
cmpl   $0x4,-0xb4(%rbp)
jge    17cc <matLUdcmp+0x6bc>
movsd  -0xd0(%rbp),%xmm0
ff
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xb4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
jmp    177d <matLUdcmp+0x66d>
jmp    17d1 <matLUdcmp+0x6c1>
jmp    17d6 <matLUdcmp+0x6c6>
mov    -0xbc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xbc(%rbp)
jmp    12dd <matLUdcmp+0x1cd>
movl   $0x0,-0xb4(%rbp)
cmpl   $0x4,-0xb4(%rbp)
jge    1881 <matLUdcmp+0x771>
movl   $0x0,-0xbc(%rbp)
cmpl   $0x4,-0xbc(%rbp)
jge    1868 <matLUdcmp+0x758>
movslq -0xb4(%rbp),%rcx
lea    -0xb0(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0xbc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0xb4(%rbp),%ecx
shl    $0x2,%ecx
add    -0xbc(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xbc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xbc(%rbp)
jmp    180b <matLUdcmp+0x6fb>
jmp    186d <matLUdcmp+0x75d>
mov    -0xb4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xb4(%rbp)
jmp    17f4 <matLUdcmp+0x6e4>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x100,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matLUdcmp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
mov    %rcx,%rbx
mov    %rdx,%r15
mov    %rsi,%r14
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rdx)
movups (%rdi),%xmm0
movups 0x10(%rdi),%xmm1
movupd 0x20(%rdi),%xmm2
movupd 0x30(%rdi),%xmm3
movaps %xmm0,0x30(%rsp)
movups 0x70(%rdi),%xmm0
movaps %xmm0,0xa0(%rsp)
movups 0x60(%rdi),%xmm0
movaps %xmm0,0x90(%rsp)
movups 0x50(%rdi),%xmm0
movaps %xmm0,0x80(%rsp)
movups 0x40(%rdi),%xmm0
movaps %xmm0,0x70(%rsp)
movapd %xmm3,0x60(%rsp)
movapd %xmm2,0x50(%rsp)
movaps %xmm1,0x40(%rsp)
movsd  0x30(%rsp),%xmm0
call   1030 <fabs@plt>
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
maxsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0x38(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x40(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x48(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
mov    $0xffffffff,%ebp
ucomisd 0xe14(%rip),%xmm0
jne    1204 <matLUdcmp+0xf4>
jnp    175d <matLUdcmp+0x64d>
movsd  0xdf4(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0x20(%rsp)
movsd  0x50(%rsp),%xmm0
call   1030 <fabs@plt>
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
maxsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0x58(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x60(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x68(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
ucomisd 0xd87(%rip),%xmm0
jne    1291 <matLUdcmp+0x181>
jnp    175d <matLUdcmp+0x64d>
movsd  0xd67(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0x24(%rsp)
movsd  0x70(%rsp),%xmm0
call   1030 <fabs@plt>
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
maxsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0x78(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x80(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x88(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
ucomisd 0xcf4(%rip),%xmm0
jne    1324 <matLUdcmp+0x214>
jnp    175d <matLUdcmp+0x64d>
movsd  0xcd4(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0x28(%rsp)
movsd  0x90(%rsp),%xmm0
call   1030 <fabs@plt>
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
maxsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0x98(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0xa0(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0xa8(%rsp),%xmm0
call   1030 <fabs@plt>
cvtss2sd %xmm0,%xmm0
maxsd  (%rsp),%xmm0
ucomisd 0xc5b(%rip),%xmm0
jne    13bd <matLUdcmp+0x2ad>
jnp    175d <matLUdcmp+0x64d>
mov    %r15,0xc0(%rsp)
mov    %rbx,0xb8(%rsp)
movsd  0xc2b(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0x2c(%rsp)
lea    0x70(%rsp),%rdx
lea    0x38(%rsp),%r15
lea    0x50(%rsp),%rsi
mov    $0x3,%edi
xor    %ebx,%ebx
mov    %rsi,%r13
xor    %ebp,%ebp
mov    %r14,0xc8(%rsp)
jmp    1436 <matLUdcmp+0x326>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rbp
add    $0x8,%rdx
add    $0xffffffffffffffe0,%rbx
add    $0x20,%r15
add    $0x8,%r13
add    $0xffffffffffffffff,%rdi
add    $0x28,%rsi
cmp    $0x4,%rbp
je     16fb <matLUdcmp+0x5eb>
mov    %rdi,0xd8(%rsp)
mov    %rsi,0xe0(%rsp)
test   %rbp,%rbp
je     14c0 <matLUdcmp+0x3b0>
cmp    $0x1,%rbp
je     14c0 <matLUdcmp+0x3b0>
movsd  0x50(%rsp),%xmm0
mulsd  0x30(%rsp,%rbp,8),%xmm0
movsd  0x50(%rsp,%rbp,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x50(%rsp,%rbp,8)
cmp    $0x2,%rbp
je     14c0 <matLUdcmp+0x3b0>
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdx,%rax,1),%xmm0
movsd  0x70(%rsp,%rax,1),%xmm1
movsd  0x78(%rsp,%rax,1),%xmm2
mulsd  0x30(%rsp,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
mulsd  0x50(%rsp,%rbp,8),%xmm2
subsd  %xmm2,%xmm0
movsd  %xmm0,(%rdx,%rax,1)
lea    (%rbx,%rax,1),%rcx
add    $0x20,%rcx
add    $0x20,%rax
cmp    $0xffffffffffffffc0,%rcx
jne    1480 <matLUdcmp+0x370>
nopl   0x0(%rax)
mov    %rbx,0xd0(%rsp)
mov    %rdx,0xf0(%rsp)
mov    %rbp,%r14
and    $0xfffffffffffffffe,%r14
xorpd  %xmm1,%xmm1
mov    %r15,0xe8(%rsp)
mov    %rbp,%rbx
jmp    154b <matLUdcmp+0x43b>
nopl   0x0(%rax,%rax,1)
lea    (%rax,%rbp,8),%rax
movsd  %xmm0,(%rax)
movss  0x20(%rsp,%rbx,4),%xmm1
movss  %xmm1,0x1c(%rsp)
call   1030 <fabs@plt>
mulss  0x1c(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movapd (%rsp),%xmm3
ucomisd %xmm3,%xmm0
cmovae %ebx,%r12d
movapd %xmm3,%xmm1
cmpnlesd %xmm0,%xmm1
movapd %xmm1,%xmm2
andnpd %xmm0,%xmm2
andpd  %xmm3,%xmm1
orpd   %xmm2,%xmm1
add    $0x1,%rbx
add    $0x20,%r15
cmp    $0x4,%rbx
je     15f0 <matLUdcmp+0x4e0>
mov    %rbx,%rax
shl    $0x5,%rax
add    %rsp,%rax
add    $0x30,%rax
movsd  (%rax,%rbp,8),%xmm0
test   %rbp,%rbp
movapd %xmm1,(%rsp)
je     14f0 <matLUdcmp+0x3e0>
cmp    $0x1,%rbp
jne    1580 <matLUdcmp+0x470>
xor    %ecx,%ecx
jmp    15bb <matLUdcmp+0x4ab>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r13,%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x8(%r15,%rcx,8),%xmm1
mulsd  -0x20(%rdx),%xmm1
movsd  (%r15,%rcx,8),%xmm2
subsd  %xmm1,%xmm0
mulsd  (%rdx),%xmm2
subsd  %xmm2,%xmm0
add    $0x2,%rcx
add    $0x40,%rdx
cmp    %rcx,%r14
jne    1590 <matLUdcmp+0x480>
test   $0x1,%bpl
je     14f0 <matLUdcmp+0x3e0>
movsd  (%rax,%rcx,8),%xmm1
shl    $0x5,%rcx
add    %rsp,%rcx
add    $0x30,%rcx
mulsd  (%rcx,%rbp,8),%xmm1
subsd  %xmm1,%xmm0
jmp    14f0 <matLUdcmp+0x3e0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r12d,%eax
cmp    %rax,%rbp
je     1655 <matLUdcmp+0x545>
movslq %r12d,%rax
mov    %rax,%rcx
shl    $0x5,%rcx
mov    %rbp,%rdx
shl    $0x5,%rdx
movaps 0x30(%rsp,%rcx,1),%xmm0
movaps 0x30(%rsp,%rdx,1),%xmm1
movaps %xmm1,0x30(%rsp,%rcx,1)
movaps %xmm0,0x30(%rsp,%rdx,1)
movaps 0x40(%rsp,%rcx,1),%xmm0
movapd 0x40(%rsp,%rdx,1),%xmm1
movapd %xmm1,0x40(%rsp,%rcx,1)
movaps %xmm0,0x40(%rsp,%rdx,1)
mov    0xc0(%rsp),%rcx
movsd  (%rcx),%xmm0
xorps  0x9da(%rip),%xmm0
movlps %xmm0,(%rcx)
movss  0x20(%rsp,%rbp,4),%xmm0
movss  %xmm0,0x20(%rsp,%rax,4)
mov    0xc8(%rsp),%r14
mov    %r12d,(%r14,%rbp,4)
mov    %rbp,%rax
shl    $0x5,%rax
add    %rsp,%rax
add    $0x30,%rax
movsd  (%rax,%rbp,8),%xmm1
ucomisd 0x994(%rip),%xmm1
movsd  0x97c(%rip),%xmm0
mov    0xf0(%rsp),%rdx
mov    0xe8(%rsp),%r15
mov    0xe0(%rsp),%rsi
mov    0xd8(%rsp),%rdi
mov    0xd0(%rsp),%rbx
jne    16c9 <matLUdcmp+0x5b9>
jp     16c9 <matLUdcmp+0x5b9>
lea    (%rax,%rbp,8),%rax
movabs $0x3bc79ca10c924223,%rcx
9c c7 3b
mov    %rcx,(%rax)
movsd  0x93f(%rip),%xmm1
cmp    $0x2,%rbp
ja     1410 <matLUdcmp+0x300>
divsd  %xmm1,%xmm0
mov    %rsi,%rax
mov    %rdi,%rcx
nopl   (%rax)
movsd  (%rax),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax)
add    $0x20,%rax
add    $0xffffffffffffffff,%rcx
jne    16e0 <matLUdcmp+0x5d0>
jmp    1410 <matLUdcmp+0x300>
movaps 0xa0(%rsp),%xmm0
mov    0xb8(%rsp),%rax
movups %xmm0,0x70(%rax)
movaps 0x90(%rsp),%xmm0
movups %xmm0,0x60(%rax)
movaps 0x80(%rsp),%xmm0
movups %xmm0,0x50(%rax)
movaps 0x70(%rsp),%xmm0
movups %xmm0,0x40(%rax)
movapd 0x30(%rsp),%xmm0
movapd 0x40(%rsp),%xmm1
movapd 0x50(%rsp),%xmm2
movapd 0x60(%rsp),%xmm3
movupd %xmm3,0x30(%rax)
movupd %xmm2,0x20(%rax)
movupd %xmm1,0x10(%rax)
movupd %xmm0,(%rax)
xor    %ebp,%ebp
mov    %ebp,%eax
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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