double paw_Corrector_In(int i, double y[], double f[]) {
  double tmp;
  tmp = y[i] - 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>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x4(%rbp),%rcx
movsd  0xeda(%rip),%xmm2
mulsd  (%rax,%rcx,8),%xmm2
movsd  0xec5(%rip),%xmm3
mulsd  %xmm3,%xmm0
addsd  %xmm2,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  0xe9b(%rip),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  %xmm2,%xmm0
movsd  0xe6d(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

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

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

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