#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ CalcFloorCeil(int, double const, double *, double *,
                                 double *);
scalar_t__ CreateMatrixD(int, int);
int /*<<< orphan*/ DeleteMatrixD(double **);
scalar_t__ calloc(size_t, int);
double exp(double);
double const fabs(double);
int /*<<< orphan*/ free(double *);
int getLogFactorial(int);
double ldexp(double, int const);
int log(double);
double Pomeranz(int n, double x) {
  const double EPS = 1.0e-15;
  const int ENO = 350;
  const double RENO = ldexp(1.0, ENO);
  int coreno;
  const double t = n * x;
  double w, sum, minsum;
  int i, j, k, s;
  int r1, r2;
  int jlow, jup, klow, kup, kup0;
  double *A;
  double *Atflo;
  double *Atcei;
  double **V;
  double **H;
  A = (double *)calloc((size_t)(2 * n + 3), sizeof(double));
  Atflo = (double *)calloc((size_t)(2 * n + 3), sizeof(double));
  Atcei = (double *)calloc((size_t)(2 * n + 3), sizeof(double));
  V = (double **)CreateMatrixD(2, n + 2);
  H = (double **)CreateMatrixD(4, n + 2);
  CalcFloorCeil(n, t, A, Atflo, Atcei);
  for (j = 1; j <= n + 1; j++)
    V[0][j] = 0;
  for (j = 2; j <= n + 1; j++)
    V[1][j] = 0;
  V[1][1] = RENO;
  coreno = 1;
  H[0][0] = 1;
  w = 2.0 * A[2] / n;
  for (j = 1; j <= n + 1; j++)
    H[0][j] = w * H[0][j - 1] / j;
  H[1][0] = 1;
  w = (1.0 - 2.0 * A[2]) / n;
  for (j = 1; j <= n + 1; j++)
    H[1][j] = w * H[1][j - 1] / j;
  H[2][0] = 1;
  w = A[2] / n;
  for (j = 1; j <= n + 1; j++)
    H[2][j] = w * H[2][j - 1] / j;
  H[3][0] = 1;
  for (j = 1; j <= n + 1; j++)
    H[3][j] = 0;
  r1 = 0;
  r2 = 1;
  for (i = 2; i <= 2 * n + 2; i++) {
    jlow = 2 + (int)Atflo[i];
    if (jlow < 1)
      jlow = 1;
    jup = (int)Atcei[i];
    if (jup > n + 1)
      jup = n + 1;
    klow = 2 + (int)Atflo[i - 1];
    if (klow < 1)
      klow = 1;
    kup0 = (int)Atcei[i - 1];
    w = (A[i] - A[i - 1]) / n;
    s = -1;
    for (j = 0; j < 4; j++) {
      if (fabs(w - H[j][1]) <= EPS) {
        s = j;
        break;
      }
    }
    minsum = RENO;
    r1 = (r1 + 1) & 1;
    r2 = (r2 + 1) & 1;
    for (j = jlow; j <= jup; j++) {
      kup = kup0;
      if (kup > j)
        kup = j;
      sum = 0;
      for (k = kup; k >= klow; k--)
        sum += V[r1][k] * H[s][j - k];
      V[r2][j] = sum;
      if (sum < minsum)
        minsum = sum;
    }
    if (minsum < 1.0e-280) {
      for (j = jlow; j <= jup; j++)
        V[r2][j] *= RENO;
      coreno++;
    }
  }
  sum = V[r2][n + 1];
  free(A);
  free(Atflo);
  free(Atcei);
  DeleteMatrixD(H);
  DeleteMatrixD(V);
  w = getLogFactorial(n) - coreno * ENO * 0.69314718055994530941 + log(sum);
  if (w >= 0.)
    return 1.;
  return exp(w);
}


x86 assembly (clang, -O0):
<Pomeranz>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    %edi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  0xe4d(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x15e,-0x24(%rbp)
movsd  0xe31(%rip),%xmm0
mov    $0x15e,%edi
call   1040 <ldexp@plt>
movsd  %xmm0,-0x30(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
add    $0x3,%eax
movslq %eax,%rdi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x90(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
add    $0x3,%eax
movslq %eax,%rdi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0x98(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
add    $0x3,%eax
movslq %eax,%rdi
mov    $0x8,%esi
call   1070 <calloc@plt>
mov    %rax,-0xa0(%rbp)
mov    -0xc(%rbp),%esi
add    $0x2,%esi
mov    $0x2,%edi
call   1090 <CreateMatrixD@plt>
mov    %rax,-0xa8(%rbp)
mov    -0xc(%rbp),%esi
add    $0x2,%esi
mov    $0x4,%edi
call   1090 <CreateMatrixD@plt>
mov    %rax,-0xb0(%rbp)
mov    -0xc(%rbp),%edi
movsd  -0x40(%rbp),%xmm0
mov    -0x90(%rbp),%rsi
mov    -0x98(%rbp),%rdx
mov    -0xa0(%rbp),%rcx
call   10a0 <CalcFloorCeil@plt>
movl   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     12d0 <Pomeranz+0x130>
mov    -0xa8(%rbp),%rax
mov    (%rax),%rax
movslq -0x60(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    129b <Pomeranz+0xfb>
movl   $0x2,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     130d <Pomeranz+0x16d>
mov    -0xa8(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x60(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    12d7 <Pomeranz+0x137>
movsd  -0x30(%rbp),%xmm0
mov    -0xa8(%rbp),%rax
mov    0x8(%rax),%rax
movsd  %xmm0,0x8(%rax)
movl   $0x1,-0x34(%rbp)
mov    -0xb0(%rbp),%rax
mov    (%rax),%rax
movsd  0xcc5(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x90(%rbp),%rax
movsd  0xcc2(%rip),%xmm0
mulsd  0x10(%rax),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     13c0 <Pomeranz+0x220>
movsd  -0x48(%rbp),%xmm0
mov    -0xb0(%rbp),%rax
mov    (%rax),%rax
mov    -0x60(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0xb0(%rbp),%rax
mov    (%rax),%rax
movslq -0x60(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    1368 <Pomeranz+0x1c8>
mov    -0xb0(%rbp),%rax
mov    0x8(%rax),%rax
movsd  0xc2d(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x90(%rbp),%rax
movsd  0x10(%rax),%xmm1
addsd  %xmm1,%xmm1
movsd  0xc11(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     1462 <Pomeranz+0x2c2>
movsd  -0x48(%rbp),%xmm0
mov    -0xb0(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x60(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0xb0(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x60(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    1408 <Pomeranz+0x268>
mov    -0xb0(%rbp),%rax
mov    0x10(%rax),%rax
movsd  0xb8b(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x90(%rbp),%rax
movsd  0x10(%rax),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     14f4 <Pomeranz+0x354>
movsd  -0x48(%rbp),%xmm0
mov    -0xb0(%rbp),%rax
mov    0x10(%rax),%rax
mov    -0x60(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0xb0(%rbp),%rax
mov    0x10(%rax),%rax
movslq -0x60(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    149a <Pomeranz+0x2fa>
mov    -0xb0(%rbp),%rax
mov    0x18(%rax),%rax
movsd  0xaf9(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     1548 <Pomeranz+0x3a8>
mov    -0xb0(%rbp),%rax
mov    0x18(%rax),%rax
movslq -0x60(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    1512 <Pomeranz+0x372>
movl   $0x0,-0x6c(%rbp)
movl   $0x1,-0x70(%rbp)
movl   $0x2,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
mov    -0xc(%rbp),%ecx
shl    $0x1,%ecx
add    $0x2,%ecx
cmp    %ecx,%eax
jg     1825 <Pomeranz+0x685>
mov    -0x98(%rbp),%rax
movslq -0x5c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%eax
add    $0x2,%eax
mov    %eax,-0x74(%rbp)
cmpl   $0x1,-0x74(%rbp)
jge    1598 <Pomeranz+0x3f8>
movl   $0x1,-0x74(%rbp)
mov    -0xa0(%rbp),%rax
movslq -0x5c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%eax
mov    %eax,-0x78(%rbp)
mov    -0x78(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jle    15c5 <Pomeranz+0x425>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x78(%rbp)
mov    -0x98(%rbp),%rax
mov    -0x5c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cvttsd2si (%rax,%rcx,8),%eax
add    $0x2,%eax
mov    %eax,-0x7c(%rbp)
cmpl   $0x1,-0x7c(%rbp)
jge    15f1 <Pomeranz+0x451>
movl   $0x1,-0x7c(%rbp)
mov    -0xa0(%rbp),%rax
mov    -0x5c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cvttsd2si (%rax,%rcx,8),%eax
mov    %eax,-0x84(%rbp)
mov    -0x90(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x90(%rbp),%rax
mov    -0x5c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0xffffffff,-0x68(%rbp)
movl   $0x0,-0x60(%rbp)
cmpl   $0x4,-0x60(%rbp)
jge    16a8 <Pomeranz+0x508>
movsd  -0x48(%rbp),%xmm0
mov    -0xb0(%rbp),%rax
movslq -0x60(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
subsd  0x8(%rax),%xmm0
call   1080 <fabs@plt>
movaps %xmm0,%xmm1
movsd  0x988(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     1695 <Pomeranz+0x4f5>
mov    -0x60(%rbp),%eax
mov    %eax,-0x68(%rbp)
jmp    16a8 <Pomeranz+0x508>
jmp    169a <Pomeranz+0x4fa>
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    164d <Pomeranz+0x4ad>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
and    $0x1,%eax
mov    %eax,-0x6c(%rbp)
mov    -0x70(%rbp),%eax
add    $0x1,%eax
and    $0x1,%eax
mov    %eax,-0x70(%rbp)
mov    -0x74(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
cmp    -0x78(%rbp),%eax
jg     17b4 <Pomeranz+0x614>
mov    -0x84(%rbp),%eax
mov    %eax,-0x80(%rbp)
mov    -0x80(%rbp),%eax
cmp    -0x60(%rbp),%eax
jle    16f7 <Pomeranz+0x557>
mov    -0x60(%rbp),%eax
mov    %eax,-0x80(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x80(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x7c(%rbp),%eax
jl     1766 <Pomeranz+0x5c6>
mov    -0xa8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xb0(%rbp),%rax
movslq -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x60(%rbp),%ecx
sub    -0x64(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x64(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x64(%rbp)
jmp    1705 <Pomeranz+0x565>
movsd  -0x50(%rbp),%xmm0
mov    -0xa8(%rbp),%rax
movslq -0x70(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x60(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    17a1 <Pomeranz+0x601>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    17a6 <Pomeranz+0x606>
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    16d0 <Pomeranz+0x530>
movsd  0x864(%rip),%xmm0
ucomisd -0x58(%rbp),%xmm0
jbe    1812 <Pomeranz+0x672>
mov    -0x74(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
cmp    -0x78(%rbp),%eax
jg     1809 <Pomeranz+0x669>
movsd  -0x30(%rbp),%xmm0
mov    -0xa8(%rbp),%rax
movslq -0x70(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x60(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    17cd <Pomeranz+0x62d>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1817 <Pomeranz+0x677>
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    155d <Pomeranz+0x3bd>
mov    -0xa8(%rbp),%rax
movslq -0x70(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x90(%rbp),%rdi
call   1030 <free@plt>
mov    -0x98(%rbp),%rdi
call   1030 <free@plt>
mov    -0xa0(%rbp),%rdi
call   1030 <free@plt>
mov    -0xb0(%rbp),%rdi
call   1060 <DeleteMatrixD@plt>
mov    -0xa8(%rbp),%rdi
call   1060 <DeleteMatrixD@plt>
mov    -0xc(%rbp),%edi
call   1050 <getLogFactorial@plt>
cvtsi2sd %eax,%xmm1
imul   $0x15e,-0x34(%rbp),%eax
cvtsi2sd %eax,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0x75f(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
call   10c0 <log@plt>
movsd  -0xb8(%rbp),%xmm0
ff
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     190c <Pomeranz+0x76c>
movsd  0x6fe(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    191b <Pomeranz+0x77b>
movsd  -0x48(%rbp),%xmm0
call   10b0 <exp@plt>
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Pomeranz>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
movsd  %xmm0,(%rsp)
mov    %edi,%ebp
movsd  0xe40(%rip),%xmm0
mov    $0x15e,%edi
call   1040 <ldexp@plt>
cvtsi2sd %ebp,%xmm1
movaps %xmm0,0x70(%rsp)
movsd  %xmm1,0x30(%rsp)
movsd  (%rsp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%rsp)
lea    0x3(%rbp,%rbp,1),%eax
movslq %eax,%rbx
mov    %rbx,%rdi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,%r14
mov    %rbx,%rdi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,%r15
mov    %rbx,%rdi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,%r13
lea    0x2(%rbp),%r12d
mov    $0x2,%edi
mov    %r12d,%esi
call   1090 <CreateMatrixD@plt>
mov    %rax,%rbx
mov    $0x4,%edi
mov    %r12d,%esi
call   1090 <CreateMatrixD@plt>
mov    %rax,0x10(%rsp)
mov    %ebp,%edi
movsd  (%rsp),%xmm0
mov    %r14,%rsi
mov    %r15,%rdx
mov    %r13,0x28(%rsp)
mov    %r13,%rcx
call   10a0 <CalcFloorCeil@plt>
lea    0x1(%rbp),%r9d
test   %ebp,%ebp
mov    %rbp,0x38(%rsp)
mov    %rbx,(%rsp)
mov    %r9d,0xc(%rsp)
js     12db <Pomeranz+0x13b>
mov    (%rbx),%rdi
add    $0x8,%rdi
mov    %r9d,%edx
shl    $0x3,%rdx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0xc(%rsp),%r9d
mov    0x38(%rsp),%rbp
test   %ebp,%ebp
je     12db <Pomeranz+0x13b>
mov    0x8(%rbx),%rbx
lea    0x10(%rbx),%rdi
mov    %ebp,%edx
shl    $0x3,%rdx
xor    %esi,%esi
call   1060 <memset@plt>
mov    %rbx,%r10
mov    0xc(%rsp),%r9d
mov    0x38(%rsp),%rbp
movapd 0x70(%rsp),%xmm0
movsd  %xmm0,0x8(%rbx)
mov    0x10(%rsp),%r8
mov    (%r8),%rsi
movabs $0x3ff0000000000000,%rbx
mov    %rbx,(%rsi)
jmp    1308 <Pomeranz+0x168>
mov    0x8(%rbx),%r10
movapd 0x70(%rsp),%xmm0
movsd  %xmm0,0x8(%r10)
mov    0x10(%rsp),%r8
mov    (%r8),%rsi
movabs $0x3ff0000000000000,%rbx
mov    %rbx,(%rsi)
test   %ebp,%ebp
js     1915 <Pomeranz+0x775>
movsd  0x30(%rsp),%xmm3
movsd  0xcea(%rip),%xmm2
movsd  0x10(%r14),%xmm0
addsd  %xmm0,%xmm0
divsd  %xmm3,%xmm0
add    $0xffffffffffffffff,%r12
test   %ebp,%ebp
je     1901 <Pomeranz+0x761>
mov    %r12,%rcx
and    $0xfffffffffffffffe,%rcx
mov    $0x1,%eax
movsd  0xcbc(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mulsd  %xmm0,%xmm2
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm2
movsd  %xmm2,(%rsi,%rax,8)
mulsd  %xmm0,%xmm2
lea    0x1(%rax),%rdx
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
divsd  %xmm1,%xmm2
movsd  %xmm2,0x8(%rsi,%rax,8)
add    $0x2,%rax
cmp    %rcx,%rdx
jne    1350 <Pomeranz+0x1b0>
test   $0x1,%r12b
je     13a0 <Pomeranz+0x200>
mulsd  %xmm2,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rsi,%rax,8)
mov    0x8(%r8),%r12
mov    %rbx,(%r12)
test   %ebp,%ebp
js     191c <Pomeranz+0x77c>
movsd  0x10(%r14),%xmm2
addsd  %xmm2,%xmm2
movsd  0xc3e(%rip),%xmm1
movapd %xmm1,%xmm0
subsd  %xmm2,%xmm0
divsd  %xmm3,%xmm0
cmp    $0x2,%r9d
mov    $0x1,%eax
cmovge %r9d,%eax
mov    %rax,%rcx
add    $0xffffffffffffffff,%rcx
mov    %rsi,0x20(%rsp)
je     1949 <Pomeranz+0x7a9>
mov    %eax,%esi
and    $0xfffffffe,%esi
mov    $0x1,%edx
movsd  0xc01(%rip),%xmm1
nop
mulsd  %xmm0,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,(%r12,%rdx,8)
mulsd  %xmm0,%xmm1
lea    0x1(%rdx),%rdi
xorps  %xmm2,%xmm2
cvtsi2sd %edi,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,0x8(%r12,%rdx,8)
add    $0x2,%rdx
cmp    %rsi,%rdi
jne    1400 <Pomeranz+0x260>
test   $0x1,%al
je     1451 <Pomeranz+0x2b1>
mulsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rdx,8)
mov    0x10(%r8),%rdi
mov    %rbx,(%rdi)
test   %ebp,%ebp
js     1923 <Pomeranz+0x783>
add    $0x1,%rax
movsd  0x10(%r14),%xmm0
divsd  %xmm3,%xmm0
add    $0xffffffffffffffff,%rax
test   %rcx,%rcx
je     195b <Pomeranz+0x7bb>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
mov    $0x1,%ecx
movsd  0xb71(%rip),%xmm1
nop
mulsd  %xmm0,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
mulsd  %xmm0,%xmm1
lea    0x1(%rcx),%rsi
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,0x8(%rdi,%rcx,8)
add    $0x2,%rcx
cmp    %rdx,%rsi
jne    1490 <Pomeranz+0x2f0>
test   $0x1,%al
je     14de <Pomeranz+0x33e>
mulsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rdi,%rcx,8)
mov    0x18(%r8),%rax
mov    %rbx,(%rax)
test   %ebp,%ebp
js     1975 <Pomeranz+0x7d5>
mov    %rdi,0x88(%rsp)
lea    0x0(%rbp,%rbp,1),%ebx
mov    %rax,0x80(%rsp)
lea    0x8(%rax),%rdi
cmp    $0x2,%r9d
movl   $0x1,0x1c(%rsp)
mov    $0x1,%edx
cmovge %r9d,%edx
shl    $0x3,%rdx
xor    %esi,%esi
call   1060 <memset@plt>
mov    0x20(%rsp),%r8
mov    0xc(%rsp),%r9d
movapd 0x70(%rsp),%xmm1
movsd  0x30(%rsp),%xmm3
add    $0x2,%ebx
cmp    $0x3,%ebx
mov    $0x2,%eax
cmovge %ebx,%eax
xor    %ebx,%ebx
or     $0x1,%eax
mov    %rax,0xa0(%rsp)
movapd %xmm1,%xmm0
unpcklpd %xmm1,%xmm0
mov    $0x2,%ecx
movapd 0xab4(%rip),%xmm9
movsd  0xa94(%rip),%xmm2
movsd  0xa93(%rip),%xmm8
mov    $0x1,%edi
mov    0x28(%rsp),%r11
mov    %r14,0x60(%rsp)
mov    %r15,0x58(%rsp)
mov    %r12,0x50(%rsp)
jmp    15bb <Pomeranz+0x41b>
nopl   0x0(%rax,%rax,1)
addl   $0x1,0x1c(%rsp)
mov    0x18(%rsp),%edi
add    $0x1,%rcx
cmp    0xa0(%rsp),%rcx
je     1934 <Pomeranz+0x794>
cvttsd2si (%r15,%rcx,8),%r10d
test   %r10d,%r10d
mov    $0xffffffff,%eax
cmovs  %eax,%r10d
cvttsd2si (%r11,%rcx,8),%edx
cmp    %edx,%r9d
cmovl  %r9d,%edx
mov    %rdx,0x68(%rsp)
cvttsd2si -0x8(%r15,%rcx,8),%ebp
test   %ebp,%ebp
cmovs  %eax,%ebp
movsd  (%r14,%rcx,8),%xmm4
subsd  -0x8(%r14,%rcx,8),%xmm4
divsd  %xmm3,%xmm4
movapd %xmm4,%xmm5
subsd  0x8(%r8),%xmm5
lea    0x2(%r10),%eax
andpd  %xmm9,%xmm5
ucomisd %xmm5,%xmm2
jae    1680 <Pomeranz+0x4e0>
movapd %xmm4,%xmm5
subsd  0x8(%r12),%xmm5
andpd  %xmm9,%xmm5
mov    $0x1,%esi
ucomisd %xmm5,%xmm2
jae    1682 <Pomeranz+0x4e2>
movapd %xmm4,%xmm5
mov    %eax,%r13d
mov    0x88(%rsp),%rax
subsd  0x8(%rax),%xmm5
mov    %r13d,%eax
andpd  %xmm9,%xmm5
mov    $0x2,%esi
ucomisd %xmm5,%xmm2
jae    1682 <Pomeranz+0x4e2>
mov    0x80(%rsp),%rax
subsd  0x8(%rax),%xmm4
mov    %r13d,%eax
andpd  %xmm9,%xmm4
mov    $0x3,%esi
ucomisd %xmm4,%xmm2
jae    1682 <Pomeranz+0x4e2>
mov    $0xffffffffffffffff,%rsi
jmp    1682 <Pomeranz+0x4e2>
xor    %esi,%esi
not    %ebx
and    $0x1,%ebx
not    %edi
and    $0x1,%edi
mov    0x68(%rsp),%rdx
cmp    %edx,%eax
mov    %edi,0x18(%rsp)
jle    16b0 <Pomeranz+0x510>
ucomisd %xmm1,%xmm8
ja     15a0 <Pomeranz+0x400>
jmp    15a5 <Pomeranz+0x405>
nopl   0x0(%rax)
cvttsd2si -0x8(%r11,%rcx,8),%r14d
add    $0x2,%ebp
mov    %rbx,0x40(%rsp)
mov    %ebx,%r9d
mov    %edi,%ebx
mov    (%rsp),%rdi
mov    %rbx,0x98(%rsp)
mov    (%rdi,%rbx,8),%r11
mov    %ebp,%r8d
mov    %eax,0x4c(%rsp)
movslq %eax,%rdi
movslq %r14d,%r15
lea    0x1(%rdx),%eax
mov    %rdi,%rdx
mov    %rdi,0x90(%rsp)
mov    %rdi,%r14
movapd %xmm1,%xmm4
jmp    1717 <Pomeranz+0x577>
nopl   0x0(%rax,%rax,1)
movsd  %xmm5,(%r11,%r14,8)
minsd  %xmm4,%xmm5
add    $0x1,%r14
movapd %xmm5,%xmm4
cmp    %r14d,%eax
je     1780 <Pomeranz+0x5e0>
cmp    %r14d,%r15d
mov    %r14d,%edi
cmovl  %r15d,%edi
cmp    %r15,%r14
mov    %r15d,%ebx
cmovl  %r14d,%ebx
xorpd  %xmm5,%xmm5
cmp    %ebp,%ebx
jl     1700 <Pomeranz+0x560>
movslq %edi,%r12
mov    %r14,%r13
sub    %r12,%r13
add    $0x1,%r12
mov    (%rsp),%rdi
mov    (%rdi,%r9,8),%rdi
shl    $0x3,%r13
mov    0x10(%rsp),%rbx
add    (%rbx,%rsi,8),%r13
xorpd  %xmm5,%xmm5
nopl   0x0(%rax)
movsd  -0x8(%rdi,%r12,8),%xmm6
mulsd  0x0(%r13),%xmm6
addsd  %xmm6,%xmm5
add    $0xffffffffffffffff,%r12
add    $0x8,%r13
cmp    %r8,%r12
jg     1760 <Pomeranz+0x5c0>
jmp    1700 <Pomeranz+0x560>
ucomisd %xmm5,%xmm8
mov    0xc(%rsp),%r9d
jbe    18d0 <Pomeranz+0x730>
mov    0x4c(%rsp),%edi
mov    0x68(%rsp),%r13
cmp    %r13d,%edi
mov    0x60(%rsp),%r14
mov    0x58(%rsp),%r15
mov    0x28(%rsp),%r11
mov    0x20(%rsp),%r8
mov    0x50(%rsp),%r12
mov    0x40(%rsp),%rbx
jg     15a0 <Pomeranz+0x400>
mov    (%rsp),%rdx
mov    0x98(%rsp),%rsi
mov    (%rdx,%rsi,8),%rsi
movslq %edi,%rbp
sub    %r10d,%r13d
add    $0xfffffffe,%r13d
cmp    $0x3,%r13d
jb     18ad <Pomeranz+0x70d>
add    $0x1,%r13
mov    %r13,%r8
and    $0xfffffffffffffffc,%r8
lea    -0x4(%r8),%rdx
mov    %rdx,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rdx,%rdx
je     18f3 <Pomeranz+0x753>
mov    %r9,%rbx
and    $0xfffffffffffffffe,%rbx
mov    0x90(%rsp),%rdx
lea    (%rsi,%rdx,8),%rdi
add    $0x30,%rdi
xor    %edx,%edx
movupd -0x30(%rdi,%rdx,8),%xmm4
movupd -0x20(%rdi,%rdx,8),%xmm5
movupd -0x10(%rdi,%rdx,8),%xmm6
movupd (%rdi,%rdx,8),%xmm7
mulpd  %xmm0,%xmm4
mulpd  %xmm0,%xmm5
movupd %xmm4,-0x30(%rdi,%rdx,8)
movupd %xmm5,-0x20(%rdi,%rdx,8)
mulpd  %xmm0,%xmm6
mulpd  %xmm0,%xmm7
movupd %xmm6,-0x10(%rdi,%rdx,8)
movupd %xmm7,(%rdi,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rbx
jne    1820 <Pomeranz+0x680>
test   $0x1,%r9b
je     188f <Pomeranz+0x6ef>
add    %rbp,%rdx
movupd (%rsi,%rdx,8),%xmm4
movupd 0x10(%rsi,%rdx,8),%xmm5
mulpd  %xmm0,%xmm4
mulpd  %xmm0,%xmm5
movupd %xmm4,(%rsi,%rdx,8)
movupd %xmm5,0x10(%rsi,%rdx,8)
cmp    %r8,%r13
mov    0xc(%rsp),%r9d
mov    %r8,%rdx
mov    0x20(%rsp),%r8
mov    0x40(%rsp),%rbx
je     15a0 <Pomeranz+0x400>
add    %rdx,%rbp
mov    %eax,%eax
nop
movsd  (%rsi,%rbp,8),%xmm4
mulsd  %xmm1,%xmm4
movsd  %xmm4,(%rsi,%rbp,8)
add    $0x1,%rbp
cmp    %ebp,%eax
jne    18b0 <Pomeranz+0x710>
jmp    15a0 <Pomeranz+0x400>
nopl   0x0(%rax,%rax,1)
mov    0x60(%rsp),%r14
mov    0x58(%rsp),%r15
mov    0x28(%rsp),%r11
mov    0x20(%rsp),%r8
mov    0x50(%rsp),%r12
mov    0x40(%rsp),%rbx
jmp    15a5 <Pomeranz+0x405>
xor    %edx,%edx
test   $0x1,%r9b
jne    186e <Pomeranz+0x6ce>
jmp    188f <Pomeranz+0x6ef>
mov    $0x1,%eax
test   $0x1,%r12b
jne    138c <Pomeranz+0x1ec>
jmp    13a0 <Pomeranz+0x200>
mov    0x8(%r8),%rax
mov    %rbx,(%rax)
mov    0x10(%r8),%rax
mov    %rbx,(%rax)
mov    0x18(%r8),%rax
mov    %rbx,(%rax)
mov    $0x1,%ebx
mov    %r8,%rbp
jmp    197f <Pomeranz+0x7df>
mov    %edi,%eax
mov    (%rsp),%rcx
mov    (%rcx,%rax,8),%r10
mov    0x10(%rsp),%rbp
mov    0x1c(%rsp),%ebx
jmp    197f <Pomeranz+0x7df>
mov    $0x1,%edx
test   $0x1,%al
jne    143c <Pomeranz+0x29c>
jmp    1451 <Pomeranz+0x2b1>
mov    $0x1,%ecx
movsd  0x698(%rip),%xmm1
test   $0x1,%al
jne    14ca <Pomeranz+0x32a>
jmp    14de <Pomeranz+0x33e>
mov    $0x1,%ebx
mov    0x10(%rsp),%rbp
movslq %r9d,%rax
movsd  (%r10,%rax,8),%xmm0
movsd  %xmm0,0x10(%rsp)
mov    %r14,%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    0x28(%rsp),%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
call   1070 <DeleteMatrixD@plt>
mov    (%rsp),%rdi
call   1070 <DeleteMatrixD@plt>
mov    0x38(%rsp),%rdi
call   1050 <getLogFactorial@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
imul   $0x15e,%ebx,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
mulsd  0x639(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0x10(%rsp),%xmm0
call   10c0 <log@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
addsd  (%rsp),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jae    1a1f <Pomeranz+0x87f>
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    10b0 <exp@plt>
movsd  0x5d9(%rip),%xmm0
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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