#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 */
double pow(double, int);
void interp_quad(double t0, double t1, double y, double y1, double *dts,
                 double *values, int ndts) {
  double DT = t1 - t0;
  double a, b;
  a = (y1 * DT - y) / pow(DT, 2);
  b = 2 * y / DT - y1;
  for (int v = 1; v <= ndts; v++)
    values[v] = a * ((dts[v] - t0) * (dts[v] - t0)) + b * (dts[v] - t0);
  return;
}


x86 assembly (clang, -O0):
<interp_quad>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
mov    %rdi,-0x28(%rbp)
mov    %rsi,-0x30(%rbp)
mov    %edx,-0x34(%rbp)
movsd  -0x10(%rbp),%xmm0
subsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x18(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0xe5f(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x1,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     1233 <interp_quad+0x123>
movsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
subsd  -0x8(%rbp),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
subsd  -0x8(%rbp),%xmm1
mulsd  %xmm1,%xmm2
movsd  -0x50(%rbp),%xmm1
mov    -0x28(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm3
subsd  -0x8(%rbp),%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    11bc <interp_quad+0xac>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<interp_quad>:
push   %rbp
push   %r14
push   %rbx
sub    $0x30,%rsp
mov    %edx,%ebp
mov    %rsi,%r14
mov    %rdi,%rbx
movsd  %xmm3,0x8(%rsp)
movapd %xmm2,0x20(%rsp)
movapd %xmm0,0x10(%rsp)
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movapd %xmm1,%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movapd 0x10(%rsp),%xmm9
movapd 0x20(%rsp),%xmm1
test   %ebp,%ebp
jle    12e5 <interp_quad+0x1d5>
movsd  0x8(%rsp),%xmm2
movapd %xmm2,%xmm8
movsd  (%rsp),%xmm3
mulsd  %xmm3,%xmm8
subsd  %xmm1,%xmm8
divsd  %xmm0,%xmm8
addsd  %xmm1,%xmm1
divsd  %xmm3,%xmm1
subsd  %xmm2,%xmm1
mov    %ebp,%edx
lea    0x1(%rdx),%rax
mov    $0x1,%ecx
cmp    $0x4,%ebp
jb     124c <interp_quad+0x13c>
lea    0x8(%r14),%rsi
lea    (%rbx,%rax,8),%rdi
cmp    %rdi,%rsi
jae    11bb <interp_quad+0xab>
lea    (%r14,%rax,8),%rsi
lea    0x8(%rbx),%rdi
cmp    %rsi,%rdi
jb     124c <interp_quad+0x13c>
mov    %edx,%esi
and    $0xfffffffc,%esi
lea    0x1(%rsi),%rcx
movapd %xmm9,%xmm0
unpcklpd %xmm9,%xmm0
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
movapd %xmm8,%xmm3
unpcklpd %xmm8,%xmm3
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd 0x8(%rbx,%rdi,8),%xmm4
movupd 0x18(%rbx,%rdi,8),%xmm5
subpd  %xmm0,%xmm4
subpd  %xmm0,%xmm5
movapd %xmm2,%xmm6
mulpd  %xmm4,%xmm6
mulpd  %xmm4,%xmm4
movapd %xmm2,%xmm7
mulpd  %xmm5,%xmm7
mulpd  %xmm5,%xmm5
mulpd  %xmm3,%xmm4
addpd  %xmm6,%xmm4
mulpd  %xmm3,%xmm5
addpd  %xmm7,%xmm5
movupd %xmm4,0x8(%r14,%rdi,8)
movupd %xmm5,0x18(%r14,%rdi,8)
add    $0x4,%rdi
cmp    %rdi,%rsi
jne    11f0 <interp_quad+0xe0>
cmp    %rsi,%rdx
je     12e5 <interp_quad+0x1d5>
mov    %rcx,%rdx
not    %rdx
test   $0x1,%al
jne    127f <interp_quad+0x16f>
movsd  (%rbx,%rcx,8),%xmm0
subsd  %xmm9,%xmm0
movapd %xmm1,%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm0
mulsd  %xmm8,%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,(%r14,%rcx,8)
add    $0x1,%rcx
add    %rax,%rdx
je     12e5 <interp_quad+0x1d5>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rbx,%rcx,8),%xmm0
subsd  %xmm9,%xmm0
movapd %xmm1,%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm0
mulsd  %xmm8,%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,(%r14,%rcx,8)
movsd  0x8(%rbx,%rcx,8),%xmm0
subsd  %xmm9,%xmm0
movapd %xmm1,%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm0
mulsd  %xmm8,%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,0x8(%r14,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rax
jne    1290 <interp_quad+0x180>
add    $0x30,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define void @interp_quad(double %0, double %1, double %2, double %3, double* %4, double* %5, i32 %6) {
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca double*
%13 = alloca double*
%14 = alloca i32
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca i32
store double %0, double* %8
store double %1, double* %9
store double %2, double* %10
store double %3, double* %11
store double* %4, double** %12
store double* %5, double** %13
store i32 %6, i32* %14
%19 = load double, double* %9
%20 = load double, double* %8
%21 = fsub double %19, %20
store double %21, double* %15
%22 = load double, double* %11
%23 = load double, double* %15
%24 = load double, double* %10
%25 = fneg double %24
%26 = call double @llvm.fmuladd.f64(double %22, double %23, double %25)
%27 = load double, double* %15
%28 = call double @pow(double %27, i32 2)
%29 = fdiv double %26, %28
store double %29, double* %16
%30 = load double, double* %10
%31 = fmul double 2.000000e+00, %30
%32 = load double, double* %15
%33 = fdiv double %31, %32
%34 = load double, double* %11
%35 = fsub double %33, %34
store double %35, double* %17
store i32 1, i32* %18
br label %36
36: ; preds = %71, %7
%37 = load i32, i32* %18
%38 = load i32, i32* %14
%39 = icmp sle i32 %37, %38
br i1 %39, label %40, label %74
40: ; preds = %36
%41 = load double, double* %16
%42 = load double*, double** %12
%43 = load i32, i32* %18
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %42, i64 %44
%46 = load double, double* %45
%47 = load double, double* %8
%48 = fsub double %46, %47
%49 = load double*, double** %12
%50 = load i32, i32* %18
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
%53 = load double, double* %52
%54 = load double, double* %8
%55 = fsub double %53, %54
%56 = fmul double %48, %55
%57 = load double, double* %17
%58 = load double*, double** %12
%59 = load i32, i32* %18
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %58, i64 %60
%62 = load double, double* %61
%63 = load double, double* %8
%64 = fsub double %62, %63
%65 = fmul double %57, %64
%66 = call double @llvm.fmuladd.f64(double %41, double %56, double %65)
%67 = load double*, double** %13
%68 = load i32, i32* %18
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %67, i64 %69
store double %66, double* %70
br label %71
71: ; preds = %40
%72 = load i32, i32* %18
%73 = add nsw i32 %72, 1
store i32 %73, i32* %18
br label %36
74: ; preds = %36
ret void
}

LLVM IR (-O2):
define void @interp_quad(double %0, double %1, double %2, double %3, double* nocapture %4, double* nocapture writeonly %5, i32 %6) {
%8 = fsub double %1, %0
%9 = fneg double %2
%10 = tail call double @llvm.fmuladd.f64(double %3, double %8, double %9)
%11 = tail call double @pow(double %8, i32 2)
%12 = fdiv double %10, %11
%13 = fmul double %2, 2.000000e+00
%14 = fdiv double %13, %8
%15 = fsub double %14, %3
%16 = icmp slt i32 %6, 1
br i1 %16, label %88, label %17
17: ; preds = %7
%18 = add nuw i32 %6, 1
%19 = zext i32 %18 to i64
%20 = add nsw i64 %19, -1
%21 = icmp ult i64 %20, 4
br i1 %21, label %70, label %22
22: ; preds = %17
%23 = getelementptr double, double* %5, i64 1
%24 = getelementptr double, double* %5, i64 %19
%25 = getelementptr double, double* %4, i64 1
%26 = getelementptr double, double* %4, i64 %19
%27 = icmp ult double* %23, %26
%28 = icmp ult double* %25, %24
%29 = and i1 %27, %28
br i1 %29, label %70, label %30
30: ; preds = %22
%31 = and i64 %20, -4
%32 = or i64 %31, 1
%33 = insertelement <2 x double> poison, double %0, i64 0
%34 = shufflevector <2 x double> %33, <2 x double> poison, <2 x i32> zeroinitializer
%35 = insertelement <2 x double> poison, double %0, i64 0
%36 = shufflevector <2 x double> %35, <2 x double> poison, <2 x i32> zeroinitializer
%37 = insertelement <2 x double> poison, double %15, i64 0
%38 = shufflevector <2 x double> %37, <2 x double> poison, <2 x i32> zeroinitializer
%39 = insertelement <2 x double> poison, double %15, i64 0
%40 = shufflevector <2 x double> %39, <2 x double> poison, <2 x i32> zeroinitializer
%41 = insertelement <2 x double> poison, double %12, i64 0
%42 = shufflevector <2 x double> %41, <2 x double> poison, <2 x i32> zeroinitializer
%43 = insertelement <2 x double> poison, double %12, i64 0
%44 = shufflevector <2 x double> %43, <2 x double> poison, <2 x i32> zeroinitializer
br label %45
45: ; preds = %45, %30
%46 = phi i64 [ 0, %30 ], [ %66, %45 ]
%47 = or i64 %46, 1
%48 = getelementptr inbounds double, double* %4, i64 %47
%49 = bitcast double* %48 to <2 x double>*
%50 = load <2 x double>, <2 x double>* %49
%51 = getelementptr inbounds double, double* %48, i64 2
%52 = bitcast double* %51 to <2 x double>*
%53 = load <2 x double>, <2 x double>* %52
%54 = fsub <2 x double> %50, %34
%55 = fsub <2 x double> %53, %36
%56 = fmul <2 x double> %54, %54
%57 = fmul <2 x double> %55, %55
%58 = fmul <2 x double> %38, %54
%59 = fmul <2 x double> %40, %55
%60 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %42, <2 x double> %56, <2 x double> %58)
%61 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %44, <2 x double> %57, <2 x double> %59)
%62 = getelementptr inbounds double, double* %5, i64 %47
%63 = bitcast double* %62 to <2 x double>*
store <2 x double> %60, <2 x double>* %63
%64 = getelementptr inbounds double, double* %62, i64 2
%65 = bitcast double* %64 to <2 x double>*
store <2 x double> %61, <2 x double>* %65
%66 = add nuw i64 %46, 4
%67 = icmp eq i64 %66, %31
br i1 %67, label %68, label %45
68: ; preds = %45
%69 = icmp eq i64 %20, %31
br i1 %69, label %88, label %70
70: ; preds = %22, %17, %68
%71 = phi i64 [ 1, %22 ], [ 1, %17 ], [ %32, %68 ]
%72 = xor i64 %71, -1
%73 = and i64 %19, 1
%74 = icmp eq i64 %73, 0
br i1 %74, label %75, label %84
75: ; preds = %70
%76 = getelementptr inbounds double, double* %4, i64 %71
%77 = load double, double* %76
%78 = fsub double %77, %0
%79 = fmul double %78, %78
%80 = fmul double %15, %78
%81 = tail call double @llvm.fmuladd.f64(double %12, double %79, double %80)
%82 = getelementptr inbounds double, double* %5, i64 %71
store double %81, double* %82
%83 = add nuw nsw i64 %71, 1
br label %84
84: ; preds = %75, %70
%85 = phi i64 [ %71, %70 ], [ %83, %75 ]
%86 = sub nsw i64 0, %19
%87 = icmp eq i64 %72, %86
br i1 %87, label %88, label %89
88: ; preds = %84, %89, %68, %7
ret void
89: ; preds = %84, %89
%90 = phi i64 [ %106, %89 ], [ %85, %84 ]
%91 = getelementptr inbounds double, double* %4, i64 %90
%92 = load double, double* %91
%93 = fsub double %92, %0
%94 = fmul double %93, %93
%95 = fmul double %15, %93
%96 = tail call double @llvm.fmuladd.f64(double %12, double %94, double %95)
%97 = getelementptr inbounds double, double* %5, i64 %90
store double %96, double* %97
%98 = add nuw nsw i64 %90, 1
%99 = getelementptr inbounds double, double* %4, i64 %98
%100 = load double, double* %99
%101 = fsub double %100, %0
%102 = fmul double %101, %101
%103 = fmul double %15, %101
%104 = tail call double @llvm.fmuladd.f64(double %12, double %102, double %103)
%105 = getelementptr inbounds double, double* %5, i64 %98
store double %104, double* %105
%106 = add nuw nsw i64 %90, 2
%107 = icmp eq i64 %106, %19
br i1 %107, label %88, label %89
}