#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*/ printf(char *);
double splint_monotonic(int N, double *x, double *y, double *yprime,
                        double xt) {
  int klo, khi, k;
  double h, b, a, c, d, delx, s, retval, del;
  if (xt > x[N - 1]) {
    klo = N - 2;
    khi = N - 1;
  } else if (xt < x[0]) {
    klo = 0;
    khi = 1;
  } else {
    klo = 0;
    khi = N - 1;
    while (khi - klo > 1) {
      k = (khi + klo) >> 1;
      if (x[k] > xt)
        khi = k;
      else
        klo = k;
    }
  }
  h = x[khi] - x[klo];
  if (h == 0.0)
    printf("Bad XA input to routine SPLINT\n");
  s = (y[khi] - y[klo]) / h;
  a = (yprime[klo] + yprime[khi] - 2.0 * s) / h / h;
  b = (3.0 * s - 2.0 * yprime[klo] - yprime[khi]) / h;
  c = yprime[klo];
  d = y[klo];
  del = (xt - x[klo]);
  retval = a * del * del * del + b * del * del + c * del + d;
  return retval;
}


x86 assembly (clang, -O0):
<splint_monotonic>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1163 <splint_monotonic+0x53>
mov    -0x4(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11f0 <splint_monotonic+0xe0>
movsd  -0x28(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    118d <splint_monotonic+0x7d>
movl   $0x0,-0x2c(%rbp)
movl   $0x1,-0x30(%rbp)
jmp    11eb <splint_monotonic+0xdb>
movl   $0x0,-0x2c(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
sub    -0x2c(%rbp),%eax
cmp    $0x1,%eax
jle    11e6 <splint_monotonic+0xd6>
mov    -0x30(%rbp),%eax
add    -0x2c(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    11db <splint_monotonic+0xcb>
mov    -0x34(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    11e1 <splint_monotonic+0xd1>
mov    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    119d <splint_monotonic+0x8d>
jmp    11eb <splint_monotonic+0xdb>
jmp    11f0 <splint_monotonic+0xe0>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1233 <splint_monotonic+0x123>
jp     1233 <splint_monotonic+0x123>
lea    0xde2(%rip),%rdi
call   1030 <printf@plt>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  -0x70(%rbp),%xmm1
addsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
divsd  -0x40(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x70(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  0xd66(%rip),%xmm1
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xd3a(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x50(%rbp),%xmm1
mulsd  -0x80(%rbp),%xmm1
mulsd  -0x80(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<splint_monotonic>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %rcx,%r14
mov    %rdx,%r15
mov    %rsi,%r13
lea    -0x1(%rdi),%eax
movslq %edi,%rcx
ucomisd -0x8(%rsi,%rcx,8),%xmm0
jbe    113b <splint_monotonic+0x2b>
add    $0xfffffffe,%edi
mov    %edi,%ecx
jmp    1186 <splint_monotonic+0x76>
movsd  0x0(%r13),%xmm1
xor    %ecx,%ecx
ucomisd %xmm0,%xmm1
jbe    1150 <splint_monotonic+0x40>
mov    $0x1,%eax
jmp    1186 <splint_monotonic+0x76>
cmp    $0x3,%edi
jl     1186 <splint_monotonic+0x76>
xor    %edx,%edx
jmp    116b <splint_monotonic+0x5b>
nopl   0x0(%rax)
mov    %eax,%esi
sub    %ecx,%esi
mov    %ecx,%edx
cmp    $0x1,%esi
jle    1186 <splint_monotonic+0x76>
lea    (%rdx,%rax,1),%ecx
sar    %ecx
movslq %ecx,%rsi
movsd  0x0(%r13,%rsi,8),%xmm1
ucomisd %xmm0,%xmm1
jbe    1160 <splint_monotonic+0x50>
mov    %ecx,%eax
mov    %edx,%ecx
jmp    1160 <splint_monotonic+0x50>
movslq %eax,%r12
movsd  0x0(%r13,%r12,8),%xmm7
movslq %ecx,%rbx
movsd  0x0(%r13,%rbx,8),%xmm1
subsd  %xmm1,%xmm7
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm7
jne    11d3 <splint_monotonic+0xc3>
jp     11d3 <splint_monotonic+0xc3>
lea    0xe57(%rip),%rdi
movsd  %xmm0,0x8(%rsp)
movsd  %xmm7,(%rsp)
call   1030 <puts@plt>
movsd  (%rsp),%xmm7
movsd  0x8(%rsp),%xmm0
movsd  0x0(%r13,%rbx,8),%xmm1
movsd  (%r15,%r12,8),%xmm2
movsd  (%r15,%rbx,8),%xmm8
subsd  %xmm8,%xmm2
divsd  %xmm7,%xmm2
movsd  (%r14,%rbx,8),%xmm3
movsd  (%r14,%r12,8),%xmm5
movapd %xmm3,%xmm4
addsd  %xmm5,%xmm4
movapd %xmm2,%xmm6
addsd  %xmm2,%xmm6
subsd  %xmm6,%xmm4
divsd  %xmm7,%xmm4
divsd  %xmm7,%xmm4
mulsd  0xde8(%rip),%xmm2
movapd %xmm3,%xmm6
addsd  %xmm3,%xmm6
subsd  %xmm6,%xmm2
subsd  %xmm5,%xmm2
divsd  %xmm7,%xmm2
subsd  %xmm1,%xmm0
mulsd  %xmm0,%xmm4
mulsd  %xmm0,%xmm4
mulsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm4
addsd  %xmm2,%xmm4
mulsd  %xmm3,%xmm0
addsd  %xmm4,%xmm0
addsd  %xmm8,%xmm0
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
@.str = [32 x i8] c"Bad XA input to routine SPLINT\0A\00"
define double @splint_monotonic(i32 %0, double* %1, double* %2, double* %3, double %4) {
%6 = alloca i32
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca double
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
store i32 %0, i32* %6
store double* %1, double** %7
store double* %2, double** %8
store double* %3, double** %9
store double %4, double* %10
%23 = load double, double* %10
%24 = load double*, double** %7
%25 = load i32, i32* %6
%26 = sub nsw i32 %25, 1
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %24, i64 %27
%29 = load double, double* %28
%30 = fcmp ogt double %23, %29
br i1 %30, label %31, label %36
31: ; preds = %5
%32 = load i32, i32* %6
%33 = sub nsw i32 %32, 2
store i32 %33, i32* %11
%34 = load i32, i32* %6
%35 = sub nsw i32 %34, 1
store i32 %35, i32* %12
br label %70
36: ; preds = %5
%37 = load double, double* %10
%38 = load double*, double** %7
%39 = getelementptr inbounds double, double* %38, i64 0
%40 = load double, double* %39
%41 = fcmp olt double %37, %40
br i1 %41, label %42, label %43
42: ; preds = %36
store i32 0, i32* %11
store i32 1, i32* %12
br label %69
43: ; preds = %36
store i32 0, i32* %11
%44 = load i32, i32* %6
%45 = sub nsw i32 %44, 1
store i32 %45, i32* %12
br label %46
46: ; preds = %67, %43
%47 = load i32, i32* %12
%48 = load i32, i32* %11
%49 = sub nsw i32 %47, %48
%50 = icmp sgt i32 %49, 1
br i1 %50, label %51, label %68
51: ; preds = %46
%52 = load i32, i32* %12
%53 = load i32, i32* %11
%54 = add nsw i32 %52, %53
%55 = ashr i32 %54, 1
store i32 %55, i32* %13
%56 = load double*, double** %7
%57 = load i32, i32* %13
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %56, i64 %58
%60 = load double, double* %59
%61 = load double, double* %10
%62 = fcmp ogt double %60, %61
br i1 %62, label %63, label %65
63: ; preds = %51
%64 = load i32, i32* %13
store i32 %64, i32* %12
br label %67
65: ; preds = %51
%66 = load i32, i32* %13
store i32 %66, i32* %11
br label %67
67: ; preds = %65, %63
br label %46
68: ; preds = %46
br label %69
69: ; preds = %68, %42
br label %70
70: ; preds = %69, %31
%71 = load double*, double** %7
%72 = load i32, i32* %12
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %71, i64 %73
%75 = load double, double* %74
%76 = load double*, double** %7
%77 = load i32, i32* %11
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %76, i64 %78
%80 = load double, double* %79
%81 = fsub double %75, %80
store double %81, double* %14
%82 = load double, double* %14
%83 = fcmp oeq double %82, 0.000000e+00
br i1 %83, label %84, label %86
84: ; preds = %70
%85 = call i32 @printf(i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str, i64 0, i64 0))
br label %86
86: ; preds = %84, %70
%87 = load double*, double** %8
%88 = load i32, i32* %12
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
%92 = load double*, double** %8
%93 = load i32, i32* %11
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double, double* %92, i64 %94
%96 = load double, double* %95
%97 = fsub double %91, %96
%98 = load double, double* %14
%99 = fdiv double %97, %98
store double %99, double* %20
%100 = load double*, double** %9
%101 = load i32, i32* %11
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %100, i64 %102
%104 = load double, double* %103
%105 = load double*, double** %9
%106 = load i32, i32* %12
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds double, double* %105, i64 %107
%109 = load double, double* %108
%110 = fadd double %104, %109
%111 = load double, double* %20
%112 = call double @llvm.fmuladd.f64(double -2.000000e+00, double %111, double %110)
%113 = load double, double* %14
%114 = fdiv double %112, %113
%115 = load double, double* %14
%116 = fdiv double %114, %115
store double %116, double* %16
%117 = load double, double* %20
%118 = load double*, double** %9
%119 = load i32, i32* %11
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds double, double* %118, i64 %120
%122 = load double, double* %121
%123 = fmul double 2.000000e+00, %122
%124 = fneg double %123
%125 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %117, double %124)
%126 = load double*, double** %9
%127 = load i32, i32* %12
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds double, double* %126, i64 %128
%130 = load double, double* %129
%131 = fsub double %125, %130
%132 = load double, double* %14
%133 = fdiv double %131, %132
store double %133, double* %15
%134 = load double*, double** %9
%135 = load i32, i32* %11
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds double, double* %134, i64 %136
%138 = load double, double* %137
store double %138, double* %17
%139 = load double*, double** %8
%140 = load i32, i32* %11
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds double, double* %139, i64 %141
%143 = load double, double* %142
store double %143, double* %18
%144 = load double, double* %10
%145 = load double*, double** %7
%146 = load i32, i32* %11
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %145, i64 %147
%149 = load double, double* %148
%150 = fsub double %144, %149
store double %150, double* %22
%151 = load double, double* %16
%152 = load double, double* %22
%153 = fmul double %151, %152
%154 = load double, double* %22
%155 = fmul double %153, %154
%156 = load double, double* %22
%157 = load double, double* %15
%158 = load double, double* %22
%159 = fmul double %157, %158
%160 = load double, double* %22
%161 = fmul double %159, %160
%162 = call double @llvm.fmuladd.f64(double %155, double %156, double %161)
%163 = load double, double* %17
%164 = load double, double* %22
%165 = call double @llvm.fmuladd.f64(double %163, double %164, double %162)
%166 = load double, double* %18
%167 = fadd double %165, %166
store double %167, double* %21
%168 = load double, double* %21
ret double %168
}

LLVM IR (-O2):
@str = [31 x i8] c"Bad XA input to routine SPLINT\00"
define double @splint_monotonic(i32 %0, double* nocapture %1, double* nocapture %2, double* nocapture %3, double %4) {
%6 = add nsw i32 %0, -1
%7 = sext i32 %6 to i64
%8 = getelementptr inbounds double, double* %1, i64 %7
%9 = load double, double* %8
%10 = fcmp olt double %9, %4
br i1 %10, label %11, label %13
11: ; preds = %5
%12 = add nsw i32 %0, -2
br label %31
13: ; preds = %5
%14 = load double, double* %1
%15 = fcmp ogt double %14, %4
br i1 %15, label %31, label %16
16: ; preds = %13
%17 = icmp sgt i32 %0, 2
br i1 %17, label %18, label %31
18: ; preds = %16, %18
%19 = phi i32 [ %28, %18 ], [ 0, %16 ]
%20 = phi i32 [ %27, %18 ], [ %6, %16 ]
%21 = add nsw i32 %19, %20
%22 = ashr i32 %21, 1
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %1, i64 %23
%25 = load double, double* %24
%26 = fcmp ogt double %25, %4
%27 = select i1 %26, i32 %22, i32 %20
%28 = select i1 %26, i32 %19, i32 %22
%29 = sub nsw i32 %27, %28
%30 = icmp sgt i32 %29, 1
br i1 %30, label %18, label %31
31: ; preds = %18, %16, %13, %11
%32 = phi i32 [ %6, %11 ], [ 1, %13 ], [ %6, %16 ], [ %27, %18 ]
%33 = phi i32 [ %12, %11 ], [ 0, %13 ], [ 0, %16 ], [ %28, %18 ]
%34 = sext i32 %32 to i64
%35 = getelementptr inbounds double, double* %1, i64 %34
%36 = load double, double* %35
%37 = sext i32 %33 to i64
%38 = getelementptr inbounds double, double* %1, i64 %37
%39 = load double, double* %38
%40 = fsub double %36, %39
%41 = fcmp oeq double %40, 0.000000e+00
br i1 %41, label %42, label %45
42: ; preds = %31
%43 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([31 x i8], [31 x i8]* @str, i64 0, i64 0))
%44 = load double, double* %38
br label %45
45: ; preds = %42, %31
%46 = phi double [ %44, %42 ], [ %39, %31 ]
%47 = getelementptr inbounds double, double* %2, i64 %34
%48 = load double, double* %47
%49 = getelementptr inbounds double, double* %2, i64 %37
%50 = load double, double* %49
%51 = fsub double %48, %50
%52 = fdiv double %51, %40
%53 = getelementptr inbounds double, double* %3, i64 %37
%54 = load double, double* %53
%55 = getelementptr inbounds double, double* %3, i64 %34
%56 = load double, double* %55
%57 = fadd double %54, %56
%58 = tail call double @llvm.fmuladd.f64(double %52, double -2.000000e+00, double %57)
%59 = fdiv double %58, %40
%60 = fdiv double %59, %40
%61 = fmul double %54, -2.000000e+00
%62 = tail call double @llvm.fmuladd.f64(double %52, double 3.000000e+00, double %61)
%63 = fsub double %62, %56
%64 = fdiv double %63, %40
%65 = fsub double %4, %46
%66 = fmul double %65, %60
%67 = fmul double %65, %66
%68 = fmul double %65, %64
%69 = fmul double %65, %68
%70 = tail call double @llvm.fmuladd.f64(double %67, double %65, double %69)
%71 = tail call double @llvm.fmuladd.f64(double %54, double %65, double %70)
%72 = fadd double %50, %71
ret double %72
}