#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 doublemax(int, double *);
double doublemin(int, double *);
double nthd(int ndat, double *dat, double n) {
  double a, temp;
  int i, nmax;
  double s;
  int nsame;
  if (n >= ndat - 1)
    return dat[ndat - 1];
  if (ndat == 1 || n < 0)
    return dat[0];
  if (ndat <= 2) {
    if (dat[0] < dat[1])
      return (dat[0] * (1. - n) + n * dat[1]);
    else
      return (dat[1] * (1. - n) + n * dat[0]);
  }
  a = dat[0];
  nmax = ndat;
  nsame = 1;
  for (i = 1; i < nmax; i++) {
    if (dat[i] > a) {
      temp = dat[i];
      dat[i] = dat[nmax - 1];
      dat[nmax - 1] = temp;
      nmax--;
      i--;
    } else if (dat[i] == a) {
      temp = dat[i];
      dat[i] = dat[nsame];
      dat[nsame] = temp;
      nsame++;
    }
  }
  if (nmax == ndat) {
    if (nmax - nsame <= n) {
      return a;
    } else if (nmax - nsame - 1 > n) {
      return nthd(nmax - nsame, dat + nsame, n);
    } else {
      s = n - (nmax - nsame - 1);
      return (doublemax(nmax - nsame, dat + nsame) * (1. - s) + s * dat[0]);
    }
  }
  if (n < nmax - nsame - 1) {
    return nthd(nmax - nsame, dat + nsame, n);
  } else if (n < nmax - nsame) {
    s = n - (nmax - nsame - 1);
    return (doublemax(nmax - nsame, dat + nsame) * (1. - s) + s * dat[0]);
  } else if (n <= nmax - 1) {
    return a;
  } else if (n < nmax) {
    s = n - (nmax - 1);
    return (dat[0] * (1. - s) + s * doublemin(ndat - nmax, dat + nmax));
  } else {
    return nthd(ndat - nmax, dat + nmax, n - nmax);
  }
}


x86 assembly (clang, -O0):
<nthd>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
ucomisd %xmm1,%xmm0
jb     1179 <nthd+0x49>
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
cmpl   $0x1,-0xc(%rbp)
je     1191 <nthd+0x61>
xorps  %xmm0,%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    11a3 <nthd+0x73>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
cmpl   $0x2,-0xc(%rbp)
jg     1232 <nthd+0x102>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    11fd <nthd+0xcd>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xe28(%rip),%xmm2
subsd  -0x20(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mulsd  0x8(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  0xdf2(%rip),%xmm2
subsd  -0x20(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x1,-0x44(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1353 <nthd+0x223>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    12d6 <nthd+0x1a6>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    1340 <nthd+0x210>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jne    133b <nthd+0x20b>
jp     133b <nthd+0x20b>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1340 <nthd+0x210>
jmp    1345 <nthd+0x215>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1253 <nthd+0x123>
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jne    142d <nthd+0x2fd>
mov    -0x38(%rbp),%eax
sub    -0x44(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     1387 <nthd+0x257>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
mov    -0x38(%rbp),%eax
sub    -0x44(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    13c8 <nthd+0x298>
mov    -0x38(%rbp),%edi
sub    -0x44(%rbp),%edi
mov    -0x18(%rbp),%rsi
movslq -0x44(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
movsd  -0x20(%rbp),%xmm0
call   1050 <nthd@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
movsd  -0x20(%rbp),%xmm0
mov    -0x38(%rbp),%eax
sub    -0x44(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%edi
sub    -0x44(%rbp),%edi
mov    -0x18(%rbp),%rsi
movslq -0x44(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
call   1030 <doublemax@plt>
cvtsi2sd %eax,%xmm0
movsd  0xbf7(%rip),%xmm2
subsd  -0x40(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
movsd  -0x20(%rbp),%xmm1
mov    -0x38(%rbp),%eax
sub    -0x44(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
ucomisd %xmm1,%xmm0
jbe    1472 <nthd+0x342>
mov    -0x38(%rbp),%edi
sub    -0x44(%rbp),%edi
mov    -0x18(%rbp),%rsi
movslq -0x44(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
movsd  -0x20(%rbp),%xmm0
call   1050 <nthd@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
movsd  -0x20(%rbp),%xmm1
mov    -0x38(%rbp),%eax
sub    -0x44(%rbp),%eax
cvtsi2sd %eax,%xmm0
ucomisd %xmm1,%xmm0
jbe    14f0 <nthd+0x3c0>
movsd  -0x20(%rbp),%xmm0
mov    -0x38(%rbp),%eax
sub    -0x44(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%edi
sub    -0x44(%rbp),%edi
mov    -0x18(%rbp),%rsi
movslq -0x44(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
call   1030 <doublemax@plt>
cvtsi2sd %eax,%xmm0
movsd  0xb34(%rip),%xmm2
subsd  -0x40(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
movsd  -0x20(%rbp),%xmm1
mov    -0x38(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
ucomisd %xmm1,%xmm0
jb     1518 <nthd+0x3e8>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
movsd  -0x20(%rbp),%xmm1
cvtsi2sdl -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    15af <nthd+0x47f>
movsd  -0x20(%rbp),%xmm0
mov    -0x38(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xaa7(%rip),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0xc(%rbp),%edi
sub    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
movslq -0x38(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
call   1040 <doublemin@plt>
movsd  -0x60(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm2
movaps %xmm0,%xmm3
movsd  -0x50(%rbp),%xmm0
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    15dc <nthd+0x4ac>
mov    -0xc(%rbp),%edi
sub    -0x38(%rbp),%edi
mov    -0x18(%rbp),%rsi
movslq -0x38(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x38(%rbp),%xmm1
subsd  %xmm1,%xmm0
call   1050 <nthd@plt>
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<nthd>:
push   %rbx
sub    $0x10,%rsp
movapd %xmm0,%xmm4
mov    %rsi,%rbx
mov    %edi,%r8d
lea    -0x1(%r8),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
ucomisd %xmm0,%xmm4
jae    128d <nthd+0x16d>
xorpd  %xmm1,%xmm1
mov    %r8d,%edi
jmp    1161 <nthd+0x41>
nopl   (%rax)
movslq %edx,%rcx
lea    (%rbx,%rcx,8),%rbx
ucomisd %xmm2,%xmm4
jae    128d <nthd+0x16d>
cmp    $0x1,%edi
je     129a <nthd+0x17a>
ucomisd %xmm4,%xmm1
ja     129a <nthd+0x17a>
movsd  (%rbx),%xmm0
cmp    $0x3,%edi
jl     12a0 <nthd+0x180>
mov    $0x1,%edx
mov    %edi,%esi
mov    $0x1,%eax
jmp    11b3 <nthd+0x93>
nop
movslq %esi,%rax
add    $0xffffffff,%esi
movsd  -0x8(%rbx,%rax,8),%xmm3
movsd  %xmm3,(%rbx,%rcx,8)
movsd  %xmm2,-0x8(%rbx,%rax,8)
add    $0xffffffff,%ecx
mov    %ecx,%eax
add    $0x1,%eax
cmp    %esi,%eax
jge    11e0 <nthd+0xc0>
movslq %eax,%rcx
movsd  (%rbx,%rcx,8),%xmm2
ucomisd %xmm0,%xmm2
ja     1190 <nthd+0x70>
jne    11ac <nthd+0x8c>
jp     11ac <nthd+0x8c>
movslq %edx,%rdx
movsd  (%rbx,%rdx,8),%xmm3
movsd  %xmm3,(%rbx,%rcx,8)
movsd  %xmm2,(%rbx,%rdx,8)
add    $0x1,%edx
jmp    11ac <nthd+0x8c>
nopl   0x0(%rax)
mov    %edi,%r8d
sub    %esi,%r8d
jne    1220 <nthd+0x100>
sub    %edx,%edi
xorps  %xmm2,%xmm2
cvtsi2sd %edi,%xmm2
ucomisd %xmm2,%xmm4
jae    1294 <nthd+0x174>
lea    -0x1(%rdi),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
ucomisd %xmm4,%xmm2
ja     1150 <nthd+0x30>
jmp    12d7 <nthd+0x1b7>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %esi,%ecx
sub    %edx,%ecx
lea    -0x1(%rcx),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
ucomisd %xmm4,%xmm2
jbe    123b <nthd+0x11b>
mov    %ecx,%edi
jmp    1150 <nthd+0x30>
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
ucomisd %xmm4,%xmm3
ja     1312 <nthd+0x1f2>
lea    -0x1(%rsi),%eax
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
ucomisd %xmm4,%xmm2
jae    1294 <nthd+0x174>
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
ucomisd %xmm4,%xmm0
ja     1326 <nthd+0x206>
movslq %esi,%rax
lea    (%rbx,%rax,8),%rbx
subsd  %xmm0,%xmm4
lea    -0x1(%r8),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
ucomisd %xmm0,%xmm4
jb     1148 <nthd+0x28>
cltq
movsd  (%rbx,%rax,8),%xmm0
add    $0x10,%rsp
pop    %rbx
ret
movsd  (%rbx),%xmm0
jmp    1294 <nthd+0x174>
movsd  0x8(%rbx),%xmm1
ucomisd %xmm0,%xmm1
movsd  0xd4f(%rip),%xmm2
subsd  %xmm4,%xmm2
jbe    12c5 <nthd+0x1a5>
mulsd  %xmm1,%xmm4
mulsd  %xmm2,%xmm0
addsd  %xmm4,%xmm0
jmp    1294 <nthd+0x174>
mulsd  %xmm0,%xmm4
mulsd  %xmm2,%xmm1
addsd  %xmm4,%xmm1
movapd %xmm1,%xmm0
jmp    1294 <nthd+0x174>
subsd  %xmm2,%xmm4
movsd  %xmm4,(%rsp)
movslq %edx,%rax
lea    (%rbx,%rax,8),%rsi
call   1030 <doublemax@plt>
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  0xd05(%rip),%xmm1
movsd  (%rsp),%xmm0
subsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm1
mulsd  (%rbx),%xmm0
addsd  %xmm1,%xmm0
jmp    1294 <nthd+0x174>
subsd  %xmm2,%xmm4
movsd  %xmm4,(%rsp)
movslq %edx,%rax
lea    (%rbx,%rax,8),%rsi
mov    %ecx,%edi
jmp    12e7 <nthd+0x1c7>
subsd  %xmm2,%xmm4
movsd  %xmm4,(%rsp)
movsd  0xcc9(%rip),%xmm0
subsd  %xmm4,%xmm0
sub    %esi,%edi
movslq %esi,%rax
lea    (%rbx,%rax,8),%rsi
mulsd  (%rbx),%xmm0
movsd  %xmm0,0x8(%rsp)
call   1040 <doublemin@plt>
mulsd  (%rsp),%xmm0
addsd  0x8(%rsp),%xmm0
jmp    1294 <nthd+0x174>

LLVM IR (-O0):
define double @nthd(i32 %0, double* %1, double %2) {
%4 = alloca double
%5 = alloca i32
%6 = alloca double*
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca i32
store i32 %0, i32* %5
store double* %1, double** %6
store double %2, double* %7
%14 = load double, double* %7
%15 = load i32, i32* %5
%16 = sub nsw i32 %15, 1
%17 = sitofp i32 %16 to double
%18 = fcmp oge double %14, %17
br i1 %18, label %19, label %26
19: ; preds = %3
%20 = load double*, double** %6
%21 = load i32, i32* %5
%22 = sub nsw i32 %21, 1
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %20, i64 %23
%25 = load double, double* %24
store double %25, double* %4
br label %303
26: ; preds = %3
%27 = load i32, i32* %5
%28 = icmp eq i32 %27, 1
br i1 %28, label %32, label %29
29: ; preds = %26
%30 = load double, double* %7
%31 = fcmp olt double %30, 0.000000e+00
br i1 %31, label %32, label %36
32: ; preds = %29, %26
%33 = load double*, double** %6
%34 = getelementptr inbounds double, double* %33, i64 0
%35 = load double, double* %34
store double %35, double* %4
br label %303
36: ; preds = %29
%37 = load i32, i32* %5
%38 = icmp sle i32 %37, 2
br i1 %38, label %39, label %71
39: ; preds = %36
%40 = load double*, double** %6
%41 = getelementptr inbounds double, double* %40, i64 0
%42 = load double, double* %41
%43 = load double*, double** %6
%44 = getelementptr inbounds double, double* %43, i64 1
%45 = load double, double* %44
%46 = fcmp olt double %42, %45
br i1 %46, label %47, label %59
47: ; preds = %39
%48 = load double*, double** %6
%49 = getelementptr inbounds double, double* %48, i64 0
%50 = load double, double* %49
%51 = load double, double* %7
%52 = fsub double 1.000000e+00, %51
%53 = load double, double* %7
%54 = load double*, double** %6
%55 = getelementptr inbounds double, double* %54, i64 1
%56 = load double, double* %55
%57 = fmul double %53, %56
%58 = call double @llvm.fmuladd.f64(double %50, double %52, double %57)
store double %58, double* %4
br label %303
59: ; preds = %39
%60 = load double*, double** %6
%61 = getelementptr inbounds double, double* %60, i64 1
%62 = load double, double* %61
%63 = load double, double* %7
%64 = fsub double 1.000000e+00, %63
%65 = load double, double* %7
%66 = load double*, double** %6
%67 = getelementptr inbounds double, double* %66, i64 0
%68 = load double, double* %67
%69 = fmul double %65, %68
%70 = call double @llvm.fmuladd.f64(double %62, double %64, double %69)
store double %70, double* %4
br label %303
71: ; preds = %36
%72 = load double*, double** %6
%73 = getelementptr inbounds double, double* %72, i64 0
%74 = load double, double* %73
store double %74, double* %8
%75 = load i32, i32* %5
store i32 %75, i32* %11
store i32 1, i32* %13
store i32 1, i32* %10
br label %76
76: ; preds = %146, %71
%77 = load i32, i32* %10
%78 = load i32, i32* %11
%79 = icmp slt i32 %77, %78
br i1 %79, label %80, label %149
80: ; preds = %76
%81 = load double*, double** %6
%82 = load i32, i32* %10
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds double, double* %81, i64 %83
%85 = load double, double* %84
%86 = load double, double* %8
%87 = fcmp ogt double %85, %86
br i1 %87, label %88, label %114
88: ; preds = %80
%89 = load double*, double** %6
%90 = load i32, i32* %10
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %89, i64 %91
%93 = load double, double* %92
store double %93, double* %9
%94 = load double*, double** %6
%95 = load i32, i32* %11
%96 = sub nsw i32 %95, 1
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %94, i64 %97
%99 = load double, double* %98
%100 = load double*, double** %6
%101 = load i32, i32* %10
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %100, i64 %102
store double %99, double* %103
%104 = load double, double* %9
%105 = load double*, double** %6
%106 = load i32, i32* %11
%107 = sub nsw i32 %106, 1
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds double, double* %105, i64 %108
store double %104, double* %109
%110 = load i32, i32* %11
%111 = add nsw i32 %110, -1
store i32 %111, i32* %11
%112 = load i32, i32* %10
%113 = add nsw i32 %112, -1
store i32 %113, i32* %10
br label %145
114: ; preds = %80
%115 = load double*, double** %6
%116 = load i32, i32* %10
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds double, double* %115, i64 %117
%119 = load double, double* %118
%120 = load double, double* %8
%121 = fcmp oeq double %119, %120
br i1 %121, label %122, label %144
122: ; preds = %114
%123 = load double*, double** %6
%124 = load i32, i32* %10
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds double, double* %123, i64 %125
%127 = load double, double* %126
store double %127, double* %9
%128 = load double*, double** %6
%129 = load i32, i32* %13
%130 = sext i32 %129 to i64
%131 = getelementptr inbounds double, double* %128, i64 %130
%132 = load double, double* %131
%133 = load double*, double** %6
%134 = load i32, i32* %10
%135 = sext i32 %134 to i64
%136 = getelementptr inbounds double, double* %133, i64 %135
store double %132, double* %136
%137 = load double, double* %9
%138 = load double*, double** %6
%139 = load i32, i32* %13
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double, double* %138, i64 %140
store double %137, double* %141
%142 = load i32, i32* %13
%143 = add nsw i32 %142, 1
store i32 %143, i32* %13
br label %144
144: ; preds = %122, %114
br label %145
145: ; preds = %144, %88
br label %146
146: ; preds = %145
%147 = load i32, i32* %10
%148 = add nsw i32 %147, 1
store i32 %148, i32* %10
br label %76
149: ; preds = %76
%150 = load i32, i32* %11
%151 = load i32, i32* %5
%152 = icmp eq i32 %150, %151
br i1 %152, label %153, label %205
153: ; preds = %149
%154 = load i32, i32* %11
%155 = load i32, i32* %13
%156 = sub nsw i32 %154, %155
%157 = sitofp i32 %156 to double
%158 = load double, double* %7
%159 = fcmp ole double %157, %158
br i1 %159, label %160, label %162
160: ; preds = %153
%161 = load double, double* %8
store double %161, double* %4
br label %303
162: ; preds = %153
%163 = load i32, i32* %11
%164 = load i32, i32* %13
%165 = sub nsw i32 %163, %164
%166 = sub nsw i32 %165, 1
%167 = sitofp i32 %166 to double
%168 = load double, double* %7
%169 = fcmp ogt double %167, %168
br i1 %169, label %170, label %180
170: ; preds = %162
%171 = load i32, i32* %11
%172 = load i32, i32* %13
%173 = sub nsw i32 %171, %172
%174 = load double*, double** %6
%175 = load i32, i32* %13
%176 = sext i32 %175 to i64
%177 = getelementptr inbounds double, double* %174, i64 %176
%178 = load double, double* %7
%179 = call double @nthd(i32 %173, double* %177, double %178)
store double %179, double* %4
br label %303
180: ; preds = %162
%181 = load double, double* %7
%182 = load i32, i32* %11
%183 = load i32, i32* %13
%184 = sub nsw i32 %182, %183
%185 = sub nsw i32 %184, 1
%186 = sitofp i32 %185 to double
%187 = fsub double %181, %186
store double %187, double* %12
%188 = load i32, i32* %11
%189 = load i32, i32* %13
%190 = sub nsw i32 %188, %189
%191 = load double*, double** %6
%192 = load i32, i32* %13
%193 = sext i32 %192 to i64
%194 = getelementptr inbounds double, double* %191, i64 %193
%195 = call i32 @doublemax(i32 %190, double* %194)
%196 = sitofp i32 %195 to double
%197 = load double, double* %12
%198 = fsub double 1.000000e+00, %197
%199 = load double, double* %12
%200 = load double*, double** %6
%201 = getelementptr inbounds double, double* %200, i64 0
%202 = load double, double* %201
%203 = fmul double %199, %202
%204 = call double @llvm.fmuladd.f64(double %196, double %198, double %203)
store double %204, double* %4
br label %303
205: ; preds = %149
%206 = load double, double* %7
%207 = load i32, i32* %11
%208 = load i32, i32* %13
%209 = sub nsw i32 %207, %208
%210 = sub nsw i32 %209, 1
%211 = sitofp i32 %210 to double
%212 = fcmp olt double %206, %211
br i1 %212, label %213, label %223
213: ; preds = %205
%214 = load i32, i32* %11
%215 = load i32, i32* %13
%216 = sub nsw i32 %214, %215
%217 = load double*, double** %6
%218 = load i32, i32* %13
%219 = sext i32 %218 to i64
%220 = getelementptr inbounds double, double* %217, i64 %219
%221 = load double, double* %7
%222 = call double @nthd(i32 %216, double* %220, double %221)
store double %222, double* %4
br label %303
223: ; preds = %205
%224 = load double, double* %7
%225 = load i32, i32* %11
%226 = load i32, i32* %13
%227 = sub nsw i32 %225, %226
%228 = sitofp i32 %227 to double
%229 = fcmp olt double %224, %228
br i1 %229, label %230, label %255
230: ; preds = %223
%231 = load double, double* %7
%232 = load i32, i32* %11
%233 = load i32, i32* %13
%234 = sub nsw i32 %232, %233
%235 = sub nsw i32 %234, 1
%236 = sitofp i32 %235 to double
%237 = fsub double %231, %236
store double %237, double* %12
%238 = load i32, i32* %11
%239 = load i32, i32* %13
%240 = sub nsw i32 %238, %239
%241 = load double*, double** %6
%242 = load i32, i32* %13
%243 = sext i32 %242 to i64
%244 = getelementptr inbounds double, double* %241, i64 %243
%245 = call i32 @doublemax(i32 %240, double* %244)
%246 = sitofp i32 %245 to double
%247 = load double, double* %12
%248 = fsub double 1.000000e+00, %247
%249 = load double, double* %12
%250 = load double*, double** %6
%251 = getelementptr inbounds double, double* %250, i64 0
%252 = load double, double* %251
%253 = fmul double %249, %252
%254 = call double @llvm.fmuladd.f64(double %246, double %248, double %253)
store double %254, double* %4
br label %303
255: ; preds = %223
%256 = load double, double* %7
%257 = load i32, i32* %11
%258 = sub nsw i32 %257, 1
%259 = sitofp i32 %258 to double
%260 = fcmp ole double %256, %259
br i1 %260, label %261, label %263
261: ; preds = %255
%262 = load double, double* %8
store double %262, double* %4
br label %303
263: ; preds = %255
%264 = load double, double* %7
%265 = load i32, i32* %11
%266 = sitofp i32 %265 to double
%267 = fcmp olt double %264, %266
br i1 %267, label %268, label %290
268: ; preds = %263
%269 = load double, double* %7
%270 = load i32, i32* %11
%271 = sub nsw i32 %270, 1
%272 = sitofp i32 %271 to double
%273 = fsub double %269, %272
store double %273, double* %12
%274 = load double*, double** %6
%275 = getelementptr inbounds double, double* %274, i64 0
%276 = load double, double* %275
%277 = load double, double* %12
%278 = fsub double 1.000000e+00, %277
%279 = load double, double* %12
%280 = load i32, i32* %5
%281 = load i32, i32* %11
%282 = sub nsw i32 %280, %281
%283 = load double*, double** %6
%284 = load i32, i32* %11
%285 = sext i32 %284 to i64
%286 = getelementptr inbounds double, double* %283, i64 %285
%287 = call double @doublemin(i32 %282, double* %286)
%288 = fmul double %279, %287
%289 = call double @llvm.fmuladd.f64(double %276, double %278, double %288)
store double %289, double* %4
br label %303
290: ; preds = %263
%291 = load i32, i32* %5
%292 = load i32, i32* %11
%293 = sub nsw i32 %291, %292
%294 = load double*, double** %6
%295 = load i32, i32* %11
%296 = sext i32 %295 to i64
%297 = getelementptr inbounds double, double* %294, i64 %296
%298 = load double, double* %7
%299 = load i32, i32* %11
%300 = sitofp i32 %299 to double
%301 = fsub double %298, %300
%302 = call double @nthd(i32 %293, double* %297, double %301)
store double %302, double* %4
br label %303
303: ; preds = %290, %268, %261, %230, %213, %180, %170, %160, %59, %47, %32, %19
%304 = load double, double* %4
ret double %304
}

LLVM IR (-O2):
define double @nthd(i32 %0, double* %1, double %2) {
%4 = add nsw i32 %0, -1
%5 = sitofp i32 %4 to double
%6 = fcmp ugt double %5, %2
br i1 %6, label %7, label %12
7: ; preds = %3, %128
%8 = phi double [ %132, %128 ], [ %2, %3 ]
%9 = phi double* [ %131, %128 ], [ %1, %3 ]
%10 = phi i32 [ %129, %128 ], [ %0, %3 ]
%11 = fcmp olt double %8, 0.000000e+00
br label %18
12: ; preds = %128, %76, %3
%13 = phi double* [ %1, %3 ], [ %81, %76 ], [ %131, %128 ]
%14 = phi i32 [ %4, %3 ], [ %78, %76 ], [ %133, %128 ]
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %13, i64 %15
%17 = load double, double* %16
br label %136
18: ; preds = %7, %76
%19 = phi double* [ %9, %7 ], [ %81, %76 ]
%20 = phi i32 [ %10, %7 ], [ %79, %76 ]
%21 = icmp eq i32 %20, 1
%22 = or i1 %21, %11
br i1 %22, label %23, label %25
23: ; preds = %18
%24 = load double, double* %19
br label %136
25: ; preds = %18
%26 = icmp slt i32 %20, 3
%27 = load double, double* %19
br i1 %26, label %28, label %39
28: ; preds = %25
%29 = getelementptr inbounds double, double* %19, i64 1
%30 = load double, double* %29
%31 = fcmp olt double %27, %30
%32 = fsub double 1.000000e+00, %8
br i1 %31, label %33, label %36
33: ; preds = %28
%34 = fmul double %8, %30
%35 = tail call double @llvm.fmuladd.f64(double %27, double %32, double %34)
br label %136
36: ; preds = %28
%37 = fmul double %8, %27
%38 = tail call double @llvm.fmuladd.f64(double %30, double %32, double %37)
br label %136
39: ; preds = %25, %60
%40 = phi i32 [ %63, %60 ], [ 1, %25 ]
%41 = phi i32 [ %62, %60 ], [ %20, %25 ]
%42 = phi i32 [ %64, %60 ], [ 1, %25 ]
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %19, i64 %43
%45 = load double, double* %44
%46 = fcmp ogt double %45, %27
br i1 %46, label %47, label %53
47: ; preds = %39
%48 = add nsw i32 %41, -1
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %19, i64 %49
%51 = load double, double* %50
store double %51, double* %44
store double %45, double* %50
%52 = add nsw i32 %42, -1
br label %60
53: ; preds = %39
%54 = fcmp oeq double %45, %27
br i1 %54, label %55, label %60
55: ; preds = %53
%56 = sext i32 %40 to i64
%57 = getelementptr inbounds double, double* %19, i64 %56
%58 = load double, double* %57
store double %58, double* %44
store double %45, double* %57
%59 = add nsw i32 %40, 1
br label %60
60: ; preds = %47, %55, %53
%61 = phi i32 [ %52, %47 ], [ %42, %55 ], [ %42, %53 ]
%62 = phi i32 [ %48, %47 ], [ %41, %55 ], [ %41, %53 ]
%63 = phi i32 [ %40, %47 ], [ %59, %55 ], [ %40, %53 ]
%64 = add nsw i32 %61, 1
%65 = icmp slt i32 %64, %62
br i1 %65, label %39, label %66
66: ; preds = %60
%67 = icmp eq i32 %20, %62
br i1 %67, label %68, label %93
68: ; preds = %66
%69 = sub nsw i32 %20, %63
%70 = sitofp i32 %69 to double
%71 = fcmp ult double %8, %70
br i1 %71, label %72, label %136
72: ; preds = %68
%73 = add nsw i32 %69, -1
%74 = sitofp i32 %73 to double
%75 = fcmp olt double %8, %74
br i1 %75, label %76, label %83
76: ; preds = %93, %72
%77 = phi double [ %74, %72 ], [ %96, %93 ]
%78 = phi i32 [ %73, %72 ], [ %95, %93 ]
%79 = phi i32 [ %69, %72 ], [ %94, %93 ]
%80 = sext i32 %63 to i64
%81 = getelementptr inbounds double, double* %19, i64 %80
%82 = fcmp ult double %8, %77
br i1 %82, label %18, label %12
83: ; preds = %72
%84 = fsub double %8, %74
%85 = sext i32 %63 to i64
%86 = getelementptr inbounds double, double* %19, i64 %85
%87 = tail call i32 @doublemax(i32 %69, double* %86)
%88 = sitofp i32 %87 to double
%89 = fsub double 1.000000e+00, %84
%90 = load double, double* %19
%91 = fmul double %84, %90
%92 = tail call double @llvm.fmuladd.f64(double %88, double %89, double %91)
br label %136
93: ; preds = %66
%94 = sub nsw i32 %62, %63
%95 = add nsw i32 %94, -1
%96 = sitofp i32 %95 to double
%97 = fcmp olt double %8, %96
br i1 %97, label %76, label %98
98: ; preds = %93
%99 = sitofp i32 %94 to double
%100 = fcmp olt double %8, %99
br i1 %100, label %101, label %111
101: ; preds = %98
%102 = fsub double %8, %96
%103 = sext i32 %63 to i64
%104 = getelementptr inbounds double, double* %19, i64 %103
%105 = tail call i32 @doublemax(i32 %94, double* %104)
%106 = sitofp i32 %105 to double
%107 = fsub double 1.000000e+00, %102
%108 = load double, double* %19
%109 = fmul double %102, %108
%110 = tail call double @llvm.fmuladd.f64(double %106, double %107, double %109)
br label %136
111: ; preds = %98
%112 = add nsw i32 %62, -1
%113 = sitofp i32 %112 to double
%114 = fcmp ugt double %8, %113
br i1 %114, label %115, label %136
115: ; preds = %111
%116 = sitofp i32 %62 to double
%117 = fcmp olt double %8, %116
br i1 %117, label %118, label %128
118: ; preds = %115
%119 = fsub double %8, %113
%120 = load double, double* %19
%121 = fsub double 1.000000e+00, %119
%122 = sub nsw i32 %20, %62
%123 = sext i32 %62 to i64
%124 = getelementptr inbounds double, double* %19, i64 %123
%125 = tail call double @doublemin(i32 %122, double* %124)
%126 = fmul double %119, %125
%127 = tail call double @llvm.fmuladd.f64(double %120, double %121, double %126)
br label %136
128: ; preds = %115
%129 = sub nsw i32 %20, %62
%130 = sext i32 %62 to i64
%131 = getelementptr inbounds double, double* %19, i64 %130
%132 = fsub double %8, %116
%133 = add nsw i32 %129, -1
%134 = sitofp i32 %133 to double
%135 = fcmp ult double %132, %134
br i1 %135, label %7, label %12
136: ; preds = %111, %68, %118, %101, %83, %36, %33, %23, %12
%137 = phi double [ %17, %12 ], [ %24, %23 ], [ %35, %33 ], [ %38, %36 ], [ %92, %83 ], [ %110, %101 ], [ %127, %118 ], [ %27, %68 ], [ %27, %111 ]
ret double %137
}