#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 DBL_MAX;
double RAND_MAX;
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
int /*<<< orphan*/ mir_xnyn(double **, int, int, int, int, double *, double *,
                            double *);
double rand();
double sqrt(double);
int /*<<< orphan*/ srand(int);
double estimate_mi(double *inp, int N, int ncols, int dimx, int dimy, int K,
                   double addnoise) {
  int i;
  double *scal;
  double *min;
  double *max;
  double *psi;
  int d;
  double mir = 0;
  double s, me;
  int BOX1;
  addnoise = -1;
  double **x;
  x = (double **)malloc(sizeof(*x) * ncols);
  for (int i = 0; i < ncols; ++i) {
    x[i] = (double *)malloc(sizeof(*x[i]) * N);
  }
  for (int i = 0; i < ncols; i++)
    for (int j = 0; j < N; j++)
      x[i][j] = inp[(j * ncols) + i];
  scal = (double *)malloc((dimx + dimy) * sizeof(double));
  min = (double *)malloc((dimx + dimy) * sizeof(double));
  max = (double *)malloc((dimx + dimy) * sizeof(double));
  for (d = 0; d < dimx + dimy; d++) {
    min[d] = DBL_MAX / 2;
    max[d] = -DBL_MAX / 2;
  }
  if (addnoise) {
    srand((dimx + dimy) * N * K * (int)(x[(dimx + dimy) / 2][N / 10]));
    if (addnoise == -1)
      for (d = 0; d < dimx + dimy; d++)
        for (i = 0; i < N; i++)
          x[d][i] += (1.0 * rand() / RAND_MAX) * 1e-8;
    else
      for (d = 0; d < dimx + dimy; d++)
        for (i = 0; i < N; i++)
          x[d][i] += (1.0 * rand() / RAND_MAX) * addnoise;
  }
  for (d = 0; d < dimx + dimy; d++) {
    me = s = 0;
    for (i = 0; i < N; i++)
      me += x[d][i];
    me /= N;
    for (i = 0; i < N; i++)
      s += (x[d][i] - me) * (x[d][i] - me);
    s /= (N - 1);
    s = sqrt(s);
    if (s == 0) {
      ;
    }
    for (i = 0; i < N; i++) {
      x[d][i] = (x[d][i] - me) / s;
      if (x[d][i] < min[d])
        min[d] = x[d][i];
      if (x[d][i] > max[d])
        max[d] = x[d][i];
    }
    for (i = 0; i < N; i++)
      x[d][i] = x[d][i] - min[d];
  }
  psi = (double *)malloc((N + 1) * sizeof(double));
  psi[1] = -(double).57721566490153;
  for (i = 1; i < N; i++)
    psi[i + 1] = psi[i] + 1 / (double)i;
  BOX1 = N - 5;
  for (d = 0; d < dimx + dimy; d++)
    scal[d] = BOX1 / (max[d] - min[d]);
  mir_xnyn(x, dimx, dimy, N, K, psi, scal, &mir);
  free(scal);
  free(min);
  free(max);
  free(psi);
  for (i = 0; i < ncols; ++i)
    free(x[i]);
  return mir;
}


x86 assembly (clang, -O0):
<estimate_mi>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9d,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xe6b(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movslq -0x10(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x80(%rbp)
movl   $0x0,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11fb <estimate_mi+0x9b>
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,%rdx
mov    -0x80(%rbp),%rax
movslq -0x84(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    11b7 <estimate_mi+0x57>
movl   $0x0,-0x88(%rbp)
mov    -0x88(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1291 <estimate_mi+0x131>
movl   $0x0,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1278 <estimate_mi+0x118>
mov    -0x8(%rbp),%rax
mov    -0x8c(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x88(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x80(%rbp),%rax
movslq -0x88(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    121e <estimate_mi+0xbe>
jmp    127d <estimate_mi+0x11d>
mov    -0x88(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x88(%rbp)
jmp    1205 <estimate_mi+0xa5>
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x38(%rbp)
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x40(%rbp)
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    1342 <estimate_mi+0x1e2>
mov    0x2ce3(%rip),%rax
mov    (%rax),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cvtsi2sd %eax,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
xor    %eax,%eax
mov    0x2cbf(%rip),%rcx
sub    (%rcx),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cvtsi2sd %eax,%xmm0
mov    -0x48(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    12dd <estimate_mi+0x17d>
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    135f <estimate_mi+0x1ff>
jp     135f <estimate_mi+0x1ff>
jmp    1518 <estimate_mi+0x3b8>
mov    -0x14(%rbp),%edi
add    -0x18(%rbp),%edi
imul   -0xc(%rbp),%edi
imul   -0x1c(%rbp),%edi
mov    -0x80(%rbp),%rax
mov    %rax,-0xa0(%rbp)
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0xa0(%rbp),%rax
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x98(%rbp)
mov    -0xc(%rbp),%eax
mov    $0xa,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x98(%rbp),%rax
movslq %ecx,%rcx
cvttsd2si (%rax,%rcx,8),%eax
imul   %eax,%edi
call   1050 <srand@plt>
movsd  -0x28(%rbp),%xmm0
movsd  0xc32(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    147a <estimate_mi+0x31a>
jp     147a <estimate_mi+0x31a>
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    1475 <estimate_mi+0x315>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1462 <estimate_mi+0x302>
mov    $0x0,%al
call   1080 <rand@plt>
movaps %xmm0,%xmm1
movsd  0xbed(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    0x2bc2(%rip),%rax
divsd  (%rax),%xmm0
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xbd1(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    13fd <estimate_mi+0x29d>
jmp    1467 <estimate_mi+0x307>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    13e5 <estimate_mi+0x285>
jmp    1513 <estimate_mi+0x3b3>
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    150e <estimate_mi+0x3ae>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    14fb <estimate_mi+0x39b>
mov    $0x0,%al
call   1080 <rand@plt>
movaps %xmm0,%xmm1
movsd  0xb51(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    0x2b26(%rip),%rax
divsd  (%rax),%xmm0
movsd  -0x28(%rbp),%xmm2
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1499 <estimate_mi+0x339>
jmp    1500 <estimate_mi+0x3a0>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1481 <estimate_mi+0x321>
jmp    1513 <estimate_mi+0x3b3>
jmp    1518 <estimate_mi+0x3b8>
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    179b <estimate_mi+0x63b>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1580 <estimate_mi+0x420>
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1547 <estimate_mi+0x3e7>
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  -0x70(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    15fa <estimate_mi+0x49a>
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x70(%rbp),%xmm0
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
subsd  -0x70(%rbp),%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    159a <estimate_mi+0x43a>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  -0x68(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
call   1070 <sqrt@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    163e <estimate_mi+0x4de>
jp     163e <estimate_mi+0x4de>
jmp    163e <estimate_mi+0x4de>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1730 <estimate_mi+0x5d0>
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x70(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x40(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    16d3 <estimate_mi+0x573>
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x54(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    171d <estimate_mi+0x5bd>
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1722 <estimate_mi+0x5c2>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1645 <estimate_mi+0x4e5>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1788 <estimate_mi+0x628>
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x54(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x80(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1737 <estimate_mi+0x5d7>
jmp    178d <estimate_mi+0x62d>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    151f <estimate_mi+0x3bf>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
movsd  0x852(%rip),%xmm0
movsd  %xmm0,0x8(%rax)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1818 <estimate_mi+0x6b8>
mov    -0x50(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm2
movsd  0x818(%rip),%xmm1
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x50(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    17ca <estimate_mi+0x66a>
mov    -0xc(%rbp),%eax
sub    $0x5,%eax
mov    %eax,-0x74(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    1877 <estimate_mi+0x717>
cvtsi2sdl -0x74(%rbp),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x40(%rbp),%rax
movslq -0x54(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1828 <estimate_mi+0x6c8>
mov    -0x80(%rbp),%rdi
mov    -0x14(%rbp),%esi
mov    -0x18(%rbp),%edx
mov    -0xc(%rbp),%ecx
mov    -0x1c(%rbp),%r8d
mov    -0x50(%rbp),%r9
mov    -0x38(%rbp),%r10
lea    -0x60(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
call   1040 <mir_xnyn@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    18f8 <estimate_mi+0x798>
mov    -0x80(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
call   1030 <free@plt>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    18cd <estimate_mi+0x76d>
movsd  -0x60(%rbp),%xmm0
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<estimate_mi>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %r9d,0x1c(%rsp)
mov    %r8,0x40(%rsp)
mov    %rcx,0x38(%rsp)
mov    %edx,%ebx
mov    %esi,%ebp
mov    %rdi,0x8(%rsp)
movq   $0x0,0x68(%rsp)
lea    0x0(,%rbx,8),%edi
call   1060 <malloc@plt>
mov    %rax,%r15
mov    %ebx,%eax
mov    %rax,0x10(%rsp)
mov    %rbp,%rax
mov    %rbp,(%rsp)
mov    %ebp,%r14d
mov    %rbx,0x30(%rsp)
test   %ebx,%ebx
jle    130a <estimate_mi+0x1aa>
mov    (%rsp),%rax
lea    0x0(,%rax,8),%r12d
xor    %ebx,%ebx
mov    0x10(%rsp),%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r12d,%edi
call   1060 <malloc@plt>
mov    %rax,(%r15,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%rbp
jne    11e0 <estimate_mi+0x80>
cmpl   $0x0,0x30(%rsp)
jle    130a <estimate_mi+0x1aa>
lea    -0x1(%r14),%r8
mov    %r14d,%ecx
and    $0x3,%ecx
mov    %r14d,%r11d
and    $0xfffffffc,%r11d
mov    0x10(%rsp),%rax
mov    %rax,%r13
shl    $0x4,%r13
mov    %rax,%rdi
shl    $0x5,%rdi
lea    0x0(,%rax,8),%rbp
lea    0x0(,%rbp,2),%rax
add    %rbp,%rax
xor    %r10d,%r10d
mov    0x8(%rsp),%r9
mov    0x10(%rsp),%r12
jmp    1263 <estimate_mi+0x103>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r10
add    $0x8,%r9
cmp    0x10(%rsp),%r10
je     130a <estimate_mi+0x1aa>
cmpl   $0x0,(%rsp)
jle    1250 <estimate_mi+0xf0>
mov    (%r15,%r10,8),%rbx
cmp    $0x3,%r8
jae    1280 <estimate_mi+0x120>
xor    %edx,%edx
jmp    12c8 <estimate_mi+0x168>
nopw   0x0(%rax,%rax,1)
mov    %r9,%rsi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rsi),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
movsd  (%rsi,%r12,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rdx,8)
movsd  (%rsi,%r13,1),%xmm0
movsd  %xmm0,0x10(%rbx,%rdx,8)
movsd  (%rsi,%rax,1),%xmm0
movsd  %xmm0,0x18(%rbx,%rdx,8)
add    $0x4,%rdx
add    %rdi,%rsi
cmp    %rdx,%r11
jne    1290 <estimate_mi+0x130>
test   %rcx,%rcx
je     1250 <estimate_mi+0xf0>
lea    (%rbx,%rdx,8),%rsi
imul   0x10(%rsp),%rdx
add    %r10,%rdx
mov    0x8(%rsp),%rbx
lea    (%rbx,%rdx,8),%rdx
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdx),%xmm0
movsd  %xmm0,(%rsi,%rbx,8)
add    $0x1,%rbx
add    %rbp,%rdx
cmp    %rbx,%rcx
jne    12f0 <estimate_mi+0x190>
jmp    1250 <estimate_mi+0xf0>
mov    0x40(%rsp),%rax
mov    0x38(%rsp),%rcx
lea    (%rax,%rcx,1),%r13d
lea    0x0(,%r13,8),%ebp
mov    %ebp,%edi
call   1060 <malloc@plt>
mov    %rax,0x88(%rsp)
mov    %ebp,%edi
call   1060 <malloc@plt>
mov    %rax,%r12
mov    %ebp,%edi
call   1060 <malloc@plt>
mov    %r13d,%ebx
test   %r13d,%r13d
jle    14f8 <estimate_mi+0x398>
mov    0x2c85(%rip),%rcx
mov    (%rcx),%edx
mov    %edx,%ecx
shr    $0x1f,%ecx
add    %edx,%ecx
sar    %ecx
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
neg    %ecx
cvtsi2sd %ecx,%xmm1
cmp    $0x4,%r13d
jb     1387 <estimate_mi+0x227>
lea    (%rax,%rbx,8),%rcx
cmp    %rcx,%r12
jae    1400 <estimate_mi+0x2a0>
lea    (%r12,%rbx,8),%rcx
cmp    %rcx,%rax
jae    1400 <estimate_mi+0x2a0>
xor    %ecx,%ecx
mov    %rcx,%rsi
not    %rsi
add    %rbx,%rsi
mov    %rbx,%rdx
and    $0x3,%rdx
je     13b5 <estimate_mi+0x255>
nopl   0x0(%rax,%rax,1)
movsd  %xmm0,(%r12,%rcx,8)
movsd  %xmm1,(%rax,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rdx
jne    13a0 <estimate_mi+0x240>
cmp    $0x3,%rsi
jb     14f8 <estimate_mi+0x398>
nop
movsd  %xmm0,(%r12,%rcx,8)
movsd  %xmm1,(%rax,%rcx,8)
movsd  %xmm0,0x8(%r12,%rcx,8)
movsd  %xmm1,0x8(%rax,%rcx,8)
movsd  %xmm0,0x10(%r12,%rcx,8)
movsd  %xmm1,0x10(%rax,%rcx,8)
movsd  %xmm0,0x18(%r12,%rcx,8)
movsd  %xmm1,0x18(%rax,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rbx
jne    13c0 <estimate_mi+0x260>
jmp    14f8 <estimate_mi+0x398>
mov    %ebx,%ecx
and    $0xfffffffc,%ecx
movapd %xmm0,%xmm2
unpcklpd %xmm0,%xmm2
movapd %xmm1,%xmm3
unpcklpd %xmm1,%xmm3
lea    -0x4(%rcx),%rsi
mov    %rsi,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %edx,%ebp
and    $0x3,%ebp
cmp    $0xc,%rsi
jae    1436 <estimate_mi+0x2d6>
xor    %esi,%esi
jmp    14b0 <estimate_mi+0x350>
and    $0xfffffffffffffffc,%rdx
xor    %esi,%esi
nopl   0x0(%rax)
movupd %xmm2,(%r12,%rsi,8)
movupd %xmm2,0x10(%r12,%rsi,8)
movupd %xmm3,(%rax,%rsi,8)
movupd %xmm3,0x10(%rax,%rsi,8)
movupd %xmm2,0x20(%r12,%rsi,8)
movupd %xmm2,0x30(%r12,%rsi,8)
movupd %xmm3,0x20(%rax,%rsi,8)
movupd %xmm3,0x30(%rax,%rsi,8)
movupd %xmm2,0x40(%r12,%rsi,8)
movupd %xmm2,0x50(%r12,%rsi,8)
movupd %xmm3,0x40(%rax,%rsi,8)
movupd %xmm3,0x50(%rax,%rsi,8)
movupd %xmm2,0x60(%r12,%rsi,8)
movupd %xmm2,0x70(%r12,%rsi,8)
movupd %xmm3,0x60(%rax,%rsi,8)
movupd %xmm3,0x70(%rax,%rsi,8)
add    $0x10,%rsi
add    $0xfffffffffffffffc,%rdx
jne    1440 <estimate_mi+0x2e0>
test   %rbp,%rbp
je     14ef <estimate_mi+0x38f>
lea    (%rax,%rsi,8),%rdx
add    $0x10,%rdx
lea    (%r12,%rsi,8),%rsi
add    $0x10,%rsi
shl    $0x5,%rbp
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movupd %xmm2,-0x10(%rsi,%rdi,1)
movupd %xmm2,(%rsi,%rdi,1)
movupd %xmm3,-0x10(%rdx,%rdi,1)
movupd %xmm3,(%rdx,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rbp
jne    14d0 <estimate_mi+0x370>
cmp    %rbx,%rcx
jne    1389 <estimate_mi+0x229>
mov    %r13d,%ebp
mov    (%rsp),%rdx
imul   %edx,%ebp
imul   0x1c(%rsp),%ebp
mov    %r13d,%ecx
shr    $0x1f,%ecx
add    %r13d,%ecx
sar    %ecx
movslq %ecx,%rcx
mov    (%r15,%rcx,8),%rcx
movslq %edx,%rdx
imul   $0x66666667,%rdx,%rdx
mov    %rdx,%rsi
shr    $0x3f,%rsi
sar    $0x22,%rdx
add    %esi,%edx
movslq %edx,%rdx
cvttsd2si (%rcx,%rdx,8),%edi
imul   %ebp,%edi
mov    %rax,0x20(%rsp)
mov    %rbx,0x8(%rsp)
call   1050 <srand@plt>
mov    %r13,0x28(%rsp)
test   %r13d,%r13d
jle    1a33 <estimate_mi+0x8d3>
mov    0x2a89(%rip),%r13
xor    %ebp,%ebp
jmp    157e <estimate_mi+0x41e>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rbp
mov    0x8(%rsp),%r8
cmp    %r8,%rbp
je     15c2 <estimate_mi+0x462>
cmpl   $0x0,(%rsp)
jle    1570 <estimate_mi+0x410>
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
xor    %eax,%eax
call   1080 <rand@plt>
movsd  0xa61(%rip),%xmm1
divsd  0x0(%r13),%xmm0
mov    (%r15,%rbp,8),%rax
mulsd  %xmm1,%xmm0
addsd  (%rax,%rbx,8),%xmm0
movsd  %xmm0,(%rax,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%r14
jne    1590 <estimate_mi+0x430>
jmp    1570 <estimate_mi+0x410>
cmpl   $0x0,0x28(%rsp)
mov    0x20(%rsp),%rdi
jle    1a33 <estimate_mi+0x8d3>
mov    (%rsp),%rax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0xa0(%rsp)
add    $0xffffffff,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xorpd  %xmm5,%xmm5
xorpd  %xmm1,%xmm1
movsd  %xmm0,0x98(%rsp)
divsd  %xmm0,%xmm1
movsd  %xmm1,0x60(%rsp)
lea    -0x1(%r14),%r9
mov    %r14d,%r10d
and    $0xfffffffc,%r10d
lea    -0x4(%r10),%rax
mov    %rax,0x80(%rsp)
shr    $0x2,%rax
add    $0x1,%rax
mov    %r14d,%ebx
and    $0x7,%ebx
mov    %r14d,%ebp
and    $0xfffffff8,%ebp
mov    %r14d,%r11d
and    $0x3,%r11d
movabs $0x7ffffffffffffffc,%rdx
ff ff 7f
or     $0x2,%rdx
mov    %rax,%rcx
mov    %rax,0x78(%rsp)
and    %rax,%rdx
mov    %rdx,0x70(%rsp)
xor    %r13d,%r13d
mov    %r9,0x58(%rsp)
mov    %r10,0x50(%rsp)
mov    %r11,0x48(%rsp)
jmp    167d <estimate_mi+0x51d>
nop
add    $0x1,%r13
cmp    %r8,%r13
je     1a33 <estimate_mi+0x8d3>
cmpl   $0x0,(%rsp)
jle    17d0 <estimate_mi+0x670>
mov    (%r15,%r13,8),%rax
xorpd  %xmm3,%xmm3
xor    %ecx,%ecx
cmp    $0x7,%r9
jb     16d8 <estimate_mi+0x578>
nopw   0x0(%rax,%rax,1)
addsd  (%rax,%rcx,8),%xmm3
addsd  0x8(%rax,%rcx,8),%xmm3
addsd  0x10(%rax,%rcx,8),%xmm3
addsd  0x18(%rax,%rcx,8),%xmm3
addsd  0x20(%rax,%rcx,8),%xmm3
addsd  0x28(%rax,%rcx,8),%xmm3
addsd  0x30(%rax,%rcx,8),%xmm3
addsd  0x38(%rax,%rcx,8),%xmm3
add    $0x8,%rcx
cmp    %rcx,%rbp
jne    16a0 <estimate_mi+0x540>
test   %rbx,%rbx
je     16fe <estimate_mi+0x59e>
lea    (%rax,%rcx,8),%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rax,%rcx,8),%xmm3
add    $0x1,%rcx
cmp    %rcx,%rbx
jne    16f0 <estimate_mi+0x590>
cmpl   $0x0,(%rsp)
jle    17d0 <estimate_mi+0x670>
divsd  0xa0(%rsp),%xmm3
mov    (%r15,%r13,8),%rax
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
cmp    $0x3,%r9
jb     1780 <estimate_mi+0x620>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rcx,8),%xmm1
movsd  0x8(%rax,%rcx,8),%xmm2
subsd  %xmm3,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rax,%rcx,8),%xmm1
subsd  %xmm3,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
movsd  0x18(%rax,%rcx,8),%xmm0
subsd  %xmm3,%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rcx
cmp    %rcx,%r10
jne    1730 <estimate_mi+0x5d0>
test   %r11,%r11
je     17aa <estimate_mi+0x64a>
lea    (%rax,%rcx,8),%rax
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rcx,8),%xmm1
subsd  %xmm3,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rcx
cmp    %rcx,%r11
jne    1790 <estimate_mi+0x630>
divsd  0x98(%rsp),%xmm0
ucomisd %xmm5,%xmm0
jb     180d <estimate_mi+0x6ad>
sqrtsd %xmm0,%xmm0
cmpl   $0x0,(%rsp)
jg     184b <estimate_mi+0x6eb>
jmp    1670 <estimate_mi+0x510>
nopl   0x0(%rax)
movsd  0x60(%rsp),%xmm0
ucomisd %xmm5,%xmm0
jae    1670 <estimate_mi+0x510>
movsd  0x60(%rsp),%xmm0
call   1070 <sqrt@plt>
mov    0x48(%rsp),%r11
mov    0x50(%rsp),%r10
mov    0x58(%rsp),%r9
xorpd  %xmm5,%xmm5
mov    0x8(%rsp),%r8
mov    0x20(%rsp),%rdi
jmp    1670 <estimate_mi+0x510>
movsd  %xmm3,0x90(%rsp)
call   1070 <sqrt@plt>
movsd  0x90(%rsp),%xmm3
mov    0x48(%rsp),%r11
mov    0x50(%rsp),%r10
mov    0x58(%rsp),%r9
xorpd  %xmm5,%xmm5
mov    0x8(%rsp),%r8
mov    0x20(%rsp),%rdi
cmpl   $0x0,(%rsp)
jle    1670 <estimate_mi+0x510>
lea    (%r12,%r13,8),%rcx
lea    (%r12,%r13,8),%rdx
add    $0x8,%rdx
mov    (%r15,%r13,8),%rax
xor    %esi,%esi
jmp    1869 <estimate_mi+0x709>
nop
add    $0x1,%rsi
cmp    %rsi,%r14
je     18b0 <estimate_mi+0x750>
movsd  (%rax,%rsi,8),%xmm1
subsd  %xmm3,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rsi,8)
movsd  (%r12,%r13,8),%xmm2
ucomisd %xmm1,%xmm2
jbe    1892 <estimate_mi+0x732>
movsd  %xmm1,(%r12,%r13,8)
movsd  (%rax,%rsi,8),%xmm1
ucomisd (%rdi,%r13,8),%xmm1
jbe    1860 <estimate_mi+0x700>
movsd  %xmm1,(%rdi,%r13,8)
jmp    1860 <estimate_mi+0x700>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cmpl   $0x0,(%rsp)
jle    1670 <estimate_mi+0x510>
mov    (%r15,%r13,8),%rax
cmpl   $0x4,(%rsp)
jb     18da <estimate_mi+0x77a>
cmp    %rdx,%rax
jae    1974 <estimate_mi+0x814>
lea    (%rax,%r14,8),%rdx
cmp    %rdx,%rcx
jae    1974 <estimate_mi+0x814>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %r14,%rdx
test   $0x3,%r14b
je     190a <estimate_mi+0x7aa>
mov    %r11,%rsi
xchg   %ax,%ax
movsd  (%rax,%rcx,8),%xmm0
subsd  (%r12,%r13,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    18f0 <estimate_mi+0x790>
cmp    $0x3,%rdx
jb     1670 <estimate_mi+0x510>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rax,%rcx,8),%xmm0
subsd  (%r12,%r13,8),%xmm0
movsd  0x8(%rax,%rcx,8),%xmm1
movsd  %xmm0,(%rax,%rcx,8)
subsd  (%r12,%r13,8),%xmm1
movsd  %xmm1,0x8(%rax,%rcx,8)
movsd  0x10(%rax,%rcx,8),%xmm0
subsd  (%r12,%r13,8),%xmm0
movsd  %xmm0,0x10(%rax,%rcx,8)
movsd  0x18(%rax,%rcx,8),%xmm0
subsd  (%r12,%r13,8),%xmm0
movsd  %xmm0,0x18(%rax,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%r14
jne    1920 <estimate_mi+0x7c0>
jmp    1670 <estimate_mi+0x510>
cmpq   $0x0,0x80(%rsp)
je     1a28 <estimate_mi+0x8c8>
movsd  (%r12,%r13,8),%xmm0
unpcklpd %xmm0,%xmm0
mov    0x70(%rsp),%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%rax,%rcx,8),%xmm1
movupd 0x10(%rax,%rcx,8),%xmm2
movupd 0x20(%rax,%rcx,8),%xmm3
movupd 0x30(%rax,%rcx,8),%xmm4
subpd  %xmm0,%xmm1
subpd  %xmm0,%xmm2
movupd %xmm1,(%rax,%rcx,8)
movupd %xmm2,0x10(%rax,%rcx,8)
subpd  %xmm0,%xmm3
subpd  %xmm0,%xmm4
movupd %xmm3,0x20(%rax,%rcx,8)
movupd %xmm4,0x30(%rax,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rdx
jne    19a0 <estimate_mi+0x840>
testb  $0x1,0x78(%rsp)
je     1a17 <estimate_mi+0x8b7>
movupd (%rax,%rcx,8),%xmm0
movupd 0x10(%rax,%rcx,8),%xmm1
movsd  (%r12,%r13,8),%xmm2
unpcklpd %xmm2,%xmm2
subpd  %xmm2,%xmm0
subpd  %xmm2,%xmm1
movupd %xmm0,(%rax,%rcx,8)
movupd %xmm1,0x10(%rax,%rcx,8)
mov    %r10,%rcx
cmp    %r14,%r10
je     1670 <estimate_mi+0x510>
jmp    18dc <estimate_mi+0x77c>
xor    %ecx,%ecx
testb  $0x1,0x78(%rsp)
jne    19ef <estimate_mi+0x88f>
jmp    1a17 <estimate_mi+0x8b7>
mov    (%rsp),%rbx
lea    0x8(,%rbx,8),%edi
call   1060 <malloc@plt>
mov    0x20(%rsp),%rbp
mov    %rax,%r9
movabs $0xbfe2788cfc6fb5ff,%rax
78 e2 bf
mov    %rax,0x8(%r9)
cmp    $0x2,%ebx
mov    0x10(%rsp),%r13
jl     1b05 <estimate_mi+0x9a5>
add    $0xffffffffffffffff,%r14
cmpl   $0x2,(%rsp)
jne    1a80 <estimate_mi+0x920>
mov    $0x1,%eax
movsd  0x58a(%rip),%xmm0
jmp    1ae1 <estimate_mi+0x981>
mov    %r14,%rcx
and    $0xfffffffffffffffe,%rcx
mov    $0x1,%eax
movsd  0x574(%rip),%xmm0
movsd  0x574(%rip),%xmm1
nopl   0x0(%rax)
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movapd %xmm1,%xmm3
divsd  %xmm2,%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm3,0x8(%r9,%rax,8)
lea    0x1(%rax),%rdx
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
movapd %xmm1,%xmm0
divsd  %xmm2,%xmm0
addsd  %xmm3,%xmm0
movsd  %xmm0,0x10(%r9,%rax,8)
add    $0x2,%rax
cmp    %rcx,%rdx
jne    1aa0 <estimate_mi+0x940>
test   $0x1,%r14b
je     1b05 <estimate_mi+0x9a5>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x51a(%rip),%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,0x8(%r9,%rax,8)
mov    0x28(%rsp),%rcx
test   %ecx,%ecx
mov    0x88(%rsp),%rbx
mov    0x8(%rsp),%rdi
jle    1c8d <estimate_mi+0xb2d>
mov    (%rsp),%rax
add    $0xfffffffb,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
cmp    $0x1,%ecx
jne    1bb2 <estimate_mi+0xa52>
xor    %eax,%eax
mov    %rax,%rcx
not    %rcx
test   $0x1,%dil
je     1b61 <estimate_mi+0xa01>
movsd  0x0(%rbp,%rax,8),%xmm1
subsd  (%r12,%rax,8),%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rax,8)
or     $0x1,%rax
add    %rdi,%rcx
je     1c8d <estimate_mi+0xb2d>
nopw   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rax,8),%xmm1
subsd  (%r12,%rax,8),%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rax,8)
movsd  0x8(%rbp,%rax,8),%xmm1
subsd  0x8(%r12,%rax,8),%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,0x8(%rbx,%rax,8)
add    $0x2,%rax
cmp    %rax,%rdi
jne    1b70 <estimate_mi+0xa10>
jmp    1c8d <estimate_mi+0xb2d>
lea    (%rbx,%rdi,8),%rax
lea    0x0(,%rdi,8),%rcx
add    %rbp,%rcx
lea    (%r12,%rdi,8),%rdx
cmp    %rcx,%rbx
setb   %sil
cmp    %rax,%rbp
setb   %r8b
cmp    %rdx,%rbx
setb   %cl
cmp    %rax,%r12
setb   %dl
xor    %eax,%eax
test   %r8b,%sil
jne    1b38 <estimate_mi+0x9d8>
and    %dl,%cl
jne    1b38 <estimate_mi+0x9d8>
mov    %edi,%eax
and    $0xfffffffe,%eax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rax),%rdx
mov    %rdx,%rcx
shr    %rcx
add    $0x1,%rcx
test   %rdx,%rdx
je     1d1a <estimate_mi+0xbba>
mov    %rcx,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %edx,%edx
nop
movupd 0x0(%rbp,%rdx,8),%xmm2
movupd (%r12,%rdx,8),%xmm3
subpd  %xmm3,%xmm2
movapd %xmm1,%xmm3
divpd  %xmm2,%xmm3
movupd %xmm3,(%rbx,%rdx,8)
movupd 0x10(%rbp,%rdx,8),%xmm2
movupd 0x10(%r12,%rdx,8),%xmm3
subpd  %xmm3,%xmm2
movapd %xmm1,%xmm3
divpd  %xmm2,%xmm3
movupd %xmm3,0x10(%rbx,%rdx,8)
add    $0x4,%rdx
add    $0xfffffffffffffffe,%rsi
jne    1c20 <estimate_mi+0xac0>
test   $0x1,%cl
je     1c84 <estimate_mi+0xb24>
movupd 0x0(%rbp,%rdx,8),%xmm2
movupd (%r12,%rdx,8),%xmm3
subpd  %xmm3,%xmm2
divpd  %xmm2,%xmm1
movupd %xmm1,(%rbx,%rdx,8)
cmp    %rdi,%rax
jne    1b38 <estimate_mi+0x9d8>
lea    0x68(%rsp),%rax
mov    %r15,%rdi
mov    0x38(%rsp),%rsi
mov    0x40(%rsp),%rdx
mov    (%rsp),%rcx
mov    %rbp,%r14
mov    0x1c(%rsp),%r8d
mov    %r9,%rbp
push   %rax
push   %rbx
call   1040 <mir_xnyn@plt>
add    $0x10,%rsp
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
cmpl   $0x0,0x30(%rsp)
jle    1d02 <estimate_mi+0xba2>
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%r15,%rbx,8),%rdi
call   1030 <free@plt>
add    $0x1,%rbx
cmp    %rbx,%r13
jne    1cf0 <estimate_mi+0xb90>
movsd  0x68(%rsp),%xmm0
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edx,%edx
test   $0x1,%cl
jne    1c6b <estimate_mi+0xb0b>
jmp    1c84 <estimate_mi+0xb24>

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

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