#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 nllkexchord(int np, double *th, double *fnv) {
  extern double **x;
  extern int *y, nn, nc, dmax, ncl, norc;
  extern int *dvec, *dstart;
  double *pa, *pb, rh[20][20];
  double pr, nlk, eps, *b0, *b, tem, rr;
  int i, j, k, jc, nsim, d, startp;
  double pnorms(double);
  double exchmvn(int, double[], double[], double, double);
  b0 = (double *)malloc((norc + 1) * sizeof(double));
  b = (double *)malloc((nc + 1) * sizeof(double));
  pa = (double *)malloc((dmax + 1) * sizeof(double));
  pb = (double *)malloc((dmax + 1) * sizeof(double));
  nlk = 0.;
  b0[0] = 0.0;
  eps = 1.e-6;
  for (i = 1; i <= norc - 1; i++)
    b0[i] = th[i];
  for (i = 1; i <= nc; i++)
    b[i] = th[norc - 1 + i];
  rr = th[np];
  if (rr < 0.0 || rr >= 1.) {
    *fnv = 1.e10;
    return;
  }
  for (j = 1; j <= dmax; j++)
    for (k = 1; k <= dmax; k++)
      if (j == k)
        rh[j][k] = 1.;
      else
        rh[j][k] = rr;
  for (i = 1; i <= ncl; i++) {
    d = dvec[i];
    startp = dstart[i];
    if (d <= 7)
      nsim = 0;
    else
      nsim = 2000;
    for (j = 1; j <= d; j++) {
      b0[0] = 0.0;
      for (jc = 1, tem = 0.0; jc <= nc; jc++) {
        tem += b[jc] * x[startp + j - 1][jc];
      }
      k = y[startp + j - 1];
      if (k == 1) {
        pa[j] = -6.;
        pb[j] = tem + b0[k];
      } else if (k == norc) {
        pa[j] = tem + b0[k - 1];
        pb[j] = 6.;
      } else {
        pa[j] = tem + b0[k - 1];
        pb[j] = tem + b0[k];
      }
    }
    if (d >= 2)
      pr = exchmvn(d, pa, pb, rr, eps);
    else if (k >= 1)
      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(b0);
  free(b);
  free(pa);
  free(pb);
  return;
}


x86 assembly (clang, -O0):
<nllkexchord>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd10,%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   1070 <malloc@plt>
mov    %rax,-0xcd0(%rbp)
mov    0x2e23(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0xcd8(%rbp)
mov    0x2e23(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x20(%rbp)
mov    0x2e06(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xcc0(%rbp)
ff
mov    -0xcd0(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movsd  0xdef(%rip),%xmm0
movsd  %xmm0,-0xcc8(%rbp)
ff
movl   $0x1,-0xcec(%rbp)
mov    -0xcec(%rbp),%eax
mov    0x2d80(%rip),%rcx
mov    (%rcx),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1274 <nllkexchord+0x114>
mov    -0x10(%rbp),%rax
movslq -0xcec(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xcd0(%rbp),%rax
movslq -0xcec(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xcec(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcec(%rbp)
jmp    1223 <nllkexchord+0xc3>
movl   $0x1,-0xcec(%rbp)
mov    -0xcec(%rbp),%eax
mov    0x2d35(%rip),%rcx
cmp    (%rcx),%eax
jg     12d8 <nllkexchord+0x178>
mov    -0x10(%rbp),%rax
mov    0x2d12(%rip),%rcx
mov    (%rcx),%ecx
sub    $0x1,%ecx
add    -0xcec(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xcd8(%rbp),%rax
movslq -0xcec(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xcec(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcec(%rbp)
jmp    127e <nllkexchord+0x11e>
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xce8(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0xce8(%rbp),%xmm0
ff
ja     1318 <nllkexchord+0x1b8>
movsd  -0xce8(%rbp),%xmm0
ff
movsd  0xcfa(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     132d <nllkexchord+0x1cd>
mov    -0x18(%rbp),%rax
movsd  0xd14(%rip),%xmm0
movsd  %xmm0,(%rax)
jmp    184b <nllkexchord+0x6eb>
movl   $0x1,-0xcf0(%rbp)
mov    -0xcf0(%rbp),%eax
mov    0x2c9c(%rip),%rcx
cmp    (%rcx),%eax
jg     140c <nllkexchord+0x2ac>
movl   $0x1,-0xcf4(%rbp)
mov    -0xcf4(%rbp),%eax
mov    0x2c7d(%rip),%rcx
cmp    (%rcx),%eax
jg     13f3 <nllkexchord+0x293>
mov    -0xcf0(%rbp),%eax
cmp    -0xcf4(%rbp),%eax
jne    13ae <nllkexchord+0x24e>
movslq -0xcf0(%rbp),%rcx
lea    -0xcb0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0xcf4(%rbp),%rcx
movsd  0xc64(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    13da <nllkexchord+0x27a>
movsd  -0xce8(%rbp),%xmm0
ff
movslq -0xcf0(%rbp),%rcx
lea    -0xcb0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0xcf4(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    13df <nllkexchord+0x27f>
mov    -0xcf4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcf4(%rbp)
jmp    1356 <nllkexchord+0x1f6>
jmp    13f8 <nllkexchord+0x298>
mov    -0xcf0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcf0(%rbp)
jmp    1337 <nllkexchord+0x1d7>
movl   $0x1,-0xcec(%rbp)
mov    -0xcec(%rbp),%eax
mov    0x2bc5(%rip),%rcx
cmp    (%rcx),%eax
jg     1811 <nllkexchord+0x6b1>
mov    0x2b9e(%rip),%rax
mov    (%rax),%rax
movslq -0xcec(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xd00(%rbp)
mov    0x2b7c(%rip),%rax
mov    (%rax),%rax
movslq -0xcec(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xd04(%rbp)
cmpl   $0x7,-0xd00(%rbp)
jg     147b <nllkexchord+0x31b>
movl   $0x0,-0xcfc(%rbp)
jmp    1485 <nllkexchord+0x325>
movl   $0x7d0,-0xcfc(%rbp)
movl   $0x1,-0xcf0(%rbp)
mov    -0xcf0(%rbp),%eax
cmp    -0xd00(%rbp),%eax
jg     1699 <nllkexchord+0x539>
mov    -0xcd0(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0xcf8(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xce0(%rbp)
ff
mov    -0xcf8(%rbp),%eax
mov    0x2aef(%rip),%rcx
cmp    (%rcx),%eax
jg     1544 <nllkexchord+0x3e4>
mov    -0xcd8(%rbp),%rax
movslq -0xcf8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2ab5(%rip),%rax
mov    (%rax),%rax
mov    -0xd04(%rbp),%ecx
add    -0xcf0(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xcf8(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xce0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xce0(%rbp)
ff
mov    -0xcf8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcf8(%rbp)
jmp    14c4 <nllkexchord+0x364>
mov    0x2a8d(%rip),%rax
mov    (%rax),%rax
mov    -0xd04(%rbp),%ecx
add    -0xcf0(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xcf4(%rbp)
cmpl   $0x1,-0xcf4(%rbp)
jne    15be <nllkexchord+0x45e>
mov    -0x20(%rbp),%rax
movslq -0xcf0(%rbp),%rcx
movsd  0xaa7(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xce0(%rbp),%xmm0
ff
mov    -0xcd0(%rbp),%rax
movslq -0xcf4(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0xcf0(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1680 <nllkexchord+0x520>
mov    -0xcf4(%rbp),%eax
mov    0x29e5(%rip),%rcx
cmp    (%rcx),%eax
jne    1620 <nllkexchord+0x4c0>
movsd  -0xce0(%rbp),%xmm0
ff
mov    -0xcd0(%rbp),%rax
mov    -0xcf4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0xcf0(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
movslq -0xcf0(%rbp),%rcx
movsd  0xa12(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    167b <nllkexchord+0x51b>
movsd  -0xce0(%rbp),%xmm0
ff
mov    -0xcd0(%rbp),%rax
mov    -0xcf4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0xcf0(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xce0(%rbp),%xmm0
ff
mov    -0xcd0(%rbp),%rax
movslq -0xcf4(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0xcf0(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1680 <nllkexchord+0x520>
jmp    1685 <nllkexchord+0x525>
mov    -0xcf0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcf0(%rbp)
jmp    148f <nllkexchord+0x32f>
cmpl   $0x2,-0xd00(%rbp)
jl     16d6 <nllkexchord+0x576>
mov    -0xd00(%rbp),%edi
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
movsd  -0xce8(%rbp),%xmm0
ff
movsd  -0xcc8(%rbp),%xmm1
ff
call   1060 <exchmvn@plt>
movsd  %xmm0,-0xcb8(%rbp)
ff
jmp    1723 <nllkexchord+0x5c3>
cmpl   $0x1,-0xcf4(%rbp)
jl     171e <nllkexchord+0x5be>
mov    -0x28(%rbp),%rax
movsd  0x8(%rax),%xmm0
call   1040 <pnorms@plt>
movsd  %xmm0,-0xd10(%rbp)
ff
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
call   1040 <pnorms@plt>
movaps %xmm0,%xmm1
movsd  -0xd10(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xcb8(%rbp)
ff
jmp    1723 <nllkexchord+0x5c3>
movsd  0x8e5(%rip),%xmm0
ucomisd -0xcb8(%rbp),%xmm0
ff
jb     17b6 <nllkexchord+0x656>
mov    -0x18(%rbp),%rax
movsd  0x8db(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0xcf0(%rbp)
mov    -0xcf0(%rbp),%eax
cmp    -0xd00(%rbp),%eax
jg     17b1 <nllkexchord+0x651>
mov    -0xcf0(%rbp),%esi
mov    -0x20(%rbp),%rax
movslq -0xcf0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xcf0(%rbp),%edx
mov    -0x28(%rbp),%rax
movslq -0xcf0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
lea    0x8a8(%rip),%rdi
call   1050 <printf@plt>
mov    -0xcf0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcf0(%rbp)
jmp    1753 <nllkexchord+0x5f3>
jmp    184b <nllkexchord+0x6eb>
xorps  %xmm0,%xmm0
ucomisd -0xcb8(%rbp),%xmm0
ff
jb     17d7 <nllkexchord+0x677>
movsd  0x849(%rip),%xmm0
movsd  %xmm0,-0xcb8(%rbp)
ff
movsd  -0xcb8(%rbp),%xmm0
ff
call   1080 <log@plt>
cvtsi2sd %rax,%xmm1
movsd  -0xcc0(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xcc0(%rbp)
ff
mov    -0xcec(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcec(%rbp)
jmp    1416 <nllkexchord+0x2b6>
movsd  -0xcc0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0xcd0(%rbp),%rdi
call   1030 <free@plt>
mov    -0xcd8(%rbp),%rdi
call   1030 <free@plt>
mov    -0x20(%rbp),%rdi
call   1030 <free@plt>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
add    $0xd10,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<nllkexchord>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %rdx,0x18(%rsp)
mov    %rsi,%rbp
mov    %edi,0x10(%rsp)
mov    0x2e2f(%rip),%r13
mov    0x0(%r13),%eax
lea    0x8(,%rax,8),%edi
call   1070 <malloc@plt>
mov    %rax,%r12
mov    0x2e25(%rip),%rax
mov    (%rax),%eax
lea    0x8(,%rax,8),%edi
call   1070 <malloc@plt>
mov    %rax,%rbx
mov    0x2e2d(%rip),%r15
mov    (%r15),%eax
lea    0x8(,%rax,8),%edi
call   1070 <malloc@plt>
mov    %rax,%r14
mov    (%r15),%eax
lea    0x8(,%rax,8),%edi
call   1070 <malloc@plt>
mov    %r12,%rdx
mov    %rax,%r12
movq   $0x0,(%rdx)
mov    0x0(%r13),%r9d
cmp    $0x2,%r9d
mov    %rdx,0x8(%rsp)
jl     1398 <nllkexchord+0x238>
lea    -0x1(%r9),%r8
mov    $0x1,%ecx
cmp    $0x4,%r8
jb     1318 <nllkexchord+0x1b8>
add    $0x8,%rdx
lea    0x0(,%r9,8),%rsi
add    %rbp,%rsi
cmp    %rsi,%rdx
jae    1231 <nllkexchord+0xd1>
mov    0x8(%rsp),%rax
lea    (%rax,%r9,8),%rdx
lea    0x8(%rbp),%rsi
cmp    %rdx,%rsi
jb     1318 <nllkexchord+0x1b8>
mov    %r8,%rcx
and    $0xfffffffffffffffc,%rcx
lea    -0x4(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x3,%edx
cmp    $0xc,%rdi
jae    125b <nllkexchord+0xfb>
xor    %edi,%edi
mov    0x8(%rsp),%rax
jmp    12cc <nllkexchord+0x16c>
and    $0xfffffffffffffffc,%rsi
xor    %edi,%edi
mov    0x8(%rsp),%rax
cs nopw 0x0(%rax,%rax,1)
movups 0x8(%rbp,%rdi,8),%xmm0
movups 0x18(%rbp,%rdi,8),%xmm1
movups %xmm0,0x8(%rax,%rdi,8)
movups %xmm1,0x18(%rax,%rdi,8)
movups 0x28(%rbp,%rdi,8),%xmm0
movups 0x38(%rbp,%rdi,8),%xmm1
movups %xmm0,0x28(%rax,%rdi,8)
movups %xmm1,0x38(%rax,%rdi,8)
movups 0x48(%rbp,%rdi,8),%xmm0
movups 0x58(%rbp,%rdi,8),%xmm1
movups %xmm0,0x48(%rax,%rdi,8)
movups %xmm1,0x58(%rax,%rdi,8)
movups 0x68(%rbp,%rdi,8),%xmm0
movupd 0x78(%rbp,%rdi,8),%xmm1
movups %xmm0,0x68(%rax,%rdi,8)
movupd %xmm1,0x78(%rax,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1270 <nllkexchord+0x110>
test   %rdx,%rdx
je     130b <nllkexchord+0x1ab>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    0x0(,%rdi,8),%rsi
or     $0x8,%rsi
movups 0x0(%rbp,%rsi,1),%xmm0
movupd 0x10(%rbp,%rsi,1),%xmm1
movups %xmm0,(%rax,%rsi,1)
movupd %xmm1,0x10(%rax,%rsi,1)
add    $0x4,%rdi
add    $0xffffffffffffffff,%rdx
jne    12e0 <nllkexchord+0x180>
cmp    %rcx,%r8
je     1398 <nllkexchord+0x238>
or     $0x1,%rcx
mov    %r9d,%esi
sub    %ecx,%esi
mov    %rcx,%rdx
not    %rdx
add    %r9,%rdx
and    $0x3,%rsi
je     1922 <nllkexchord+0x7c2>
mov    0x8(%rsp),%rax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  0x0(%rbp,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    1340 <nllkexchord+0x1e0>
cmp    $0x3,%rdx
jb     1398 <nllkexchord+0x238>
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x8(%rbp,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rax,%rcx,8)
movsd  0x10(%rbp,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rax,%rcx,8)
movsd  0x18(%rbp,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rax,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%r9
jne    1360 <nllkexchord+0x200>
mov    0x2c21(%rip),%rax
mov    (%rax),%r10d
test   %r10d,%r10d
jle    1527 <nllkexchord+0x3c7>
add    $0xffffffff,%r9d
movslq %r9d,%r8
lea    0x1(%r10),%rax
mov    $0x1,%ecx
cmp    $0x4,%r10
jb     149e <nllkexchord+0x33e>
lea    0x8(%rbx),%rdx
lea    (%r8,%rax,1),%rsi
lea    0x0(,%rsi,8),%rsi
add    %rbp,%rsi
cmp    %rsi,%rdx
jae    13f5 <nllkexchord+0x295>
lea    (%rbx,%rax,8),%rdx
lea    0x8(,%r8,8),%rsi
add    %rbp,%rsi
cmp    %rdx,%rsi
jb     149e <nllkexchord+0x33e>
mov    %r10d,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rdx
mov    %rdx,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rdx,%rdx
je     1986 <nllkexchord+0x826>
movabs $0x7ffffffffffffffc,%rdi
ff ff 7f
or     $0x2,%rdi
and    %r9,%rdi
lea    0x0(,%r8,8),%rsi
add    %rbp,%rsi
mov    $0x1,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movups (%rsi,%rdx,8),%xmm0
movups 0x10(%rsi,%rdx,8),%xmm1
movups %xmm0,(%rbx,%rdx,8)
movups %xmm1,0x10(%rbx,%rdx,8)
movups 0x20(%rsi,%rdx,8),%xmm0
movupd 0x30(%rsi,%rdx,8),%xmm1
movups %xmm0,0x20(%rbx,%rdx,8)
movupd %xmm1,0x30(%rbx,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1440 <nllkexchord+0x2e0>
test   $0x1,%r9b
je     1491 <nllkexchord+0x331>
lea    (%rdx,%r8,1),%rsi
movups 0x0(%rbp,%rsi,8),%xmm0
movupd 0x10(%rbp,%rsi,8),%xmm1
movups %xmm0,(%rbx,%rdx,8)
movupd %xmm1,0x10(%rbx,%rdx,8)
cmp    %rcx,%r10
je     1527 <nllkexchord+0x3c7>
or     $0x1,%rcx
mov    %eax,%esi
sub    %ecx,%esi
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
and    $0x3,%rsi
je     14d4 <nllkexchord+0x374>
lea    0x0(,%r8,8),%rdi
add    %rbp,%rdi
nopl   0x0(%rax)
movsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    14c0 <nllkexchord+0x360>
cmp    $0x3,%rdx
jb     1527 <nllkexchord+0x3c7>
lea    0x18(,%r8,8),%rdx
add    %rbp,%rdx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x18(%rdx,%rcx,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
movsd  -0x10(%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rcx,8)
movsd  -0x8(%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rbx,%rcx,8)
movsd  (%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rbx,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    14f0 <nllkexchord+0x390>
movslq 0x10(%rsp),%rax
movsd  0x0(%rbp,%rax,8),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm0,%xmm1
ja     18c9 <nllkexchord+0x769>
ucomisd 0xab8(%rip),%xmm0
jae    18c9 <nllkexchord+0x769>
mov    0x2a93(%rip),%rax
cmpl   $0x0,(%rax)
mov    0x8(%rsp),%rdi
jle    18ea <nllkexchord+0x78a>
mov    0x2a76(%rip),%rax
mov    (%rax),%ebp
add    $0x1,%ebp
mov    $0x1,%r15d
movsd  0xa8b(%rip),%xmm6
xorpd  %xmm5,%xmm5
mov    %r14,0x10(%rsp)
mov    %r12,0x28(%rsp)
movsd  %xmm0,0x48(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x2a29(%rip),%rax
mov    (%rax),%rax
mov    (%rax,%r15,4),%r13d
test   %r13d,%r13d
movsd  %xmm5,0x58(%rsp)
jle    1805 <nllkexchord+0x6a5>
mov    0x2a04(%rip),%rax
mov    (%rax),%rax
mov    %r15,0x20(%rsp)
mov    (%rax,%r15,4),%r9d
mov    0x29e9(%rip),%rax
mov    (%rax),%r8d
mov    0x29c7(%rip),%rax
mov    (%rax),%rax
mov    %rax,0x70(%rsp)
add    $0xffffffff,%r9d
mov    0x29e4(%rip),%rax
mov    (%rax),%r11
mov    0x29b2(%rip),%rax
movslq (%rax),%rax
mov    %rax,0x60(%rsp)
lea    0x1(%r8),%eax
mov    %r13,0x50(%rsp)
add    $0x1,%r13
lea    -0x1(%rax),%r10
add    $0xfffffffffffffffe,%rax
mov    %rax,0x68(%rsp)
mov    %r10d,%esi
and    $0x3,%esi
and    $0xfffffffffffffffc,%r10
mov    $0x1,%r15d
jmp    1675 <nllkexchord+0x515>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x10(%rsp),%r14
movabs $0xc018000000000000,%rax
mov    %rax,(%r14,%r15,8)
mov    0x8(%rsp),%rdi
addsd  0x8(%rdi),%xmm4
mov    0x28(%rsp),%r12
movsd  %xmm4,(%r12,%r15,8)
add    $0x1,%r15
cmp    %r13,%r15
je     17d0 <nllkexchord+0x670>
lea    (%r9,%r15,1),%eax
movq   $0x0,(%rdi)
movslq %eax,%rbp
test   %r8d,%r8d
jle    16c0 <nllkexchord+0x560>
mov    0x70(%rsp),%rax
mov    (%rax,%rbp,8),%r14
cmpq   $0x3,0x68(%rsp)
jae    16d0 <nllkexchord+0x570>
xorpd  %xmm4,%xmm4
mov    $0x1,%r12d
test   %rsi,%rsi
jne    173a <nllkexchord+0x5da>
jmp    1767 <nllkexchord+0x607>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorpd  %xmm4,%xmm4
jmp    1767 <nllkexchord+0x607>
nopl   0x0(%rax)
xorpd  %xmm4,%xmm4
xor    %r12d,%r12d
nopw   0x0(%rax,%rax,1)
movsd  0x8(%rbx,%r12,8),%xmm1
movsd  0x10(%rbx,%r12,8),%xmm2
mulsd  0x8(%r14,%r12,8),%xmm1
mulsd  0x10(%r14,%r12,8),%xmm2
addsd  %xmm4,%xmm1
movsd  0x18(%rbx,%r12,8),%xmm3
mulsd  0x18(%r14,%r12,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x20(%rbx,%r12,8),%xmm4
mulsd  0x20(%r14,%r12,8),%xmm4
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm4
add    $0x4,%r12
cmp    %r12,%r10
jne    16e0 <nllkexchord+0x580>
add    $0x1,%r12
test   %rsi,%rsi
je     1767 <nllkexchord+0x607>
lea    (%r14,%r12,8),%rcx
lea    (%rbx,%r12,8),%rdx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rdx,%rax,8),%xmm1
mulsd  (%rcx,%rax,8),%xmm1
addsd  %xmm1,%xmm4
add    $0x1,%rax
cmp    %rax,%rsi
jne    1750 <nllkexchord+0x5f0>
movslq (%r11,%rbp,4),%rbp
cmp    $0x1,%rbp
je     1640 <nllkexchord+0x4e0>
mov    0x60(%rsp),%rax
cmp    %eax,%ebp
mov    0x10(%rsp),%r14
mov    0x28(%rsp),%r12
mov    0x8(%rsp),%rdi
jne    17b0 <nllkexchord+0x650>
addsd  -0x8(%rdi,%rax,8),%xmm4
movsd  %xmm4,(%r14,%r15,8)
movapd %xmm6,%xmm4
jmp    1662 <nllkexchord+0x502>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  -0x8(%rdi,%rbp,8),%xmm1
addsd  %xmm4,%xmm1
movsd  %xmm1,(%r14,%r15,8)
addsd  (%rdi,%rbp,8),%xmm4
jmp    1662 <nllkexchord+0x502>
nopw   0x0(%rax,%rax,1)
mov    0x50(%rsp),%r13
cmp    $0x2,%r13d
jl     1800 <nllkexchord+0x6a0>
mov    %r13d,%edi
mov    %r14,%rsi
mov    %r12,%rdx
movsd  0x824(%rip),%xmm1
call   1060 <exchmvn@plt>
movapd %xmm0,%xmm2
mov    0x20(%rsp),%r15
jmp    1846 <nllkexchord+0x6e6>
nopl   0x0(%rax)
mov    0x20(%rsp),%r15
test   %ebp,%ebp
jle    1840 <nllkexchord+0x6e0>
movsd  0x8(%r12),%xmm0
call   1040 <pnorms@plt>
movaps %xmm0,0x30(%rsp)
movsd  0x8(%r14),%xmm0
call   1040 <pnorms@plt>
movapd 0x30(%rsp),%xmm2
subsd  %xmm0,%xmm2
jmp    1846 <nllkexchord+0x6e6>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd 0x30(%rsp),%xmm2
movsd  0x7ca(%rip),%xmm0
ucomisd %xmm2,%xmm0
jae    1936 <nllkexchord+0x7d6>
movapd %xmm2,%xmm0
cmpnlesd 0x7c3(%rip),%xmm0
andpd  %xmm0,%xmm2
movsd  0x7af(%rip),%xmm1
andnpd %xmm1,%xmm0
orpd   %xmm0,%xmm2
movapd %xmm2,0x30(%rsp)
movapd %xmm2,%xmm0
call   1080 <log@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  0x58(%rsp),%xmm5
subsd  %xmm0,%xmm5
lea    0x1(%r15),%rax
mov    0x2743(%rip),%rcx
movslq (%rcx),%rcx
cmp    %rcx,%r15
mov    %rax,%r15
mov    0x8(%rsp),%rdi
movsd  0x48(%rsp),%xmm0
movsd  0x747(%rip),%xmm6
jl     15a0 <nllkexchord+0x440>
jmp    18ee <nllkexchord+0x78e>
movabs $0x4202a05f20000000,%rax
a0 02 42
mov    0x18(%rsp),%rcx
mov    %rax,(%rcx)
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xorpd  %xmm5,%xmm5
mov    0x18(%rsp),%rax
movsd  %xmm5,(%rax)
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %r12,%rdi
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>
mov    0x8(%rsp),%rax
cmp    $0x3,%rdx
jae    1360 <nllkexchord+0x200>
jmp    1398 <nllkexchord+0x238>
movabs $0x41cdcd6500000000,%rax
cd cd 41
mov    0x18(%rsp),%rcx
mov    %rax,(%rcx)
test   %r13d,%r13d
jle    18db <nllkexchord+0x77b>
add    $0x1,%r13
mov    $0x1,%ebx
lea    0x6d3(%rip),%r15
nopl   (%rax)
movsd  (%r14,%rbx,8),%xmm0
movsd  (%r12,%rbx,8),%xmm1
mov    %r15,%rdi
mov    %ebx,%esi
mov    %ebx,%edx
call   1050 <printf@plt>
add    $0x1,%rbx
cmp    %rbx,%r13
jne    1960 <nllkexchord+0x800>
jmp    18db <nllkexchord+0x77b>
mov    $0x1,%edx
test   $0x1,%r9b
jne    1478 <nllkexchord+0x318>
jmp    1491 <nllkexchord+0x331>

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

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