#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*/ exit(int);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
int /*<<< orphan*/ stderr;
void interp(size_t inpLen, size_t outpLen, const double x[outpLen],
            const double xp[inpLen], const double fp[inpLen],
            double r[outpLen]) {
  int i, j;
  double a, b;
  j = -1;
  for (i = 0; i < outpLen; i++) {
    if (x[i] <= xp[0]) {
      r[i] = fp[0];
      continue;
    }
    if (x[i] >= xp[inpLen - 1]) {
      r[i] = fp[inpLen - 1];
      continue;
    }
    if (x[i] > xp[j + 1]) {
      while (x[i] > xp[j + 1]) {
        j++;
        if (j > inpLen) {
          fprintf(stderr, "Nonincreasing interpolation values! Quitting!\n");
          exit(1);
        }
      }
      a = (fp[j + 1] - fp[j]) / (xp[j + 1] - xp[j]);
      b = fp[j] - a * xp[j];
    }
    r[i] = a * x[i] + b;
  }
}


x86 assembly (clang, -O0):
<interp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0xffffffff,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
movslq -0x34(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    132a <interp+0x20a>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jb     1195 <interp+0x75>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    131c <interp+0x1fc>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rcx
sub    $0x1,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     11dc <interp+0xbc>
mov    -0x28(%rbp),%rax
mov    -0x8(%rbp),%rcx
sub    $0x1,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    131c <interp+0x1fc>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    12f0 <interp+0x1d0>
jmp    1206 <interp+0xe6>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1266 <interp+0x146>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
movslq -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jbe    1261 <interp+0x141>
mov    0x2d97(%rip),%rax
mov    (%rax),%edi
lea    0xdae(%rip),%rsi
call   1030 <fprintf@plt>
mov    $0x1,%edi
call   1040 <exit@plt>
jmp    1206 <interp+0xe6>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    114e <interp+0x2e>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<interp>:
push   %r15
push   %r14
push   %rbx
test   %rsi,%rsi
je     1202 <interp+0xe2>
mov    $0xffffffff,%r11d
xor    %r14d,%r14d
movabs $0x100000000,%r10
jmp    1168 <interp+0x48>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%r8),%xmm0
movsd  %xmm0,(%r9,%r14,8)
add    $0x1,%r14
cmp    %rsi,%r14
je     1202 <interp+0xe2>
movsd  (%rdx,%r14,8),%xmm0
movsd  (%rcx),%xmm3
ucomisd %xmm0,%xmm3
jae    1150 <interp+0x30>
ucomisd -0x8(%rcx,%rdi,8),%xmm0
jae    11c0 <interp+0xa0>
movslq %r11d,%r15
ucomisd 0x8(%rcx,%r15,8),%xmm0
jbe    11f5 <interp+0xd5>
mov    %r15,%rbx
shl    $0x20,%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    0x1(%r15),%rax
movsd  0x8(%rcx,%r15,8),%xmm3
ucomisd %xmm3,%xmm0
jbe    11c9 <interp+0xa9>
add    %r10,%rbx
mov    %rax,%r15
cmp    %rdi,%rax
jbe    11a0 <interp+0x80>
jmp    1208 <interp+0xe8>
xchg   %ax,%ax
movsd  -0x8(%r8,%rdi,8),%xmm0
jmp    1155 <interp+0x35>
movsd  (%r8,%rax,8),%xmm1
sar    $0x1d,%rbx
movsd  (%r8,%rbx,1),%xmm2
subsd  %xmm2,%xmm1
movsd  (%rcx,%rbx,1),%xmm4
subsd  %xmm4,%xmm3
divsd  %xmm3,%xmm1
mulsd  %xmm1,%xmm4
subsd  %xmm4,%xmm2
mov    %r15d,%r11d
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
jmp    1155 <interp+0x35>
pop    %rbx
pop    %r14
pop    %r15
ret
mov    0x2dd1(%rip),%rax
mov    (%rax),%edi
lea    0xde8(%rip),%rsi
call   1030 <fprintf@plt>
mov    $0x1,%edi
call   1040 <exit@plt>

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

LLVM IR (-O2):
@stderr = global i32 0
@.str = [47 x i8] c"Nonincreasing interpolation values! Quitting!\0A\00"
define void @interp(i64 %0, i64 %1, double* nocapture %2, double* nocapture %3, double* nocapture %4, double* nocapture writeonly %5) {
%7 = icmp eq i64 %1, 0
br i1 %7, label %76, label %8
8: ; preds = %6
%9 = add i64 %0, -1
%10 = getelementptr inbounds double, double* %3, i64 %9
%11 = getelementptr inbounds double, double* %4, i64 %9
br label %12
12: ; preds = %8, %68
%13 = phi i64 [ 0, %8 ], [ %74, %68 ]
%14 = phi double [ undef, %8 ], [ %72, %68 ]
%15 = phi double [ undef, %8 ], [ %71, %68 ]
%16 = phi i32 [ -1, %8 ], [ %70, %68 ]
%17 = getelementptr inbounds double, double* %2, i64 %13
%18 = load double, double* %17
%19 = load double, double* %3
%20 = fcmp ugt double %18, %19
br i1 %20, label %23, label %21
21: ; preds = %12
%22 = load double, double* %4
br label %68
23: ; preds = %12
%24 = load double, double* %10
%25 = fcmp ult double %18, %24
br i1 %25, label %28, label %26
26: ; preds = %23
%27 = load double, double* %11
br label %68
28: ; preds = %23
%29 = add nsw i32 %16, 1
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %3, i64 %30
%32 = load double, double* %31
%33 = fcmp ogt double %18, %32
br i1 %33, label %34, label %63
34: ; preds = %28
%35 = sext i32 %16 to i64
br label %36
36: ; preds = %34, %42
%37 = phi i64 [ %35, %34 ], [ %38, %42 ]
%38 = add i64 %37, 1
%39 = getelementptr inbounds double, double* %3, i64 %38
%40 = load double, double* %39
%41 = fcmp ogt double %18, %40
br i1 %41, label %42, label %48
42: ; preds = %36
%43 = icmp ugt i64 %38, %0
br i1 %43, label %44, label %36
44: ; preds = %42
%45 = load i32, i32* @stderr
%46 = tail call i32 @fprintf(i32 %45, i8* getelementptr inbounds ([47 x i8], [47 x i8]* @.str, i64 0, i64 0))
%47 = tail call i32 @exit(i32 1)
unreachable
48: ; preds = %36
%49 = trunc i64 %37 to i32
%50 = getelementptr inbounds double, double* %4, i64 %38
%51 = load double, double* %50
%52 = shl i64 %37, 32
%53 = ashr exact i64 %52, 32
%54 = getelementptr inbounds double, double* %4, i64 %53
%55 = load double, double* %54
%56 = fsub double %51, %55
%57 = getelementptr inbounds double, double* %3, i64 %53
%58 = load double, double* %57
%59 = fsub double %40, %58
%60 = fdiv double %56, %59
%61 = fneg double %60
%62 = tail call double @llvm.fmuladd.f64(double %61, double %58, double %55)
br label %63
63: ; preds = %48, %28
%64 = phi i32 [ %49, %48 ], [ %16, %28 ]
%65 = phi double [ %60, %48 ], [ %15, %28 ]
%66 = phi double [ %62, %48 ], [ %14, %28 ]
%67 = tail call double @llvm.fmuladd.f64(double %65, double %18, double %66)
br label %68
68: ; preds = %63, %26, %21
%69 = phi double [ %67, %63 ], [ %27, %26 ], [ %22, %21 ]
%70 = phi i32 [ %64, %63 ], [ %16, %26 ], [ %16, %21 ]
%71 = phi double [ %65, %63 ], [ %15, %26 ], [ %15, %21 ]
%72 = phi double [ %66, %63 ], [ %14, %26 ], [ %14, %21 ]
%73 = getelementptr inbounds double, double* %5, i64 %13
store double %69, double* %73
%74 = add nuw i64 %13, 1
%75 = icmp eq i64 %74, %1
br i1 %75, label %76, label %12
76: ; preds = %68, %6
ret void
}