void filter_r3(double *x, double *y, double *t, int len) {
  double init_val, pow_val;
  int i, lookahead;
  lookahead = (len < (8) ? len : (8));
  init_val = x[0];
  pow_val = -(0.3333333333333333333333333333333);
  for (i = 1; i <= lookahead; i++) {
    init_val += pow_val * x[i - 1];
    pow_val *= -(0.3333333333333333333333333333333);
  }
  y[0] = init_val;
  for (i = 1; i < len; i++)
    y[i] = x[i - 1] - (0.3333333333333333333333333333333) * y[i - 1];
  init_val = x[len - 1];
  pow_val = -(0.3333333333333333333333333333333);
  for (i = 1; i <= lookahead; i++) {
    init_val += pow_val * x[len - i];
    pow_val *= -(0.3333333333333333333333333333333);
  }
  t[len - 1] = init_val;
  for (i = len - 2; i >= 0; i--)
    t[i] = x[i] - (0.3333333333333333333333333333333) * t[i + 1];
  for (i = 0; i < len - 1; i++)
    y[i] =
        (-8.0 * t[i] - 8.0 / 9.0 * y[i] + x[i + 1] + 35.0 / 3.0 * x[i]) / 6.0;
  y[len - 1] = (-8.0 * t[len - 1] - 8.0 / 9.0 * y[len - 1] + x[len - 1] +
                35.0 / 3.0 * x[len - 1]) /
               6.0;
}


x86 assembly (clang, -O0):
<filter_r3>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
cmpl   $0x8,-0x1c(%rbp)
jge    1128 <filter_r3+0x28>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
jmp    1135 <filter_r3+0x35>
mov    $0x8,%eax
mov    %eax,-0x3c(%rbp)
jmp    1135 <filter_r3+0x35>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xeb0(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     11b1 <filter_r3+0xb1>
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xe67(%rip),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    115c <filter_r3+0x5c>
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1220 <filter_r3+0x120>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe03(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11c5 <filter_r3+0xc5>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xdc1(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     12a0 <filter_r3+0x1a0>
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xd78(%rip),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    124b <filter_r3+0x14b>
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jl     1314 <filter_r3+0x214>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xd0f(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    12c0 <filter_r3+0x1c0>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13cb <filter_r3+0x2cb>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  0xcd7(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xcab(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xc74(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xc5c(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    131b <filter_r3+0x21b>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  0xc2e(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xc02(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xbc6(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xbae(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<filter_r3>:
push   %r14
push   %rbx
cmp    $0x8,%ecx
mov    $0x8,%r9d
cmovl  %ecx,%r9d
movsd  (%rdi),%xmm1
test   %r9d,%r9d
jle    116b <filter_r3+0x6b>
lea    0x1(%r9),%eax
movsd  0xedb(%rip),%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
cmp    $0x2,%eax
je     121e <filter_r3+0x11e>
lea    -0x2(%rax),%r10
add    $0xfffffffffffffffd,%rax
mov    %r10d,%r8d
and    $0x3,%r8d
cmp    $0x3,%rax
jae    1174 <filter_r3+0x74>
mov    $0x1,%ebx
movsd  0xea8(%rip),%xmm1
mov    $0x2,%eax
test   %r8,%r8
jne    11f2 <filter_r3+0xf2>
jmp    121e <filter_r3+0x11e>
movapd %xmm1,%xmm0
jmp    121e <filter_r3+0x11e>
and    $0xfffffffffffffffc,%r10
neg    %r10
movsd  0xe7d(%rip),%xmm2
mov    $0x2,%eax
movapd %xmm2,%xmm1
nopl   0x0(%rax)
mulsd  %xmm2,%xmm1
movsd  -0x8(%rdi,%rax,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm0,%xmm3
mulsd  %xmm2,%xmm1
movsd  (%rdi,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm1
movsd  0x8(%rdi,%rax,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm0,%xmm3
mulsd  %xmm2,%xmm1
movsd  0x10(%rdi,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
lea    (%r10,%rax,1),%rbx
add    $0x4,%rbx
add    $0x4,%rax
cmp    $0x2,%rbx
jne    1190 <filter_r3+0x90>
lea    -0x1(%rax),%rbx
test   %r8,%r8
je     121e <filter_r3+0x11e>
movsd  0xe06(%rip),%xmm2
nopw   0x0(%rax,%rax,1)
mulsd  %xmm2,%xmm1
movsd  (%rdi,%rbx,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm3,%xmm0
mov    %rax,%rbx
add    $0x1,%rax
add    $0xffffffffffffffff,%r8
jne    1200 <filter_r3+0x100>
movsd  %xmm0,(%rsi)
cmp    $0x2,%ecx
jl     1299 <filter_r3+0x199>
mov    %ecx,%r8d
add    $0xffffffffffffffff,%r8
mov    $0x1,%ebx
cmp    $0x2,%ecx
je     1280 <filter_r3+0x180>
mov    %r8,%r10
and    $0xfffffffffffffffe,%r10
neg    %r10
mov    $0x1,%ebx
movsd  0xdb1(%rip),%xmm1
nop
mulsd  %xmm1,%xmm0
addsd  -0x8(%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rsi,%rbx,8)
mulsd  %xmm1,%xmm0
addsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,0x8(%rsi,%rbx,8)
lea    (%r10,%rbx,1),%rax
add    $0x2,%rax
add    $0x2,%rbx
cmp    $0x1,%rax
jne    1250 <filter_r3+0x150>
test   $0x1,%r8b
je     1299 <filter_r3+0x199>
mulsd  0xd72(%rip),%xmm0
addsd  -0x8(%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rsi,%rbx,8)
movslq %ecx,%r10
lea    -0x1(%r10),%r8
movsd  -0x8(%rdi,%r10,8),%xmm0
test   %r9d,%r9d
jle    139a <filter_r3+0x29a>
mov    %r9d,%r11d
lea    -0x1(%r11),%rax
mov    %r11d,%r9d
and    $0x3,%r9d
cmp    $0x3,%rax
jae    12df <filter_r3+0x1df>
movsd  0xd34(%rip),%xmm1
mov    $0x1,%eax
test   %r9,%r9
jne    1365 <filter_r3+0x265>
jmp    139a <filter_r3+0x29a>
and    $0xfffffffc,%r11d
lea    (%rdi,%r10,8),%r14
add    $0xfffffffffffffff8,%r14
movsd  0xd0d(%rip),%xmm2
mov    $0xffffffffffffffff,%rax
movapd %xmm2,%xmm1
xchg   %ax,%ax
movsd  0x8(%r14,%rax,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm0,%xmm3
mulsd  %xmm2,%xmm1
movsd  (%r14,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm1
movsd  -0x8(%r14,%rax,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm0,%xmm3
mulsd  %xmm2,%xmm1
movsd  -0x10(%r14,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm1
lea    (%r11,%rax,1),%rbx
add    $0xfffffffffffffffc,%rbx
add    $0xfffffffffffffffc,%rax
cmp    $0xffffffffffffffff,%rbx
jne    1300 <filter_r3+0x200>
neg    %rax
test   %r9,%r9
je     139a <filter_r3+0x29a>
sub    %rax,%r10
lea    (%rdi,%r10,8),%rax
neg    %r9
xor    %ebx,%ebx
movsd  0xc87(%rip),%xmm2
nopl   0x0(%rax)
movsd  (%rax,%rbx,8),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm1
add    $0xffffffffffffffff,%rbx
cmp    %rbx,%r9
jne    1380 <filter_r3+0x280>
movsd  %xmm0,(%rdx,%r8,8)
cmp    $0x2,%ecx
jl     153e <filter_r3+0x43e>
lea    -0x2(%rcx),%ebx
mov    %rbx,%rax
test   $0x1,%bl
jne    13d0 <filter_r3+0x2d0>
movsd  0x8(%rdx,%rbx,8),%xmm0
mulsd  0xc3e(%rip),%xmm0
addsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rdx,%rbx,8)
lea    -0x1(%rbx),%rax
test   %rbx,%rbx
je     140e <filter_r3+0x30e>
movsd  0xc23(%rip),%xmm0
nopl   (%rax)
movsd  0x8(%rdx,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdi,%rax,8),%xmm1
movsd  %xmm1,(%rdx,%rax,8)
mulsd  %xmm0,%xmm1
addsd  -0x8(%rdi,%rax,8),%xmm1
movsd  %xmm1,-0x8(%rdx,%rax,8)
add    $0xfffffffffffffffe,%rax
cmp    $0xffffffffffffffff,%rax
jne    13e0 <filter_r3+0x2e0>
cmp    $0x2,%ecx
jl     153e <filter_r3+0x43e>
mov    %r8d,%r14d
cmp    $0x2,%r8d
jae    1427 <filter_r3+0x327>
xor    %eax,%eax
jmp    14d4 <filter_r3+0x3d4>
lea    (%rsi,%r14,8),%rcx
lea    (%rdx,%r14,8),%rbx
cmp    %rsi,%rbx
seta   %r10b
lea    (%rdi,%r14,8),%rbx
add    $0x8,%rbx
cmp    %rdx,%rcx
seta   %r11b
cmp    %rsi,%rbx
seta   %bl
cmp    %rdi,%rcx
seta   %r9b
xor    %eax,%eax
test   %r11b,%r10b
jne    14d4 <filter_r3+0x3d4>
and    %r9b,%bl
jne    14d4 <filter_r3+0x3d4>
mov    %r14d,%eax
and    $0xfffffffe,%eax
xor    %ebx,%ebx
movapd 0xbc2(%rip),%xmm0
movapd 0xbca(%rip),%xmm1
movapd 0xbd2(%rip),%xmm2
movapd 0xbda(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
movupd (%rdx,%rbx,8),%xmm4
movupd (%rsi,%rbx,8),%xmm5
mulpd  %xmm0,%xmm5
mulpd  %xmm1,%xmm4
addpd  %xmm5,%xmm4
movupd (%rdi,%rbx,8),%xmm5
movupd 0x8(%rdi,%rbx,8),%xmm6
addpd  %xmm4,%xmm6
mulpd  %xmm2,%xmm5
addpd  %xmm6,%xmm5
divpd  %xmm3,%xmm5
movupd %xmm5,(%rsi,%rbx,8)
add    $0x2,%rbx
cmp    %rbx,%rax
jne    1490 <filter_r3+0x390>
cmp    %r14,%rax
je     153e <filter_r3+0x43e>
movsd  0xb2c(%rip),%xmm0
movsd  0xb2c(%rip),%xmm1
movsd  0xb2c(%rip),%xmm2
movsd  0xb2c(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%rax,8),%xmm4
mulsd  %xmm0,%xmm4
movsd  (%rdx,%rax,8),%xmm5
mulsd  %xmm1,%xmm5
addsd  %xmm4,%xmm5
addsd  0x8(%rdi,%rax,8),%xmm5
movsd  (%rdi,%rax,8),%xmm4
mulsd  %xmm2,%xmm4
addsd  %xmm5,%xmm4
divsd  %xmm3,%xmm4
movsd  %xmm4,(%rsi,%rax,8)
lea    0x1(%rax),%rcx
mov    %rcx,%rax
cmp    %rcx,%r14
jne    1500 <filter_r3+0x400>
movsd  (%rsi,%r8,8),%xmm0
mulsd  0xabc(%rip),%xmm0
movsd  (%rdx,%r8,8),%xmm1
mulsd  0xab6(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  (%rdi,%r8,8),%xmm0
addsd  %xmm0,%xmm1
mulsd  0xaa8(%rip),%xmm0
addsd  %xmm1,%xmm0
divsd  0xaa4(%rip),%xmm0
movsd  %xmm0,(%rsi,%r8,8)
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @filter_r3(double* %0, double* %1, double* %2, i32 %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca double
%10 = alloca double
%11 = alloca i32
%12 = alloca i32
store double* %0, double** %5
store double* %1, double** %6
store double* %2, double** %7
store i32 %3, i32* %8
%13 = load i32, i32* %8
%14 = icmp slt i32 %13, 8
br i1 %14, label %15, label %17
15: ; preds = %4
%16 = load i32, i32* %8
br label %18
17: ; preds = %4
br label %18
18: ; preds = %17, %15
%19 = phi i32 [ %16, %15 ], [ 8, %17 ]
store i32 %19, i32* %12
%20 = load double*, double** %5
%21 = getelementptr inbounds double, double* %20, i64 0
%22 = load double, double* %21
store double %22, double* %9
store double 0xBFD5555555555555, double* %10
store i32 1, i32* %11
br label %23
23: ; preds = %39, %18
%24 = load i32, i32* %11
%25 = load i32, i32* %12
%26 = icmp sle i32 %24, %25
br i1 %26, label %27, label %42
27: ; preds = %23
%28 = load double, double* %10
%29 = load double*, double** %5
%30 = load i32, i32* %11
%31 = sub nsw i32 %30, 1
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %29, i64 %32
%34 = load double, double* %33
%35 = load double, double* %9
%36 = call double @llvm.fmuladd.f64(double %28, double %34, double %35)
store double %36, double* %9
%37 = load double, double* %10
%38 = fmul double %37, 0xBFD5555555555555
store double %38, double* %10
br label %39
39: ; preds = %27
%40 = load i32, i32* %11
%41 = add nsw i32 %40, 1
store i32 %41, i32* %11
br label %23
42: ; preds = %23
%43 = load double, double* %9
%44 = load double*, double** %6
%45 = getelementptr inbounds double, double* %44, i64 0
store double %43, double* %45
store i32 1, i32* %11
br label %46
46: ; preds = %68, %42
%47 = load i32, i32* %11
%48 = load i32, i32* %8
%49 = icmp slt i32 %47, %48
br i1 %49, label %50, label %71
50: ; preds = %46
%51 = load double*, double** %5
%52 = load i32, i32* %11
%53 = sub nsw i32 %52, 1
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %51, i64 %54
%56 = load double, double* %55
%57 = load double*, double** %6
%58 = load i32, i32* %11
%59 = sub nsw i32 %58, 1
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %57, i64 %60
%62 = load double, double* %61
%63 = call double @llvm.fmuladd.f64(double 0xBFD5555555555555, double %62, double %56)
%64 = load double*, double** %6
%65 = load i32, i32* %11
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %64, i64 %66
store double %63, double* %67
br label %68
68: ; preds = %50
%69 = load i32, i32* %11
%70 = add nsw i32 %69, 1
store i32 %70, i32* %11
br label %46
71: ; preds = %46
%72 = load double*, double** %5
%73 = load i32, i32* %8
%74 = sub nsw i32 %73, 1
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %72, i64 %75
%77 = load double, double* %76
store double %77, double* %9
store double 0xBFD5555555555555, double* %10
store i32 1, i32* %11
br label %78
78: ; preds = %95, %71
%79 = load i32, i32* %11
%80 = load i32, i32* %12
%81 = icmp sle i32 %79, %80
br i1 %81, label %82, label %98
82: ; preds = %78
%83 = load double, double* %10
%84 = load double*, double** %5
%85 = load i32, i32* %8
%86 = load i32, i32* %11
%87 = sub nsw i32 %85, %86
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %84, i64 %88
%90 = load double, double* %89
%91 = load double, double* %9
%92 = call double @llvm.fmuladd.f64(double %83, double %90, double %91)
store double %92, double* %9
%93 = load double, double* %10
%94 = fmul double %93, 0xBFD5555555555555
store double %94, double* %10
br label %95
95: ; preds = %82
%96 = load i32, i32* %11
%97 = add nsw i32 %96, 1
store i32 %97, i32* %11
br label %78
98: ; preds = %78
%99 = load double, double* %9
%100 = load double*, double** %7
%101 = load i32, i32* %8
%102 = sub nsw i32 %101, 1
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %100, i64 %103
store double %99, double* %104
%105 = load i32, i32* %8
%106 = sub nsw i32 %105, 2
store i32 %106, i32* %11
br label %107
107: ; preds = %127, %98
%108 = load i32, i32* %11
%109 = icmp sge i32 %108, 0
br i1 %109, label %110, label %130
110: ; preds = %107
%111 = load double*, double** %5
%112 = load i32, i32* %11
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %111, i64 %113
%115 = load double, double* %114
%116 = load double*, double** %7
%117 = load i32, i32* %11
%118 = add nsw i32 %117, 1
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds double, double* %116, i64 %119
%121 = load double, double* %120
%122 = call double @llvm.fmuladd.f64(double 0xBFD5555555555555, double %121, double %115)
%123 = load double*, double** %7
%124 = load i32, i32* %11
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds double, double* %123, i64 %125
store double %122, double* %126
br label %127
127: ; preds = %110
%128 = load i32, i32* %11
%129 = add nsw i32 %128, -1
store i32 %129, i32* %11
br label %107
130: ; preds = %107
store i32 0, i32* %11
br label %131
131: ; preds = %168, %130
%132 = load i32, i32* %11
%133 = load i32, i32* %8
%134 = sub nsw i32 %133, 1
%135 = icmp slt i32 %132, %134
br i1 %135, label %136, label %171
136: ; preds = %131
%137 = load double*, double** %7
%138 = load i32, i32* %11
%139 = sext i32 %138 to i64
%140 = getelementptr inbounds double, double* %137, i64 %139
%141 = load double, double* %140
%142 = load double*, double** %6
%143 = load i32, i32* %11
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds double, double* %142, i64 %144
%146 = load double, double* %145
%147 = fmul double 0x3FEC71C71C71C71C, %146
%148 = fneg double %147
%149 = call double @llvm.fmuladd.f64(double -8.000000e+00, double %141, double %148)
%150 = load double*, double** %5
%151 = load i32, i32* %11
%152 = add nsw i32 %151, 1
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds double, double* %150, i64 %153
%155 = load double, double* %154
%156 = fadd double %149, %155
%157 = load double*, double** %5
%158 = load i32, i32* %11
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds double, double* %157, i64 %159
%161 = load double, double* %160
%162 = call double @llvm.fmuladd.f64(double 0x4027555555555555, double %161, double %156)
%163 = fdiv double %162, 6.000000e+00
%164 = load double*, double** %6
%165 = load i32, i32* %11
%166 = sext i32 %165 to i64
%167 = getelementptr inbounds double, double* %164, i64 %166
store double %163, double* %167
br label %168
168: ; preds = %136
%169 = load i32, i32* %11
%170 = add nsw i32 %169, 1
store i32 %170, i32* %11
br label %131
171: ; preds = %131
%172 = load double*, double** %7
%173 = load i32, i32* %8
%174 = sub nsw i32 %173, 1
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds double, double* %172, i64 %175
%177 = load double, double* %176
%178 = load double*, double** %6
%179 = load i32, i32* %8
%180 = sub nsw i32 %179, 1
%181 = sext i32 %180 to i64
%182 = getelementptr inbounds double, double* %178, i64 %181
%183 = load double, double* %182
%184 = fmul double 0x3FEC71C71C71C71C, %183
%185 = fneg double %184
%186 = call double @llvm.fmuladd.f64(double -8.000000e+00, double %177, double %185)
%187 = load double*, double** %5
%188 = load i32, i32* %8
%189 = sub nsw i32 %188, 1
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %187, i64 %190
%192 = load double, double* %191
%193 = fadd double %186, %192
%194 = load double*, double** %5
%195 = load i32, i32* %8
%196 = sub nsw i32 %195, 1
%197 = sext i32 %196 to i64
%198 = getelementptr inbounds double, double* %194, i64 %197
%199 = load double, double* %198
%200 = call double @llvm.fmuladd.f64(double 0x4027555555555555, double %199, double %193)
%201 = fdiv double %200, 6.000000e+00
%202 = load double*, double** %6
%203 = load i32, i32* %8
%204 = sub nsw i32 %203, 1
%205 = sext i32 %204 to i64
%206 = getelementptr inbounds double, double* %202, i64 %205
store double %201, double* %206
ret void
}

LLVM IR (-O2):
define void @filter_r3(double* nocapture %0, double* nocapture %1, double* nocapture %2, i32 %3) {
%5 = icmp slt i32 %3, 8
%6 = select i1 %5, i32 %3, i32 8
%7 = load double, double* %0
%8 = icmp slt i32 %6, 1
br i1 %8, label %69, label %9
9: ; preds = %4
%10 = add nuw nsw i32 %6, 1
%11 = zext i32 %10 to i64
%12 = tail call double @llvm.fmuladd.f64(double %7, double 0xBFD5555555555555, double %7)
%13 = icmp eq i32 %10, 2
br i1 %13, label %69, label %14
14: ; preds = %9
%15 = add nsw i64 %11, -2
%16 = add nsw i64 %11, -3
%17 = and i64 %15, 3
%18 = icmp ult i64 %16, 3
br i1 %18, label %49, label %19
19: ; preds = %14
%20 = and i64 %15, -4
br label %21
21: ; preds = %21, %19
%22 = phi i64 [ 2, %19 ], [ %46, %21 ]
%23 = phi double [ %12, %19 ], [ %45, %21 ]
%24 = phi double [ 0xBFD5555555555555, %19 ], [ %42, %21 ]
%25 = phi i64 [ 1, %19 ], [ %41, %21 ]
%26 = phi i64 [ 0, %19 ], [ %47, %21 ]
%27 = fmul double %24, 0xBFD5555555555555
%28 = getelementptr inbounds double, double* %0, i64 %25
%29 = load double, double* %28
%30 = tail call double @llvm.fmuladd.f64(double %27, double %29, double %23)
%31 = or i64 %22, 1
%32 = fmul double %27, 0xBFD5555555555555
%33 = getelementptr inbounds double, double* %0, i64 %22
%34 = load double, double* %33
%35 = tail call double @llvm.fmuladd.f64(double %32, double %34, double %30)
%36 = add nuw nsw i64 %22, 2
%37 = fmul double %32, 0xBFD5555555555555
%38 = getelementptr inbounds double, double* %0, i64 %31
%39 = load double, double* %38
%40 = tail call double @llvm.fmuladd.f64(double %37, double %39, double %35)
%41 = add nuw nsw i64 %22, 3
%42 = fmul double %37, 0xBFD5555555555555
%43 = getelementptr inbounds double, double* %0, i64 %36
%44 = load double, double* %43
%45 = tail call double @llvm.fmuladd.f64(double %42, double %44, double %40)
%46 = add nuw nsw i64 %22, 4
%47 = add i64 %26, 4
%48 = icmp eq i64 %47, %20
br i1 %48, label %49, label %21
49: ; preds = %21, %14
%50 = phi double [ undef, %14 ], [ %45, %21 ]
%51 = phi i64 [ 2, %14 ], [ %46, %21 ]
%52 = phi double [ %12, %14 ], [ %45, %21 ]
%53 = phi double [ 0xBFD5555555555555, %14 ], [ %42, %21 ]
%54 = phi i64 [ 1, %14 ], [ %41, %21 ]
%55 = icmp eq i64 %17, 0
br i1 %55, label %69, label %56
56: ; preds = %49, %56
%57 = phi i64 [ %66, %56 ], [ %51, %49 ]
%58 = phi double [ %65, %56 ], [ %52, %49 ]
%59 = phi double [ %62, %56 ], [ %53, %49 ]
%60 = phi i64 [ %57, %56 ], [ %54, %49 ]
%61 = phi i64 [ %67, %56 ], [ 0, %49 ]
%62 = fmul double %59, 0xBFD5555555555555
%63 = getelementptr inbounds double, double* %0, i64 %60
%64 = load double, double* %63
%65 = tail call double @llvm.fmuladd.f64(double %62, double %64, double %58)
%66 = add nuw nsw i64 %57, 1
%67 = add i64 %61, 1
%68 = icmp eq i64 %67, %17
br i1 %68, label %69, label %56
69: ; preds = %49, %56, %9, %4
%70 = phi double [ %7, %4 ], [ %12, %9 ], [ %50, %49 ], [ %65, %56 ]
store double %70, double* %1
%71 = icmp sgt i32 %3, 1
br i1 %71, label %72, label %106
72: ; preds = %69
%73 = zext i32 %3 to i64
%74 = add nsw i64 %73, -1
%75 = and i64 %74, 1
%76 = icmp eq i32 %3, 2
br i1 %76, label %96, label %77
77: ; preds = %72
%78 = and i64 %74, -2
br label %79
79: ; preds = %79, %77
%80 = phi double [ %70, %77 ], [ %91, %79 ]
%81 = phi i64 [ 1, %77 ], [ %93, %79 ]
%82 = phi i64 [ 0, %77 ], [ %94, %79 ]
%83 = add nsw i64 %81, -1
%84 = getelementptr inbounds double, double* %0, i64 %83
%85 = load double, double* %84
%86 = tail call double @llvm.fmuladd.f64(double %80, double 0xBFD5555555555555, double %85)
%87 = getelementptr inbounds double, double* %1, i64 %81
store double %86, double* %87
%88 = add nuw nsw i64 %81, 1
%89 = getelementptr inbounds double, double* %0, i64 %81
%90 = load double, double* %89
%91 = tail call double @llvm.fmuladd.f64(double %86, double 0xBFD5555555555555, double %90)
%92 = getelementptr inbounds double, double* %1, i64 %88
store double %91, double* %92
%93 = add nuw nsw i64 %81, 2
%94 = add i64 %82, 2
%95 = icmp eq i64 %94, %78
br i1 %95, label %96, label %79
96: ; preds = %79, %72
%97 = phi double [ %70, %72 ], [ %91, %79 ]
%98 = phi i64 [ 1, %72 ], [ %93, %79 ]
%99 = icmp eq i64 %75, 0
br i1 %99, label %106, label %100
100: ; preds = %96
%101 = add nsw i64 %98, -1
%102 = getelementptr inbounds double, double* %0, i64 %101
%103 = load double, double* %102
%104 = tail call double @llvm.fmuladd.f64(double %97, double 0xBFD5555555555555, double %103)
%105 = getelementptr inbounds double, double* %1, i64 %98
store double %104, double* %105
br label %106
106: ; preds = %100, %96, %69
%107 = add nsw i32 %3, -1
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds double, double* %0, i64 %108
%110 = load double, double* %109
br i1 %8, label %171, label %111
111: ; preds = %106
%112 = sext i32 %3 to i64
%113 = add nuw nsw i32 %6, 1
%114 = zext i32 %113 to i64
%115 = add nsw i64 %114, -1
%116 = add nsw i64 %114, -2
%117 = and i64 %115, 3
%118 = icmp ult i64 %116, 3
br i1 %118, label %152, label %119
119: ; preds = %111
%120 = and i64 %115, -4
br label %121
121: ; preds = %121, %119
%122 = phi i64 [ 1, %119 ], [ %149, %121 ]
%123 = phi double [ %110, %119 ], [ %147, %121 ]
%124 = phi double [ 0xBFD5555555555555, %119 ], [ %148, %121 ]
%125 = phi i64 [ 0, %119 ], [ %150, %121 ]
%126 = sub nsw i64 %112, %122
%127 = getelementptr inbounds double, double* %0, i64 %126
%128 = load double, double* %127
%129 = tail call double @llvm.fmuladd.f64(double %124, double %128, double %123)
%130 = fmul double %124, 0xBFD5555555555555
%131 = xor i64 %122, -1
%132 = add i64 %131, %112
%133 = getelementptr inbounds double, double* %0, i64 %132
%134 = load double, double* %133
%135 = tail call double @llvm.fmuladd.f64(double %130, double %134, double %129)
%136 = fmul double %130, 0xBFD5555555555555
%137 = add nuw nsw i64 %122, 2
%138 = sub nsw i64 %112, %137
%139 = getelementptr inbounds double, double* %0, i64 %138
%140 = load double, double* %139
%141 = tail call double @llvm.fmuladd.f64(double %136, double %140, double %135)
%142 = fmul double %136, 0xBFD5555555555555
%143 = add nuw nsw i64 %122, 3
%144 = sub nsw i64 %112, %143
%145 = getelementptr inbounds double, double* %0, i64 %144
%146 = load double, double* %145
%147 = tail call double @llvm.fmuladd.f64(double %142, double %146, double %141)
%148 = fmul double %142, 0xBFD5555555555555
%149 = add nuw nsw i64 %122, 4
%150 = add nuw i64 %125, 4
%151 = icmp eq i64 %150, %120
br i1 %151, label %152, label %121
152: ; preds = %121, %111
%153 = phi double [ undef, %111 ], [ %147, %121 ]
%154 = phi i64 [ 1, %111 ], [ %149, %121 ]
%155 = phi double [ %110, %111 ], [ %147, %121 ]
%156 = phi double [ 0xBFD5555555555555, %111 ], [ %148, %121 ]
%157 = icmp eq i64 %117, 0
br i1 %157, label %171, label %158
158: ; preds = %152, %158
%159 = phi i64 [ %168, %158 ], [ %154, %152 ]
%160 = phi double [ %166, %158 ], [ %155, %152 ]
%161 = phi double [ %167, %158 ], [ %156, %152 ]
%162 = phi i64 [ %169, %158 ], [ 0, %152 ]
%163 = sub nsw i64 %112, %159
%164 = getelementptr inbounds double, double* %0, i64 %163
%165 = load double, double* %164
%166 = tail call double @llvm.fmuladd.f64(double %161, double %165, double %160)
%167 = fmul double %161, 0xBFD5555555555555
%168 = add nuw nsw i64 %159, 1
%169 = add i64 %162, 1
%170 = icmp eq i64 %169, %117
br i1 %170, label %171, label %158
171: ; preds = %152, %158, %106
%172 = phi double [ %110, %106 ], [ %153, %152 ], [ %166, %158 ]
%173 = getelementptr inbounds double, double* %2, i64 %108
store double %172, double* %173
br i1 %71, label %174, label %269
174: ; preds = %171
%175 = add nsw i32 %3, -2
%176 = zext i32 %175 to i64
%177 = and i64 %176, 1
%178 = icmp eq i64 %177, 0
br i1 %178, label %179, label %188
179: ; preds = %174
%180 = getelementptr inbounds double, double* %0, i64 %176
%181 = load double, double* %180
%182 = add nuw nsw i64 %176, 1
%183 = getelementptr inbounds double, double* %2, i64 %182
%184 = load double, double* %183
%185 = tail call double @llvm.fmuladd.f64(double %184, double 0xBFD5555555555555, double %181)
%186 = getelementptr inbounds double, double* %2, i64 %176
store double %185, double* %186
%187 = add nsw i64 %176, -1
br label %188
188: ; preds = %179, %174
%189 = phi i64 [ %176, %174 ], [ %187, %179 ]
%190 = icmp eq i32 %175, 0
br i1 %190, label %191, label %236
191: ; preds = %236, %188
br i1 %71, label %192, label %269
192: ; preds = %191
%193 = zext i32 %107 to i64
%194 = icmp ult i32 %107, 2
br i1 %194, label %234, label %195
195: ; preds = %192
%196 = getelementptr double, double* %1, i64 %193
%197 = getelementptr double, double* %2, i64 %193
%198 = add nuw nsw i64 %193, 1
%199 = getelementptr double, double* %0, i64 %198
%200 = icmp ugt double* %197, %1
%201 = icmp ugt double* %196, %2
%202 = and i1 %200, %201
%203 = icmp ugt double* %199, %1
%204 = icmp ugt double* %196, %0
%205 = and i1 %203, %204
%206 = or i1 %202, %205
br i1 %206, label %234, label %207
207: ; preds = %195
%208 = and i64 %193, 4294967294
br label %209
209: ; preds = %209, %207
%210 = phi i64 [ 0, %207 ], [ %230, %209 ]
%211 = getelementptr inbounds double, double* %2, i64 %210
%212 = bitcast double* %211 to <2 x double>*
%213 = load <2 x double>, <2 x double>* %212
%214 = getelementptr inbounds double, double* %1, i64 %210
%215 = bitcast double* %214 to <2 x double>*
%216 = load <2 x double>, <2 x double>* %215
%217 = fmul <2 x double> %216, <double 0xBFEC71C71C71C71C, double 0xBFEC71C71C71C71C>
%218 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %213, <2 x double> <double -8.000000e+00, double -8.000000e+00>, <2 x double> %217)
%219 = or i64 %210, 1
%220 = getelementptr inbounds double, double* %0, i64 %219
%221 = bitcast double* %220 to <2 x double>*
%222 = load <2 x double>, <2 x double>* %221
%223 = fadd <2 x double> %222, %218
%224 = getelementptr inbounds double, double* %0, i64 %210
%225 = bitcast double* %224 to <2 x double>*
%226 = load <2 x double>, <2 x double>* %225
%227 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %226, <2 x double> <double 0x4027555555555555, double 0x4027555555555555>, <2 x double> %223)
%228 = fdiv <2 x double> %227, <double 6.000000e+00, double 6.000000e+00>
%229 = bitcast double* %214 to <2 x double>*
store <2 x double> %228, <2 x double>* %229
%230 = add nuw i64 %210, 2
%231 = icmp eq i64 %230, %208
br i1 %231, label %232, label %209
232: ; preds = %209
%233 = icmp eq i64 %208, %193
br i1 %233, label %269, label %234
234: ; preds = %195, %192, %232
%235 = phi i64 [ 0, %195 ], [ 0, %192 ], [ %208, %232 ]
br label %252
236: ; preds = %188, %236
%237 = phi i64 [ %250, %236 ], [ %189, %188 ]
%238 = getelementptr inbounds double, double* %0, i64 %237
%239 = load double, double* %238
%240 = add nuw nsw i64 %237, 1
%241 = getelementptr inbounds double, double* %2, i64 %240
%242 = load double, double* %241
%243 = tail call double @llvm.fmuladd.f64(double %242, double 0xBFD5555555555555, double %239)
%244 = getelementptr inbounds double, double* %2, i64 %237
store double %243, double* %244
%245 = add nsw i64 %237, -1
%246 = getelementptr inbounds double, double* %0, i64 %245
%247 = load double, double* %246
%248 = tail call double @llvm.fmuladd.f64(double %243, double 0xBFD5555555555555, double %247)
%249 = getelementptr inbounds double, double* %2, i64 %245
store double %248, double* %249
%250 = add nsw i64 %237, -2
%251 = icmp eq i64 %245, 0
br i1 %251, label %191, label %236
252: ; preds = %234, %252
%253 = phi i64 [ %260, %252 ], [ %235, %234 ]
%254 = getelementptr inbounds double, double* %2, i64 %253
%255 = load double, double* %254
%256 = getelementptr inbounds double, double* %1, i64 %253
%257 = load double, double* %256
%258 = fmul double %257, 0xBFEC71C71C71C71C
%259 = tail call double @llvm.fmuladd.f64(double %255, double -8.000000e+00, double %258)
%260 = add nuw nsw i64 %253, 1
%261 = getelementptr inbounds double, double* %0, i64 %260
%262 = load double, double* %261
%263 = fadd double %262, %259
%264 = getelementptr inbounds double, double* %0, i64 %253
%265 = load double, double* %264
%266 = tail call double @llvm.fmuladd.f64(double %265, double 0x4027555555555555, double %263)
%267 = fdiv double %266, 6.000000e+00
store double %267, double* %256
%268 = icmp eq i64 %260, %193
br i1 %268, label %269, label %252
269: ; preds = %252, %232, %171, %191
%270 = load double, double* %173
%271 = getelementptr inbounds double, double* %1, i64 %108
%272 = load double, double* %271
%273 = fmul double %272, 0xBFEC71C71C71C71C
%274 = tail call double @llvm.fmuladd.f64(double %270, double -8.000000e+00, double %273)
%275 = load double, double* %109
%276 = fadd double %275, %274
%277 = tail call double @llvm.fmuladd.f64(double %275, double 0x4027555555555555, double %276)
%278 = fdiv double %277, 6.000000e+00
store double %278, double* %271
ret void
}