double paw_Corrector_In_F(int i, double f[]) {
  double tmp;
  tmp = -4.16666666667e-2 *
        (9.0 * f[i - 1] + 19.0 * f[i] - 5.0 * f[i + 1] + 1.0 * f[i + 2]);
  return tmp;
}


x86 assembly (clang, -O0):
<paw_Corrector_In_F>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  0xeeb(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  0xed6(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xeac(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  0xe7e(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<paw_Corrector_In_F>:
movslq %edi,%rax
movsd  (%rsi,%rax,8),%xmm0
mulsd  0xef0(%rip),%xmm0
movsd  -0x8(%rsi,%rax,8),%xmm1
mulsd  0xeea(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  0x8(%rsi,%rax,8),%xmm0
mulsd  0xee0(%rip),%xmm0
addsd  %xmm1,%xmm0
addsd  0x10(%rsi,%rax,8),%xmm0
mulsd  0xed6(%rip),%xmm0
ret

LLVM IR (-O0):
define double @paw_Corrector_In_F(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, 1
%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 = sext i32 %13 to i64
%15 = getelementptr inbounds double, double* %12, i64 %14
%16 = load double, double* %15
%17 = fmul double 1.900000e+01, %16
%18 = call double @llvm.fmuladd.f64(double 9.000000e+00, double %11, double %17)
%19 = load double*, double** %4
%20 = load i32, i32* %3
%21 = add nsw i32 %20, 1
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %19, i64 %22
%24 = load double, double* %23
%25 = call double @llvm.fmuladd.f64(double -5.000000e+00, double %24, double %18)
%26 = load double*, double** %4
%27 = load i32, i32* %3
%28 = add nsw i32 %27, 2
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %26, i64 %29
%31 = load double, double* %30
%32 = call double @llvm.fmuladd.f64(double 1.000000e+00, double %31, double %25)
%33 = fmul double 0xBFA5555555556819, %32
store double %33, double* %5
%34 = load double, double* %5
ret double %34
}

LLVM IR (-O2):
define double @paw_Corrector_In_F(i32 %0, double* nocapture %1) {
%3 = add nsw i32 %0, -1
%4 = sext i32 %3 to i64
%5 = getelementptr inbounds double, double* %1, i64 %4
%6 = load double, double* %5
%7 = sext i32 %0 to i64
%8 = getelementptr inbounds double, double* %1, i64 %7
%9 = load double, double* %8
%10 = fmul double %9, 1.900000e+01
%11 = tail call double @llvm.fmuladd.f64(double %6, double 9.000000e+00, double %10)
%12 = add nsw i32 %0, 1
%13 = sext i32 %12 to i64
%14 = getelementptr inbounds double, double* %1, i64 %13
%15 = load double, double* %14
%16 = tail call double @llvm.fmuladd.f64(double %15, double -5.000000e+00, double %11)
%17 = add nsw i32 %0, 2
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds double, double* %1, i64 %18
%20 = load double, double* %19
%21 = fadd double %20, %16
%22 = fmul double %21, 0xBFA5555555556819
ret double %22
}