#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 */
float modff(float, float *);
int cubic_Intrpl(float *ip, int ipsz, float *op, int opsz) {
  int j;
  float map;
  int i;
  float iI, iF;
  op[0] = ip[0];
  for (j = 1; j < opsz - 1; j++) {
    map = ((float)j) * (ipsz - 1) / (opsz - 1);
    iF = modff(map, &iI);
    i = (int)iI;
    if (i == 0 || i == ipsz - 2)
      op[j] = (1 - iF) * ip[i] + iF * ip[i + 1];
    else
      op[j] = 0.5 * (ip[i - 1] * (-iF * iF * iF + 2 * iF * iF - iF) +
                     ip[i] * (3 * iF * iF * iF - 5 * iF * iF + 2) +
                     ip[i + 1] * (-3 * iF * iF * iF + 4 * iF * iF + iF) +
                     ip[i + 2] * (iF * iF * iF - iF * iF));
  }
  op[opsz - 1] = ip[ipsz - 1];
  return 0;
}


x86 assembly (clang, -O0):
<cubic_Intrpl>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,(%rax)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1371 <cubic_Intrpl+0x261>
cvtsi2ssl -0x20(%rbp),%xmm0
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
lea    -0x2c(%rbp),%rdi
call   1030 <modff@plt>
movss  %xmm0,-0x30(%rbp)
cvttss2si -0x2c(%rbp),%eax
mov    %eax,-0x28(%rbp)
cmpl   $0x0,-0x28(%rbp)
je     11aa <cubic_Intrpl+0x9a>
mov    -0x28(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jne    11f5 <cubic_Intrpl+0xe5>
movss  0xe6a(%rip),%xmm0
subss  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x30(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    135e <cubic_Intrpl+0x24e>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x30(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm2
mulss  -0x30(%rbp),%xmm2
movss  -0x30(%rbp),%xmm3
movss  0xde5(%rip),%xmm1
mulss  -0x30(%rbp),%xmm1
mulss  -0x30(%rbp),%xmm1
mulss  %xmm3,%xmm2
addss  %xmm1,%xmm2
subss  -0x30(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0xdc1(%rip),%xmm3
mulss  -0x30(%rbp),%xmm3
mulss  -0x30(%rbp),%xmm3
movss  -0x30(%rbp),%xmm5
movss  0xda6(%rip),%xmm4
mulss  -0x30(%rbp),%xmm4
mulss  -0x30(%rbp),%xmm4
movd   %xmm4,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm4
mulss  %xmm5,%xmm3
addss  %xmm4,%xmm3
movss  0xd7b(%rip),%xmm4
addss  %xmm4,%xmm3
mulss  %xmm3,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0xd4d(%rip),%xmm2
mulss  -0x30(%rbp),%xmm2
mulss  -0x30(%rbp),%xmm2
movss  -0x30(%rbp),%xmm4
movss  0xd32(%rip),%xmm3
mulss  -0x30(%rbp),%xmm3
mulss  -0x30(%rbp),%xmm3
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
addss  -0x30(%rbp),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x30(%rbp),%xmm2
mulss  -0x30(%rbp),%xmm2
movss  -0x30(%rbp),%xmm4
movss  -0x30(%rbp),%xmm3
mulss  -0x30(%rbp),%xmm3
movd   %xmm3,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm3
mulss  %xmm4,%xmm2
addss  %xmm3,%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xcb7(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1363 <cubic_Intrpl+0x253>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    113d <cubic_Intrpl+0x2d>
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
xor    %eax,%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cubic_Intrpl>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rdx,%r14
mov    %esi,%r15d
mov    %rdi,%rbx
movss  (%rdi),%xmm0
movss  %xmm0,(%rdx)
lea    -0x1(%rcx),%eax
mov    %eax,0x8(%rsp)
lea    -0x1(%r15),%eax
mov    %eax,0x4(%rsp)
cmp    $0x3,%ecx
jl     128c <cubic_Intrpl+0x17c>
xorps  %xmm0,%xmm0
cvtsi2ssl 0x4(%rsp),%xmm0
movss  %xmm0,0x10(%rsp)
mov    0x8(%rsp),%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  %xmm0,0xc(%rsp)
add    $0xfffffffe,%r15d
mov    %eax,%r13d
mov    $0x1,%ebp
lea    0x14(%rsp),%r12
jmp    11b2 <cubic_Intrpl+0xa2>
nopw   0x0(%rax,%rax,1)
movss  0xe90(%rip),%xmm1
subss  %xmm0,%xmm1
movss  0x4(%rbx,%rax,4),%xmm2
mulss  %xmm0,%xmm2
mulss  (%rbx,%rax,4),%xmm1
addss  %xmm2,%xmm1
movss  %xmm1,(%r14,%rbp,4)
add    $0x1,%rbp
cmp    %rbp,%r13
je     128c <cubic_Intrpl+0x17c>
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  0x10(%rsp),%xmm0
divss  0xc(%rsp),%xmm0
mov    %r12,%rdi
call   1030 <modff@plt>
cvttss2si 0x14(%rsp),%ecx
movslq %ecx,%rax
test   %ecx,%ecx
je     1180 <cubic_Intrpl+0x70>
cmp    %ecx,%r15d
je     1180 <cubic_Intrpl+0x70>
movaps %xmm0,%xmm2
addss  %xmm0,%xmm2
mulss  %xmm0,%xmm2
movaps %xmm0,%xmm3
mulss  %xmm0,%xmm3
movaps %xmm3,%xmm1
mulss  %xmm0,%xmm1
subss  %xmm1,%xmm2
subss  %xmm0,%xmm2
movaps %xmm0,%xmm4
mulss  0xdf5(%rip),%xmm4
mulss  %xmm0,%xmm4
movaps %xmm0,%xmm5
mulss  0xdea(%rip),%xmm5
mulss  %xmm0,%xmm5
mulss  %xmm0,%xmm4
addss  %xmm5,%xmm4
addss  0xdda(%rip),%xmm4
mulss  (%rbx,%rax,4),%xmm4
mulss  -0x4(%rbx,%rax,4),%xmm2
addss  %xmm4,%xmm2
movaps %xmm0,%xmm4
mulss  0xdc4(%rip),%xmm4
mulss  %xmm0,%xmm4
movaps %xmm0,%xmm5
mulss  0xdb9(%rip),%xmm5
mulss  %xmm0,%xmm5
mulss  %xmm0,%xmm4
addss  %xmm5,%xmm4
addss  %xmm0,%xmm4
mulss  0x4(%rbx,%rax,4),%xmm4
addss  %xmm2,%xmm4
subss  %xmm3,%xmm1
mulss  0x8(%rbx,%rax,4),%xmm1
addss  %xmm4,%xmm1
mulss  0xd8d(%rip),%xmm1
jmp    119f <cubic_Intrpl+0x8f>
movslq 0x4(%rsp),%rax
movss  (%rbx,%rax,4),%xmm0
movslq 0x8(%rsp),%rax
movss  %xmm0,(%r14,%rax,4)
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @cubic_Intrpl(float* nocapture %0, i32 %1, float* nocapture writeonly %2, i32 %3) {
%5 = alloca float
%6 = bitcast float* %5 to i8*
%7 = load float, float* %0
store float %7, float* %2
%8 = add i32 %3, -1
%9 = icmp sgt i32 %3, 2
%10 = add nsw i32 %1, -1
br i1 %9, label %11, label %86
11: ; preds = %4
%12 = sitofp i32 %10 to float
%13 = sitofp i32 %8 to float
%14 = add nsw i32 %1, -2
%15 = zext i32 %8 to i64
br label %16
16: ; preds = %11, %81
%17 = phi i64 [ 1, %11 ], [ %84, %81 ]
%18 = trunc i64 %17 to i32
%19 = sitofp i32 %18 to float
%20 = fmul float %12, %19
%21 = fdiv float %20, %13
%22 = call float @modff(float %21, float* %5)
%23 = load float, float* %5
%24 = fptosi float %23 to i32
%25 = icmp eq i32 %24, 0
%26 = icmp eq i32 %14, %24
%27 = select i1 %25, i1 true, i1 %26
br i1 %27, label %28, label %39
28: ; preds = %16
%29 = fsub float 1.000000e+00, %22
%30 = sext i32 %24 to i64
%31 = getelementptr inbounds float, float* %0, i64 %30
%32 = load float, float* %31
%33 = add nsw i32 %24, 1
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds float, float* %0, i64 %34
%36 = load float, float* %35
%37 = fmul float %22, %36
%38 = tail call float @llvm.fmuladd.f32(float %29, float %32, float %37)
br label %81
39: ; preds = %16
%40 = add nsw i32 %24, -1
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float, float* %0, i64 %41
%43 = load float, float* %42
%44 = fneg float %22
%45 = fmul float %22, %44
%46 = fmul float %22, 2.000000e+00
%47 = fmul float %22, %46
%48 = tail call float @llvm.fmuladd.f32(float %45, float %22, float %47)
%49 = fsub float %48, %22
%50 = sext i32 %24 to i64
%51 = getelementptr inbounds float, float* %0, i64 %50
%52 = load float, float* %51
%53 = fmul float %22, 3.000000e+00
%54 = fmul float %22, %53
%55 = fmul float %22, -5.000000e+00
%56 = fmul float %22, %55
%57 = tail call float @llvm.fmuladd.f32(float %54, float %22, float %56)
%58 = fadd float %57, 2.000000e+00
%59 = fmul float %58, %52
%60 = tail call float @llvm.fmuladd.f32(float %43, float %49, float %59)
%61 = add nsw i32 %24, 1
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %0, i64 %62
%64 = load float, float* %63
%65 = fmul float %22, -3.000000e+00
%66 = fmul float %22, %65
%67 = fmul float %22, 4.000000e+00
%68 = fmul float %22, %67
%69 = tail call float @llvm.fmuladd.f32(float %66, float %22, float %68)
%70 = fadd float %22, %69
%71 = tail call float @llvm.fmuladd.f32(float %64, float %70, float %60)
%72 = add nsw i32 %24, 2
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds float, float* %0, i64 %73
%75 = load float, float* %74
%76 = fmul float %22, %22
%77 = fneg float %76
%78 = tail call float @llvm.fmuladd.f32(float %76, float %22, float %77)
%79 = tail call float @llvm.fmuladd.f32(float %75, float %78, float %71)
%80 = fmul float %79, 5.000000e-01
br label %81
81: ; preds = %28, %39
%82 = phi float [ %38, %28 ], [ %80, %39 ]
%83 = getelementptr inbounds float, float* %2, i64 %17
store float %82, float* %83
%84 = add nuw nsw i64 %17, 1
%85 = icmp eq i64 %84, %15
br i1 %85, label %86, label %16
86: ; preds = %81, %4
%87 = sext i32 %10 to i64
%88 = getelementptr inbounds float, float* %0, i64 %87
%89 = load float, float* %88
%90 = sext i32 %8 to i64
%91 = getelementptr inbounds float, float* %2, i64 %90
store float %89, float* %91
ret i32 0
}