#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 fabs(double);
void lm_decompose_LUP(double *const A, const int n, const double Tol,
                      int *const P, int *const failure) {
  int i, j, k, imax;
  double maxA, absA, tmp;
  for (i = 0; i < n; i++)
    P[i] = i;
  for (i = 0; i < n; i++) {
    maxA = 0.0;
    imax = i;
    for (k = i; k < n; k++)
      if ((absA = fabs(A[k * n + i])) > maxA) {
        maxA = absA;
        imax = k;
      }
    if (maxA < Tol) {
      *failure = 21;
      return;
    }
    if (imax != i) {
      j = P[i];
      P[i] = P[imax];
      P[imax] = j;
      for (k = 0; k < n; ++k) {
        tmp = A[i * n + k];
        A[i * n + k] = A[imax * n + k];
        A[imax * n + k] = tmp;
      }
    }
    for (j = i + 1; j < n; j++) {
      A[j * n + i] /= A[i * n + i];
      for (k = i + 1; k < n; k++)
        A[j * n + k] -= A[j * n + i] * A[i * n + k];
    }
  }
}


x86 assembly (clang, -O0):
<lm_decompose_LUP>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1147 <lm_decompose_LUP+0x47>
mov    -0x2c(%rbp),%edx
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    111f <lm_decompose_LUP+0x1f>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13b0 <lm_decompose_LUP+0x2b0>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11ce <lm_decompose_LUP+0xce>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe69(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
ucomisd -0x40(%rbp),%xmm0
jbe    11bb <lm_decompose_LUP+0xbb>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    11c0 <lm_decompose_LUP+0xc0>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    116e <lm_decompose_LUP+0x6e>
movsd  -0x40(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11f1 <lm_decompose_LUP+0xf1>
mov    -0x28(%rbp),%rax
movl   $0x15,(%rax)
jmp    13b0 <lm_decompose_LUP+0x2b0>
mov    -0x38(%rbp),%eax
cmp    -0x2c(%rbp),%eax
je     12b7 <lm_decompose_LUP+0x1b7>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x30(%rbp),%edx
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12b2 <lm_decompose_LUP+0x1b2>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1236 <lm_decompose_LUP+0x136>
jmp    12b7 <lm_decompose_LUP+0x1b7>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    139d <lm_decompose_LUP+0x29d>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    138a <lm_decompose_LUP+0x28a>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%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    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    130a <lm_decompose_LUP+0x20a>
jmp    138f <lm_decompose_LUP+0x28f>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12c0 <lm_decompose_LUP+0x1c0>
jmp    13a2 <lm_decompose_LUP+0x2a2>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    114e <lm_decompose_LUP+0x4e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<lm_decompose_LUP>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xc0,%rsp
mov    %rcx,(%rsp)
mov    %rdx,-0x60(%rsp)
mov    %esi,-0x7c(%rsp)
test   %esi,%esi
jle    1935 <lm_decompose_LUP+0x835>
mov    -0x7c(%rsp),%ecx
mov    %ecx,%eax
cmp    $0x8,%ecx
jae    113d <lm_decompose_LUP+0x3d>
xor    %ecx,%ecx
mov    -0x60(%rsp),%rbx
jmp    12a0 <lm_decompose_LUP+0x1a0>
mov    %eax,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rbp
mov    %rbp,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x3,%edx
cmp    $0x18,%rbp
mov    -0x60(%rsp),%rbx
jae    1170 <lm_decompose_LUP+0x70>
movdqa 0xe97(%rip),%xmm1
xor    %ebp,%ebp
jmp    1249 <lm_decompose_LUP+0x149>
and    $0xfffffffffffffffc,%rsi
movdqa 0xe84(%rip),%xmm1
xor    %ebp,%ebp
movdqa 0xe89(%rip),%xmm8
movdqa 0xe90(%rip),%xmm9
movdqa 0xe97(%rip),%xmm10
movdqa 0xe9e(%rip),%xmm11
movdqa 0xea6(%rip),%xmm6
movdqa 0xeae(%rip),%xmm7
movdqa 0xeb6(%rip),%xmm2
movdqa 0xebe(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movdqa %xmm1,%xmm4
paddd  %xmm8,%xmm4
movdqu %xmm1,(%rbx,%rbp,4)
movdqu %xmm4,0x10(%rbx,%rbp,4)
movdqa %xmm1,%xmm4
paddd  %xmm9,%xmm4
movdqa %xmm1,%xmm5
paddd  %xmm10,%xmm5
movdqu %xmm4,0x20(%rbx,%rbp,4)
movdqu %xmm5,0x30(%rbx,%rbp,4)
movdqa %xmm1,%xmm4
paddd  %xmm11,%xmm4
movdqa %xmm1,%xmm5
paddd  %xmm6,%xmm5
movdqu %xmm4,0x40(%rbx,%rbp,4)
movdqu %xmm5,0x50(%rbx,%rbp,4)
movdqa %xmm1,%xmm4
paddd  %xmm7,%xmm4
movdqa %xmm1,%xmm5
paddd  %xmm2,%xmm5
movdqu %xmm4,0x60(%rbx,%rbp,4)
movdqu %xmm5,0x70(%rbx,%rbp,4)
add    $0x20,%rbp
paddd  %xmm3,%xmm1
add    $0xfffffffffffffffc,%rsi
jne    11d0 <lm_decompose_LUP+0xd0>
test   %rdx,%rdx
je     1290 <lm_decompose_LUP+0x190>
lea    (%rbx,%rbp,4),%rsi
add    $0x10,%rsi
shl    $0x5,%rdx
xor    %ebp,%ebp
movdqa 0xdac(%rip),%xmm2
movdqa 0xdb4(%rip),%xmm3
nopl   0x0(%rax)
movdqa %xmm1,%xmm4
paddd  %xmm2,%xmm4
movdqu %xmm1,-0x10(%rsi,%rbp,1)
movdqu %xmm4,(%rsi,%rbp,1)
paddd  %xmm3,%xmm1
add    $0x20,%rbp
cmp    %rbp,%rdx
jne    1270 <lm_decompose_LUP+0x170>
cmp    %rax,%rcx
je     12ac <lm_decompose_LUP+0x1ac>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %ecx,(%rbx,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    12a0 <lm_decompose_LUP+0x1a0>
cmpl   $0x0,-0x7c(%rsp)
jle    1935 <lm_decompose_LUP+0x835>
lea    0x2(%rax),%rcx
mov    %rcx,0x88(%rsp)
lea    0x8(,%rax,8),%rcx
mov    %rcx,0x38(%rsp)
lea    (%rax,%rax,1),%rcx
mov    %rcx,0x30(%rsp)
lea    0x1(%rax),%rcx
mov    %rcx,0x80(%rsp)
lea    0x10(%rdi),%rcx
mov    %rcx,0x28(%rsp)
lea    -0x2(%rax),%r8
lea    (%rdi,%rax,8),%rdx
mov    %rdx,0x10(%rsp)
mov    %eax,%edx
and    $0xfffffffc,%edx
add    $0xfffffffffffffffc,%rdx
mov    %rdx,-0x8(%rsp)
mov    %rdx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %eax,%edx
and    $0xfffffffc,%edx
mov    %rdx,-0x18(%rsp)
mov    %rsi,%rdx
mov    %rsi,-0x10(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x20(%rsp)
mov    %rax,%rdx
neg    %rdx
mov    %rdx,-0x30(%rsp)
lea    0x0(,%rax,8),%rbx
lea    0x8(,%rax,8),%rdx
mov    %rdx,0x18(%rsp)
lea    0x30(%rdi),%rdx
movapd 0xd35(%rip),%xmm1
movapd 0xd3d(%rip),%xmm2
lea    0x8(%rdi),%rcx
lea    0x20(%rdi,%rax,8),%rbp
mov    %rbp,-0x40(%rsp)
lea    (%rdi,%rax,8),%rbp
mov    %rbp,-0x48(%rsp)
mov    $0x1,%ebp
mov    %rbp,-0x50(%rsp)
xor    %r9d,%r9d
mov    %rdi,%r14
mov    %r8,0x20(%rsp)
mov    %rcx,0x8(%rsp)
mov    %rcx,%rsi
mov    %rdx,-0x28(%rsp)
mov    %rdx,%r12
mov    %rdi,%r15
mov    %rbx,-0x70(%rsp)
jmp    13f2 <lm_decompose_LUP+0x2f2>
nopl   0x0(%rax)
addq   $0x1,-0x50(%rsp)
mov    0x18(%rsp),%rcx
add    %rcx,%r15
mov    0x48(%rsp),%r12
add    %rbx,%r12
mov    0x50(%rsp),%rsi
add    %rbx,%rsi
mov    0x58(%rsp),%r8
add    $0xffffffffffffffff,%r8
add    %rcx,-0x40(%rsp)
add    %rbx,-0x48(%rsp)
add    %rbx,%r14
mov    0x40(%rsp),%r9
cmp    %rax,%r9
je     1935 <lm_decompose_LUP+0x835>
mov    %r8,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x70(%rsp)
mov    %r9,%rdx
imul   %rax,%rdx
lea    (%rdi,%rdx,8),%r10
lea    (%rdx,%rax,1),%rcx
lea    (%rdi,%rcx,8),%rcx
mov    %rcx,-0x38(%rsp)
mov    0x20(%rsp),%rcx
sub    %r9,%rcx
mov    %rcx,-0x58(%rsp)
mov    0x38(%rsp),%rbp
imul   %r9,%rbp
add    0x30(%rsp),%rdx
mov    %rdx,0xa0(%rsp)
mov    0x28(%rsp),%rcx
mov    %rbp,0x98(%rsp)
add    %rbp,%rcx
mov    %rcx,0x90(%rsp)
pxor   %xmm3,%xmm3
mov    %r15,%rdx
mov    %r9,%rbp
mov    %r9,%rcx
mov    %r9,-0x78(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movdqa %xmm3,%xmm4
movsd  (%rdx),%xmm3
andpd  %xmm1,%xmm3
ucomisd %xmm4,%xmm3
cmova  %ebp,%ecx
maxsd  %xmm4,%xmm3
add    $0x1,%rbp
add    %rbx,%rdx
cmp    %rbp,%rax
jne    1470 <lm_decompose_LUP+0x370>
ucomisd %xmm3,%xmm0
ja     192b <lm_decompose_LUP+0x82b>
mov    %r8,0x58(%rsp)
mov    %ecx,%edx
cmp    %rdx,-0x78(%rsp)
je     1670 <lm_decompose_LUP+0x570>
mov    -0x78(%rsp),%rbx
mov    %rbx,%r9
imul   %rax,%r9
movslq %ecx,%r11
mov    -0x7c(%rsp),%r13d
imul   %r13d,%ecx
mov    -0x60(%rsp),%rdx
mov    (%rdx,%rbx,4),%r8d
mov    (%rdx,%r11,4),%ebp
mov    %ebp,(%rdx,%rbx,4)
mov    %r8d,(%rdx,%r11,4)
movslq %ecx,%r8
cmp    $0x4,%r13d
jb     1505 <lm_decompose_LUP+0x405>
mov    0x10(%rsp),%rcx
lea    (%rcx,%r8,8),%rcx
cmp    %rcx,%r10
jae    158a <lm_decompose_LUP+0x48a>
lea    (%rdi,%r8,8),%rcx
cmp    -0x38(%rsp),%rcx
jae    158a <lm_decompose_LUP+0x48a>
xor    %edx,%edx
mov    -0x70(%rsp),%rbx
mov    %rdx,%rcx
test   $0x1,%al
je     1537 <lm_decompose_LUP+0x437>
add    %rdx,%r9
movsd  (%rdi,%r9,8),%xmm3
lea    (%rdx,%r8,1),%rcx
movsd  (%rdi,%rcx,8),%xmm4
movsd  %xmm4,(%rdi,%r9,8)
movsd  %xmm3,(%rdi,%rcx,8)
mov    %rdx,%rcx
or     $0x1,%rcx
not    %rdx
cmp    -0x30(%rsp),%rdx
je     1670 <lm_decompose_LUP+0x570>
mov    0x8(%rsp),%rdx
lea    (%rdx,%r8,8),%rdx
xchg   %ax,%ax
movsd  -0x8(%rsi,%rcx,8),%xmm3
movsd  -0x8(%rdx,%rcx,8),%xmm4
movsd  %xmm4,-0x8(%rsi,%rcx,8)
movsd  %xmm3,-0x8(%rdx,%rcx,8)
movsd  (%rsi,%rcx,8),%xmm3
movsd  (%rdx,%rcx,8),%xmm4
movsd  %xmm4,(%rsi,%rcx,8)
movsd  %xmm3,(%rdx,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rax
jne    1550 <lm_decompose_LUP+0x450>
jmp    1670 <lm_decompose_LUP+0x570>
cmpq   $0x0,-0x8(%rsp)
mov    -0x70(%rsp),%rbx
je     1919 <lm_decompose_LUP+0x819>
mov    -0x28(%rsp),%rcx
lea    (%rcx,%r8,8),%rcx
mov    -0x20(%rsp),%rdx
xor    %ebp,%ebp
nopl   0x0(%rax,%rax,1)
movups -0x30(%r12,%rbp,8),%xmm3
movups -0x20(%r12,%rbp,8),%xmm4
movups -0x30(%rcx,%rbp,8),%xmm5
movups -0x20(%rcx,%rbp,8),%xmm6
movups %xmm5,-0x30(%r12,%rbp,8)
movups %xmm6,-0x20(%r12,%rbp,8)
movups %xmm3,-0x30(%rcx,%rbp,8)
movups %xmm4,-0x20(%rcx,%rbp,8)
movupd -0x10(%r12,%rbp,8),%xmm3
movupd (%r12,%rbp,8),%xmm4
movdqu -0x10(%rcx,%rbp,8),%xmm5
movdqu (%rcx,%rbp,8),%xmm6
movdqu %xmm5,-0x10(%r12,%rbp,8)
movdqu %xmm6,(%r12,%rbp,8)
movupd %xmm3,-0x10(%rcx,%rbp,8)
movupd %xmm4,(%rcx,%rbp,8)
add    $0x8,%rbp
add    $0xfffffffffffffffe,%rdx
jne    15b0 <lm_decompose_LUP+0x4b0>
testb  $0x1,-0x10(%rsp)
je     1650 <lm_decompose_LUP+0x550>
lea    (%r9,%rbp,1),%rcx
movupd (%rdi,%rcx,8),%xmm3
movupd 0x10(%rdi,%rcx,8),%xmm4
add    %r8,%rbp
movdqu (%rdi,%rbp,8),%xmm5
movdqu 0x10(%rdi,%rbp,8),%xmm6
movdqu %xmm5,(%rdi,%rcx,8)
movdqu %xmm6,0x10(%rdi,%rcx,8)
movupd %xmm3,(%rdi,%rbp,8)
movupd %xmm4,0x10(%rdi,%rbp,8)
mov    -0x18(%rsp),%rcx
mov    %rcx,%rdx
cmp    %rax,%rcx
jne    150c <lm_decompose_LUP+0x40c>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r12,0x48(%rsp)
mov    %rsi,0x50(%rsp)
mov    -0x78(%rsp),%rcx
add    $0x1,%rcx
mov    %rcx,0x40(%rsp)
cmp    %rax,%rcx
jae    13b0 <lm_decompose_LUP+0x2b0>
mov    -0x78(%rsp),%rcx
mov    %rcx,%rdx
imul   %rax,%rdx
add    %rdx,%rcx
mov    %rcx,0xb8(%rsp)
mov    -0x50(%rsp),%r9
mov    %rdx,0x78(%rsp)
lea    (%r9,%rdx,1),%rcx
mov    %rcx,0xb0(%rsp)
lea    0x1(%r9),%rcx
mov    %rcx,-0x68(%rsp)
mov    -0x58(%rsp),%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x68(%rsp)
add    %r9,%rcx
add    $0x1,%rcx
mov    %rcx,0x60(%rsp)
mov    -0x48(%rsp),%r11
mov    -0x40(%rsp),%r13
xor    %r10d,%r10d
jmp    1707 <lm_decompose_LUP+0x607>
add    $0x1,%r9
add    $0x1,%r10
add    %rbx,%r13
add    %rbx,%r11
cmp    %rax,%r9
je     13b0 <lm_decompose_LUP+0x2b0>
mov    %r9,%rbp
imul   %rax,%rbp
mov    -0x78(%rsp),%rcx
lea    (%rcx,%rbp,1),%r8
movsd  (%rdi,%r8,8),%xmm3
mov    0xb8(%rsp),%rcx
divsd  (%rdi,%rcx,8),%xmm3
mov    -0x50(%rsp),%rcx
add    %rbp,%rcx
movsd  %xmm3,(%rdi,%r8,8)
mov    0xb0(%rsp),%rdx
mulsd  (%rdi,%rdx,8),%xmm3
movsd  (%rdi,%rcx,8),%xmm4
subsd  %xmm3,%xmm4
movsd  %xmm4,(%rdi,%rcx,8)
cmp    %rax,-0x68(%rsp)
je     16f0 <lm_decompose_LUP+0x5f0>
mov    -0x68(%rsp),%r12
cmpq   $0x4,-0x58(%rsp)
jb     1880 <lm_decompose_LUP+0x780>
mov    %r10,%rcx
imul   %rax,%rcx
mov    0x88(%rsp),%rdx
add    %rcx,%rdx
lea    (%rdi,%rdx,8),%rdx
mov    0x98(%rsp),%r12
lea    (%rdx,%r12,1),%rsi
mov    %rsi,0xa8(%rsp)
mov    0xa0(%rsp),%rsi
add    %rcx,%rsi
lea    (%rcx,%rax,1),%rbx
lea    (%rdi,%rbx,8),%rbx
add    %r12,%rbx
add    0x80(%rsp),%rcx
lea    (%rdi,%rsi,8),%rsi
lea    (%rdi,%rcx,8),%rcx
cmp    %rcx,%rdx
setb   %dl
cmp    %rsi,%rbx
setb   %bl
mov    0xa8(%rsp),%rcx
cmp    -0x38(%rsp),%rcx
setb   %cl
cmp    %rsi,0x90(%rsp)
setb   %sil
mov    -0x68(%rsp),%r12
test   %bl,%dl
mov    -0x70(%rsp),%rbx
jne    1880 <lm_decompose_LUP+0x780>
mov    -0x68(%rsp),%r12
and    %sil,%cl
jne    1880 <lm_decompose_LUP+0x780>
movsd  (%rdi,%r8,8),%xmm3
unpcklpd %xmm3,%xmm3
xorpd  %xmm2,%xmm3
xor    %ecx,%ecx
mov    0x70(%rsp),%rdx
nopl   0x0(%rax,%rax,1)
movupd 0x10(%r15,%rcx,8),%xmm4
movupd 0x20(%r15,%rcx,8),%xmm5
movupd -0x10(%r13,%rcx,8),%xmm6
movupd 0x0(%r13,%rcx,8),%xmm7
mulpd  %xmm3,%xmm4
addpd  %xmm6,%xmm4
mulpd  %xmm3,%xmm5
addpd  %xmm7,%xmm5
movupd %xmm4,-0x10(%r13,%rcx,8)
movupd %xmm5,0x0(%r13,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1820 <lm_decompose_LUP+0x720>
mov    0x60(%rsp),%r12
mov    0x68(%rsp),%rcx
cmp    %rcx,-0x58(%rsp)
je     16f0 <lm_decompose_LUP+0x5f0>
nopl   0x0(%rax,%rax,1)
mov    %eax,%edx
sub    %r12d,%edx
mov    %r12,%rcx
test   $0x1,%dl
je     18b6 <lm_decompose_LUP+0x7b6>
movsd  (%rdi,%r8,8),%xmm3
mov    0x78(%rsp),%rcx
add    %r12,%rcx
add    %r12,%rbp
mulsd  (%rdi,%rcx,8),%xmm3
movsd  (%rdi,%rbp,8),%xmm4
subsd  %xmm3,%xmm4
movsd  %xmm4,(%rdi,%rbp,8)
lea    0x1(%r12),%rcx
not    %r12
cmp    -0x30(%rsp),%r12
je     16f0 <lm_decompose_LUP+0x5f0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rdi,%r8,8),%xmm3
mulsd  (%r14,%rcx,8),%xmm3
movsd  (%r11,%rcx,8),%xmm4
subsd  %xmm3,%xmm4
movsd  0x8(%r11,%rcx,8),%xmm3
movsd  %xmm4,(%r11,%rcx,8)
movsd  (%rdi,%r8,8),%xmm4
mulsd  0x8(%r14,%rcx,8),%xmm4
subsd  %xmm4,%xmm3
movsd  %xmm3,0x8(%r11,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rax
jne    18d0 <lm_decompose_LUP+0x7d0>
jmp    16f0 <lm_decompose_LUP+0x5f0>
xor    %ebp,%ebp
testb  $0x1,-0x10(%rsp)
jne    161d <lm_decompose_LUP+0x51d>
jmp    1650 <lm_decompose_LUP+0x550>
mov    (%rsp),%rax
movl   $0x15,(%rax)
add    $0xc0,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @lm_decompose_LUP(double* %0, i32 %1, double %2, i32* %3, i32* %4) {
%6 = alloca double*
%7 = alloca i32
%8 = alloca double
%9 = alloca i32*
%10 = alloca i32*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca double
%16 = alloca double
%17 = alloca double
store double* %0, double** %6
store i32 %1, i32* %7
store double %2, double* %8
store i32* %3, i32** %9
store i32* %4, i32** %10
store i32 0, i32* %11
br label %18
18: ; preds = %28, %5
%19 = load i32, i32* %11
%20 = load i32, i32* %7
%21 = icmp slt i32 %19, %20
br i1 %21, label %22, label %31
22: ; preds = %18
%23 = load i32, i32* %11
%24 = load i32*, i32** %9
%25 = load i32, i32* %11
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds i32, i32* %24, i64 %26
store i32 %23, i32* %27
br label %28
28: ; preds = %22
%29 = load i32, i32* %11
%30 = add nsw i32 %29, 1
store i32 %30, i32* %11
br label %18
31: ; preds = %18
store i32 0, i32* %11
br label %32
32: ; preds = %208, %31
%33 = load i32, i32* %11
%34 = load i32, i32* %7
%35 = icmp slt i32 %33, %34
br i1 %35, label %36, label %211
36: ; preds = %32
store double 0.000000e+00, double* %15
%37 = load i32, i32* %11
store i32 %37, i32* %14
%38 = load i32, i32* %11
store i32 %38, i32* %13
br label %39
39: ; preds = %60, %36
%40 = load i32, i32* %13
%41 = load i32, i32* %7
%42 = icmp slt i32 %40, %41
br i1 %42, label %43, label %63
43: ; preds = %39
%44 = load double*, double** %6
%45 = load i32, i32* %13
%46 = load i32, i32* %7
%47 = mul nsw i32 %45, %46
%48 = load i32, i32* %11
%49 = add nsw i32 %47, %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %44, i64 %50
%52 = load double, double* %51
%53 = call double @llvm.fabs.f64(double %52)
store double %53, double* %16
%54 = load double, double* %15
%55 = fcmp ogt double %53, %54
br i1 %55, label %56, label %59
56: ; preds = %43
%57 = load double, double* %16
store double %57, double* %15
%58 = load i32, i32* %13
store i32 %58, i32* %14
br label %59
59: ; preds = %56, %43
br label %60
60: ; preds = %59
%61 = load i32, i32* %13
%62 = add nsw i32 %61, 1
store i32 %62, i32* %13
br label %39
63: ; preds = %39
%64 = load double, double* %15
%65 = load double, double* %8
%66 = fcmp olt double %64, %65
br i1 %66, label %67, label %69
67: ; preds = %63
%68 = load i32*, i32** %10
store i32 21, i32* %68
br label %211
69: ; preds = %63
%70 = load i32, i32* %14
%71 = load i32, i32* %11
%72 = icmp ne i32 %70, %71
br i1 %72, label %73, label %137
73: ; preds = %69
%74 = load i32*, i32** %9
%75 = load i32, i32* %11
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds i32, i32* %74, i64 %76
%78 = load i32, i32* %77
store i32 %78, i32* %12
%79 = load i32*, i32** %9
%80 = load i32, i32* %14
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds i32, i32* %79, i64 %81
%83 = load i32, i32* %82
%84 = load i32*, i32** %9
%85 = load i32, i32* %11
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds i32, i32* %84, i64 %86
store i32 %83, i32* %87
%88 = load i32, i32* %12
%89 = load i32*, i32** %9
%90 = load i32, i32* %14
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds i32, i32* %89, i64 %91
store i32 %88, i32* %92
store i32 0, i32* %13
br label %93
93: ; preds = %133, %73
%94 = load i32, i32* %13
%95 = load i32, i32* %7
%96 = icmp slt i32 %94, %95
br i1 %96, label %97, label %136
97: ; preds = %93
%98 = load double*, double** %6
%99 = load i32, i32* %11
%100 = load i32, i32* %7
%101 = mul nsw i32 %99, %100
%102 = load i32, i32* %13
%103 = add nsw i32 %101, %102
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds double, double* %98, i64 %104
%106 = load double, double* %105
store double %106, double* %17
%107 = load double*, double** %6
%108 = load i32, i32* %14
%109 = load i32, i32* %7
%110 = mul nsw i32 %108, %109
%111 = load i32, i32* %13
%112 = add nsw i32 %110, %111
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %107, i64 %113
%115 = load double, double* %114
%116 = load double*, double** %6
%117 = load i32, i32* %11
%118 = load i32, i32* %7
%119 = mul nsw i32 %117, %118
%120 = load i32, i32* %13
%121 = add nsw i32 %119, %120
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds double, double* %116, i64 %122
store double %115, double* %123
%124 = load double, double* %17
%125 = load double*, double** %6
%126 = load i32, i32* %14
%127 = load i32, i32* %7
%128 = mul nsw i32 %126, %127
%129 = load i32, i32* %13
%130 = add nsw i32 %128, %129
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds double, double* %125, i64 %131
store double %124, double* %132
br label %133
133: ; preds = %97
%134 = load i32, i32* %13
%135 = add nsw i32 %134, 1
store i32 %135, i32* %13
br label %93
136: ; preds = %93
br label %137
137: ; preds = %136, %69
%138 = load i32, i32* %11
%139 = add nsw i32 %138, 1
store i32 %139, i32* %12
br label %140
140: ; preds = %204, %137
%141 = load i32, i32* %12
%142 = load i32, i32* %7
%143 = icmp slt i32 %141, %142
br i1 %143, label %144, label %207
144: ; preds = %140
%145 = load double*, double** %6
%146 = load i32, i32* %11
%147 = load i32, i32* %7
%148 = mul nsw i32 %146, %147
%149 = load i32, i32* %11
%150 = add nsw i32 %148, %149
%151 = sext i32 %150 to i64
%152 = getelementptr inbounds double, double* %145, i64 %151
%153 = load double, double* %152
%154 = load double*, double** %6
%155 = load i32, i32* %12
%156 = load i32, i32* %7
%157 = mul nsw i32 %155, %156
%158 = load i32, i32* %11
%159 = add nsw i32 %157, %158
%160 = sext i32 %159 to i64
%161 = getelementptr inbounds double, double* %154, i64 %160
%162 = load double, double* %161
%163 = fdiv double %162, %153
store double %163, double* %161
%164 = load i32, i32* %11
%165 = add nsw i32 %164, 1
store i32 %165, i32* %13
br label %166
166: ; preds = %200, %144
%167 = load i32, i32* %13
%168 = load i32, i32* %7
%169 = icmp slt i32 %167, %168
br i1 %169, label %170, label %203
170: ; preds = %166
%171 = load double*, double** %6
%172 = load i32, i32* %12
%173 = load i32, i32* %7
%174 = mul nsw i32 %172, %173
%175 = load i32, i32* %11
%176 = add nsw i32 %174, %175
%177 = sext i32 %176 to i64
%178 = getelementptr inbounds double, double* %171, i64 %177
%179 = load double, double* %178
%180 = load double*, double** %6
%181 = load i32, i32* %11
%182 = load i32, i32* %7
%183 = mul nsw i32 %181, %182
%184 = load i32, i32* %13
%185 = add nsw i32 %183, %184
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds double, double* %180, i64 %186
%188 = load double, double* %187
%189 = load double*, double** %6
%190 = load i32, i32* %12
%191 = load i32, i32* %7
%192 = mul nsw i32 %190, %191
%193 = load i32, i32* %13
%194 = add nsw i32 %192, %193
%195 = sext i32 %194 to i64
%196 = getelementptr inbounds double, double* %189, i64 %195
%197 = load double, double* %196
%198 = fneg double %179
%199 = call double @llvm.fmuladd.f64(double %198, double %188, double %197)
store double %199, double* %196
br label %200
200: ; preds = %170
%201 = load i32, i32* %13
%202 = add nsw i32 %201, 1
store i32 %202, i32* %13
br label %166
203: ; preds = %166
br label %204
204: ; preds = %203
%205 = load i32, i32* %12
%206 = add nsw i32 %205, 1
store i32 %206, i32* %12
br label %140
207: ; preds = %140
br label %208
208: ; preds = %207
%209 = load i32, i32* %11
%210 = add nsw i32 %209, 1
store i32 %210, i32* %11
br label %32
211: ; preds = %67, %32
ret void
}

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