#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);
void ERI_Interpolate_Linear_Complex(double out[2], 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) {
    out[0] = sy[0];
    out[1] = sy[1];
    return;
  }
  if (i == ns) {
    out[0] = sy[2 * ns - 2];
    out[1] = sy[2 * ns - 1];
    return;
  }
  x0 = sx[i - 1];
  x1 = sx[i];
  y0 = sy[2 * i - 2];
  y1 = sy[2 * i + 0];
  out[0] = y0 + (y1 - y0) * (x - x0) / (x1 - x0);
  y0 = sy[2 * i - 1];
  y1 = sy[2 * i + 1];
  out[1] = y0 + (y1 - y0) * (x - x0) / (x1 - x0);
}


x86 assembly (clang, -O0):
<ERI_Interpolate_Linear_Complex>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-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    1172 <ERI_Interpolate_Linear_Complex+0x62>
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    1182 <ERI_Interpolate_Linear_Complex+0x72>
jmp    1195 <ERI_Interpolate_Linear_Complex+0x85>
jmp    1187 <ERI_Interpolate_Linear_Complex+0x77>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1145 <ERI_Interpolate_Linear_Complex+0x35>
cmpl   $0x0,-0x28(%rbp)
jne    11c6 <ERI_Interpolate_Linear_Complex+0xb6>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x8(%rax)
jmp    1312 <ERI_Interpolate_Linear_Complex+0x202>
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jne    1212 <ERI_Interpolate_Linear_Complex+0x102>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x8(%rax)
jmp    1312 <ERI_Interpolate_Linear_Complex+0x202>
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
shl    $0x1,%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%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
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%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
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x8(%rax)
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ERI_Interpolate_Linear_Complex>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %ecx,%ebp
mov    %rdx,%r15
mov    %rsi,%rbx
movsd  %xmm0,(%rsp)
mov    %rdi,%r14
xor    %edi,%edi
test   %ecx,%ecx
setg   %dil
call   1030 <assert@plt>
movsd  (%rsp),%xmm0
test   %ebp,%ebp
jle    1184 <ERI_Interpolate_Linear_Complex+0x74>
mov    %ebp,%ecx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
ucomisd (%rbx,%rax,8),%xmm0
jbe    1162 <ERI_Interpolate_Linear_Complex+0x52>
add    $0x1,%rax
cmp    %rax,%rcx
jne    1150 <ERI_Interpolate_Linear_Complex+0x40>
mov    %ecx,%eax
test   %eax,%eax
je     1184 <ERI_Interpolate_Linear_Complex+0x74>
cmp    %ebp,%eax
jne    1196 <ERI_Interpolate_Linear_Complex+0x86>
add    %ebp,%ebp
movslq %ebp,%rax
movsd  -0x10(%r15,%rax,8),%xmm0
movsd  %xmm0,(%r14)
movsd  -0x8(%r15,%rax,8),%xmm0
jmp    11fb <ERI_Interpolate_Linear_Complex+0xeb>
movsd  (%r15),%xmm0
movsd  %xmm0,(%r14)
movsd  0x8(%r15),%xmm0
jmp    11fb <ERI_Interpolate_Linear_Complex+0xeb>
movslq %eax,%rcx
movapd %xmm0,%xmm4
movsd  -0x8(%rbx,%rcx,8),%xmm0
mov    %eax,%ecx
movsd  (%rbx,%rcx,8),%xmm1
lea    (%rax,%rax,1),%ecx
movslq %ecx,%rdx
movsd  -0x10(%r15,%rdx,8),%xmm2
movsd  (%r15,%rcx,8),%xmm3
subsd  %xmm2,%xmm3
subsd  %xmm0,%xmm4
mulsd  %xmm4,%xmm3
subsd  %xmm0,%xmm1
divsd  %xmm1,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%r14)
movsd  -0x8(%r15,%rdx,8),%xmm2
lea    0x1(%rax,%rax,1),%eax
movsd  (%r15,%rax,8),%xmm0
subsd  %xmm2,%xmm0
mulsd  %xmm4,%xmm0
divsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,0x8(%r14)
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @ERI_Interpolate_Linear_Complex(double* nocapture writeonly %0, double %1, double* nocapture %2, double* nocapture %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
%7 = zext i1 %6 to i32
%8 = tail call i32 @assert(i32 %7)
br i1 %6, label %9, label %24
9: ; preds = %5
%10 = zext i32 %4 to i64
br label %11
11: ; preds = %9, %16
%12 = phi i64 [ 0, %9 ], [ %17, %16 ]
%13 = getelementptr inbounds double, double* %2, i64 %12
%14 = load double, double* %13
%15 = fcmp olt double %14, %1
br i1 %15, label %16, label %19
16: ; preds = %11
%17 = add nuw nsw i64 %12, 1
%18 = icmp eq i64 %17, %10
br i1 %18, label %21, label %11
19: ; preds = %11
%20 = trunc i64 %12 to i32
br label %21
21: ; preds = %16, %19
%22 = phi i32 [ %20, %19 ], [ %4, %16 ]
%23 = icmp eq i32 %22, 0
br i1 %23, label %24, label %28
24: ; preds = %5, %21
%25 = load double, double* %3
store double %25, double* %0
%26 = getelementptr inbounds double, double* %3, i64 1
%27 = load double, double* %26
br label %74
28: ; preds = %21
%29 = icmp eq i32 %22, %4
br i1 %29, label %30, label %40
30: ; preds = %28
%31 = shl nsw i32 %4, 1
%32 = add nsw i32 %31, -2
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %3, i64 %33
%35 = load double, double* %34
store double %35, double* %0
%36 = add nsw i32 %31, -1
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %3, i64 %37
%39 = load double, double* %38
br label %74
40: ; preds = %28
%41 = add nsw i32 %22, -1
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %2, i64 %42
%44 = load double, double* %43
%45 = zext i32 %22 to i64
%46 = getelementptr inbounds double, double* %2, i64 %45
%47 = load double, double* %46
%48 = shl nuw nsw i32 %22, 1
%49 = add nsw i32 %48, -2
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %3, i64 %50
%52 = load double, double* %51
%53 = zext i32 %48 to i64
%54 = getelementptr inbounds double, double* %3, i64 %53
%55 = load double, double* %54
%56 = fsub double %55, %52
%57 = fsub double %1, %44
%58 = fmul double %57, %56
%59 = fsub double %47, %44
%60 = fdiv double %58, %59
%61 = fadd double %52, %60
store double %61, double* %0
%62 = add nsw i32 %48, -1
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %3, i64 %63
%65 = load double, double* %64
%66 = or i32 %48, 1
%67 = zext i32 %66 to i64
%68 = getelementptr inbounds double, double* %3, i64 %67
%69 = load double, double* %68
%70 = fsub double %69, %65
%71 = fmul double %57, %70
%72 = fdiv double %71, %59
%73 = fadd double %65, %72
br label %74
74: ; preds = %40, %30, %24
%75 = phi double [ %73, %40 ], [ %39, %30 ], [ %27, %24 ]
%76 = getelementptr inbounds double, double* %0, i64 1
store double %75, double* %76
ret void
}