#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*/ fprintf(int /*<<< orphan*/, char *);
int /*<<< orphan*/ stderr;
float interpfloat(float *arrx, float *arry, int npts, float x) {
  int n;
  int nLo = 0, nHi = npts - 1;
  if (arrx[1] > arrx[0]) {
    for (n = 1; n < npts - 1; n++)
      if (arrx[n + 1] <= arrx[n]) {
        fprintf(stderr,
                "Error: interpfloats() x is not monotonically increasesing");
        return 0;
      }
  } else if (arrx[1] < arrx[0]) {
    for (n = 1; n < npts - 1; n++)
      if (arrx[n + 1] >= arrx[n]) {
        fprintf(stderr,
                "Error: interpfloats() x is not monotonically decreasing");
        return 0;
      }
  } else {
    fprintf(stderr, "Error: interpfloats() x is not monotonically decreasing "
                    "or increasesing");
    return 0;
  }
  do {
    n = (int)(.5 * (nLo + nHi));
    if (x > arrx[n])
      nLo = n;
    else
      nHi = n;
  } while ((nHi - nLo) > 1);
  return ((arrx[nHi] - x) * arry[nLo] + (x - arrx[nLo]) * arry[nHi]) /
         (arrx[nHi] - arrx[nLo]);
}


x86 assembly (clang, -O0):
<interpfloat>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
ucomiss (%rax),%xmm0
jbe    11c8 <interpfloat+0xb8>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11c3 <interpfloat+0xb3>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss %xmm1,%xmm0
jb     11b0 <interpfloat+0xa0>
mov    0x2e4b(%rip),%rax
mov    (%rax),%edi
lea    0xe6a(%rip),%rsi
call   1030 <fprintf@plt>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1346 <interpfloat+0x236>
jmp    11b5 <interpfloat+0xa5>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1155 <interpfloat+0x45>
jmp    127f <interpfloat+0x16f>
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm1
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
ucomiss %xmm1,%xmm0
jbe    1258 <interpfloat+0x148>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1253 <interpfloat+0x143>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
jb     1240 <interpfloat+0x130>
mov    0x2dbb(%rip),%rax
mov    (%rax),%edi
lea    0xe14(%rip),%rsi
call   1030 <fprintf@plt>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1346 <interpfloat+0x236>
jmp    1245 <interpfloat+0x135>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11e9 <interpfloat+0xd9>
jmp    127a <interpfloat+0x16a>
mov    0x2d81(%rip),%rax
mov    (%rax),%edi
lea    0xe12(%rip),%rsi
call   1030 <fprintf@plt>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1346 <interpfloat+0x236>
jmp    127f <interpfloat+0x16f>
jmp    1284 <interpfloat+0x174>
mov    -0x28(%rbp),%eax
add    -0x2c(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xd6a(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x24(%rbp)
movss  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
jbe    12c3 <interpfloat+0x1b3>
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
jmp    12c9 <interpfloat+0x1b9>
mov    -0x24(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    12ce <interpfloat+0x1be>
mov    -0x2c(%rbp),%eax
sub    -0x28(%rbp),%eax
cmp    $0x1,%eax
jg     1284 <interpfloat+0x174>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x20(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x20(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<interpfloat>:
push   %rax
add    $0xffffffff,%edx
movss  (%rdi),%xmm2
movss  0x4(%rdi),%xmm1
ucomiss %xmm2,%xmm1
jbe    116c <interpfloat+0x5c>
cmp    $0x2,%edx
mov    $0x1,%eax
cmovge %edx,%eax
add    $0xffffffffffffffff,%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
cmp    %rcx,%rax
je     11bc <interpfloat+0xac>
movss  0x8(%rdi,%rcx,4),%xmm2
add    $0x1,%rcx
ucomiss %xmm2,%xmm1
movaps %xmm2,%xmm1
jb     1140 <interpfloat+0x30>
mov    0x2e82(%rip),%rax
mov    (%rax),%edi
lea    0xea1(%rip),%rsi
jmp    123c <interpfloat+0x12c>
ucomiss %xmm1,%xmm2
jbe    122c <interpfloat+0x11c>
cmp    $0x2,%edx
mov    $0x1,%eax
cmovge %edx,%eax
add    $0xffffffffffffffff,%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
cmp    %rcx,%rax
je     11bc <interpfloat+0xac>
movss  0x8(%rdi,%rcx,4),%xmm2
add    $0x1,%rcx
ucomiss %xmm1,%xmm2
movaps %xmm2,%xmm1
jb     1190 <interpfloat+0x80>
mov    0x2e32(%rip),%rax
mov    (%rax),%edi
lea    0xe8b(%rip),%rsi
jmp    123c <interpfloat+0x12c>
xor    %eax,%eax
movsd  0xe3a(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
mov    %edx,%ecx
lea    (%rcx,%rax,1),%edx
xorps  %xmm2,%xmm2
cvtsi2sd %edx,%xmm2
mulsd  %xmm1,%xmm2
cvttsd2si %xmm2,%edx
movslq %edx,%rdx
ucomiss (%rdi,%rdx,4),%xmm0
cmova  %edx,%eax
cmova  %ecx,%edx
mov    %edx,%ecx
sub    %eax,%ecx
cmp    $0x1,%ecx
jg     11d0 <interpfloat+0xc0>
movslq %edx,%rcx
movss  (%rdi,%rcx,4),%xmm2
movaps %xmm2,%xmm1
subss  %xmm0,%xmm1
cltq
movss  (%rdi,%rax,4),%xmm3
subss  %xmm3,%xmm0
mulss  (%rsi,%rcx,4),%xmm0
mulss  (%rsi,%rax,4),%xmm1
addss  %xmm0,%xmm1
subss  %xmm3,%xmm2
divss  %xmm2,%xmm1
jmp    1244 <interpfloat+0x134>
mov    0x2dad(%rip),%rax
mov    (%rax),%edi
lea    0xe3e(%rip),%rsi
call   1030 <fprintf@plt>
xorps  %xmm1,%xmm1
movaps %xmm1,%xmm0
pop    %rax
ret

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

LLVM IR (-O2):
@stderr = global i32 0
@.str = [58 x i8] c"Error: interpfloats() x is not monotonically increasesing\00"
@.str.1 = [56 x i8] c"Error: interpfloats() x is not monotonically decreasing\00"
@.str.2 = [72 x i8] c"Error: interpfloats() x is not monotonically decreasing or increasesing\00"
define float @interpfloat(float* nocapture %0, float* nocapture %1, i32 %2, float %3) {
%5 = add i32 %2, -1
%6 = getelementptr inbounds float, float* %0, i64 1
%7 = load float, float* %6
%8 = load float, float* %0
%9 = fcmp ogt float %7, %8
br i1 %9, label %10, label %26
10: ; preds = %4
%11 = call i32 @llvm.smax.i32(i32 %5, i32 1)
%12 = zext i32 %11 to i64
br label %13
13: ; preds = %10, %18
%14 = phi float [ %7, %10 ], [ %21, %18 ]
%15 = phi i64 [ 1, %10 ], [ %19, %18 ]
%16 = icmp eq i64 %15, %12
br i1 %16, label %17, label %18
17: ; preds = %31, %13
br label %46
18: ; preds = %13
%19 = add nuw nsw i64 %15, 1
%20 = getelementptr inbounds float, float* %0, i64 %19
%21 = load float, float* %20
%22 = fcmp ugt float %21, %14
br i1 %22, label %13, label %23
23: ; preds = %18
%24 = load i32, i32* @stderr
%25 = tail call i32 @fprintf(i32 %24, i8* getelementptr inbounds ([58 x i8], [58 x i8]* @.str, i64 0, i64 0))
br label %78
26: ; preds = %4
%27 = fcmp olt float %7, %8
br i1 %27, label %28, label %43
28: ; preds = %26
%29 = call i32 @llvm.smax.i32(i32 %5, i32 1)
%30 = zext i32 %29 to i64
br label %31
31: ; preds = %28, %35
%32 = phi float [ %7, %28 ], [ %38, %35 ]
%33 = phi i64 [ 1, %28 ], [ %36, %35 ]
%34 = icmp eq i64 %33, %30
br i1 %34, label %17, label %35
35: ; preds = %31
%36 = add nuw nsw i64 %33, 1
%37 = getelementptr inbounds float, float* %0, i64 %36
%38 = load float, float* %37
%39 = fcmp ult float %38, %32
br i1 %39, label %31, label %40
40: ; preds = %35
%41 = load i32, i32* @stderr
%42 = tail call i32 @fprintf(i32 %41, i8* getelementptr inbounds ([56 x i8], [56 x i8]* @.str.1, i64 0, i64 0))
br label %78
43: ; preds = %26
%44 = load i32, i32* @stderr
%45 = tail call i32 @fprintf(i32 %44, i8* getelementptr inbounds ([72 x i8], [72 x i8]* @.str.2, i64 0, i64 0))
br label %78
46: ; preds = %17, %46
%47 = phi i32 [ %57, %46 ], [ 0, %17 ]
%48 = phi i32 [ %58, %46 ], [ %5, %17 ]
%49 = add nsw i32 %48, %47
%50 = sitofp i32 %49 to double
%51 = fmul double %50, 5.000000e-01
%52 = fptosi double %51 to i32
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds float, float* %0, i64 %53
%55 = load float, float* %54
%56 = fcmp olt float %55, %3
%57 = select i1 %56, i32 %52, i32 %47
%58 = select i1 %56, i32 %48, i32 %52
%59 = sub nsw i32 %58, %57
%60 = icmp sgt i32 %59, 1
br i1 %60, label %46, label %61
61: ; preds = %46
%62 = sext i32 %58 to i64
%63 = getelementptr inbounds float, float* %0, i64 %62
%64 = load float, float* %63
%65 = fsub float %64, %3
%66 = sext i32 %57 to i64
%67 = getelementptr inbounds float, float* %1, i64 %66
%68 = load float, float* %67
%69 = getelementptr inbounds float, float* %0, i64 %66
%70 = load float, float* %69
%71 = fsub float %3, %70
%72 = getelementptr inbounds float, float* %1, i64 %62
%73 = load float, float* %72
%74 = fmul float %71, %73
%75 = tail call float @llvm.fmuladd.f32(float %65, float %68, float %74)
%76 = fsub float %64, %70
%77 = fdiv float %75, %76
br label %78
78: ; preds = %61, %43, %40, %23
%79 = phi float [ 0.000000e+00, %23 ], [ %77, %61 ], [ 0.000000e+00, %40 ], [ 0.000000e+00, %43 ]
ret float %79
}