#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 */
size_t abs(int);
int /*<<< orphan*/ free(double *);
scalar_t__ log(double);
scalar_t__ malloc(int);
int /*<<< orphan*/ printf(char *, int, double, int, double);
void nllkarord(int np, double *th, double *fnv) {
  extern int *y, *dvec, *dstart, nn, nc, ncl, dmax, norc;
  extern double **x;
  double rh[20][20], *pa, *pb;
  double pr, nlk, eps, sd;
  double *b0, *b, rho, *rr, tem;
  int i, j, k, jk, jc;
  int nsim, ifail, d, startp;
  double pnorms(double);
  void mvn1(int, double[], double[], double[][20], int, double, double *,
            double *, int *);
  pa = (double *)malloc((dmax + 1) * sizeof(double));
  pb = (double *)malloc((dmax + 1) * sizeof(double));
  b0 = (double *)malloc((norc + 1) * sizeof(double));
  b = (double *)malloc((nc + 1) * sizeof(double));
  rr = (double *)malloc((dmax + 1) * sizeof(double));
  nlk = 0.;
  eps = 1.e-6;
  for (i = 1; i < norc; i++)
    b0[i] = th[i];
  for (i = 1; i <= nc; i++)
    b[i] = th[norc - 1 + i];
  rho = th[np];
  if (rho < 0.0 || rho >= 1.) {
    *fnv = 1.e10;
    return;
  }
  for (jk = 1, rr[0] = 1.; jk < dmax; jk++)
    rr[jk] = rr[jk - 1] * rho;
  for (j = 1; j <= dmax; j++) {
    for (k = 1; k <= dmax; k++)
      rh[j][k] = rr[abs(k - j)];
  }
  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++) {
      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) {
      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(b0);
  free(b);
  free(pa);
  free(pb);
  free(rr);
  return;
}


x86 assembly (clang, -O0):
<nllkarord>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    0x2e53(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0xca8(%rbp)
mov    0x2e33(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0xcb0(%rbp)
mov    0x2de3(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0xcd8(%rbp)
mov    0x2dd3(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0xce0(%rbp)
mov    0x2dd3(%rip),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   1080 <malloc@plt>
mov    %rax,-0xcf0(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xcc0(%rbp)
ff
movsd  0xdc7(%rip),%xmm0
movsd  %xmm0,-0xcc8(%rbp)
ff
movl   $0x1,-0xcfc(%rbp)
mov    -0xcfc(%rbp),%eax
mov    0x2d58(%rip),%rcx
cmp    (%rcx),%eax
jge    1297 <nllkarord+0x127>
mov    -0x10(%rbp),%rax
movslq -0xcfc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xcd8(%rbp),%rax
movslq -0xcfc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xcfc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcfc(%rbp)
jmp    124b <nllkarord+0xdb>
movl   $0x1,-0xcfc(%rbp)
mov    -0xcfc(%rbp),%eax
mov    0x2d12(%rip),%rcx
cmp    (%rcx),%eax
jg     12fb <nllkarord+0x18b>
mov    -0x10(%rbp),%rax
mov    0x2cef(%rip),%rcx
mov    (%rcx),%ecx
sub    $0x1,%ecx
add    -0xcfc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xce0(%rbp),%rax
movslq -0xcfc(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xcfc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcfc(%rbp)
jmp    12a1 <nllkarord+0x131>
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     133b <nllkarord+0x1cb>
movsd  -0xce8(%rbp),%xmm0
ff
movsd  0xcd7(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1350 <nllkarord+0x1e0>
mov    -0x18(%rbp),%rax
movsd  0xcf1(%rip),%xmm0
movsd  %xmm0,(%rax)
jmp    1949 <nllkarord+0x7d9>
movl   $0x1,-0xd08(%rbp)
mov    -0xcf0(%rbp),%rax
movsd  0xc9f(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0xd08(%rbp),%eax
mov    0x2c66(%rip),%rcx
cmp    (%rcx),%eax
jge    13c9 <nllkarord+0x259>
mov    -0xcf0(%rbp),%rax
mov    -0xd08(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0xce8(%rbp),%xmm0
ff
mov    -0xcf0(%rbp),%rax
movslq -0xd08(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xd08(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd08(%rbp)
jmp    136d <nllkarord+0x1fd>
movl   $0x1,-0xd00(%rbp)
mov    -0xd00(%rbp),%eax
mov    0x2c00(%rip),%rcx
cmp    (%rcx),%eax
jg     1486 <nllkarord+0x316>
movl   $0x1,-0xd04(%rbp)
mov    -0xd04(%rbp),%eax
mov    0x2be1(%rip),%rcx
cmp    (%rcx),%eax
jg     146d <nllkarord+0x2fd>
mov    -0xcf0(%rbp),%rax
mov    %rax,-0xd28(%rbp)
mov    -0xd04(%rbp),%edi
sub    -0xd00(%rbp),%edi
call   1040 <abs@plt>
mov    %rax,%rcx
mov    -0xd28(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm0
movslq -0xd00(%rbp),%rcx
lea    -0xca0(%rbp),%rax
imul   $0xa0,%rcx,%rcx
add    %rcx,%rax
movslq -0xd04(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xd04(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd04(%rbp)
jmp    13f2 <nllkarord+0x282>
jmp    1472 <nllkarord+0x302>
mov    -0xd00(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd00(%rbp)
jmp    13d3 <nllkarord+0x263>
movl   $0x1,-0xcfc(%rbp)
mov    -0xcfc(%rbp),%eax
mov    0x2b4b(%rip),%rcx
cmp    (%rcx),%eax
jg     18fd <nllkarord+0x78d>
mov    0x2b24(%rip),%rax
mov    (%rax),%rax
movslq -0xcfc(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xd18(%rbp)
mov    0x2b02(%rip),%rax
mov    (%rax),%rax
movslq -0xcfc(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xd1c(%rbp)
cmpl   $0x7,-0xd18(%rbp)
jg     14f5 <nllkarord+0x385>
movl   $0x0,-0xd10(%rbp)
jmp    14ff <nllkarord+0x38f>
movl   $0x7d0,-0xd10(%rbp)
movl   $0x1,-0xd00(%rbp)
mov    -0xd00(%rbp),%eax
cmp    -0xd18(%rbp),%eax
jg     1717 <nllkarord+0x5a7>
movl   $0x1,-0xd0c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xcf8(%rbp)
ff
mov    -0xd0c(%rbp),%eax
mov    0x2a83(%rip),%rcx
cmp    (%rcx),%eax
jg     15b0 <nllkarord+0x440>
mov    -0xce0(%rbp),%rax
movslq -0xd0c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2a49(%rip),%rax
mov    (%rax),%rax
mov    -0xd1c(%rbp),%ecx
add    -0xd00(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xd0c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xcf8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xcf8(%rbp)
ff
mov    -0xd0c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd0c(%rbp)
jmp    1530 <nllkarord+0x3c0>
mov    0x2a21(%rip),%rax
mov    (%rax),%rax
mov    -0xd1c(%rbp),%ecx
add    -0xd00(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0xd04(%rbp)
cmpl   $0x1,-0xd04(%rbp)
jne    1630 <nllkarord+0x4c0>
mov    -0xca8(%rbp),%rax
movslq -0xd00(%rbp),%rcx
movsd  0xa38(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xcf8(%rbp),%xmm0
ff
mov    -0xcd8(%rbp),%rax
movslq -0xd04(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0xcb0(%rbp),%rax
movslq -0xd00(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    16fe <nllkarord+0x58e>
mov    -0xd04(%rbp),%eax
mov    0x2973(%rip),%rcx
cmp    (%rcx),%eax
jne    1698 <nllkarord+0x528>
movsd  -0xcf8(%rbp),%xmm0
ff
mov    -0xcd8(%rbp),%rax
mov    -0xd04(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0xca8(%rbp),%rax
movslq -0xd00(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xcb0(%rbp),%rax
movslq -0xd00(%rbp),%rcx
movsd  0x99a(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    16f9 <nllkarord+0x589>
movsd  -0xcf8(%rbp),%xmm0
ff
mov    -0xcd8(%rbp),%rax
mov    -0xd04(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0xca8(%rbp),%rax
movslq -0xd00(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xcf8(%rbp),%xmm0
ff
mov    -0xcd8(%rbp),%rax
movslq -0xd04(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0xcb0(%rbp),%rax
movslq -0xd00(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    16fe <nllkarord+0x58e>
jmp    1703 <nllkarord+0x593>
mov    -0xd00(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd00(%rbp)
jmp    1509 <nllkarord+0x399>
cmpl   $0x2,-0xd18(%rbp)
jl     17c8 <nllkarord+0x658>
mov    -0xd18(%rbp),%edi
mov    -0xca8(%rbp),%rsi
mov    -0xcb0(%rbp),%rdx
lea    -0xca0(%rbp),%rcx
mov    -0xd10(%rbp),%r8d
movsd  -0xcc8(%rbp),%xmm0
ff
lea    -0xcb8(%rbp),%r9
lea    -0xcd0(%rbp),%r10
lea    -0xd14(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
call   1060 <mvn1@plt>
cmpl   $0x0,-0xd14(%rbp)
jle    17c3 <nllkarord+0x653>
mov    -0x18(%rbp),%rax
movsd  0x896(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0xce0(%rbp),%rdi
call   1030 <free@plt>
mov    -0xca8(%rbp),%rdi
call   1030 <free@plt>
mov    -0xcb0(%rbp),%rdi
call   1030 <free@plt>
mov    -0xcf0(%rbp),%rdi
call   1030 <free@plt>
jmp    1949 <nllkarord+0x7d9>
jmp    1809 <nllkarord+0x699>
mov    -0xcb0(%rbp),%rax
movsd  0x8(%rax),%xmm0
call   1050 <pnorms@plt>
movsd  %xmm0,-0xd30(%rbp)
ff
mov    -0xca8(%rbp),%rax
movsd  0x8(%rax),%xmm0
call   1050 <pnorms@plt>
movaps %xmm0,%xmm1
movsd  -0xd30(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xcb8(%rbp)
ff
movsd  0x7ff(%rip),%xmm0
ucomisd -0xcb8(%rbp),%xmm0
ff
jb     18a2 <nllkarord+0x732>
mov    -0x18(%rbp),%rax
movsd  0x7f5(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x1,-0xd00(%rbp)
mov    -0xd00(%rbp),%eax
cmp    -0xd18(%rbp),%eax
jg     189d <nllkarord+0x72d>
mov    -0xd00(%rbp),%esi
mov    -0xca8(%rbp),%rax
movslq -0xd00(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0xd00(%rbp),%edx
mov    -0xcb0(%rbp),%rax
movslq -0xd00(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
lea    0x7bc(%rip),%rdi
call   1070 <printf@plt>
mov    -0xd00(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xd00(%rbp)
jmp    1839 <nllkarord+0x6c9>
jmp    1949 <nllkarord+0x7d9>
xorps  %xmm0,%xmm0
ucomisd -0xcb8(%rbp),%xmm0
ff
jb     18c3 <nllkarord+0x753>
movsd  0x75d(%rip),%xmm0
movsd  %xmm0,-0xcb8(%rbp)
ff
movsd  -0xcb8(%rbp),%xmm0
ff
call   1090 <log@plt>
cvtsi2sd %rax,%xmm1
movsd  -0xcc0(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xcc0(%rbp)
ff
mov    -0xcfc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xcfc(%rbp)
jmp    1490 <nllkarord+0x320>
movsd  -0xcc0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0xcd8(%rbp),%rdi
call   1030 <free@plt>
mov    -0xce0(%rbp),%rdi
call   1030 <free@plt>
mov    -0xca8(%rbp),%rdi
call   1030 <free@plt>
mov    -0xcb0(%rbp),%rdi
call   1030 <free@plt>
mov    -0xcf0(%rbp),%rdi
call   1030 <free@plt>
add    $0xd40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<nllkarord>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xcf8,%rsp
mov    %rdx,0x10(%rsp)
mov    %rsi,%r15
mov    %edi,0x8(%rsp)
mov    0x2e4c(%rip),%r12
mov    (%r12),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,0x20(%rsp)
mov    (%r12),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,%r13
mov    0x2ded(%rip),%rbx
mov    (%rbx),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,%r14
mov    0x2de5(%rip),%rax
mov    (%rax),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,%rbp
mov    (%r12),%eax
lea    0x8(,%rax,8),%edi
call   1080 <malloc@plt>
mov    %rax,%rsi
mov    (%rbx),%eax
cmp    $0x2,%eax
jl     139f <nllkarord+0x22f>
lea    -0x1(%rax),%r8
mov    $0x1,%ecx
cmp    $0x4,%r8
jb     1327 <nllkarord+0x1b7>
lea    0x8(%r14),%rdx
lea    (%r15,%rax,8),%rdi
cmp    %rdi,%rdx
jae    123b <nllkarord+0xcb>
lea    (%r14,%rax,8),%rdx
lea    0x8(%r15),%rdi
cmp    %rdx,%rdi
jb     1327 <nllkarord+0x1b7>
mov    %r8,%rcx
and    $0xfffffffffffffffc,%rcx
lea    -0x4(%rcx),%rdi
mov    %rdi,%rbx
shr    $0x2,%rbx
add    $0x1,%rbx
mov    %ebx,%edx
and    $0x3,%edx
cmp    $0xc,%rdi
jae    1260 <nllkarord+0xf0>
xor    %edi,%edi
jmp    12dc <nllkarord+0x16c>
and    $0xfffffffffffffffc,%rbx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movups 0x8(%r15,%rdi,8),%xmm0
movups 0x18(%r15,%rdi,8),%xmm1
movups %xmm0,0x8(%r14,%rdi,8)
movups %xmm1,0x18(%r14,%rdi,8)
movups 0x28(%r15,%rdi,8),%xmm0
movups 0x38(%r15,%rdi,8),%xmm1
movups %xmm0,0x28(%r14,%rdi,8)
movups %xmm1,0x38(%r14,%rdi,8)
movups 0x48(%r15,%rdi,8),%xmm0
movups 0x58(%r15,%rdi,8),%xmm1
movups %xmm0,0x48(%r14,%rdi,8)
movups %xmm1,0x58(%r14,%rdi,8)
movups 0x68(%r15,%rdi,8),%xmm0
movupd 0x78(%r15,%rdi,8),%xmm1
movups %xmm0,0x68(%r14,%rdi,8)
movupd %xmm1,0x78(%r14,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rbx
jne    1270 <nllkarord+0x100>
test   %rdx,%rdx
je     131e <nllkarord+0x1ae>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    0x0(,%rdi,8),%rbx
or     $0x8,%rbx
movups (%r15,%rbx,1),%xmm0
movupd 0x10(%r15,%rbx,1),%xmm1
movups %xmm0,(%r14,%rbx,1)
movupd %xmm1,0x10(%r14,%rbx,1)
add    $0x4,%rdi
add    $0xffffffffffffffff,%rdx
jne    12f0 <nllkarord+0x180>
cmp    %rcx,%r8
je     139f <nllkarord+0x22f>
or     $0x1,%rcx
mov    %eax,%edi
sub    %ecx,%edi
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
and    $0x3,%rdi
je     1356 <nllkarord+0x1e6>
nopw   0x0(%rax,%rax,1)
movsd  (%r15,%rcx,8),%xmm0
movsd  %xmm0,(%r14,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rdi
jne    1340 <nllkarord+0x1d0>
cmp    $0x3,%rdx
jb     139f <nllkarord+0x22f>
nopl   0x0(%rax)
movsd  (%r15,%rcx,8),%xmm0
movsd  %xmm0,(%r14,%rcx,8)
movsd  0x8(%r15,%rcx,8),%xmm0
movsd  %xmm0,0x8(%r14,%rcx,8)
movsd  0x10(%r15,%rcx,8),%xmm0
movsd  %xmm0,0x10(%r14,%rcx,8)
movsd  0x18(%r15,%rcx,8),%xmm0
movsd  %xmm0,0x18(%r14,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    1360 <nllkarord+0x1f0>
mov    0x2c1a(%rip),%rcx
mov    (%rcx),%r10d
test   %r10d,%r10d
jle    1528 <nllkarord+0x3b8>
add    $0xffffffff,%eax
movslq %eax,%r8
lea    0x1(%r10),%rax
mov    $0x1,%ecx
cmp    $0x4,%r10
jb     14a5 <nllkarord+0x335>
lea    0x8(%rbp),%rdx
lea    (%r8,%rax,1),%rdi
lea    (%r15,%rdi,8),%rdi
cmp    %rdi,%rdx
jae    13f8 <nllkarord+0x288>
lea    0x0(,%rax,8),%rdx
add    %rbp,%rdx
lea    (%r15,%r8,8),%rdi
add    $0x8,%rdi
cmp    %rdx,%rdi
jb     14a5 <nllkarord+0x335>
mov    %r10d,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rdx
mov    %rdx,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rdx,%rdx
je     1a93 <nllkarord+0x923>
movabs $0x7ffffffffffffffc,%rdx
ff ff 7f
or     $0x2,%rdx
and    %r9,%rdx
lea    (%r15,%r8,8),%rdi
add    $0x38,%rdi
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movups -0x30(%rdi,%rbx,8),%xmm0
movups -0x20(%rdi,%rbx,8),%xmm1
movups %xmm0,0x8(%rbp,%rbx,8)
movups %xmm1,0x18(%rbp,%rbx,8)
movups -0x10(%rdi,%rbx,8),%xmm0
movupd (%rdi,%rbx,8),%xmm1
movups %xmm0,0x28(%rbp,%rbx,8)
movupd %xmm1,0x38(%rbp,%rbx,8)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%rdx
jne    1440 <nllkarord+0x2d0>
test   $0x1,%r9b
je     1498 <nllkarord+0x328>
or     $0x1,%rbx
lea    (%rbx,%r8,1),%rdx
movups (%r15,%rdx,8),%xmm0
movupd 0x10(%r15,%rdx,8),%xmm1
movups %xmm0,0x0(%rbp,%rbx,8)
movupd %xmm1,0x10(%rbp,%rbx,8)
cmp    %rcx,%r10
je     1528 <nllkarord+0x3b8>
or     $0x1,%rcx
mov    %eax,%edx
sub    %ecx,%edx
mov    %rcx,%rbx
not    %rbx
add    %rax,%rbx
and    $0x3,%rdx
je     14d5 <nllkarord+0x365>
lea    (%r15,%r8,8),%rdi
nopl   0x0(%rax)
movsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rdx
jne    14c0 <nllkarord+0x350>
cmp    $0x3,%rbx
jb     1528 <nllkarord+0x3b8>
lea    (%r15,%r8,8),%rdx
add    $0x18,%rdx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  -0x18(%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rcx,8)
movsd  -0x10(%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rbp,%rcx,8)
movsd  -0x8(%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rbp,%rcx,8)
movsd  (%rdx,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rbp,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    14f0 <nllkarord+0x380>
movslq 0x8(%rsp),%rax
movsd  (%r15,%rax,8),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm0,%xmm1
ja     1595 <nllkarord+0x425>
ucomisd 0xabb(%rip),%xmm0
jae    1595 <nllkarord+0x425>
mov    %r13,0x8(%rsp)
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rsi)
mov    0x2a80(%rip),%rax
mov    (%rax),%eax
cmp    $0x2,%eax
jl     1612 <nllkarord+0x4a2>
lea    -0x1(%rax),%rdi
lea    -0x2(%rax),%rdx
mov    %edi,%ecx
and    $0x3,%ecx
cmp    $0x3,%rdx
jae    15ac <nllkarord+0x43c>
mov    $0x1,%edx
movsd  0xa75(%rip),%xmm1
test   %rcx,%rcx
jne    15fa <nllkarord+0x48a>
jmp    1612 <nllkarord+0x4a2>
movabs $0x4202a05f20000000,%rax
a0 02 42
mov    0x10(%rsp),%rcx
mov    %rax,(%rcx)
jmp    1a28 <nllkarord+0x8b8>
and    $0xfffffffffffffffc,%rdi
movsd  0xa48(%rip),%xmm1
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsi,%rdx,8)
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsi,%rdx,8)
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rsi,%rdx,8)
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsi,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rdi
jne    15c0 <nllkarord+0x450>
add    $0x1,%rdx
test   %rcx,%rcx
je     1612 <nllkarord+0x4a2>
lea    (%rsi,%rdx,8),%rdx
xor    %edi,%edi
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rdi,8)
add    $0x1,%rdi
cmp    %rdi,%rcx
jne    1600 <nllkarord+0x490>
test   %eax,%eax
mov    %r14,0x18(%rsp)
jle    1680 <nllkarord+0x510>
lea    0x118(%rsp),%r15
mov    $0x1,%r13d
xor    %r12d,%r12d
jmp    164a <nllkarord+0x4da>
xchg   %ax,%ax
movslq %eax,%rcx
lea    0x1(%r13),%rdx
add    $0xffffffffffffffff,%r12
add    $0xa0,%r15
cmp    %rcx,%r13
mov    %rdx,%r13
jge    1680 <nllkarord+0x510>
test   %eax,%eax
jle    1630 <nllkarord+0x4c0>
xor    %ebx,%ebx
lea    (%r12,%rbx,1),%edi
mov    %rsi,%r14
call   1040 <abs@plt>
mov    %r14,%rsi
movsd  (%r14,%rax,8),%xmm0
movsd  %xmm0,(%r15,%rbx,8)
mov    0x296e(%rip),%rax
movslq (%rax),%rax
add    $0x1,%rbx
cmp    %rax,%rbx
jl     1650 <nllkarord+0x4e0>
jmp    1630 <nllkarord+0x4c0>
mov    0x2961(%rip),%rax
cmpl   $0x0,(%rax)
mov    %rsi,0x40(%rsp)
jle    19e5 <nllkarord+0x875>
mov    $0x1,%ebx
movsd  0x966(%rip),%xmm5
xorpd  %xmm4,%xmm4
mov    0x20(%rsp),%r14
mov    0x8(%rsp),%r12
jmp    16fd <nllkarord+0x58d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm1,%xmm0
call   1090 <log@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  0x58(%rsp),%xmm4
subsd  %xmm0,%xmm4
lea    0x1(%rbx),%rax
mov    0x2902(%rip),%rcx
movslq (%rcx),%rcx
cmp    %rcx,%rbx
mov    %rax,%rbx
movsd  0x911(%rip),%xmm5
jge    19f3 <nllkarord+0x883>
mov    0x28cc(%rip),%rax
mov    (%rax),%rax
mov    (%rax,%rbx,4),%r15d
cmp    $0x8,%r15d
mov    $0x7d0,%ecx
mov    $0x0,%eax
cmovl  %eax,%ecx
test   %r15d,%r15d
movsd  %xmm4,0x58(%rsp)
jle    1980 <nllkarord+0x810>
mov    %ecx,0x38(%rsp)
mov    0x2892(%rip),%rax
mov    (%rax),%rax
mov    %rbx,0x50(%rsp)
mov    (%rax,%rbx,4),%eax
mov    0x2878(%rip),%rcx
mov    (%rcx),%r8d
mov    0x2856(%rip),%rcx
mov    (%rcx),%rcx
mov    %rcx,0x28(%rsp)
add    $0xffffffff,%eax
mov    0x2874(%rip),%rcx
mov    (%rcx),%r9
mov    0x2842(%rip),%rcx
movslq (%rcx),%r11
lea    0x1(%r8),%ecx
mov    %r15,0x48(%rsp)
lea    0x1(%r15),%r10
lea    -0x1(%rcx),%rdi
add    $0xfffffffffffffffe,%rcx
mov    %rcx,0x60(%rsp)
mov    %edi,%ebx
and    $0x3,%ebx
and    $0xfffffffffffffffc,%rdi
mov    $0x1,%r15d
jmp    17d5 <nllkarord+0x665>
nopl   0x0(%rax)
mov    0x20(%rsp),%r14
movabs $0xc018000000000000,%rcx
mov    %rcx,(%r14,%r15,8)
mov    0x18(%rsp),%rcx
addsd  0x8(%rcx),%xmm0
mov    0x8(%rsp),%r12
movsd  %xmm0,(%r12,%r15,8)
add    $0x1,%r15
cmp    %r10,%r15
je     1920 <nllkarord+0x7b0>
lea    (%rax,%r15,1),%ecx
movslq %ecx,%r14
test   %r8d,%r8d
jle    1810 <nllkarord+0x6a0>
mov    0x28(%rsp),%rcx
mov    (%rcx,%r14,8),%r13
cmpq   $0x3,0x60(%rsp)
jae    1820 <nllkarord+0x6b0>
xorpd  %xmm0,%xmm0
mov    $0x1,%r12d
test   %rbx,%rbx
jne    188a <nllkarord+0x71a>
jmp    18c7 <nllkarord+0x757>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
jmp    18c7 <nllkarord+0x757>
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
xor    %r12d,%r12d
nopw   0x0(%rax,%rax,1)
movsd  0x8(%rbp,%r12,8),%xmm1
movsd  0x10(%rbp,%r12,8),%xmm2
mulsd  0x8(%r13,%r12,8),%xmm1
mulsd  0x10(%r13,%r12,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x18(%rbp,%r12,8),%xmm3
mulsd  0x18(%r13,%r12,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x20(%rbp,%r12,8),%xmm0
mulsd  0x20(%r13,%r12,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%r12
cmp    %r12,%rdi
jne    1830 <nllkarord+0x6c0>
add    $0x1,%r12
test   %rbx,%rbx
je     18c7 <nllkarord+0x757>
lea    0x0(,%r12,8),%rsi
add    %r13,%rsi
lea    0x0(,%r12,8),%rcx
add    %rbp,%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rcx,%rdx,8),%xmm1
mulsd  (%rsi,%rdx,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%rbx
jne    18b0 <nllkarord+0x740>
movslq (%r9,%r14,4),%rcx
cmp    $0x1,%rcx
je     17a0 <nllkarord+0x630>
cmp    %r11d,%ecx
mov    0x20(%rsp),%r14
mov    0x8(%rsp),%r12
jne    1900 <nllkarord+0x790>
mov    0x18(%rsp),%rcx
addsd  -0x8(%rcx,%r11,8),%xmm0
movsd  %xmm0,(%r14,%r15,8)
movapd %xmm5,%xmm0
jmp    17c2 <nllkarord+0x652>
nop
mov    0x18(%rsp),%rdx
movsd  -0x8(%rdx,%rcx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%r15,8)
addsd  (%rdx,%rcx,8),%xmm0
jmp    17c2 <nllkarord+0x652>
nop
mov    0x48(%rsp),%r15
cmp    $0x2,%r15d
mov    0x50(%rsp),%rbx
mov    0x38(%rsp),%r8d
jl     1980 <nllkarord+0x810>
mov    %r15d,%edi
mov    %r14,%rsi
mov    %r12,%rdx
lea    0x70(%rsp),%rcx
movsd  0x6c5(%rip),%xmm0
lea    0x30(%rsp),%r9
lea    0x3c(%rsp),%rax
push   %rax
lea    0x70(%rsp),%rax
push   %rax
call   1060 <mvn1@plt>
add    $0x10,%rsp
cmpl   $0x0,0x3c(%rsp)
jg     1aa4 <nllkarord+0x934>
movsd  0x30(%rsp),%xmm1
jmp    19ad <nllkarord+0x83d>
nopl   0x0(%rax,%rax,1)
movsd  0x8(%r12),%xmm0
call   1050 <pnorms@plt>
movsd  %xmm0,0x28(%rsp)
movsd  0x8(%r14),%xmm0
call   1050 <pnorms@plt>
movsd  0x28(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x30(%rsp)
movsd  0x663(%rip),%xmm0
ucomisd %xmm1,%xmm0
jae    1a3a <nllkarord+0x8ca>
xorpd  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jb     16c0 <nllkarord+0x550>
movabs $0x3cd203af9ee75616,%rax
03 d2 3c
mov    %rax,0x30(%rsp)
movsd  0x640(%rip),%xmm1
jmp    16c0 <nllkarord+0x550>
xorpd  %xmm4,%xmm4
mov    0x20(%rsp),%r14
mov    0x8(%rsp),%r12
mov    0x10(%rsp),%rax
movsd  %xmm4,(%rax)
mov    0x18(%rsp),%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    0x40(%rsp),%rdi
call   1030 <free@plt>
add    $0xcf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
movabs $0x41cdcd6500000000,%rax
cd cd 41
mov    0x10(%rsp),%rcx
mov    %rax,(%rcx)
test   %r15d,%r15d
jle    1a28 <nllkarord+0x8b8>
add    $0x1,%r15
mov    $0x1,%ebx
lea    0x5cf(%rip),%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%r14,%rbx,8),%xmm0
movsd  (%r12,%rbx,8),%xmm1
mov    %rbp,%rdi
mov    %ebx,%esi
mov    %ebx,%edx
call   1070 <printf@plt>
add    $0x1,%rbx
cmp    %rbx,%r15
jne    1a70 <nllkarord+0x900>
jmp    1a28 <nllkarord+0x8b8>
xor    %ebx,%ebx
test   $0x1,%r9b
jne    1479 <nllkarord+0x309>
jmp    1498 <nllkarord+0x328>
movabs $0x41cdcd6500000000,%rax
cd cd 41
mov    0x10(%rsp),%rcx
mov    %rax,(%rcx)
jmp    1a06 <nllkarord+0x896>

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

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