#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 */
int /*<<< orphan*/ BubbleSort(int *, int);
int /*<<< orphan*/ ComputeFdualXb(int, double *);
int /*<<< orphan*/ ComputeH(double *, double *, int);
int /*<<< orphan*/ ComputeHs(int *, int, int, double *, double *);
int /*<<< orphan*/ TriSymGaxpy(double *, double *, double *, int, double *);
int /*<<< orphan*/ TrisolveREG(double *, double *, double *, double *, double *,
                               int);
double *calloc(int, int);
int /*<<< orphan*/ free(double *);
int /*<<< orphan*/ reconstruct(double *, int, double *);
void Project(double *y, int M, int *I, int L, double *xI, double *wI) {
  double aux_double = 0;
  double ymean = 0;
  int i = 0;
  double *h0;
  double *h1;
  double *w0;
  double *z;
  double *temp;
  double *wr;
  double *aux_vec;
  h0 = calloc(M - 1, sizeof(double));
  h1 = calloc(M - 2, sizeof(double));
  z = calloc(M - 1, sizeof(double));
  temp = calloc(L, sizeof(double));
  wr = calloc(M, sizeof(double));
  aux_vec = calloc(M, sizeof(double));
  aux_double = 0;
  for (i = 0; i < M; i++) {
    aux_double = aux_double + y[i];
  }
  ymean = (double)(aux_double / M);
  for (i = 0; i < M; i++) {
    y[i] = y[i] - ymean;
  }
  if (L > 0) {
    BubbleSort(I, L);
    ComputeH(h0, h1, M);
    ComputeFdualXb(M, y);
    w0 = y;
    TrisolveREG(h0, h1, h1, w0, z, M - 1);
    ComputeHs(I, M, L, h0, h1);
    for (i = 0; i < L; i++) {
      temp[i] = z[I[i]];
    }
    for (i = 0; i < L; i++) {
      y[i] = 0;
    }
    TriSymGaxpy(h0, h1, temp, L, wI);
    for (i = 0; i < L; i++) {
      wr[I[i]] = wI[i];
    }
    reconstruct(wr, M, aux_vec);
    for (i = 0; i < M; i++) {
      xI[i] = ymean + wr[i];
    }
  } else {
    L = 0;
    for (i = 0; i < M; i++) {
      xI[i] = ymean;
    }
  }
  free(h0);
  free(h1);
  free(z);
  free(temp);
  free(wr);
  free(aux_vec);
}


x86 assembly (clang, -O0):
<Project>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0xc(%rbp),%edi
sub    $0x1,%edi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x50(%rbp)
mov    -0xc(%rbp),%edi
sub    $0x2,%edi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x58(%rbp)
mov    -0xc(%rbp),%edi
sub    $0x1,%edi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x68(%rbp)
mov    -0x1c(%rbp),%edi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x70(%rbp)
mov    -0xc(%rbp),%edi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x78(%rbp)
mov    -0xc(%rbp),%edi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x80(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1277 <Project+0xe7>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1246 <Project+0xb6>
movsd  -0x38(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12ca <Project+0x13a>
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1291 <Project+0x101>
cmpl   $0x0,-0x1c(%rbp)
jle    1460 <Project+0x2d0>
mov    -0x18(%rbp),%rdi
mov    -0x1c(%rbp),%esi
call   10b0 <BubbleSort@plt>
mov    -0x50(%rbp),%rdi
mov    -0x58(%rbp),%rsi
mov    -0xc(%rbp),%edx
call   1050 <ComputeH@plt>
mov    -0xc(%rbp),%edi
mov    -0x8(%rbp),%rsi
call   10a0 <ComputeFdualXb@plt>
mov    -0x8(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x50(%rbp),%rdi
mov    -0x58(%rbp),%rsi
mov    -0x58(%rbp),%rdx
mov    -0x60(%rbp),%rcx
mov    -0x68(%rbp),%r8
mov    -0xc(%rbp),%r9d
sub    $0x1,%r9d
call   1090 <TrisolveREG@plt>
mov    -0x18(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x1c(%rbp),%edx
mov    -0x50(%rbp),%rcx
mov    -0x58(%rbp),%r8
call   1080 <ComputeHs@plt>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    137f <Project+0x1ef>
mov    -0x68(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x44(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x70(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1343 <Project+0x1b3>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13b0 <Project+0x220>
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1386 <Project+0x1f6>
mov    -0x50(%rbp),%rdi
mov    -0x58(%rbp),%rsi
mov    -0x70(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
mov    -0x30(%rbp),%r8
call   1060 <TriSymGaxpy@plt>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    140b <Project+0x27b>
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x78(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x44(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    13cf <Project+0x23f>
mov    -0x78(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x80(%rbp),%rdx
call   1030 <reconstruct@plt>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    145b <Project+0x2cb>
movsd  -0x40(%rbp),%xmm0
mov    -0x78(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1422 <Project+0x292>
jmp    149f <Project+0x30f>
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    149a <Project+0x30a>
movsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    146e <Project+0x2de>
jmp    149f <Project+0x30f>
mov    -0x50(%rbp),%rdi
call   1040 <free@plt>
mov    -0x58(%rbp),%rdi
call   1040 <free@plt>
mov    -0x68(%rbp),%rdi
call   1040 <free@plt>
mov    -0x70(%rbp),%rdi
call   1040 <free@plt>
mov    -0x78(%rbp),%rdi
call   1040 <free@plt>
mov    -0x80(%rbp),%rdi
call   1040 <free@plt>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Project>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r9,0x48(%rsp)
mov    %r8,%rbx
mov    %ecx,%r15d
mov    %rdx,%r14
mov    %esi,%r12d
mov    %rdi,%r13
lea    -0x1(%r12),%ebp
mov    %ebp,%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,0x30(%rsp)
lea    -0x2(%r12),%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,0x28(%rsp)
mov    %ebp,0x10(%rsp)
mov    %ebp,%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,0x8(%rsp)
mov    %r15d,0x4(%rsp)
mov    %r15d,%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,0x20(%rsp)
mov    %r12d,%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,%r15
mov    %r12d,%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,0x18(%rsp)
cvtsi2sd %r12d,%xmm0
test   %r12d,%r12d
jle    125e <Project+0xbe>
mov    %r12d,%eax
lea    -0x1(%rax),%rdx
mov    %eax,%ecx
and    $0x7,%ecx
cmp    $0x7,%rdx
jae    126b <Project+0xcb>
xorpd  %xmm5,%xmm5
xor    %edx,%edx
jmp    12c1 <Project+0x121>
xorpd  %xmm5,%xmm5
divsd  %xmm0,%xmm5
jmp    13eb <Project+0x24b>
mov    %eax,%esi
and    $0xfffffff8,%esi
xorpd  %xmm5,%xmm5
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
addsd  0x0(%r13,%rdx,8),%xmm5
addsd  0x8(%r13,%rdx,8),%xmm5
addsd  0x10(%r13,%rdx,8),%xmm5
addsd  0x18(%r13,%rdx,8),%xmm5
addsd  0x20(%r13,%rdx,8),%xmm5
addsd  0x28(%r13,%rdx,8),%xmm5
addsd  0x30(%r13,%rdx,8),%xmm5
addsd  0x38(%r13,%rdx,8),%xmm5
add    $0x8,%rdx
cmp    %rdx,%rsi
jne    1280 <Project+0xe0>
test   %rcx,%rcx
je     12ee <Project+0x14e>
lea    0x0(,%rdx,8),%rdx
add    %r13,%rdx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rdx,%rsi,8),%xmm5
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    12e0 <Project+0x140>
divsd  %xmm0,%xmm5
test   %r12d,%r12d
jle    13eb <Project+0x24b>
cmp    $0x4,%r12d
jae    1308 <Project+0x168>
xor    %ecx,%ecx
jmp    13d0 <Project+0x230>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm5,%xmm0
unpcklpd %xmm5,%xmm0
lea    -0x4(%rcx),%rdx
mov    %rdx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
test   %rdx,%rdx
je     174b <Project+0x5ab>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movupd 0x0(%r13,%rdx,8),%xmm1
movupd 0x10(%r13,%rdx,8),%xmm2
movupd 0x20(%r13,%rdx,8),%xmm3
movupd 0x30(%r13,%rdx,8),%xmm4
subpd  %xmm0,%xmm1
subpd  %xmm0,%xmm2
movupd %xmm1,0x0(%r13,%rdx,8)
movupd %xmm2,0x10(%r13,%rdx,8)
subpd  %xmm0,%xmm3
subpd  %xmm0,%xmm4
movupd %xmm3,0x20(%r13,%rdx,8)
movupd %xmm4,0x30(%r13,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1340 <Project+0x1a0>
test   $0x1,%sil
je     13bc <Project+0x21c>
movupd 0x0(%r13,%rdx,8),%xmm1
movupd 0x10(%r13,%rdx,8),%xmm2
subpd  %xmm0,%xmm1
subpd  %xmm0,%xmm2
movupd %xmm1,0x0(%r13,%rdx,8)
movupd %xmm2,0x10(%r13,%rdx,8)
cmp    %rax,%rcx
je     13eb <Project+0x24b>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rcx,8),%xmm0
subsd  %xmm5,%xmm0
movsd  %xmm0,0x0(%r13,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    13d0 <Project+0x230>
mov    0x4(%rsp),%esi
test   %esi,%esi
jle    148b <Project+0x2eb>
mov    %rbx,0x40(%rsp)
movapd %xmm5,0x50(%rsp)
mov    %r14,%rdi
call   10c0 <BubbleSort@plt>
mov    0x30(%rsp),%rbx
mov    %rbx,%rdi
mov    %r12,%rbp
mov    0x28(%rsp),%r12
mov    %r12,%rsi
mov    %ebp,%edx
call   1050 <ComputeH@plt>
mov    %ebp,%edi
mov    %r13,%rsi
call   10b0 <ComputeFdualXb@plt>
mov    %rbx,%rdi
mov    %r12,%rsi
mov    %r12,%rdx
mov    %r13,%rcx
mov    0x8(%rsp),%r8
mov    0x10(%rsp),%r9d
call   10a0 <TrisolveREG@plt>
mov    %r14,%rdi
mov    %rbp,0x10(%rsp)
mov    %ebp,%esi
mov    0x4(%rsp),%ebp
mov    %ebp,%edx
mov    %rbx,%rcx
mov    %r12,%r8
call   1090 <ComputeHs@plt>
mov    %ebp,%r12d
lea    -0x1(%r12),%rcx
mov    %r12d,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
mov    %rcx,0x38(%rsp)
jae    14ad <Project+0x30d>
xor    %ecx,%ecx
mov    0x8(%rsp),%rbx
mov    0x20(%rsp),%rbp
jmp    1508 <Project+0x368>
test   %r12d,%r12d
mov    0x18(%rsp),%r13
jle    184e <Project+0x6ae>
mov    %r12d,%eax
cmp    $0x4,%r12d
jae    171d <Project+0x57d>
xor    %ecx,%ecx
jmp    1840 <Project+0x6a0>
mov    %r12d,%edx
and    $0xfffffffc,%edx
xor    %ecx,%ecx
mov    0x8(%rsp),%rbx
mov    0x20(%rsp),%rbp
nop
movslq (%r14,%rcx,4),%rsi
movsd  (%rbx,%rsi,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rcx,8)
movslq 0x4(%r14,%rcx,4),%rsi
movsd  (%rbx,%rsi,8),%xmm0
movsd  %xmm0,0x8(%rbp,%rcx,8)
movslq 0x8(%r14,%rcx,4),%rsi
movsd  (%rbx,%rsi,8),%xmm0
movsd  %xmm0,0x10(%rbp,%rcx,8)
movslq 0xc(%r14,%rcx,4),%rsi
movsd  (%rbx,%rsi,8),%xmm0
movsd  %xmm0,0x18(%rbp,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    14c0 <Project+0x320>
test   %rax,%rax
je     1537 <Project+0x397>
lea    0x0(,%rcx,8),%rdx
add    %rbp,%rdx
lea    (%r14,%rcx,4),%rcx
xor    %esi,%esi
xchg   %ax,%ax
movslq (%rcx,%rsi,4),%rdi
movsd  (%rbx,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    1520 <Project+0x380>
lea    0x0(,%r12,8),%rdx
xor    %ebx,%ebx
mov    %r13,%rdi
xor    %esi,%esi
mov    %rbp,%r13
call   1070 <memset@plt>
mov    0x30(%rsp),%rdi
mov    0x28(%rsp),%rsi
mov    %rbp,%rdx
mov    0x4(%rsp),%ecx
mov    0x48(%rsp),%rbp
mov    %rbp,%r8
call   1060 <TriSymGaxpy@plt>
mov    %r12d,%eax
and    $0x3,%eax
cmpq   $0x3,0x38(%rsp)
jae    1581 <Project+0x3e1>
mov    0x18(%rsp),%r13
jmp    15dc <Project+0x43c>
and    $0xfffffffc,%r12d
xor    %ebx,%ebx
mov    0x18(%rsp),%r13
nopl   0x0(%rax)
movsd  0x0(%rbp,%rbx,8),%xmm0
movslq (%r14,%rbx,4),%rcx
movsd  %xmm0,(%r15,%rcx,8)
movsd  0x8(%rbp,%rbx,8),%xmm0
movslq 0x4(%r14,%rbx,4),%rcx
movsd  %xmm0,(%r15,%rcx,8)
movsd  0x10(%rbp,%rbx,8),%xmm0
movslq 0x8(%r14,%rbx,4),%rcx
movsd  %xmm0,(%r15,%rcx,8)
movsd  0x18(%rbp,%rbx,8),%xmm0
movslq 0xc(%r14,%rbx,4),%rcx
movsd  %xmm0,(%r15,%rcx,8)
add    $0x4,%rbx
cmp    %rbx,%r12
jne    1590 <Project+0x3f0>
test   %rax,%rax
je     1618 <Project+0x478>
lea    (%r14,%rbx,4),%rcx
lea    0x0(,%rbx,8),%rdx
add    %rbp,%rdx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rdx,%rsi,8),%xmm0
movslq (%rcx,%rsi,4),%rdi
movsd  %xmm0,(%r15,%rdi,8)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    1600 <Project+0x460>
mov    %r15,%rdi
mov    0x10(%rsp),%rbp
mov    %ebp,%esi
mov    %r13,%rdx
call   1030 <reconstruct@plt>
test   %ebp,%ebp
mov    0x40(%rsp),%rbx
movapd 0x50(%rsp),%xmm3
jle    184e <Project+0x6ae>
mov    %ebp,%eax
cmp    $0x4,%ebp
jae    164b <Project+0x4ab>
xor    %ecx,%ecx
jmp    1700 <Project+0x560>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
lea    -0x4(%rcx),%rdx
mov    %rdx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
test   %rdx,%rdx
je     175c <Project+0x5bc>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %edx,%edx
nopl   0x0(%rax)
movupd (%r15,%rdx,8),%xmm1
movupd 0x10(%r15,%rdx,8),%xmm2
addpd  %xmm0,%xmm1
addpd  %xmm0,%xmm2
movupd %xmm1,(%rbx,%rdx,8)
movupd %xmm2,0x10(%rbx,%rdx,8)
movupd 0x20(%r15,%rdx,8),%xmm1
movupd 0x30(%r15,%rdx,8),%xmm2
addpd  %xmm0,%xmm1
addpd  %xmm0,%xmm2
movupd %xmm1,0x20(%rbx,%rdx,8)
movupd %xmm2,0x30(%rbx,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1680 <Project+0x4e0>
test   $0x1,%sil
je     16f2 <Project+0x552>
movupd (%r15,%rdx,8),%xmm1
movupd 0x10(%r15,%rdx,8),%xmm2
addpd  %xmm0,%xmm1
addpd  %xmm0,%xmm2
movupd %xmm1,(%rbx,%rdx,8)
movupd %xmm2,0x10(%rbx,%rdx,8)
cmp    %rax,%rcx
je     184e <Project+0x6ae>
nopl   0x0(%rax,%rax,1)
movsd  (%r15,%rcx,8),%xmm0
addsd  %xmm3,%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1700 <Project+0x560>
jmp    184e <Project+0x6ae>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm5,%xmm0
unpcklpd %xmm5,%xmm0
lea    -0x4(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x7,%edx
cmp    $0x1c,%rdi
jae    176a <Project+0x5ca>
xor    %edi,%edi
jmp    1805 <Project+0x665>
xor    %edx,%edx
test   $0x1,%sil
jne    1398 <Project+0x1f8>
jmp    13bc <Project+0x21c>
xor    %edx,%edx
test   $0x1,%sil
jne    16d2 <Project+0x532>
jmp    16f2 <Project+0x552>
movabs $0x7ffffffffffffffe,%rdi
ff ff 7f
add    $0xfffffffffffffffa,%rdi
and    %rdi,%rsi
xor    %edi,%edi
nopl   (%rax)
movupd %xmm0,(%rbx,%rdi,8)
movupd %xmm0,0x10(%rbx,%rdi,8)
movupd %xmm0,0x20(%rbx,%rdi,8)
movupd %xmm0,0x30(%rbx,%rdi,8)
movupd %xmm0,0x40(%rbx,%rdi,8)
movupd %xmm0,0x50(%rbx,%rdi,8)
movupd %xmm0,0x60(%rbx,%rdi,8)
movupd %xmm0,0x70(%rbx,%rdi,8)
movupd %xmm0,0x80(%rbx,%rdi,8)
movupd %xmm0,0x90(%rbx,%rdi,8)
movupd %xmm0,0xa0(%rbx,%rdi,8)
movupd %xmm0,0xb0(%rbx,%rdi,8)
movupd %xmm0,0xc0(%rbx,%rdi,8)
movupd %xmm0,0xd0(%rbx,%rdi,8)
movupd %xmm0,0xe0(%rbx,%rdi,8)
movupd %xmm0,0xf0(%rbx,%rdi,8)
add    $0x20,%rdi
add    $0xfffffffffffffff8,%rsi
jne    1780 <Project+0x5e0>
test   %rdx,%rdx
je     1834 <Project+0x694>
lea    (%rbx,%rdi,8),%rsi
add    $0x10,%rsi
shl    $0x5,%rdx
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movupd %xmm0,-0x10(%rsi,%rdi,1)
movupd %xmm0,(%rsi,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rdx
jne    1820 <Project+0x680>
cmp    %rax,%rcx
je     184e <Project+0x6ae>
nopl   0x0(%rax)
movsd  %xmm5,(%rbx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1840 <Project+0x6a0>
mov    0x30(%rsp),%rdi
call   1040 <free@plt>
mov    0x28(%rsp),%rdi
call   1040 <free@plt>
mov    0x8(%rsp),%rdi
call   1040 <free@plt>
mov    0x20(%rsp),%rdi
call   1040 <free@plt>
mov    %r15,%rdi
call   1040 <free@plt>
mov    %r13,%rdi
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free@plt>

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

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