#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*/ free(double *);
scalar_t__ log(double);
scalar_t__ malloc(int);
int /*<<< orphan*/ printf(char *, int, double, int, double);
void nllkunstr(int np, double *th, double *fnv) {
  extern int *y, nn, nc, ncl, d;
  extern double **x;
  double rh[20][20], L[20][20], *pa, *pb;
  double pr, nlk, eps, sd;
  double b0, *b, *rr, tem;
  int i, j, k, jc;
  int nsim, ifail;
  double pnorms(double);
  int pdef(double[][20], double[][20], int, double);
  void mvn1(int, double[], double[], double[][20], int, double, double *,
            double *, int *);
  pa = (double *)malloc((d + 1) * sizeof(double));
  pb = (double *)malloc((d + 1) * sizeof(double));
  b = (double *)malloc((nc + 1) * sizeof(double));
  rr = (double *)malloc((d + 1) * sizeof(double));
  nlk = 0.;
  eps = 1.e-6;
  b0 = th[1];
  for (j = 1; j <= nc; j++)
    b[j] = th[j + 1];
  for (j = nc + 2; j <= np; j++) {
    if (th[j] >= 1. || th[j] <= -1.) {
      *fnv = 1.e9;
      free(b);
      free(pa);
      free(pb);
      free(rr);
      return;
    }
  }
  for (i = 1, k = 1; i <= d; i++) {
    for (j = i; j <= d; j++) {
      if (i == j)
        rh[i][j] = 1;
      else {
        rh[i][j] = th[nc + 1 + k];
        k++;
        if (rh[i][j] > 1 || rh[i][j] < -1) {
          *fnv = 1.e9;
          free(b);
          free(pa);
          free(pb);
          free(rr);
          return;
        }
      }
    }
    for (j = 1; j <= i - 1; j++)
      rh[i][j] = rh[j][i];
  }
  if (!pdef(rh, L, d, 1.e-10)) {
    *fnv = 1.e9;
    free(b);
    free(pa);
    free(pb);
    free(rr);
    return;
  }
  for (i = 1; i <= ncl; i++) {
    if (d <= 7)
      nsim = 0;
    else
      nsim = 2000;
    for (j = 1; j <= d; j++) {
      for (jc = 1, tem = b0; jc <= nc; jc++)
        tem += b[jc] * x[d * (i - 1) + j][jc];
      if (y[d * (i - 1) + j] == 1) {
        pa[j] = -6.;
        pb[j] = tem;
      } else {
        pa[j] = tem;
        pb[j] = 6.;
      }
    }
    if (d >= 2) {
      mvn1(d, pa, pb, rh, nsim, eps, &pr, &sd, &ifail);
      if (ifail > 0) {
        *fnv = 1.e9;
        free(b);
        free(pa);
        free(pb);
        free(rr);
        return;
      }
    } else
      pr = pnorms(pb[1]) - pnorms(pa[1]);
    if (pr <= -.1) {
      *fnv = 1.e9;
      for (j = 1; j <= d; j++)
        printf("pa[%d]=%6.3f, pb[%d]=%6.3f\n", j, pa[j], j, pb[j]);
      return;
    }
    if (pr <= 0.)
      pr = 1.e-15;
    nlk -= log(pr);
  }
  *fnv = nlk;
  free(b);
  free(pa);
  free(pb);
  free(rr);
  return;
}


x86 assembly (clang, -O0):
<nllkunstr>:
push   %rbp
mov    %rsp,%rbp
sub    $0x19a0,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    0x2e33(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0x1928(%rbp)
mov    0x2e13(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0x1930(%rbp)
mov    0x2e0b(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0x1960(%rbp)
mov    0x2dd3(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0x1968(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1940(%rbp)
ff
movsd  0xde7(%rip),%xmm0
movsd  %xmm0,-0x1948(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x1958(%rbp)
ff
movl   $0x1,-0x1978(%rbp)
mov    -0x1978(%rbp),%eax
mov    0x2d8f(%rip),%rcx
cmp    (%rcx),%eax
jg     128d <nllkunstr+0x11d>
mov    -0x10(%rbp),%rax
mov    -0x1978(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x1960(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1978(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1978(%rbp)
jmp    123c <nllkunstr+0xcc>
mov    0x2d44(%rip),%rax
mov    (%rax),%eax
add    $0x2,%eax
mov    %eax,-0x1978(%rbp)
mov    -0x1978(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     134c <nllkunstr+0x1dc>
mov    -0x10(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xd72(%rip),%xmm1
ucomisd %xmm1,%xmm0
jae    12ee <nllkunstr+0x17e>
mov    -0x10(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  0xd5d(%rip),%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jb     1333 <nllkunstr+0x1c3>
mov    -0x18(%rbp),%rax
movsd  0xd16(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x1960(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1928(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1930(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1968(%rbp),%rdi
call   1030 <free@plt>
jmp    19f5 <nllkunstr+0x885>
jmp    1338 <nllkunstr+0x1c8>
mov    -0x1978(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1978(%rbp)
jmp    129f <nllkunstr+0x12f>
movl   $0x1,-0x1974(%rbp)
movl   $0x1,-0x197c(%rbp)
mov    -0x1974(%rbp),%eax
mov    0x2c53(%rip),%rcx
cmp    (%rcx),%eax
jg     158b <nllkunstr+0x41b>
mov    -0x1974(%rbp),%eax
mov    %eax,-0x1978(%rbp)
mov    -0x1978(%rbp),%eax
mov    0x2c32(%rip),%rcx
cmp    (%rcx),%eax
jg     14f5 <nllkunstr+0x385>
mov    -0x1974(%rbp),%eax
cmp    -0x1978(%rbp),%eax
jne    13d9 <nllkunstr+0x269>
movslq -0x1974(%rbp),%rcx
lea    -0xca0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x1978(%rbp),%rcx
movsd  0xc69(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    14dc <nllkunstr+0x36c>
mov    -0x10(%rbp),%rax
mov    0x2bf4(%rip),%rcx
mov    (%rcx),%ecx
add    $0x1,%ecx
add    -0x197c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0x1974(%rbp),%rcx
lea    -0xca0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x1978(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x197c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x197c(%rbp)
movslq -0x1974(%rbp),%rcx
lea    -0xca0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x1978(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xbe2(%rip),%xmm1
ucomisd %xmm1,%xmm0
ja     1492 <nllkunstr+0x322>
movslq -0x1974(%rbp),%rcx
lea    -0xca0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x1978(%rbp),%rcx
movsd  0xbb9(%rip),%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    14d7 <nllkunstr+0x367>
mov    -0x18(%rbp),%rax
movsd  0xb72(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x1960(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1928(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1930(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1968(%rbp),%rdi
call   1030 <free@plt>
jmp    19f5 <nllkunstr+0x885>
jmp    14dc <nllkunstr+0x36c>
jmp    14e1 <nllkunstr+0x371>
mov    -0x1978(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1978(%rbp)
jmp    1381 <nllkunstr+0x211>
movl   $0x1,-0x1978(%rbp)
mov    -0x1978(%rbp),%eax
mov    -0x1974(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1572 <nllkunstr+0x402>
movslq -0x1978(%rbp),%rcx
lea    -0xca0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x1974(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0x1974(%rbp),%rcx
lea    -0xca0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0x1978(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1978(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1978(%rbp)
jmp    14ff <nllkunstr+0x38f>
jmp    1577 <nllkunstr+0x407>
mov    -0x1974(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1974(%rbp)
jmp    1360 <nllkunstr+0x1f0>
lea    -0xca0(%rbp),%rdi
lea    -0x1920(%rbp),%rsi
mov    0x2a20(%rip),%rax
mov    (%rax),%edx
movsd  0xa5e(%rip),%xmm0
call   1060 <pdef@plt>
cmp    $0x0,%eax
jne    15fd <nllkunstr+0x48d>
mov    -0x18(%rbp),%rax
movsd  0xa4c(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x1960(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1928(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1930(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1968(%rbp),%rdi
call   1030 <free@plt>
jmp    19f5 <nllkunstr+0x885>
movl   $0x1,-0x1974(%rbp)
mov    -0x1974(%rbp),%eax
mov    0x29d4(%rip),%rcx
cmp    (%rcx),%eax
jg     19b5 <nllkunstr+0x845>
mov    0x299d(%rip),%rax
cmpl   $0x7,(%rax)
jg     163b <nllkunstr+0x4cb>
movl   $0x0,-0x1984(%rbp)
jmp    1645 <nllkunstr+0x4d5>
movl   $0x7d0,-0x1984(%rbp)
movl   $0x1,-0x1978(%rbp)
mov    -0x1978(%rbp),%eax
mov    0x2964(%rip),%rcx
cmp    (%rcx),%eax
jg     17c6 <nllkunstr+0x656>
movl   $0x1,-0x1980(%rbp)
movsd  -0x1958(%rbp),%xmm0
ff
movsd  %xmm0,-0x1970(%rbp)
ff
mov    -0x1980(%rbp),%eax
mov    0x294d(%rip),%rcx
cmp    (%rcx),%eax
jg     170a <nllkunstr+0x59a>
mov    -0x1960(%rbp),%rax
movslq -0x1980(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x291b(%rip),%rax
mov    (%rax),%rax
mov    0x2909(%rip),%rcx
mov    (%rcx),%ecx
mov    -0x1974(%rbp),%edx
sub    $0x1,%edx
imul   %edx,%ecx
add    -0x1978(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1980(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x1970(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1970(%rbp)
ff
mov    -0x1980(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1980(%rbp)
jmp    167e <nllkunstr+0x50e>
mov    0x28cf(%rip),%rax
mov    (%rax),%rax
mov    0x28a5(%rip),%rcx
mov    (%rcx),%ecx
mov    -0x1974(%rbp),%edx
sub    $0x1,%edx
imul   %edx,%ecx
add    -0x1978(%rbp),%ecx
movslq %ecx,%rcx
cmpl   $0x1,(%rax,%rcx,4)
jne    1777 <nllkunstr+0x607>
mov    -0x1928(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  0x8de(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x1970(%rbp),%xmm0
ff
mov    -0x1930(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    17ad <nllkunstr+0x63d>
movsd  -0x1970(%rbp),%xmm0
ff
mov    -0x1928(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1930(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  0x880(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    17b2 <nllkunstr+0x642>
mov    -0x1978(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1978(%rbp)
jmp    164f <nllkunstr+0x4df>
mov    0x27f3(%rip),%rax
cmpl   $0x2,(%rax)
jl     187d <nllkunstr+0x70d>
mov    0x27e3(%rip),%rax
mov    (%rax),%edi
mov    -0x1928(%rbp),%rsi
mov    -0x1930(%rbp),%rdx
lea    -0xca0(%rbp),%rcx
mov    -0x1984(%rbp),%r8d
movsd  -0x1948(%rbp),%xmm0
ff
lea    -0x1938(%rbp),%r9
lea    -0x1950(%rbp),%r10
lea    -0x1988(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
call   1050 <mvn1@plt>
cmpl   $0x0,-0x1988(%rbp)
jle    1878 <nllkunstr+0x708>
mov    -0x18(%rbp),%rax
movsd  0x7d1(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x1960(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1928(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1930(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1968(%rbp),%rdi
call   1030 <free@plt>
jmp    19f5 <nllkunstr+0x885>
jmp    18be <nllkunstr+0x74e>
mov    -0x1930(%rbp),%rax
movsd  0x8(%rax),%xmm0
call   1040 <pnorms@plt>
movsd  %xmm0,-0x1990(%rbp)
ff
mov    -0x1928(%rbp),%rax
movsd  0x8(%rax),%xmm0
call   1040 <pnorms@plt>
movaps %xmm0,%xmm1
movsd  -0x1990(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x1938(%rbp)
ff
movsd  0x752(%rip),%xmm0
ucomisd -0x1938(%rbp),%xmm0
ff
jb     195a <nllkunstr+0x7ea>
mov    -0x18(%rbp),%rax
movsd  0x730(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0x1978(%rbp)
mov    -0x1978(%rbp),%eax
mov    0x26c5(%rip),%rcx
cmp    (%rcx),%eax
jg     1955 <nllkunstr+0x7e5>
mov    -0x1978(%rbp),%esi
mov    -0x1928(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x1978(%rbp),%edx
mov    -0x1930(%rbp),%rax
movslq -0x1978(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
lea    0x70c(%rip),%rdi
call   1070 <printf@plt>
mov    -0x1978(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1978(%rbp)
jmp    18ee <nllkunstr+0x77e>
jmp    19f5 <nllkunstr+0x885>
xorps  %xmm0,%xmm0
ucomisd -0x1938(%rbp),%xmm0
ff
jb     197b <nllkunstr+0x80b>
movsd  0x6ad(%rip),%xmm0
movsd  %xmm0,-0x1938(%rbp)
ff
movsd  -0x1938(%rbp),%xmm0
ff
call   1090 <log@plt>
cvtsi2sd %rax,%xmm1
movsd  -0x1940(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x1940(%rbp)
ff
mov    -0x1974(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1974(%rbp)
jmp    1607 <nllkunstr+0x497>
movsd  -0x1940(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x1960(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1928(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1930(%rbp),%rdi
call   1030 <free@plt>
mov    -0x1968(%rbp),%rdi
call   1030 <free@plt>
add    $0x19a0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<nllkunstr>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1978,%rsp
mov    %rdx,%r12
mov    %rsi,%rbp
mov    %edi,%r15d
mov    0x2e2f(%rip),%r14
mov    (%r14),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,0x10(%rsp)
mov    (%r14),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,0x18(%rsp)
mov    0x2e18(%rip),%r13
mov    0x0(%r13),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,%rbx
mov    (%r14),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,0x28(%rsp)
movsd  0x8(%rbp),%xmm3
mov    0x0(%r13),%r13d
test   %r13d,%r13d
jle    133b <nllkunstr+0x1cb>
lea    0x1(%r13),%rcx
mov    $0x1,%edx
cmp    $0x4,%r13
jb     12c9 <nllkunstr+0x159>
lea    0x8(%rbx),%rax
lea    0x8(,%rcx,8),%rsi
add    %rbp,%rsi
cmp    %rsi,%rax
jae    1231 <nllkunstr+0xc1>
lea    (%rbx,%rcx,8),%rax
lea    0x10(%rbp),%rsi
cmp    %rax,%rsi
jb     12c9 <nllkunstr+0x159>
mov    %r13d,%edx
and    $0xfffffffc,%edx
lea    -0x4(%rdx),%rsi
mov    %rsi,%rax
shr    $0x2,%rax
add    $0x1,%rax
test   %rsi,%rsi
je     1984 <nllkunstr+0x814>
mov    %rax,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movups 0x10(%rbp,%rsi,8),%xmm0
movups 0x20(%rbp,%rsi,8),%xmm1
movups %xmm0,0x8(%rbx,%rsi,8)
movups %xmm1,0x18(%rbx,%rsi,8)
movupd 0x30(%rbp,%rsi,8),%xmm0
movupd 0x40(%rbp,%rsi,8),%xmm1
movupd %xmm0,0x28(%rbx,%rsi,8)
movupd %xmm1,0x38(%rbx,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1260 <nllkunstr+0xf0>
test   $0x1,%al
je     12c0 <nllkunstr+0x150>
shl    $0x3,%rsi
mov    %rsi,%rax
or     $0x10,%rax
movupd 0x0(%rbp,%rax,1),%xmm0
movupd 0x10(%rbp,%rax,1),%xmm1
or     $0x8,%rsi
movupd %xmm0,(%rbx,%rsi,1)
movupd %xmm1,0x10(%rbx,%rsi,1)
cmp    %rdx,%r13
je     133b <nllkunstr+0x1cb>
or     $0x1,%rdx
mov    %ecx,%edi
sub    %edx,%edi
mov    %rdx,%rsi
not    %rsi
add    %rcx,%rsi
and    $0x3,%rdi
je     18fc <nllkunstr+0x78c>
lea    0x1(%rdx),%rax
movsd  0x8(%rbp,%rdx,8),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
mov    %rax,%rdx
add    $0xffffffffffffffff,%rdi
jne    12e0 <nllkunstr+0x170>
cmp    $0x3,%rsi
jb     133b <nllkunstr+0x1cb>
xchg   %ax,%ax
movsd  0x8(%rbp,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rax,8)
movsd  0x10(%rbp,%rax,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rax,8)
movsd  0x18(%rbp,%rax,8),%xmm0
movsd  %xmm0,0x10(%rbx,%rax,8)
movsd  0x20(%rbp,%rax,8),%xmm0
movsd  %xmm0,0x18(%rbx,%rax,8)
lea    0x4(%rax),%rdx
mov    %rdx,%rax
cmp    %rdx,%rcx
jne    1300 <nllkunstr+0x190>
lea    0x2(%r13),%eax
cmp    %r15d,%eax
jg     1391 <nllkunstr+0x221>
cltq
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
mov    %r13d,%ecx
not    %ecx
add    %r15d,%ecx
xor    %edx,%edx
movsd  0xc9d(%rip),%xmm0
movsd  0xc9d(%rip),%xmm1
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm2
ucomisd %xmm0,%xmm2
jae    18b7 <nllkunstr+0x747>
ucomisd %xmm2,%xmm1
jae    18b7 <nllkunstr+0x747>
add    $0x1,%rdx
cmp    %edx,%ecx
jne    1370 <nllkunstr+0x200>
mov    %r12,0x20(%rsp)
mov    (%r14),%eax
mov    %rax,0x8(%rsp)
test   %eax,%eax
jle    15a4 <nllkunstr+0x434>
mov    0x8(%rsp),%rdx
lea    0x1(%rdx),%rax
mov    %rax,0x48(%rsp)
lea    0x118(%rsp),%r10
lea    0x110(%rsp),%rsi
lea    0x130(%rsp),%r11
lea    0x2f8(%rsp),%rax
mov    %rax,0x38(%rsp)
lea    0x78(%rsp),%rax
mov    %rax,0x40(%rsp)
mov    $0x1,%edi
mov    $0x1,%r15d
xor    %r9d,%r9d
movsd  0xc07(%rip),%xmm0
movsd  0xc07(%rip),%xmm1
jmp    144b <nllkunstr+0x2db>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r15
add    $0x1,%r9
add    $0xa8,%r10
add    $0xffffffffffffffff,%rdx
add    $0xa0,%r11
addq   $0x8,0x38(%rsp)
mov    %r14,%rsi
add    $0xa0,%rsi
addq   $0x8,0x40(%rsp)
cmp    0x48(%rsp),%r15
je     15a4 <nllkunstr+0x434>
mov    %r9d,%r12d
and    $0x3,%r12d
mov    %r9,%r8
and    $0xfffffffffffffffc,%r8
lea    -0x1(%r9),%rax
mov    %rax,0x30(%rsp)
lea    (%r15,%r15,4),%rax
shl    $0x5,%rax
add    %rsp,%rax
add    $0x70,%rax
lea    (%rax,%r15,8),%rax
xor    %r14d,%r14d
jmp    1496 <nllkunstr+0x326>
nopw   0x0(%rax,%rax,1)
movabs $0x3ff0000000000000,%rcx
mov    %rcx,(%rax)
add    $0x1,%r14
cmp    %r14,%rdx
je     14d0 <nllkunstr+0x360>
test   %r14,%r14
je     1480 <nllkunstr+0x310>
lea    (%rdi,%r13,1),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  0x0(%rbp,%rcx,8),%xmm2
movsd  %xmm2,(%r10,%r14,8)
ucomisd %xmm0,%xmm2
ja     158d <nllkunstr+0x41d>
ucomisd %xmm2,%xmm1
ja     158d <nllkunstr+0x41d>
add    $0x1,%edi
jmp    148d <nllkunstr+0x31d>
nopw   0x0(%rax,%rax,1)
mov    %rsi,%r14
cmp    $0x2,%r15
jb     1410 <nllkunstr+0x2a0>
mov    $0x1,%eax
cmpq   $0x3,0x30(%rsp)
jb     154b <nllkunstr+0x3db>
mov    0x38(%rsp),%rsi
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x1e0(%rsi),%xmm2
ff
movsd  %xmm2,-0x18(%r11,%rax,8)
movsd  -0x140(%rsi),%xmm2
ff
movsd  %xmm2,-0x10(%r11,%rax,8)
movsd  -0xa0(%rsi),%xmm2
ff
movsd  %xmm2,-0x8(%r11,%rax,8)
movsd  (%rsi),%xmm2
movsd  %xmm2,(%r11,%rax,8)
add    $0x4,%rax
add    $0x280,%rsi
cmp    %rax,%r8
jne    1500 <nllkunstr+0x390>
add    $0x1,%rax
test   $0x3,%r9b
je     1410 <nllkunstr+0x2a0>
lea    (%r14,%rax,8),%rsi
lea    (%rax,%rax,4),%rax
shl    $0x5,%rax
add    0x40(%rsp),%rax
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movsd  (%rax),%xmm2
movsd  %xmm2,(%rsi,%rcx,8)
add    $0x1,%rcx
add    $0xa0,%rax
cmp    %rcx,%r12
jne    1570 <nllkunstr+0x400>
jmp    1410 <nllkunstr+0x2a0>
movabs $0x41cdcd6500000000,%rax
cd cd 41
mov    0x20(%rsp),%rcx
mov    %rax,(%rcx)
jmp    18c5 <nllkunstr+0x755>
movsd  %xmm3,0x60(%rsp)
lea    0x70(%rsp),%rdi
lea    0xcf0(%rsp),%rsi
movsd  0xa51(%rip),%xmm0
mov    0x8(%rsp),%rdx
call   1060 <pdef@plt>
test   %eax,%eax
je     190e <nllkunstr+0x79e>
mov    0x2a10(%rip),%rax
cmpl   $0x0,(%rax)
jle    1939 <nllkunstr+0x7c9>
mov    $0x1,%eax
mov    %rax,0x30(%rsp)
movsd  0xa25(%rip),%xmm5
movsd  0xa25(%rip),%xmm6
xorpd  %xmm0,%xmm0
mov    0x18(%rsp),%r15
mov    0x10(%rsp),%rbp
mov    0x29b0(%rip),%r12
jmp    166b <nllkunstr+0x4fb>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm1,%xmm0
call   1090 <log@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  0x40(%rsp),%xmm1
subsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
mov    0x30(%rsp),%rdx
lea    0x1(%rdx),%eax
mov    0x299a(%rip),%rcx
cmp    (%rcx),%edx
mov    %rax,0x30(%rsp)
movsd  0x9bb(%rip),%xmm5
movsd  0x9bb(%rip),%xmm6
jge    1947 <nllkunstr+0x7d7>
movsd  %xmm0,0x40(%rsp)
mov    (%r12),%edx
cmp    $0x8,%edx
mov    $0x7d0,%ecx
mov    $0x0,%eax
cmovl  %eax,%ecx
test   %edx,%edx
jle    1850 <nllkunstr+0x6e0>
mov    %ecx,0x48(%rsp)
mov    0x2940(%rip),%rax
mov    (%rax),%r11d
mov    0x2926(%rip),%rax
mov    (%rax),%rax
mov    %rax,0x8(%rsp)
mov    0x30(%rsp),%rax
add    $0xffffffff,%eax
imul   %edx,%eax
mov    0x2924(%rip),%rcx
mov    (%rcx),%r9
lea    0x1(%r11),%r10d
movslq %eax,%r8
mov    %rdx,0x38(%rsp)
lea    0x1(%rdx),%rcx
lea    -0x1(%r10),%rax
add    $0xfffffffffffffffe,%r10
mov    %eax,%esi
and    $0x3,%esi
and    $0xfffffffffffffffc,%rax
mov    $0x1,%r12d
movsd  0x60(%rsp),%xmm4
jmp    170a <nllkunstr+0x59a>
xchg   %ax,%ax
movsd  %xmm1,0x0(%rbp,%r12,8)
movsd  %xmm0,(%r15,%r12,8)
add    $0x1,%r12
cmp    %rcx,%r12
je     17f0 <nllkunstr+0x680>
lea    (%r12,%r8,1),%r15
movapd %xmm4,%xmm0
test   %r11d,%r11d
jle    17c8 <nllkunstr+0x658>
mov    0x8(%rsp),%rdx
mov    (%rdx,%r15,8),%r14
mov    $0x1,%r13d
movapd %xmm4,%xmm0
cmp    $0x3,%r10
jb     1795 <nllkunstr+0x625>
movapd %xmm4,%xmm0
xor    %r13d,%r13d
nopl   0x0(%rax,%rax,1)
movsd  0x8(%rbx,%r13,8),%xmm1
movsd  0x10(%rbx,%r13,8),%xmm2
mulsd  0x8(%r14,%r13,8),%xmm1
mulsd  0x10(%r14,%r13,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x18(%rbx,%r13,8),%xmm3
mulsd  0x18(%r14,%r13,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x20(%rbx,%r13,8),%xmm0
mulsd  0x20(%r14,%r13,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%r13
cmp    %r13,%rax
jne    1740 <nllkunstr+0x5d0>
add    $0x1,%r13
test   %rsi,%rsi
je     17c8 <nllkunstr+0x658>
lea    (%r14,%r13,8),%rdx
lea    (%rbx,%r13,8),%rbp
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  0x0(%rbp,%rdi,8),%xmm1
mulsd  (%rdx,%rdi,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdi
cmp    %rdi,%rsi
jne    17b0 <nllkunstr+0x640>
cmpl   $0x1,(%r9,%r15,4)
movapd %xmm5,%xmm1
je     17d7 <nllkunstr+0x667>
movapd %xmm0,%xmm1
mov    0x18(%rsp),%r15
mov    0x10(%rsp),%rbp
je     16f0 <nllkunstr+0x580>
movapd %xmm6,%xmm0
jmp    16f0 <nllkunstr+0x580>
mov    0x38(%rsp),%rdi
cmp    $0x2,%edi
mov    0x27c1(%rip),%r12
mov    0x48(%rsp),%r8d
jl     1850 <nllkunstr+0x6e0>
mov    %rbp,%rsi
mov    %r15,%rdx
lea    0x70(%rsp),%rcx
movsd  0x80f(%rip),%xmm0
lea    0x50(%rsp),%r9
lea    0x5c(%rsp),%rax
push   %rax
lea    0x70(%rsp),%rax
push   %rax
call   1050 <mvn1@plt>
add    $0x10,%rsp
cmpl   $0x0,0x5c(%rsp)
jg     19ea <nllkunstr+0x87a>
movsd  0x50(%rsp),%xmm1
jmp    187b <nllkunstr+0x70b>
cs nopw 0x0(%rax,%rax,1)
movsd  0x8(%r15),%xmm0
call   1040 <pnorms@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x8(%rbp),%xmm0
call   1040 <pnorms@plt>
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x50(%rsp)
movsd  0x7ad(%rip),%xmm0
ucomisd %xmm1,%xmm0
jae    1993 <nllkunstr+0x823>
xorpd  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jb     1620 <nllkunstr+0x4b0>
movabs $0x3cd203af9ee75616,%rax
03 d2 3c
mov    %rax,0x50(%rsp)
movsd  0x786(%rip),%xmm1
jmp    1620 <nllkunstr+0x4b0>
movabs $0x41cdcd6500000000,%rax
cd cd 41
mov    %rax,(%r12)
mov    %rbx,%rdi
call   1030 <free@plt>
mov    0x10(%rsp),%rdi
call   1030 <free@plt>
mov    0x18(%rsp),%rdi
call   1030 <free@plt>
mov    0x28(%rsp),%rdi
add    $0x1978,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>
mov    %rdx,%rax
cmp    $0x3,%rsi
jae    1300 <nllkunstr+0x190>
jmp    133b <nllkunstr+0x1cb>
movabs $0x41cdcd6500000000,%rax
cd cd 41
mov    0x20(%rsp),%rcx
mov    %rax,(%rcx)
mov    %rbx,%rdi
call   1030 <free@plt>
mov    0x10(%rsp),%rdi
call   1030 <free@plt>
mov    0x18(%rsp),%rdi
jmp    1963 <nllkunstr+0x7f3>
xorpd  %xmm0,%xmm0
mov    0x18(%rsp),%r15
mov    0x10(%rsp),%rbp
mov    0x20(%rsp),%rax
movsd  %xmm0,(%rax)
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    0x28(%rsp),%rdi
call   1030 <free@plt>
add    $0x1978,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %esi,%esi
test   $0x1,%al
jne    129a <nllkunstr+0x12a>
jmp    12c0 <nllkunstr+0x150>
movabs $0x41cdcd6500000000,%rax
cd cd 41
mov    0x20(%rsp),%rcx
mov    %rax,(%rcx)
cmpl   $0x0,(%r12)
jle    1972 <nllkunstr+0x802>
mov    $0x1,%ebx
lea    0x690(%rip),%r14
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rbx,8),%xmm0
movsd  (%r15,%rbx,8),%xmm1
mov    %r14,%rdi
mov    %ebx,%esi
mov    %ebx,%edx
call   1070 <printf@plt>
lea    0x1(%rbx),%rax
movslq (%r12),%rcx
cmp    %rcx,%rbx
mov    %rax,%rbx
jl     19c0 <nllkunstr+0x850>
jmp    1972 <nllkunstr+0x802>
movabs $0x41cdcd6500000000,%rax
cd cd 41
mov    0x20(%rsp),%rcx
mov    %rax,(%rcx)
jmp    1950 <nllkunstr+0x7e0>

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

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