double paw_Derivative5_4(int i, double f[]) {
  double tmp;
  tmp = 4.16666666667e-2 * (-2.0 * f[i - 3] + 12.0 * f[i - 2] -
                            36.0 * f[i - 1] + 20.0 * f[i] + 6.0 * f[i + 1]);
  return tmp;
}


x86 assembly (clang, -O0):
<paw_Derivative5_4>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  0xeee(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
addsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xebf(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe9a(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xe70(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe58(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<paw_Derivative5_4>:
movslq %edi,%rax
movsd  -0x18(%rsi,%rax,8),%xmm0
movsd  -0x10(%rsi,%rax,8),%xmm1
mulsd  0xee9(%rip),%xmm1
addsd  %xmm0,%xmm0
movsd  -0x8(%rsi,%rax,8),%xmm2
mulsd  0xedf(%rip),%xmm2
subsd  %xmm0,%xmm1
movsd  (%rsi,%rax,8),%xmm3
mulsd  0xed6(%rip),%xmm3
addsd  %xmm1,%xmm2
movsd  0x8(%rsi,%rax,8),%xmm0
mulsd  0xecc(%rip),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
mulsd  0xec4(%rip),%xmm0
ret

LLVM IR (-O0):
define double @paw_Derivative5_4(i32 %0, double* %1) {
%3 = alloca i32
%4 = alloca double*
%5 = alloca double
store i32 %0, i32* %3
store double* %1, double** %4
%6 = load double*, double** %4
%7 = load i32, i32* %3
%8 = sub nsw i32 %7, 3
%9 = sext i32 %8 to i64
%10 = getelementptr inbounds double, double* %6, i64 %9
%11 = load double, double* %10
%12 = load double*, double** %4
%13 = load i32, i32* %3
%14 = sub nsw i32 %13, 2
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %12, i64 %15
%17 = load double, double* %16
%18 = fmul double 1.200000e+01, %17
%19 = call double @llvm.fmuladd.f64(double -2.000000e+00, double %11, double %18)
%20 = load double*, double** %4
%21 = load i32, i32* %3
%22 = sub nsw i32 %21, 1
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %20, i64 %23
%25 = load double, double* %24
%26 = call double @llvm.fmuladd.f64(double -3.600000e+01, double %25, double %19)
%27 = load double*, double** %4
%28 = load i32, i32* %3
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %27, i64 %29
%31 = load double, double* %30
%32 = call double @llvm.fmuladd.f64(double 2.000000e+01, double %31, double %26)
%33 = load double*, double** %4
%34 = load i32, i32* %3
%35 = add nsw i32 %34, 1
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %33, i64 %36
%38 = load double, double* %37
%39 = call double @llvm.fmuladd.f64(double 6.000000e+00, double %38, double %32)
%40 = fmul double 0x3FA5555555556819, %39
store double %40, double* %5
%41 = load double, double* %5
ret double %41
}

LLVM IR (-O2):
define double @paw_Derivative5_4(i32 %0, double* nocapture %1) {
%3 = add nsw i32 %0, -3
%4 = sext i32 %3 to i64
%5 = getelementptr inbounds double, double* %1, i64 %4
%6 = load double, double* %5
%7 = add nsw i32 %0, -2
%8 = sext i32 %7 to i64
%9 = getelementptr inbounds double, double* %1, i64 %8
%10 = load double, double* %9
%11 = fmul double %10, 1.200000e+01
%12 = tail call double @llvm.fmuladd.f64(double %6, double -2.000000e+00, double %11)
%13 = add nsw i32 %0, -1
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds double, double* %1, i64 %14
%16 = load double, double* %15
%17 = tail call double @llvm.fmuladd.f64(double %16, double -3.600000e+01, double %12)
%18 = sext i32 %0 to i64
%19 = getelementptr inbounds double, double* %1, i64 %18
%20 = load double, double* %19
%21 = tail call double @llvm.fmuladd.f64(double %20, double 2.000000e+01, double %17)
%22 = add nsw i32 %0, 1
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %1, i64 %23
%25 = load double, double* %24
%26 = tail call double @llvm.fmuladd.f64(double %25, double 6.000000e+00, double %21)
%27 = fmul double %26, 0x3FA5555555556819
ret double %27
}