#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 array_max(double *, int);
double array_max_abs(double *, int);
double array_min(double *, int);
double cholmod(double *, int, double *, double);
scalar_t__ fabs(double);
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
double sqrt(double);
double modelhess(double *A, int N, double *dx, double *L) {
  double *U22;
  double sqrteps, maxdiag, mindiag, maxposdiag, u;
  double maxoffdiag, maxinp, maxadd;
  double maxev, minev, offrow, sdd;
  int step, i, j, k;
  sqrteps = sqrt((double)1e-15);
  U22 = (double *)malloc(sizeof(double) * N * N);
  for (i = 0; i < N; ++i) {
    step = i * N;
    for (j = 0; j < N; ++j) {
      A[step + j] /= (dx[i] * dx[j]);
    }
  }
  for (i = 0; i < N; ++i) {
    U22[i] = A[i + i * N];
  }
  maxdiag = array_max(U22, N);
  mindiag = array_min(U22, N);
  maxposdiag = 0.0;
  if (maxdiag > maxposdiag) {
    maxposdiag = maxdiag;
  }
  for (i = 0; i < N; ++i) {
    U22[i] = 0.0;
  }
  u = 0.0;
  if (mindiag <= sqrteps * maxposdiag) {
    u = 2 * (maxposdiag - mindiag) * sqrteps - mindiag;
    maxdiag += u;
  }
  k = 0;
  for (i = 0; i < N; ++i) {
    for (j = 0; j < N; ++j) {
      if (j > i) {
        U22[k] = A[i * N + j];
        k++;
      }
    }
  }
  maxoffdiag = array_max_abs(U22, k);
  if (maxoffdiag * (1 + 2 * sqrteps) > maxdiag) {
    u += (maxoffdiag - maxdiag) + 2 * sqrteps * maxoffdiag;
    maxdiag = maxoffdiag * (1 + 2 * sqrteps);
  }
  if (maxdiag == 0) {
    u = 1;
    maxdiag = 1;
  }
  if (u > 0) {
    for (i = 0; i < N; ++i) {
      A[i * N + i] += u;
    }
  }
  if (maxdiag > maxoffdiag / N) {
    maxinp = sqrt(maxdiag);
  } else {
    maxinp = sqrt(maxoffdiag / N);
  }
  maxadd = cholmod(A, N, L, maxinp);
  if (maxadd > 0) {
    maxev = minev = A[0];
    for (i = 0; i < N; ++i) {
      offrow = 0.0;
      step = i * N;
      for (j = 0; j < i; ++j) {
        offrow += fabs(A[step + j]);
      }
      for (j = i + 1; j < N; ++j) {
        offrow += fabs(A[step + j]);
      }
      if (maxev < A[step + i] + offrow) {
        maxev = A[step + i] + offrow;
      }
      if (minev > A[step + i] - offrow) {
        minev = A[step + i] - offrow;
      }
    }
    sdd = (maxev - minev) * sqrteps - minev;
    if (sdd < 0) {
      sdd = 0;
    }
    if (maxadd > sdd) {
      u = sdd;
    } else {
      u = maxadd;
    }
    for (i = 0; i < N; ++i) {
      A[i * N + i] += u;
    }
  }
  maxadd = cholmod(A, N, L, 0.0);
  for (i = 0; i < N; ++i) {
    step = i * N;
    for (j = 0; j < N; ++j) {
      A[step + j] *= (dx[i] * dx[j]);
    }
  }
  for (i = 0; i < N; ++i) {
    step = i * N;
    for (j = 0; j < i; ++j) {
      L[step + j] *= dx[i];
    }
  }
  free(U22);
  return maxadd;
}


x86 assembly (clang, -O0):
<modelhess>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  0xe5e(%rip),%xmm0
call   1090 <sqrt@plt>
movsd  %xmm0,-0x30(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
movslq -0xc(%rbp),%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0x28(%rbp)
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1277 <modelhess+0xf7>
mov    -0x90(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x8c(%rbp)
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    125e <modelhess+0xde>
mov    -0x18(%rbp),%rax
movslq -0x90(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x94(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x94(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    11fa <modelhess+0x7a>
jmp    1263 <modelhess+0xe3>
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    11d1 <modelhess+0x51>
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12d2 <modelhess+0x152>
mov    -0x8(%rbp),%rax
mov    -0x90(%rbp),%ecx
mov    -0x90(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x90(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1281 <modelhess+0x101>
mov    -0x28(%rbp),%rdi
mov    -0xc(%rbp),%esi
call   1070 <array_max@plt>
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%rdi
mov    -0xc(%rbp),%esi
call   1050 <array_min@plt>
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
ucomisd -0x48(%rbp),%xmm0
jbe    1316 <modelhess+0x196>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1356 <modelhess+0x1d6>
mov    -0x28(%rbp),%rax
movslq -0x90(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1320 <modelhess+0x1a0>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     13ca <modelhess+0x24a>
movsd  -0x48(%rbp),%xmm1
subsd  -0x40(%rbp),%xmm1
movsd  0xc7f(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x98(%rbp)
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1485 <modelhess+0x305>
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    146c <modelhess+0x2ec>
mov    -0x94(%rbp),%eax
cmp    -0x90(%rbp),%eax
jle    1453 <modelhess+0x2d3>
mov    -0x8(%rbp),%rax
mov    -0x90(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x94(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x98(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x98(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
jmp    1458 <modelhess+0x2d8>
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    13f7 <modelhess+0x277>
jmp    1471 <modelhess+0x2f1>
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    13de <modelhess+0x25e>
mov    -0x28(%rbp),%rdi
mov    -0x98(%rbp),%esi
call   10a0 <array_max_abs@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
addsd  %xmm1,%xmm1
movsd  0xb61(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    1513 <modelhess+0x393>
movsd  -0x58(%rbp),%xmm1
subsd  -0x38(%rbp),%xmm1
movsd  0xb34(%rip),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
addsd  %xmm1,%xmm1
movsd  0xb0a(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1545 <modelhess+0x3c5>
jp     1545 <modelhess+0x3c5>
movsd  0xadd(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xad0(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    15af <modelhess+0x42f>
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    15aa <modelhess+0x42a>
movsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x90(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x90(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1561 <modelhess+0x3e1>
jmp    15af <modelhess+0x42f>
movsd  -0x38(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
cvtsi2sdl -0xc(%rbp),%xmm2
divsd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    15e0 <modelhess+0x460>
movsd  -0x38(%rbp),%xmm0
call   1090 <sqrt@plt>
movsd  %xmm0,-0x60(%rbp)
jmp    15f8 <modelhess+0x478>
movsd  -0x58(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
call   1090 <sqrt@plt>
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x20(%rbp),%rdx
movsd  -0x60(%rbp),%xmm0
call   1040 <cholmod@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    18b4 <modelhess+0x734>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    17d6 <modelhess+0x656>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x90(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x8c(%rbp)
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0x90(%rbp),%eax
jge    16c3 <modelhess+0x543>
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x94(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1060 <fabs@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    1671 <modelhess+0x4f1>
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1721 <modelhess+0x5a1>
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x94(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1060 <fabs@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    16d2 <modelhess+0x552>
movsd  -0x70(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x90(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x80(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    176f <modelhess+0x5ef>
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x90(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x78(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x90(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
subsd  -0x80(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jbe    17bd <modelhess+0x63d>
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x90(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    17c2 <modelhess+0x642>
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1640 <modelhess+0x4c0>
movsd  -0x70(%rbp),%xmm0
subsd  -0x78(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x78(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0x88(%rbp),%xmm0
ff
jbe    182d <modelhess+0x6ad>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
ucomisd -0x88(%rbp),%xmm0
ff
jbe    1852 <modelhess+0x6d2>
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x50(%rbp)
jmp    185c <modelhess+0x6dc>
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    18af <modelhess+0x72f>
movsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x90(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x90(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1866 <modelhess+0x6e6>
jmp    18b4 <modelhess+0x734>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x20(%rbp),%rdx
xorps  %xmm0,%xmm0
call   1040 <cholmod@plt>
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1978 <modelhess+0x7f8>
mov    -0x90(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x8c(%rbp)
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    195f <modelhess+0x7df>
mov    -0x18(%rbp),%rax
movslq -0x90(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x94(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x94(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    18ff <modelhess+0x77f>
jmp    1964 <modelhess+0x7e4>
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    18d6 <modelhess+0x756>
movl   $0x0,-0x90(%rbp)
mov    -0x90(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1a17 <modelhess+0x897>
mov    -0x90(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x8c(%rbp)
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0x90(%rbp),%eax
jge    19fe <modelhess+0x87e>
mov    -0x18(%rbp),%rax
movslq -0x90(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x8c(%rbp),%ecx
add    -0x94(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    19ab <modelhess+0x82b>
jmp    1a03 <modelhess+0x883>
mov    -0x90(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x90(%rbp)
jmp    1982 <modelhess+0x802>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
movsd  -0x68(%rbp),%xmm0
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<modelhess>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %rcx,0x30(%rsp)
mov    %rdx,%r13
mov    %esi,%ebx
mov    %rdi,%rbp
mov    %esi,%edi
imul   %esi,%edi
shl    $0x3,%edi
call   1090 <malloc@plt>
mov    %rax,%r14
mov    %ebx,%r8d
test   %ebx,%ebx
mov    %rax,0x40(%rsp)
mov    %r8,0x10(%rsp)
mov    %rbx,0x28(%rsp)
mov    %rbp,0x20(%rsp)
jle    148f <modelhess+0x2ff>
mov    %rbp,%r9
lea    0x0(,%r8,8),%r11
add    %r13,%r11
mov    %r8d,%r14d
and    $0xfffffffe,%r14d
lea    -0x2(%r14),%rax
mov    %rax,0x8(%rsp)
mov    %rax,%rcx
shr    %rcx
add    $0x1,%rcx
lea    -0x1(%r8),%r15
mov    %rcx,%rax
mov    %rcx,0x18(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x38(%rsp)
add    $0x10,%rbp
lea    0x0(,%r8,8),%r12
lea    0x8(%r9),%rdi
xor    %eax,%eax
jmp    1243 <modelhess+0xb3>
xchg   %ax,%ax
add    $0x1,%rax
add    %r12,%rbp
add    %r12,%rdi
cmp    %r8,%rax
je     13fd <modelhess+0x26d>
mov    %rax,%rcx
imul   %r8,%rcx
cmpl   $0x2,0x28(%rsp)
jae    12e0 <modelhess+0x150>
xor    %esi,%esi
mov    %rsi,%rdx
test   $0x1,%r8b
je     1288 <modelhess+0xf8>
movsd  0x0(%r13,%rax,8),%xmm0
mulsd  0x0(%r13,%rsi,8),%xmm0
add    %rsi,%rcx
movsd  (%r9,%rcx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r9,%rcx,8)
mov    %rsi,%rdx
or     $0x1,%rdx
cmp    %rsi,%r15
je     1230 <modelhess+0xa0>
nopl   (%rax)
movsd  0x0(%r13,%rax,8),%xmm0
mulsd  0x0(%r13,%rdx,8),%xmm0
movsd  -0x8(%rdi,%rdx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  (%rdi,%rdx,8),%xmm0
movsd  %xmm1,-0x8(%rdi,%rdx,8)
movsd  0x0(%r13,%rax,8),%xmm1
mulsd  0x8(%r13,%rdx,8),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rdi,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%r8
jne    1290 <modelhess+0x100>
jmp    1230 <modelhess+0xa0>
nopl   0x0(%rax,%rax,1)
lea    (%rcx,%r8,1),%rdx
lea    (%r9,%rcx,8),%rsi
lea    (%r9,%rdx,8),%rbx
lea    0x8(,%rax,8),%rdx
add    %r13,%rdx
lea    0x0(,%rax,8),%r8
add    %r13,%r8
cmp    %r11,%rsi
setb   %r9b
cmp    %r13,%rbx
seta   %r10b
cmp    %rdx,%rsi
setb   %dl
cmp    %rbx,%r8
setb   %bl
test   %r10b,%r9b
jne    13e1 <modelhess+0x251>
mov    $0x0,%esi
and    %bl,%dl
mov    0x10(%rsp),%r8
mov    0x20(%rsp),%r9
jne    1257 <modelhess+0xc7>
cmpq   $0x0,0x8(%rsp)
je     13f2 <modelhess+0x262>
movsd  0x0(%r13,%rax,8),%xmm0
unpcklpd %xmm0,%xmm0
mov    0x38(%rsp),%rsi
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movupd 0x0(%r13,%rdx,8),%xmm1
mulpd  %xmm0,%xmm1
movupd -0x10(%rbp,%rdx,8),%xmm2
divpd  %xmm1,%xmm2
movupd 0x0(%rbp,%rdx,8),%xmm1
movupd %xmm2,-0x10(%rbp,%rdx,8)
movupd 0x10(%r13,%rdx,8),%xmm2
mulpd  %xmm0,%xmm2
divpd  %xmm2,%xmm1
movupd %xmm1,0x0(%rbp,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rsi
jne    1360 <modelhess+0x1d0>
testb  $0x1,0x18(%rsp)
je     13d0 <modelhess+0x240>
movsd  0x0(%r13,%rax,8),%xmm0
unpcklpd %xmm0,%xmm0
movupd 0x0(%r13,%rdx,8),%xmm1
mulpd  %xmm0,%xmm1
add    %rcx,%rdx
movupd (%r9,%rdx,8),%xmm0
divpd  %xmm1,%xmm0
movupd %xmm0,(%r9,%rdx,8)
mov    %r14,%rsi
cmp    %r8,%r14
je     1230 <modelhess+0xa0>
jmp    1257 <modelhess+0xc7>
xor    %esi,%esi
mov    0x10(%rsp),%r8
mov    0x20(%rsp),%r9
jmp    1257 <modelhess+0xc7>
xor    %edx,%edx
testb  $0x1,0x18(%rsp)
jne    13a7 <modelhess+0x217>
jmp    13d0 <modelhess+0x240>
mov    0x28(%rsp),%rbx
test   %ebx,%ebx
mov    0x40(%rsp),%r14
jle    148f <modelhess+0x2ff>
xor    %ecx,%ecx
test   %r15,%r15
je     146a <modelhess+0x2da>
and    $0xfffffffe,%r8d
mov    0x28(%rsp),%rbp
lea    0x0(%rbp,%rbp,1),%esi
mov    $0x1,%edi
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movslq %ebp,%rbp
lea    (%rdi,%rbp,1),%rbx
movslq %ecx,%rcx
lea    (%rax,%rcx,1),%rdx
movsd  (%r9,%rdx,8),%xmm0
movsd  %xmm0,(%r14,%rax,8)
movsd  (%r9,%rbx,8),%xmm0
movsd  %xmm0,0x8(%r14,%rax,8)
add    $0x2,%rax
add    %esi,%ecx
add    $0x2,%rdi
add    %esi,%ebp
cmp    %rax,%r8
jne    1430 <modelhess+0x2a0>
jmp    146c <modelhess+0x2dc>
xor    %eax,%eax
testb  $0x1,0x10(%rsp)
mov    0x28(%rsp),%rbx
je     148f <modelhess+0x2ff>
mov    %eax,%ecx
imul   %ebx,%ecx
movslq %ecx,%rcx
add    %rax,%rcx
movsd  (%r9,%rcx,8),%xmm0
movsd  %xmm0,(%r14,%rax,8)
mov    %r14,%rdi
mov    %ebx,%esi
call   1080 <array_max@plt>
movsd  %xmm0,0x8(%rsp)
mov    %r14,%rdi
mov    %ebx,%esi
call   1060 <array_min@plt>
xorpd  %xmm1,%xmm1
movsd  0x8(%rsp),%xmm3
maxsd  %xmm1,%xmm3
test   %ebx,%ebx
mov    %rbx,%r15
jle    14f1 <modelhess+0x361>
mov    0x10(%rsp),%rax
lea    0x0(,%rax,8),%rdx
mov    %r14,%rdi
xor    %esi,%esi
movsd  %xmm0,0x18(%rsp)
movsd  %xmm3,0x38(%rsp)
call   1040 <memset@plt>
movsd  0x38(%rsp),%xmm3
movsd  0x18(%rsp),%xmm0
xorpd  %xmm1,%xmm1
movsd  0xb07(%rip),%xmm2
mulsd  %xmm3,%xmm2
ucomisd %xmm0,%xmm2
mov    0x20(%rsp),%r9
jb     1530 <modelhess+0x3a0>
subsd  %xmm0,%xmm3
addsd  %xmm3,%xmm3
mulsd  0xae8(%rip),%xmm3
subsd  %xmm0,%xmm3
movsd  0x8(%rsp),%xmm0
addsd  %xmm3,%xmm0
movsd  %xmm0,0x8(%rsp)
movapd %xmm3,%xmm1
test   %r15d,%r15d
mov    0x10(%rsp),%rdx
movsd  %xmm1,0x18(%rsp)
jle    15d8 <modelhess+0x448>
mov    %edx,%eax
and    $0xfffffffe,%eax
lea    0x8(%r9),%rcx
lea    0x0(,%rdx,8),%r8
xor    %ebx,%ebx
xor    %esi,%esi
jmp    156c <modelhess+0x3dc>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rbx
add    %r8,%rcx
cmp    %rdx,%rbx
je     15da <modelhess+0x44a>
xor    %ebp,%ebp
cmp    $0x1,%r15d
jne    15a9 <modelhess+0x419>
test   $0x1,%dl
je     1560 <modelhess+0x3d0>
cmp    %rbx,%rbp
jbe    1560 <modelhess+0x3d0>
mov    %rbx,%rdi
imul   %rdx,%rdi
add    %rdi,%rbp
movsd  (%r9,%rbp,8),%xmm0
movslq %esi,%rsi
movsd  %xmm0,(%r14,%rsi,8)
add    $0x1,%esi
jmp    1560 <modelhess+0x3d0>
nopl   0x0(%rax)
add    $0x2,%rbp
cmp    %rbp,%rax
je     1574 <modelhess+0x3e4>
cmp    %rbx,%rbp
jbe    15c0 <modelhess+0x430>
movsd  -0x8(%rcx,%rbp,8),%xmm0
movslq %esi,%rsi
movsd  %xmm0,(%r14,%rsi,8)
add    $0x1,%esi
cmp    %rbx,%rbp
jb     15a0 <modelhess+0x410>
movsd  (%rcx,%rbp,8),%xmm0
movslq %esi,%rsi
movsd  %xmm0,(%r14,%rsi,8)
add    $0x1,%esi
jmp    15a0 <modelhess+0x410>
xor    %esi,%esi
mov    %r14,%rdi
call   10b0 <array_max_abs@plt>
movapd %xmm0,%xmm1
movsd  0xa1a(%rip),%xmm0
mulsd  %xmm1,%xmm0
movapd %xmm1,%xmm2
movsd  0x8(%rsp),%xmm5
subsd  %xmm5,%xmm2
movsd  0xa08(%rip),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm2,%xmm3
movapd %xmm5,%xmm2
cmpltsd %xmm0,%xmm2
andpd  %xmm2,%xmm3
movsd  0x9f3(%rip),%xmm4
andnpd %xmm4,%xmm2
orpd   %xmm3,%xmm2
movsd  0x18(%rsp),%xmm4
addsd  %xmm2,%xmm4
maxsd  %xmm5,%xmm0
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
movsd  0x9d5(%rip),%xmm3
jne    1653 <modelhess+0x4c3>
jp     1653 <modelhess+0x4c3>
movapd %xmm3,%xmm4
mov    0x10(%rsp),%r8
mov    0x20(%rsp),%rdi
jne    1665 <modelhess+0x4d5>
jp     1665 <modelhess+0x4d5>
movapd %xmm3,%xmm0
ucomisd %xmm2,%xmm4
jbe    1709 <modelhess+0x579>
test   %r15d,%r15d
jle    1709 <modelhess+0x579>
xor    %eax,%eax
cmp    $0x1,%r15d
je     16ed <modelhess+0x55d>
lea    (%r15,%r15,1),%r9d
and    $0xfffffffe,%r8d
neg    %r8
mov    $0x1,%ebx
xor    %esi,%esi
mov    %r15d,%ebp
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movslq %ebp,%rbp
lea    (%rbx,%rbp,1),%r10
cltq
lea    (%rdx,%rax,1),%rcx
movsd  (%rdi,%rcx,8),%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,(%rdi,%rcx,8)
movsd  (%rdi,%r10,8),%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,(%rdi,%r10,8)
add    $0x2,%rdx
add    %r9d,%eax
add    $0x2,%rbx
add    %r9d,%ebp
add    $0xfffffffffffffffe,%rsi
cmp    %rsi,%r8
jne    16a0 <modelhess+0x510>
neg    %rsi
mov    %rsi,%rax
mov    0x10(%rsp),%r8
test   $0x1,%r8b
je     1709 <modelhess+0x579>
mov    %eax,%ecx
imul   %r15d,%ecx
movslq %ecx,%rcx
add    %rax,%rcx
addsd  (%rdi,%rcx,8),%xmm4
movsd  %xmm4,(%rdi,%rcx,8)
xorps  %xmm2,%xmm2
cvtsi2sd %r15d,%xmm2
divsd  %xmm2,%xmm1
maxsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1729 <modelhess+0x599>
sqrtsd %xmm0,%xmm0
jmp    1733 <modelhess+0x5a3>
call   10a0 <sqrt@plt>
mov    0x20(%rsp),%rdi
mov    %r15d,%esi
mov    0x30(%rsp),%rdx
call   1050 <cholmod@plt>
mov    0x20(%rsp),%rdi
ucomisd 0x8db(%rip),%xmm0
jbe    1a79 <modelhess+0x8e9>
movsd  %xmm0,0x60(%rsp)
movsd  (%rdi),%xmm2
test   %r15d,%r15d
jle    19b1 <modelhess+0x821>
lea    0x18(%rdi),%r14
mov    0x10(%rsp),%rax
lea    0x0(,%rax,8),%rcx
mov    %rcx,0x70(%rsp)
lea    0x8(,%rax,8),%rcx
mov    %rcx,0x68(%rsp)
xor    %ebx,%ebx
mov    %rax,%rbp
mov    %rdi,%r15
mov    %rdi,0x48(%rsp)
movapd %xmm2,%xmm3
jmp    180b <modelhess+0x67b>
nopl   0x0(%rax)
mov    0x58(%rsp),%rbp
mov    0x50(%rsp),%r15
movslq 0x18(%rsp),%rax
mov    0x38(%rsp),%rdx
add    %rax,%rdx
mov    0x20(%rsp),%rdi
movsd  (%rdi,%rdx,8),%xmm0
movapd %xmm2,%xmm1
addsd  %xmm0,%xmm1
maxsd  0x78(%rsp),%xmm1
subsd  %xmm2,%xmm0
minsd  0x80(%rsp),%xmm0
mov    0x70(%rsp),%rax
add    %rax,%r14
add    %rax,0x48(%rsp)
add    0x68(%rsp),%r15
add    $0xffffffffffffffff,%rbp
mov    %r12,%rbx
movapd %xmm0,%xmm2
movapd %xmm1,%xmm3
mov    0x10(%rsp),%r8
cmp    %r8,%r12
je     19be <modelhess+0x82e>
mov    %ebx,%eax
mov    %ebx,0x18(%rsp)
test   %rbx,%rbx
mov    %rbp,0x58(%rsp)
mov    %r15,0x50(%rsp)
movsd  %xmm2,0x80(%rsp)
mov    %rbx,0x38(%rsp)
movsd  %xmm3,0x78(%rsp)
je     1850 <modelhess+0x6c0>
lea    -0x1(%rbx),%rax
cmp    $0x3,%rax
jae    1860 <modelhess+0x6d0>
xorpd  %xmm2,%xmm2
xor    %r12d,%r12d
jmp    18fc <modelhess+0x76c>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm2,%xmm2
jmp    194b <modelhess+0x7bb>
nopl   0x0(%rax)
mov    %rbx,%r15
and    $0xfffffffffffffffc,%r15
xorpd  %xmm2,%xmm2
xor    %r12d,%r12d
xchg   %ax,%ax
movsd  %xmm2,0x8(%rsp)
movsd  -0x18(%r14,%r12,8),%xmm0
call   1070 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  -0x10(%r14,%r12,8),%xmm0
call   1070 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  -0x8(%r14,%r12,8),%xmm0
call   1070 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  (%r14,%r12,8),%xmm0
call   1070 <fabs@plt>
xorps  %xmm2,%xmm2
cvtsi2sd %rax,%xmm2
addsd  0x8(%rsp),%xmm2
add    $0x4,%r12
cmp    %r12,%r15
jne    1870 <modelhess+0x6e0>
test   $0x3,%bl
je     194b <modelhess+0x7bb>
mov    0x18(%rsp),%eax
mov    %eax,%r15d
and    $0x3,%r15d
mov    0x48(%rsp),%rax
lea    (%rax,%r12,8),%rbx
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
movsd  %xmm2,0x8(%rsp)
movsd  (%rbx,%rbp,8),%xmm0
call   1070 <fabs@plt>
movsd  0x8(%rsp),%xmm2
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  %xmm0,%xmm2
add    $0x1,%rbp
cmp    %rbp,%r15
jne    1920 <modelhess+0x790>
mov    0x18(%rsp),%eax
imul   0x28(%rsp),%eax
mov    %eax,0x18(%rsp)
mov    0x38(%rsp),%rax
lea    0x1(%rax),%r12
cmp    0x10(%rsp),%r12
jae    17a0 <modelhess+0x610>
mov    $0x1,%ebx
mov    0x58(%rsp),%rbp
mov    0x50(%rsp),%r15
nopl   0x0(%rax,%rax,1)
movsd  %xmm2,0x8(%rsp)
movsd  (%r15,%rbx,8),%xmm0
call   1070 <fabs@plt>
movsd  0x8(%rsp),%xmm2
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  %xmm0,%xmm2
add    $0x1,%rbx
cmp    %rbx,%rbp
jne    1980 <modelhess+0x7f0>
jmp    17aa <modelhess+0x61a>
movapd %xmm2,%xmm1
movapd %xmm2,%xmm0
mov    0x10(%rsp),%r8
mov    0x28(%rsp),%r15
test   %r15d,%r15d
jle    1e75 <modelhess+0xce5>
subsd  %xmm0,%xmm1
mulsd  0x628(%rip),%xmm1
subsd  %xmm0,%xmm1
xorpd  %xmm0,%xmm0
maxsd  %xmm1,%xmm0
minsd  0x60(%rsp),%xmm0
xor    %eax,%eax
cmp    $0x1,%r15d
je     1a5d <modelhess+0x8cd>
lea    (%r15,%r15,1),%r9d
and    $0xfffffffe,%r8d
neg    %r8
mov    $0x1,%ebx
xor    %esi,%esi
mov    %r15d,%ebp
xor    %edx,%edx
nopl   0x0(%rax)
movslq %ebp,%rbp
lea    (%rbx,%rbp,1),%r10
cltq
lea    (%rdx,%rax,1),%rcx
movsd  (%rdi,%rcx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
movsd  (%rdi,%r10,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%r10,8)
add    $0x2,%rdx
add    %r9d,%eax
add    $0x2,%rbx
add    %r9d,%ebp
add    $0xfffffffffffffffe,%rsi
cmp    %rsi,%r8
jne    1a10 <modelhess+0x880>
neg    %rsi
mov    %rsi,%rax
mov    0x10(%rsp),%r8
test   $0x1,%r8b
je     1a79 <modelhess+0x8e9>
mov    %eax,%ecx
imul   %r15d,%ecx
movslq %ecx,%rcx
add    %rax,%rcx
addsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,(%rdi,%rcx,8)
xorpd  %xmm0,%xmm0
mov    %r15d,%esi
mov    0x30(%rsp),%rdx
mov    %rdi,%rbp
call   1050 <cholmod@plt>
test   %r15d,%r15d
movsd  %xmm0,0x18(%rsp)
jle    1e8c <modelhess+0xcfc>
mov    0x10(%rsp),%rbx
lea    0x0(,%rbx,8),%rax
add    %r13,%rax
mov    %rax,0x8(%rsp)
mov    %ebx,%r9d
and    $0xfffffffc,%r9d
mov    %rbx,%r10
neg    %r10
lea    0x10(%rbp),%rsi
lea    0x0(,%rbx,8),%r11
mov    %r11,%r14
and    $0xffffffffffffffe0,%r14
lea    0x8(%rbp),%rax
xor    %ecx,%ecx
jmp    1af3 <modelhess+0x963>
nopl   0x0(%rax)
add    $0x1,%rcx
add    %r11,%rsi
add    %r11,%rax
cmp    %rbx,%rcx
je     1c5f <modelhess+0xacf>
mov    %rcx,%rdx
imul   %rbx,%rdx
cmp    $0x4,%r15d
jae    1b80 <modelhess+0x9f0>
xor    %edi,%edi
mov    %rdi,%rbp
test   $0x1,%bl
je     1b37 <modelhess+0x9a7>
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  0x0(%r13,%rdi,8),%xmm1
add    %rdi,%rdx
mov    0x20(%rsp),%rbp
mulsd  0x0(%rbp,%rdx,8),%xmm1
movsd  %xmm1,0x0(%rbp,%rdx,8)
mov    %rdi,%rbp
or     $0x1,%rbp
not    %rdi
cmp    %r10,%rdi
je     1ae0 <modelhess+0x950>
nop
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  0x0(%r13,%rbp,8),%xmm1
mulsd  -0x8(%rax,%rbp,8),%xmm1
movsd  %xmm1,-0x8(%rax,%rbp,8)
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  0x8(%r13,%rbp,8),%xmm1
mulsd  (%rax,%rbp,8),%xmm1
movsd  %xmm1,(%rax,%rbp,8)
add    $0x2,%rbp
cmp    %rbp,%rbx
jne    1b40 <modelhess+0x9b0>
jmp    1ae0 <modelhess+0x950>
lea    (%rdx,%rbx,1),%rdi
mov    0x20(%rsp),%rbp
lea    0x0(,%rdx,8),%rbx
add    %rbp,%rbx
lea    0x0(,%rdi,8),%rdi
add    %rbp,%rdi
lea    0x8(,%rcx,8),%rbp
add    %r13,%rbp
lea    0x0(,%rcx,8),%r15
add    %r13,%r15
cmp    0x8(%rsp),%rbx
setb   %r12b
cmp    %r13,%rdi
seta   %r8b
cmp    %rbp,%rbx
setb   %bl
cmp    %rdi,%r15
setb   %bpl
test   %r8b,%r12b
jne    1c4e <modelhess+0xabe>
mov    $0x0,%edi
and    %bpl,%bl
mov    0x10(%rsp),%rbx
mov    0x28(%rsp),%r15
jne    1b06 <modelhess+0x976>
movsd  0x0(%r13,%rcx,8),%xmm4
unpcklpd %xmm4,%xmm4
xor    %edi,%edi
nopl   0x0(%rax)
movupd 0x0(%r13,%rdi,1),%xmm1
movupd 0x10(%r13,%rdi,1),%xmm2
mulpd  %xmm4,%xmm1
mulpd  %xmm4,%xmm2
movupd -0x10(%rsi,%rdi,1),%xmm3
mulpd  %xmm1,%xmm3
movupd (%rsi,%rdi,1),%xmm1
mulpd  %xmm2,%xmm1
movupd %xmm3,-0x10(%rsi,%rdi,1)
movupd %xmm1,(%rsi,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%r14
jne    1c00 <modelhess+0xa70>
mov    %r9,%rdi
cmp    %rbx,%r9
je     1ae0 <modelhess+0x950>
jmp    1b06 <modelhess+0x976>
xor    %edi,%edi
mov    0x10(%rsp),%rbx
mov    0x28(%rsp),%r15
jmp    1b06 <modelhess+0x976>
test   %r15d,%r15d
jle    1e8c <modelhess+0xcfc>
movabs $0x7ffffffffffffffc,%r8
ff ff 7f
lea    0x8(,%rbx,8),%r9
mov    0x30(%rsp),%rsi
lea    0x30(%rsi),%rdx
lea    0x0(,%rbx,8),%r10
xor    %ecx,%ecx
mov    %rsi,%rbp
jmp    1cb5 <modelhess+0xb25>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rcx
add    %r10,%rdx
add    %r10,%rbp
cmp    0x10(%rsp),%rcx
je     1e8c <modelhess+0xcfc>
test   %rcx,%rcx
je     1ca0 <modelhess+0xb10>
mov    %rcx,%r14
mov    0x10(%rsp),%rdi
imul   %rdi,%r14
cmp    $0x4,%rcx
jae    1ce0 <modelhess+0xb50>
xor    %edi,%edi
jmp    1e00 <modelhess+0xc70>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rcx,%rax
imul   %rdi,%rax
lea    (%rsi,%rax,8),%rax
mov    %rsi,%rbx
lea    0x8(,%rcx,8),%rsi
add    %r13,%rsi
cmp    %rsi,%rax
jae    1d22 <modelhess+0xb92>
mov    %r9,%rax
imul   %rcx,%rax
add    %rbx,%rax
lea    0x0(,%rcx,8),%rsi
add    %r13,%rsi
cmp    %rsi,%rax
jbe    1d22 <modelhess+0xb92>
xor    %edi,%edi
mov    %rbx,%rsi
jmp    1e00 <modelhess+0xc70>
mov    %rcx,%rdi
and    %r8,%rdi
lea    -0x4(%rdi),%rax
mov    %rax,%r11
shr    $0x2,%r11
add    $0x1,%r11
test   %rax,%rax
je     1e64 <modelhess+0xcd4>
mov    %rcx,%rax
and    $0xfffffffffffffffc,%rax
add    $0xfffffffffffffffc,%rax
shr    $0x2,%rax
add    $0x1,%rax
and    $0xfffffffffffffffe,%rax
movsd  0x0(%r13,%rcx,8),%xmm5
unpcklpd %xmm5,%xmm5
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd -0x30(%rdx,%rsi,8),%xmm1
movupd -0x20(%rdx,%rsi,8),%xmm2
movupd -0x10(%rdx,%rsi,8),%xmm3
movupd (%rdx,%rsi,8),%xmm4
mulpd  %xmm5,%xmm1
mulpd  %xmm5,%xmm2
movupd %xmm1,-0x30(%rdx,%rsi,8)
movupd %xmm2,-0x20(%rdx,%rsi,8)
mulpd  %xmm5,%xmm3
mulpd  %xmm5,%xmm4
movupd %xmm3,-0x10(%rdx,%rsi,8)
movupd %xmm4,(%rdx,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rax
jne    1d70 <modelhess+0xbe0>
test   $0x1,%r11b
je     1def <modelhess+0xc5f>
movsd  0x0(%r13,%rcx,8),%xmm3
unpcklpd %xmm3,%xmm3
add    %r14,%rsi
mov    0x30(%rsp),%rax
movupd (%rax,%rsi,8),%xmm1
movupd 0x10(%rax,%rsi,8),%xmm2
mulpd  %xmm3,%xmm1
mulpd  %xmm3,%xmm2
movupd %xmm1,(%rax,%rsi,8)
movupd %xmm2,0x10(%rax,%rsi,8)
cmp    %rdi,%rcx
mov    0x30(%rsp),%rsi
je     1ca0 <modelhess+0xb10>
nopl   (%rax)
mov    %rdi,%rax
or     $0x1,%rax
test   $0x1,%cl
je     1e23 <modelhess+0xc93>
movsd  0x0(%r13,%rcx,8),%xmm1
add    %r14,%rdi
mulsd  (%rsi,%rdi,8),%xmm1
movsd  %xmm1,(%rsi,%rdi,8)
mov    %rax,%rdi
cmp    %rax,%rcx
je     1ca0 <modelhess+0xb10>
nopl   0x0(%rax)
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  0x0(%rbp,%rdi,8),%xmm1
movsd  %xmm1,0x0(%rbp,%rdi,8)
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  0x8(%rbp,%rdi,8),%xmm1
movsd  %xmm1,0x8(%rbp,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rcx
jne    1e30 <modelhess+0xca0>
jmp    1ca0 <modelhess+0xb10>
xor    %esi,%esi
test   $0x1,%r11b
jne    1dbe <modelhess+0xc2e>
jmp    1def <modelhess+0xc5f>
xorpd  %xmm0,%xmm0
mov    %r15d,%esi
mov    0x30(%rsp),%rdx
call   1050 <cholmod@plt>
movsd  %xmm0,0x18(%rsp)
mov    0x40(%rsp),%rdi
call   1030 <free@plt>
movsd  0x18(%rsp),%xmm0
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @modelhess(double* %0, i32 %1, double* nocapture %2, double* %3) {
%5 = bitcast double* %3 to i8*
%6 = shl i32 %1, 3
%7 = mul i32 %6, %1
%8 = tail call i64 @malloc(i32 %7)
%9 = inttoptr i64 %8 to double*
%10 = inttoptr i64 %8 to i8*
%11 = icmp sgt i32 %1, 0
br i1 %11, label %12, label %178
12: ; preds = %4
%13 = zext i32 %1 to i64
%14 = zext i32 %1 to i64
%15 = getelementptr double, double* %2, i64 %14
%16 = and i64 %14, 4294967294
%17 = add nsw i64 %16, -2
%18 = lshr exact i64 %17, 1
%19 = add nuw i64 %18, 1
%20 = add nsw i64 %14, -1
%21 = icmp ult i32 %1, 2
%22 = and i64 %14, 4294967294
%23 = and i64 %19, 1
%24 = icmp eq i64 %17, 0
%25 = and i64 %19, -2
%26 = icmp eq i64 %23, 0
%27 = icmp eq i64 %22, %14
%28 = and i64 %14, 1
%29 = icmp eq i64 %28, 0
br label %36
30: ; preds = %143
br i1 %11, label %31, label %178
31: ; preds = %30
%32 = and i64 %14, 1
%33 = icmp eq i64 %20, 0
br i1 %33, label %167, label %34
34: ; preds = %31
%35 = and i64 %14, 4294967294
br label %146
36: ; preds = %143, %12
%37 = phi i64 [ 0, %12 ], [ %144, %143 ]
%38 = mul i64 %37, %14
%39 = getelementptr double, double* %0, i64 %38
%40 = add i64 %38, %14
%41 = getelementptr double, double* %0, i64 %40
%42 = mul nsw i64 %37, %13
%43 = getelementptr inbounds double, double* %2, i64 %37
br i1 %21, label %107, label %44
44: ; preds = %36
%45 = add nuw i64 %37, 1
%46 = getelementptr double, double* %2, i64 %45
%47 = getelementptr double, double* %2, i64 %37
%48 = icmp ult double* %39, %15
%49 = icmp ugt double* %41, %2
%50 = and i1 %48, %49
%51 = icmp ult double* %39, %46
%52 = icmp ult double* %47, %41
%53 = and i1 %51, %52
%54 = or i1 %50, %53
br i1 %54, label %107, label %55
55: ; preds = %44
br i1 %24, label %90, label %56
56: ; preds = %55
%57 = load double, double* %43
%58 = insertelement <2 x double> poison, double %57, i64 0
%59 = shufflevector <2 x double> %58, <2 x double> poison, <2 x i32> zeroinitializer
%60 = load double, double* %43
%61 = insertelement <2 x double> poison, double %60, i64 0
%62 = shufflevector <2 x double> %61, <2 x double> poison, <2 x i32> zeroinitializer
br label %63
63: ; preds = %63, %56
%64 = phi i64 [ 0, %56 ], [ %87, %63 ]
%65 = phi i64 [ 0, %56 ], [ %88, %63 ]
%66 = getelementptr inbounds double, double* %2, i64 %64
%67 = bitcast double* %66 to <2 x double>*
%68 = load <2 x double>, <2 x double>* %67
%69 = fmul <2 x double> %59, %68
%70 = add nuw nsw i64 %64, %42
%71 = getelementptr inbounds double, double* %0, i64 %70
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = fdiv <2 x double> %73, %69
%75 = bitcast double* %71 to <2 x double>*
store <2 x double> %74, <2 x double>* %75
%76 = or i64 %64, 2
%77 = getelementptr inbounds double, double* %2, i64 %76
%78 = bitcast double* %77 to <2 x double>*
%79 = load <2 x double>, <2 x double>* %78
%80 = fmul <2 x double> %62, %79
%81 = add nuw nsw i64 %76, %42
%82 = getelementptr inbounds double, double* %0, i64 %81
%83 = bitcast double* %82 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = fdiv <2 x double> %84, %80
%86 = bitcast double* %82 to <2 x double>*
store <2 x double> %85, <2 x double>* %86
%87 = add nuw i64 %64, 4
%88 = add nuw i64 %65, 2
%89 = icmp eq i64 %88, %25
br i1 %89, label %90, label %63
90: ; preds = %63, %55
%91 = phi i64 [ 0, %55 ], [ %87, %63 ]
br i1 %26, label %106, label %92
92: ; preds = %90
%93 = load double, double* %43
%94 = insertelement <2 x double> poison, double %93, i64 0
%95 = shufflevector <2 x double> %94, <2 x double> poison, <2 x i32> zeroinitializer
%96 = getelementptr inbounds double, double* %2, i64 %91
%97 = bitcast double* %96 to <2 x double>*
%98 = load <2 x double>, <2 x double>* %97
%99 = fmul <2 x double> %95, %98
%100 = add nuw nsw i64 %91, %42
%101 = getelementptr inbounds double, double* %0, i64 %100
%102 = bitcast double* %101 to <2 x double>*
%103 = load <2 x double>, <2 x double>* %102
%104 = fdiv <2 x double> %103, %99
%105 = bitcast double* %101 to <2 x double>*
store <2 x double> %104, <2 x double>* %105
br label %106
106: ; preds = %90, %92
br i1 %27, label %143, label %107
107: ; preds = %44, %36, %106
%108 = phi i64 [ 0, %44 ], [ 0, %36 ], [ %22, %106 ]
br i1 %29, label %119, label %109
109: ; preds = %107
%110 = load double, double* %43
%111 = getelementptr inbounds double, double* %2, i64 %108
%112 = load double, double* %111
%113 = fmul double %110, %112
%114 = add nuw nsw i64 %108, %42
%115 = getelementptr inbounds double, double* %0, i64 %114
%116 = load double, double* %115
%117 = fdiv double %116, %113
store double %117, double* %115
%118 = or i64 %108, 1
br label %119
119: ; preds = %109, %107
%120 = phi i64 [ %108, %107 ], [ %118, %109 ]
%121 = icmp eq i64 %20, %108
br i1 %121, label %143, label %122
122: ; preds = %119, %122
%123 = phi i64 [ %141, %122 ], [ %120, %119 ]
%124 = load double, double* %43
%125 = getelementptr inbounds double, double* %2, i64 %123
%126 = load double, double* %125
%127 = fmul double %124, %126
%128 = add nuw nsw i64 %123, %42
%129 = getelementptr inbounds double, double* %0, i64 %128
%130 = load double, double* %129
%131 = fdiv double %130, %127
store double %131, double* %129
%132 = add nuw nsw i64 %123, 1
%133 = load double, double* %43
%134 = getelementptr inbounds double, double* %2, i64 %132
%135 = load double, double* %134
%136 = fmul double %133, %135
%137 = add nuw nsw i64 %132, %42
%138 = getelementptr inbounds double, double* %0, i64 %137
%139 = load double, double* %138
%140 = fdiv double %139, %136
store double %140, double* %138
%141 = add nuw nsw i64 %123, 2
%142 = icmp eq i64 %141, %14
br i1 %142, label %143, label %122
143: ; preds = %119, %122, %106
%144 = add nuw nsw i64 %37, 1
%145 = icmp eq i64 %144, %14
br i1 %145, label %30, label %36
146: ; preds = %146, %34
%147 = phi i64 [ 0, %34 ], [ %164, %146 ]
%148 = phi i64 [ 0, %34 ], [ %165, %146 ]
%149 = trunc i64 %147 to i32
%150 = mul nsw i32 %149, %1
%151 = sext i32 %150 to i64
%152 = add nsw i64 %147, %151
%153 = getelementptr inbounds double, double* %0, i64 %152
%154 = load double, double* %153
%155 = getelementptr inbounds double, double* %9, i64 %147
store double %154, double* %155
%156 = or i64 %147, 1
%157 = trunc i64 %156 to i32
%158 = mul nsw i32 %157, %1
%159 = sext i32 %158 to i64
%160 = add nsw i64 %156, %159
%161 = getelementptr inbounds double, double* %0, i64 %160
%162 = load double, double* %161
%163 = getelementptr inbounds double, double* %9, i64 %156
store double %162, double* %163
%164 = add nuw nsw i64 %147, 2
%165 = add i64 %148, 2
%166 = icmp eq i64 %165, %35
br i1 %166, label %167, label %146
167: ; preds = %146, %31
%168 = phi i64 [ 0, %31 ], [ %164, %146 ]
%169 = icmp eq i64 %32, 0
br i1 %169, label %178, label %170
170: ; preds = %167
%171 = trunc i64 %168 to i32
%172 = mul nsw i32 %171, %1
%173 = sext i32 %172 to i64
%174 = add nsw i64 %168, %173
%175 = getelementptr inbounds double, double* %0, i64 %174
%176 = load double, double* %175
%177 = getelementptr inbounds double, double* %9, i64 %168
store double %176, double* %177
br label %178
178: ; preds = %170, %167, %4, %30
%179 = tail call double @array_max(double* %9, i32 %1)
%180 = tail call double @array_min(double* %9, i32 %1)
%181 = fcmp ogt double %179, 0.000000e+00
%182 = select i1 %181, double %179, double 0.000000e+00
br i1 %11, label %183, label %186
183: ; preds = %178
%184 = zext i32 %1 to i64
%185 = shl nuw nsw i64 %184, 3
br label %186
186: ; preds = %183, %178
%187 = fmul double %182, 0x3E60FA3389D6EB40
%188 = fcmp ugt double %180, %187
br i1 %188, label %195, label %189
189: ; preds = %186
%190 = fsub double %182, %180
%191 = fmul double %190, 2.000000e+00
%192 = fneg double %180
%193 = tail call double @llvm.fmuladd.f64(double %191, double 0x3E60FA3389D6EB40, double %192)
%194 = fadd double %179, %193
br label %195
195: ; preds = %189, %186
%196 = phi double [ %193, %189 ], [ 0.000000e+00, %186 ]
%197 = phi double [ %194, %189 ], [ %179, %186 ]
br i1 %11, label %198, label %254
198: ; preds = %195
%199 = zext i32 %1 to i64
%200 = zext i32 %1 to i64
%201 = and i64 %200, 1
%202 = icmp eq i32 %1, 1
%203 = and i64 %200, 4294967294
%204 = icmp eq i64 %201, 0
br label %205
205: ; preds = %250, %198
%206 = phi i64 [ 0, %198 ], [ %252, %250 ]
%207 = phi i32 [ 0, %198 ], [ %251, %250 ]
%208 = mul nsw i64 %206, %199
br i1 %202, label %237, label %209
209: ; preds = %205, %232
%210 = phi i64 [ %234, %232 ], [ 0, %205 ]
%211 = phi i32 [ %233, %232 ], [ %207, %205 ]
%212 = phi i64 [ %235, %232 ], [ 0, %205 ]
%213 = icmp ugt i64 %210, %206
br i1 %213, label %214, label %221
214: ; preds = %209
%215 = add nuw nsw i64 %210, %208
%216 = getelementptr inbounds double, double* %0, i64 %215
%217 = load double, double* %216
%218 = sext i32 %211 to i64
%219 = getelementptr inbounds double, double* %9, i64 %218
store double %217, double* %219
%220 = add nsw i32 %211, 1
br label %221
221: ; preds = %209, %214
%222 = phi i32 [ %220, %214 ], [ %211, %209 ]
%223 = icmp ult i64 %210, %206
br i1 %223, label %232, label %224
224: ; preds = %221
%225 = or i64 %210, 1
%226 = add nuw nsw i64 %225, %208
%227 = getelementptr inbounds double, double* %0, i64 %226
%228 = load double, double* %227
%229 = sext i32 %222 to i64
%230 = getelementptr inbounds double, double* %9, i64 %229
store double %228, double* %230
%231 = add nsw i32 %222, 1
br label %232
232: ; preds = %224, %221
%233 = phi i32 [ %231, %224 ], [ %222, %221 ]
%234 = add nuw nsw i64 %210, 2
%235 = add i64 %212, 2
%236 = icmp eq i64 %235, %203
br i1 %236, label %237, label %209
237: ; preds = %232, %205
%238 = phi i32 [ undef, %205 ], [ %233, %232 ]
%239 = phi i64 [ 0, %205 ], [ %234, %232 ]
%240 = phi i32 [ %207, %205 ], [ %233, %232 ]
br i1 %204, label %250, label %241
241: ; preds = %237
%242 = icmp ugt i64 %239, %206
br i1 %242, label %243, label %250
243: ; preds = %241
%244 = add nuw nsw i64 %239, %208
%245 = getelementptr inbounds double, double* %0, i64 %244
%246 = load double, double* %245
%247 = sext i32 %240 to i64
%248 = getelementptr inbounds double, double* %9, i64 %247
store double %246, double* %248
%249 = add nsw i32 %240, 1
br label %250
250: ; preds = %241, %243, %237
%251 = phi i32 [ %238, %237 ], [ %249, %243 ], [ %240, %241 ]
%252 = add nuw nsw i64 %206, 1
%253 = icmp eq i64 %252, %200
br i1 %253, label %254, label %205
254: ; preds = %250, %195
%255 = phi i32 [ 0, %195 ], [ %251, %250 ]
%256 = tail call double @array_max_abs(double* %9, i32 %255)
%257 = fmul double %256, 0x3FF0000010FA338A
%258 = fcmp ogt double %257, %197
%259 = fsub double %256, %197
%260 = tail call double @llvm.fmuladd.f64(double %256, double 0x3E70FA3389D6EB40, double %259)
%261 = select i1 %258, double %260, double -0.000000e+00
%262 = fadd double %196, %261
%263 = select i1 %258, double %257, double %197
%264 = fcmp oeq double %263, 0.000000e+00
%265 = select i1 %264, double 1.000000e+00, double %262
%266 = select i1 %264, double 1.000000e+00, double %263
%267 = fcmp ogt double %265, 0.000000e+00
%268 = and i1 %267, %11
br i1 %268, label %269, label %307
269: ; preds = %254
%270 = zext i32 %1 to i64
%271 = and i64 %270, 1
%272 = icmp eq i32 %1, 1
br i1 %272, label %296, label %273
273: ; preds = %269
%274 = and i64 %270, 4294967294
br label %275
275: ; preds = %275, %273
%276 = phi i64 [ 0, %273 ], [ %293, %275 ]
%277 = phi i64 [ 0, %273 ], [ %294, %275 ]
%278 = trunc i64 %276 to i32
%279 = mul nsw i32 %278, %1
%280 = sext i32 %279 to i64
%281 = add nsw i64 %276, %280
%282 = getelementptr inbounds double, double* %0, i64 %281
%283 = load double, double* %282
%284 = fadd double %265, %283
store double %284, double* %282
%285 = or i64 %276, 1
%286 = trunc i64 %285 to i32
%287 = mul nsw i32 %286, %1
%288 = sext i32 %287 to i64
%289 = add nsw i64 %285, %288
%290 = getelementptr inbounds double, double* %0, i64 %289
%291 = load double, double* %290
%292 = fadd double %265, %291
store double %292, double* %290
%293 = add nuw nsw i64 %276, 2
%294 = add i64 %277, 2
%295 = icmp eq i64 %294, %274
br i1 %295, label %296, label %275
296: ; preds = %275, %269
%297 = phi i64 [ 0, %269 ], [ %293, %275 ]
%298 = icmp eq i64 %271, 0
br i1 %298, label %307, label %299
299: ; preds = %296
%300 = trunc i64 %297 to i32
%301 = mul nsw i32 %300, %1
%302 = sext i32 %301 to i64
%303 = add nsw i64 %297, %302
%304 = getelementptr inbounds double, double* %0, i64 %303
%305 = load double, double* %304
%306 = fadd double %265, %305
store double %306, double* %304
br label %307
307: ; preds = %299, %296, %254
%308 = sitofp i32 %1 to double
%309 = fdiv double %256, %308
%310 = fcmp ogt double %266, %309
%311 = select i1 %310, double %266, double %309
%312 = tail call double @sqrt(double %311)
%313 = tail call double @cholmod(double* %0, i32 %1, double* %3, double %312)
%314 = fcmp ogt double %313, 0.000000e+00
br i1 %314, label %315, label %466
315: ; preds = %307
%316 = load double, double* %0
br i1 %11, label %317, label %416
317: ; preds = %315
%318 = zext i32 %1 to i64
%319 = zext i32 %1 to i64
br label %320
320: ; preds = %317, %402
%321 = phi i64 [ 0, %317 ], [ %389, %402 ]
%322 = phi i64 [ 1, %317 ], [ %414, %402 ]
%323 = phi double [ %316, %317 ], [ %413, %402 ]
%324 = phi double [ %316, %317 ], [ %410, %402 ]
%325 = trunc i64 %321 to i32
%326 = mul nsw i64 %321, %318
%327 = mul nsw i32 %325, %1
%328 = icmp eq i64 %321, 0
br i1 %328, label %387, label %329
329: ; preds = %320
%330 = add nsw i64 %321, -1
%331 = and i64 %321, 3
%332 = icmp ult i64 %330, 3
br i1 %332, label %369, label %333
333: ; preds = %329
%334 = and i64 %321, 9223372036854775804
br label %335
335: ; preds = %335, %333
%336 = phi i64 [ 0, %333 ], [ %366, %335 ]
%337 = phi double [ 0.000000e+00, %333 ], [ %365, %335 ]
%338 = phi i64 [ 0, %333 ], [ %367, %335 ]
%339 = add nuw nsw i64 %336, %326
%340 = getelementptr inbounds double, double* %0, i64 %339
%341 = load double, double* %340
%342 = tail call i64 @fabs(double %341)
%343 = sitofp i64 %342 to double
%344 = fadd double %337, %343
%345 = or i64 %336, 1
%346 = add nuw nsw i64 %345, %326
%347 = getelementptr inbounds double, double* %0, i64 %346
%348 = load double, double* %347
%349 = tail call i64 @fabs(double %348)
%350 = sitofp i64 %349 to double
%351 = fadd double %344, %350
%352 = or i64 %336, 2
%353 = add nuw nsw i64 %352, %326
%354 = getelementptr inbounds double, double* %0, i64 %353
%355 = load double, double* %354
%356 = tail call i64 @fabs(double %355)
%357 = sitofp i64 %356 to double
%358 = fadd double %351, %357
%359 = or i64 %336, 3
%360 = add nuw nsw i64 %359, %326
%361 = getelementptr inbounds double, double* %0, i64 %360
%362 = load double, double* %361
%363 = tail call i64 @fabs(double %362)
%364 = sitofp i64 %363 to double
%365 = fadd double %358, %364
%366 = add nuw nsw i64 %336, 4
%367 = add i64 %338, 4
%368 = icmp eq i64 %367, %334
br i1 %368, label %369, label %335
369: ; preds = %335, %329
%370 = phi double [ undef, %329 ], [ %365, %335 ]
%371 = phi i64 [ 0, %329 ], [ %366, %335 ]
%372 = phi double [ 0.000000e+00, %329 ], [ %365, %335 ]
%373 = icmp eq i64 %331, 0
br i1 %373, label %387, label %374
374: ; preds = %369, %374
%375 = phi i64 [ %384, %374 ], [ %371, %369 ]
%376 = phi double [ %383, %374 ], [ %372, %369 ]
%377 = phi i64 [ %385, %374 ], [ 0, %369 ]
%378 = add nuw nsw i64 %375, %326
%379 = getelementptr inbounds double, double* %0, i64 %378
%380 = load double, double* %379
%381 = tail call i64 @fabs(double %380)
%382 = sitofp i64 %381 to double
%383 = fadd double %376, %382
%384 = add nuw nsw i64 %375, 1
%385 = add i64 %377, 1
%386 = icmp eq i64 %385, %331
br i1 %386, label %387, label %374
387: ; preds = %369, %374, %320
%388 = phi double [ 0.000000e+00, %320 ], [ %370, %369 ], [ %383, %374 ]
%389 = add nuw nsw i64 %321, 1
%390 = icmp ult i64 %389, %318
br i1 %390, label %391, label %402
391: ; preds = %387, %391
%392 = phi i64 [ %400, %391 ], [ %322, %387 ]
%393 = phi double [ %399, %391 ], [ %388, %387 ]
%394 = add nuw nsw i64 %392, %326
%395 = getelementptr inbounds double, double* %0, i64 %394
%396 = load double, double* %395
%397 = tail call i64 @fabs(double %396)
%398 = sitofp i64 %397 to double
%399 = fadd double %393, %398
%400 = add nuw nsw i64 %392, 1
%401 = icmp eq i64 %400, %319
br i1 %401, label %402, label %391
402: ; preds = %391, %387
%403 = phi double [ %388, %387 ], [ %399, %391 ]
%404 = sext i32 %327 to i64
%405 = add nsw i64 %321, %404
%406 = getelementptr inbounds double, double* %0, i64 %405
%407 = load double, double* %406
%408 = fadd double %403, %407
%409 = fcmp olt double %324, %408
%410 = select i1 %409, double %408, double %324
%411 = fsub double %407, %403
%412 = fcmp ogt double %323, %411
%413 = select i1 %412, double %411, double %323
%414 = add nuw nsw i64 %322, 1
%415 = icmp eq i64 %389, %319
br i1 %415, label %416, label %320
416: ; preds = %402, %315
%417 = phi double [ %316, %315 ], [ %410, %402 ]
%418 = phi double [ %316, %315 ], [ %413, %402 ]
%419 = fsub double %417, %418
%420 = fneg double %418
%421 = tail call double @llvm.fmuladd.f64(double %419, double 0x3E60FA3389D6EB40, double %420)
%422 = fcmp olt double %421, 0.000000e+00
%423 = select i1 %422, double 0.000000e+00, double %421
%424 = fcmp ogt double %313, %423
%425 = select i1 %424, double %423, double %313
br i1 %11, label %428, label %426
426: ; preds = %416
%427 = tail call double @cholmod(double* %0, i32 %1, double* %3, double 0.000000e+00)
br label %696
428: ; preds = %416
%429 = zext i32 %1 to i64
%430 = and i64 %429, 1
%431 = icmp eq i32 %1, 1
br i1 %431, label %455, label %432
432: ; preds = %428
%433 = and i64 %429, 4294967294
br label %434
434: ; preds = %434, %432
%435 = phi i64 [ 0, %432 ], [ %452, %434 ]
%436 = phi i64 [ 0, %432 ], [ %453, %434 ]
%437 = trunc i64 %435 to i32
%438 = mul nsw i32 %437, %1
%439 = sext i32 %438 to i64
%440 = add nsw i64 %435, %439
%441 = getelementptr inbounds double, double* %0, i64 %440
%442 = load double, double* %441
%443 = fadd double %425, %442
store double %443, double* %441
%444 = or i64 %435, 1
%445 = trunc i64 %444 to i32
%446 = mul nsw i32 %445, %1
%447 = sext i32 %446 to i64
%448 = add nsw i64 %444, %447
%449 = getelementptr inbounds double, double* %0, i64 %448
%450 = load double, double* %449
%451 = fadd double %425, %450
store double %451, double* %449
%452 = add nuw nsw i64 %435, 2
%453 = add i64 %436, 2
%454 = icmp eq i64 %453, %433
br i1 %454, label %455, label %434
455: ; preds = %434, %428
%456 = phi i64 [ 0, %428 ], [ %452, %434 ]
%457 = icmp eq i64 %430, 0
br i1 %457, label %466, label %458
458: ; preds = %455
%459 = trunc i64 %456 to i32
%460 = mul nsw i32 %459, %1
%461 = sext i32 %460 to i64
%462 = add nsw i64 %456, %461
%463 = getelementptr inbounds double, double* %0, i64 %462
%464 = load double, double* %463
%465 = fadd double %425, %464
store double %465, double* %463
br label %466
466: ; preds = %458, %455, %307
%467 = tail call double @cholmod(double* %0, i32 %1, double* %3, double 0.000000e+00)
br i1 %11, label %468, label %696
468: ; preds = %466
%469 = zext i32 %1 to i64
%470 = zext i32 %1 to i64
%471 = getelementptr double, double* %2, i64 %470
%472 = icmp ult i32 %1, 4
%473 = and i64 %470, 4294967292
%474 = icmp eq i64 %473, %470
%475 = and i64 %470, 1
%476 = icmp eq i64 %475, 0
%477 = sub nsw i64 0, %470
br label %484
478: ; preds = %570
br i1 %11, label %479, label %696
479: ; preds = %478
%480 = zext i32 %1 to i64
%481 = zext i32 %1 to i64
%482 = shl nuw nsw i64 %470, 3
%483 = add nuw nsw i64 %482, 8
br label %573
484: ; preds = %570, %468
%485 = phi i64 [ 0, %468 ], [ %571, %570 ]
%486 = mul i64 %485, %470
%487 = getelementptr double, double* %0, i64 %486
%488 = add i64 %486, %470
%489 = getelementptr double, double* %0, i64 %488
%490 = mul nsw i64 %485, %469
%491 = getelementptr inbounds double, double* %2, i64 %485
br i1 %472, label %533, label %492
492: ; preds = %484
%493 = add nuw i64 %485, 1
%494 = getelementptr double, double* %2, i64 %493
%495 = getelementptr double, double* %2, i64 %485
%496 = icmp ult double* %487, %471
%497 = icmp ugt double* %489, %2
%498 = and i1 %496, %497
%499 = icmp ult double* %487, %494
%500 = icmp ult double* %495, %489
%501 = and i1 %499, %500
%502 = or i1 %498, %501
br i1 %502, label %533, label %503
503: ; preds = %492
%504 = load double, double* %491
%505 = insertelement <2 x double> poison, double %504, i64 0
%506 = shufflevector <2 x double> %505, <2 x double> poison, <2 x i32> zeroinitializer
%507 = insertelement <2 x double> poison, double %504, i64 0
%508 = shufflevector <2 x double> %507, <2 x double> poison, <2 x i32> zeroinitializer
br label %509
509: ; preds = %509, %503
%510 = phi i64 [ 0, %503 ], [ %530, %509 ]
%511 = getelementptr inbounds double, double* %2, i64 %510
%512 = bitcast double* %511 to <2 x double>*
%513 = load <2 x double>, <2 x double>* %512
%514 = getelementptr inbounds double, double* %511, i64 2
%515 = bitcast double* %514 to <2 x double>*
%516 = load <2 x double>, <2 x double>* %515
%517 = fmul <2 x double> %506, %513
%518 = fmul <2 x double> %508, %516
%519 = add nuw nsw i64 %510, %490
%520 = getelementptr inbounds double, double* %0, i64 %519
%521 = bitcast double* %520 to <2 x double>*
%522 = load <2 x double>, <2 x double>* %521
%523 = getelementptr inbounds double, double* %520, i64 2
%524 = bitcast double* %523 to <2 x double>*
%525 = load <2 x double>, <2 x double>* %524
%526 = fmul <2 x double> %522, %517
%527 = fmul <2 x double> %525, %518
%528 = bitcast double* %520 to <2 x double>*
store <2 x double> %526, <2 x double>* %528
%529 = bitcast double* %523 to <2 x double>*
store <2 x double> %527, <2 x double>* %529
%530 = add nuw i64 %510, 4
%531 = icmp eq i64 %530, %473
br i1 %531, label %532, label %509
532: ; preds = %509
br i1 %474, label %570, label %533
533: ; preds = %492, %484, %532
%534 = phi i64 [ 0, %492 ], [ 0, %484 ], [ %473, %532 ]
%535 = xor i64 %534, -1
br i1 %476, label %546, label %536
536: ; preds = %533
%537 = load double, double* %491
%538 = getelementptr inbounds double, double* %2, i64 %534
%539 = load double, double* %538
%540 = fmul double %537, %539
%541 = add nuw nsw i64 %534, %490
%542 = getelementptr inbounds double, double* %0, i64 %541
%543 = load double, double* %542
%544 = fmul double %543, %540
store double %544, double* %542
%545 = or i64 %534, 1
br label %546
546: ; preds = %536, %533
%547 = phi i64 [ %534, %533 ], [ %545, %536 ]
%548 = icmp eq i64 %535, %477
br i1 %548, label %570, label %549
549: ; preds = %546, %549
%550 = phi i64 [ %568, %549 ], [ %547, %546 ]
%551 = load double, double* %491
%552 = getelementptr inbounds double, double* %2, i64 %550
%553 = load double, double* %552
%554 = fmul double %551, %553
%555 = add nuw nsw i64 %550, %490
%556 = getelementptr inbounds double, double* %0, i64 %555
%557 = load double, double* %556
%558 = fmul double %557, %554
store double %558, double* %556
%559 = add nuw nsw i64 %550, 1
%560 = load double, double* %491
%561 = getelementptr inbounds double, double* %2, i64 %559
%562 = load double, double* %561
%563 = fmul double %560, %562
%564 = add nuw nsw i64 %559, %490
%565 = getelementptr inbounds double, double* %0, i64 %564
%566 = load double, double* %565
%567 = fmul double %566, %563
store double %567, double* %565
%568 = add nuw nsw i64 %550, 2
%569 = icmp eq i64 %568, %470
br i1 %569, label %570, label %549
570: ; preds = %546, %549, %532
%571 = add nuw nsw i64 %485, 1
%572 = icmp eq i64 %571, %470
br i1 %572, label %478, label %484
573: ; preds = %479, %693
%574 = phi i64 [ 0, %479 ], [ %694, %693 ]
%575 = and i64 %574, 9223372036854775804
%576 = add nsw i64 %575, -4
%577 = lshr exact i64 %576, 2
%578 = add nuw nsw i64 %577, 1
%579 = mul i64 %574, %470
%580 = getelementptr double, double* %3, i64 %579
%581 = mul i64 %483, %574
%582 = getelementptr i8, i8* %5, i64 %581
%583 = getelementptr double, double* %2, i64 %574
%584 = bitcast double* %583 to i8*
%585 = add nuw i64 %574, 1
%586 = getelementptr double, double* %2, i64 %585
%587 = mul nsw i64 %574, %480
%588 = icmp eq i64 %574, 0
br i1 %588, label %693, label %589
589: ; preds = %573
%590 = getelementptr inbounds double, double* %2, i64 %574
%591 = icmp ult i64 %574, 4
br i1 %591, label %663, label %592
592: ; preds = %589
%593 = icmp ult double* %580, %586
%594 = icmp ugt i8* %582, %584
%595 = and i1 %593, %594
br i1 %595, label %663, label %596
596: ; preds = %592
%597 = and i64 %574, 9223372036854775804
%598 = and i64 %578, 1
%599 = icmp eq i64 %576, 0
br i1 %599, label %641, label %600
600: ; preds = %596
%601 = and i64 %578, 9223372036854775806
%602 = load double, double* %590
%603 = insertelement <2 x double> poison, double %602, i64 0
%604 = shufflevector <2 x double> %603, <2 x double> poison, <2 x i32> zeroinitializer
%605 = insertelement <2 x double> poison, double %602, i64 0
%606 = shufflevector <2 x double> %605, <2 x double> poison, <2 x i32> zeroinitializer
%607 = load double, double* %590
%608 = insertelement <2 x double> poison, double %607, i64 0
%609 = shufflevector <2 x double> %608, <2 x double> poison, <2 x i32> zeroinitializer
%610 = insertelement <2 x double> poison, double %607, i64 0
%611 = shufflevector <2 x double> %610, <2 x double> poison, <2 x i32> zeroinitializer
br label %612
612: ; preds = %612, %600
%613 = phi i64 [ 0, %600 ], [ %638, %612 ]
%614 = phi i64 [ 0, %600 ], [ %639, %612 ]
%615 = add nuw nsw i64 %613, %587
%616 = getelementptr inbounds double, double* %3, i64 %615
%617 = bitcast double* %616 to <2 x double>*
%618 = load <2 x double>, <2 x double>* %617
%619 = getelementptr inbounds double, double* %616, i64 2
%620 = bitcast double* %619 to <2 x double>*
%621 = load <2 x double>, <2 x double>* %620
%622 = fmul <2 x double> %604, %618
%623 = fmul <2 x double> %606, %621
%624 = bitcast double* %616 to <2 x double>*
store <2 x double> %622, <2 x double>* %624
%625 = bitcast double* %619 to <2 x double>*
store <2 x double> %623, <2 x double>* %625
%626 = or i64 %613, 4
%627 = add nuw nsw i64 %626, %587
%628 = getelementptr inbounds double, double* %3, i64 %627
%629 = bitcast double* %628 to <2 x double>*
%630 = load <2 x double>, <2 x double>* %629
%631 = getelementptr inbounds double, double* %628, i64 2
%632 = bitcast double* %631 to <2 x double>*
%633 = load <2 x double>, <2 x double>* %632
%634 = fmul <2 x double> %609, %630
%635 = fmul <2 x double> %611, %633
%636 = bitcast double* %628 to <2 x double>*
store <2 x double> %634, <2 x double>* %636
%637 = bitcast double* %631 to <2 x double>*
store <2 x double> %635, <2 x double>* %637
%638 = add nuw i64 %613, 8
%639 = add i64 %614, 2
%640 = icmp eq i64 %639, %601
br i1 %640, label %641, label %612
641: ; preds = %612, %596
%642 = phi i64 [ 0, %596 ], [ %638, %612 ]
%643 = icmp eq i64 %598, 0
br i1 %643, label %661, label %644
644: ; preds = %641
%645 = load double, double* %590
%646 = insertelement <2 x double> poison, double %645, i64 0
%647 = shufflevector <2 x double> %646, <2 x double> poison, <2 x i32> zeroinitializer
%648 = insertelement <2 x double> poison, double %645, i64 0
%649 = shufflevector <2 x double> %648, <2 x double> poison, <2 x i32> zeroinitializer
%650 = add nuw nsw i64 %642, %587
%651 = getelementptr inbounds double, double* %3, i64 %650
%652 = bitcast double* %651 to <2 x double>*
%653 = load <2 x double>, <2 x double>* %652
%654 = getelementptr inbounds double, double* %651, i64 2
%655 = bitcast double* %654 to <2 x double>*
%656 = load <2 x double>, <2 x double>* %655
%657 = fmul <2 x double> %647, %653
%658 = fmul <2 x double> %649, %656
%659 = bitcast double* %651 to <2 x double>*
store <2 x double> %657, <2 x double>* %659
%660 = bitcast double* %654 to <2 x double>*
store <2 x double> %658, <2 x double>* %660
br label %661
661: ; preds = %641, %644
%662 = icmp eq i64 %574, %597
br i1 %662, label %693, label %663
663: ; preds = %592, %589, %661
%664 = phi i64 [ 0, %592 ], [ 0, %589 ], [ %597, %661 ]
%665 = or i64 %664, 1
%666 = and i64 %574, 1
%667 = icmp eq i64 %666, 0
br i1 %667, label %675, label %668
668: ; preds = %663
%669 = load double, double* %590
%670 = add nuw nsw i64 %664, %587
%671 = getelementptr inbounds double, double* %3, i64 %670
%672 = load double, double* %671
%673 = fmul double %669, %672
store double %673, double* %671
%674 = or i64 %664, 1
br label %675
675: ; preds = %668, %663
%676 = phi i64 [ %664, %663 ], [ %674, %668 ]
%677 = icmp eq i64 %574, %665
br i1 %677, label %693, label %678
678: ; preds = %675, %678
%679 = phi i64 [ %691, %678 ], [ %676, %675 ]
%680 = load double, double* %590
%681 = add nuw nsw i64 %679, %587
%682 = getelementptr inbounds double, double* %3, i64 %681
%683 = load double, double* %682
%684 = fmul double %680, %683
store double %684, double* %682
%685 = add nuw nsw i64 %679, 1
%686 = load double, double* %590
%687 = add nuw nsw i64 %685, %587
%688 = getelementptr inbounds double, double* %3, i64 %687
%689 = load double, double* %688
%690 = fmul double %686, %689
store double %690, double* %688
%691 = add nuw nsw i64 %679, 2
%692 = icmp eq i64 %691, %574
br i1 %692, label %693, label %678
693: ; preds = %675, %678, %661, %573
%694 = add nuw nsw i64 %574, 1
%695 = icmp eq i64 %694, %481
br i1 %695, label %696, label %573
696: ; preds = %693, %466, %426, %478
%697 = phi double [ %467, %478 ], [ %467, %466 ], [ %427, %426 ], [ %467, %693 ]
%698 = tail call i32 @free(double* %9)
ret double %697
}