#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*/ assert(int);
double ERI_Interpolate_Linear_Real(double x, const double *sx, const double *sy,
                                   int ns) {
  int i;
  double x0, x1, y0, y1;
  assert(0 < ns);
  for (i = 0; i < ns && sx[i] < x; i++)
    ;
  if (i == 0) {
    return sy[0];
  }
  if (i == ns) {
    return sy[ns - 1];
  }
  x0 = sx[i - 1];
  x1 = sx[i];
  y0 = sy[i - 1];
  y1 = sy[i];
  return y0 + (y1 - y0) * (x - x0) / (x1 - x0);
}


x86 assembly (clang, -O0):
<ERI_Interpolate_Linear_Real>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movsd  %xmm0,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %edx,-0x24(%rbp)
xor    %eax,%eax
cmp    -0x24(%rbp),%eax
setl   %al
and    $0x1,%al
movzbl %al,%edi
call   1030 <assert@plt>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%ecx
xor    %eax,%eax
cmp    -0x24(%rbp),%ecx
mov    %al,-0x49(%rbp)
jge    116e <ERI_Interpolate_Linear_Real+0x5e>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x49(%rbp)
mov    -0x49(%rbp),%al
test   $0x1,%al
jne    117e <ERI_Interpolate_Linear_Real+0x6e>
jmp    1191 <ERI_Interpolate_Linear_Real+0x81>
jmp    1183 <ERI_Interpolate_Linear_Real+0x73>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1141 <ERI_Interpolate_Linear_Real+0x31>
cmpl   $0x0,-0x28(%rbp)
jne    11ad <ERI_Interpolate_Linear_Real+0x9d>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    125b <ERI_Interpolate_Linear_Real+0x14b>
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jne    11d5 <ERI_Interpolate_Linear_Real+0xc5>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    125b <ERI_Interpolate_Linear_Real+0x14b>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
subsd  -0x40(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm2
subsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm1
movsd  -0x38(%rbp),%xmm2
subsd  -0x30(%rbp),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ERI_Interpolate_Linear_Real>:
push   %rbp
push   %r14
push   %rbx
sub    $0x10,%rsp
mov    %edx,%ebp
mov    %rsi,%r14
mov    %rdi,%rbx
movsd  %xmm0,0x8(%rsp)
xor    %edi,%edi
test   %edx,%edx
setg   %dil
call   1030 <assert@plt>
test   %ebp,%ebp
jle    1176 <ERI_Interpolate_Linear_Real+0x66>
mov    %ebp,%ecx
xor    %eax,%eax
movsd  0x8(%rsp),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
ucomisd (%rbx,%rax,8),%xmm0
jbe    1162 <ERI_Interpolate_Linear_Real+0x52>
add    $0x1,%rax
cmp    %rax,%rcx
jne    1150 <ERI_Interpolate_Linear_Real+0x40>
mov    %ecx,%eax
test   %eax,%eax
je     1176 <ERI_Interpolate_Linear_Real+0x66>
cmp    %ebp,%eax
jne    117d <ERI_Interpolate_Linear_Real+0x6d>
movslq %ebp,%rax
movsd  -0x8(%r14,%rax,8),%xmm0
jmp    11b8 <ERI_Interpolate_Linear_Real+0xa8>
movsd  (%r14),%xmm0
jmp    11b8 <ERI_Interpolate_Linear_Real+0xa8>
movslq %eax,%rcx
movsd  -0x8(%rbx,%rcx,8),%xmm4
mov    %eax,%eax
movsd  (%rbx,%rax,8),%xmm1
movsd  -0x8(%r14,%rcx,8),%xmm2
movsd  (%r14,%rax,8),%xmm3
subsd  %xmm2,%xmm3
movsd  0x8(%rsp),%xmm0
subsd  %xmm4,%xmm0
mulsd  %xmm3,%xmm0
subsd  %xmm4,%xmm1
divsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
add    $0x10,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define double @ERI_Interpolate_Linear_Real(double %0, double* %1, double* %2, i32 %3) {
%5 = alloca double
%6 = alloca double
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
store double %0, double* %6
store double* %1, double** %7
store double* %2, double** %8
store i32 %3, i32* %9
%15 = load i32, i32* %9
%16 = icmp slt i32 0, %15
%17 = zext i1 %16 to i32
%18 = call i32 @assert(i32 %17)
store i32 0, i32* %10
br label %19
19: ; preds = %34, %4
%20 = load i32, i32* %10
%21 = load i32, i32* %9
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %31
23: ; preds = %19
%24 = load double*, double** %7
%25 = load i32, i32* %10
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
%28 = load double, double* %27
%29 = load double, double* %6
%30 = fcmp olt double %28, %29
br label %31
31: ; preds = %23, %19
%32 = phi i1 [ false, %19 ], [ %30, %23 ]
br i1 %32, label %33, label %37
33: ; preds = %31
br label %34
34: ; preds = %33
%35 = load i32, i32* %10
%36 = add nsw i32 %35, 1
store i32 %36, i32* %10
br label %19
37: ; preds = %31
%38 = load i32, i32* %10
%39 = icmp eq i32 %38, 0
br i1 %39, label %40, label %44
40: ; preds = %37
%41 = load double*, double** %8
%42 = getelementptr inbounds double, double* %41, i64 0
%43 = load double, double* %42
store double %43, double* %5
br label %91
44: ; preds = %37
%45 = load i32, i32* %10
%46 = load i32, i32* %9
%47 = icmp eq i32 %45, %46
br i1 %47, label %48, label %55
48: ; preds = %44
%49 = load double*, double** %8
%50 = load i32, i32* %9
%51 = sub nsw i32 %50, 1
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %49, i64 %52
%54 = load double, double* %53
store double %54, double* %5
br label %91
55: ; preds = %44
%56 = load double*, double** %7
%57 = load i32, i32* %10
%58 = sub nsw i32 %57, 1
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %56, i64 %59
%61 = load double, double* %60
store double %61, double* %11
%62 = load double*, double** %7
%63 = load i32, i32* %10
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double, double* %62, i64 %64
%66 = load double, double* %65
store double %66, double* %12
%67 = load double*, double** %8
%68 = load i32, i32* %10
%69 = sub nsw i32 %68, 1
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double, double* %67, i64 %70
%72 = load double, double* %71
store double %72, double* %13
%73 = load double*, double** %8
%74 = load i32, i32* %10
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %73, i64 %75
%77 = load double, double* %76
store double %77, double* %14
%78 = load double, double* %13
%79 = load double, double* %14
%80 = load double, double* %13
%81 = fsub double %79, %80
%82 = load double, double* %6
%83 = load double, double* %11
%84 = fsub double %82, %83
%85 = fmul double %81, %84
%86 = load double, double* %12
%87 = load double, double* %11
%88 = fsub double %86, %87
%89 = fdiv double %85, %88
%90 = fadd double %78, %89
store double %90, double* %5
br label %91
91: ; preds = %55, %48, %40
%92 = load double, double* %5
ret double %92
}

LLVM IR (-O2):
define double @ERI_Interpolate_Linear_Real(double %0, double* nocapture %1, double* nocapture %2, i32 %3) {
%5 = icmp sgt i32 %3, 0
%6 = zext i1 %5 to i32
%7 = tail call i32 @assert(i32 %6)
br i1 %5, label %8, label %23
8: ; preds = %4
%9 = zext i32 %3 to i64
br label %10
10: ; preds = %8, %15
%11 = phi i64 [ 0, %8 ], [ %16, %15 ]
%12 = getelementptr inbounds double, double* %1, i64 %11
%13 = load double, double* %12
%14 = fcmp olt double %13, %0
br i1 %14, label %15, label %18
15: ; preds = %10
%16 = add nuw nsw i64 %11, 1
%17 = icmp eq i64 %16, %9
br i1 %17, label %20, label %10
18: ; preds = %10
%19 = trunc i64 %11 to i32
br label %20
20: ; preds = %15, %18
%21 = phi i32 [ %19, %18 ], [ %3, %15 ]
%22 = icmp eq i32 %21, 0
br i1 %22, label %23, label %25
23: ; preds = %4, %20
%24 = load double, double* %2
br label %50
25: ; preds = %20
%26 = icmp eq i32 %21, %3
br i1 %26, label %27, label %32
27: ; preds = %25
%28 = add nsw i32 %3, -1
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %2, i64 %29
%31 = load double, double* %30
br label %50
32: ; preds = %25
%33 = add nsw i32 %21, -1
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %1, i64 %34
%36 = load double, double* %35
%37 = zext i32 %21 to i64
%38 = getelementptr inbounds double, double* %1, i64 %37
%39 = load double, double* %38
%40 = getelementptr inbounds double, double* %2, i64 %34
%41 = load double, double* %40
%42 = getelementptr inbounds double, double* %2, i64 %37
%43 = load double, double* %42
%44 = fsub double %43, %41
%45 = fsub double %0, %36
%46 = fmul double %45, %44
%47 = fsub double %39, %36
%48 = fdiv double %46, %47
%49 = fadd double %41, %48
br label %50
50: ; preds = %32, %27, %23
%51 = phi double [ %24, %23 ], [ %31, %27 ], [ %49, %32 ]
ret double %51
}