user: What is the c source code of the assembly code below:

```asm
<filter_r3>:
  endbr64
  mov    $0x8,%r10d
  push   %rbx
  movsd  (%rdi),%xmm0
  mov    %ecx,%r9d
  cmp    %r10d,%ecx
  cmovle %ecx,%r10d
  test   %ecx,%ecx
  jle    L100
  movsd  D26(%rip),%xmm3
  movapd %xmm0,%xmm1
  mov    $0x2,%eax
  movapd %xmm3,%xmm2
  jmp    L23
  nopl   0x0(%rax,%rax,1)
L109:
  movsd  -0x8(%rdi,%rax,8),%xmm1
  add    $0x1,%rax
L23:
  mulsd  %xmm2,%xmm1
  mulsd  %xmm3,%xmm2
  addsd  %xmm1,%xmm0
  cmp    %eax,%r10d
  jge    L109
  movsd  %xmm0,(%rsi)
  cmp    $0x1,%r9d
  je     L35
  lea    -0x2(%r9),%ecx
  movsd  D78(%rip),%xmm1
  xor    %eax,%eax
  mov    %rcx,%r11
  shl    $0x3,%rcx
  jmp    L110
L61:
  add    $0x8,%rax
L110:
  movapd %xmm0,%xmm2
  movsd  (%rdi,%rax,1),%xmm0
  mulsd  %xmm1,%xmm2
  subsd  %xmm2,%xmm0
  movsd  %xmm0,0x8(%rsi,%rax,1)
  cmp    %rax,%rcx
  jne    L61
  movslq %r9d,%rax
  shl    $0x3,%rax
  lea    -0x8(%rax),%rcx
  lea    (%rdi,%rcx,1),%r8
  movsd  (%r8),%xmm1
L128:
  add    %rdi,%rax
  mov    $0x1,%ebx
  movapd %xmm3,%xmm0
L80:
  movsd  -0x8(%rax),%xmm2
  add    $0x1,%ebx
  sub    $0x8,%rax
  mulsd  %xmm0,%xmm2
  mulsd  %xmm3,%xmm0
  addsd  %xmm2,%xmm1
  cmp    %ebx,%r10d
  jge    L80
L60:
  lea    (%rdx,%rcx,1),%r10
  movslq %r11d,%rax
  movsd  %xmm1,(%r10)
  movsd  D78(%rip),%xmm1
  test   %r11d,%r11d
  js     L52
  cs nopw 0x0(%rax,%rax,1)
L33:
  movsd  0x8(%rdx,%rax,8),%xmm2
  movsd  (%rdi,%rax,8),%xmm0
  mulsd  %xmm1,%xmm2
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  sub    $0x1,%rax
  test   %eax,%eax
  jns    L33
L52:
  cmp    $0x1,%r9d
  jle    L14
  movsd  D22(%rip),%xmm1
  movsd  D92(%rip),%xmm3
  xor    %eax,%eax
  movsd  D88(%rip),%xmm5
  movsd  D80(%rip),%xmm4
  nopl   0x0(%rax)
L115:
  movsd  (%rdx,%rax,8),%xmm0
  movsd  (%rsi,%rax,8),%xmm2
  mov    %rax,%r9
  mulsd  %xmm3,%xmm2
  mulsd  %xmm1,%xmm0
  subsd  %xmm2,%xmm0
  movsd  (%rdi,%rax,8),%xmm2
  addsd  0x8(%rdi,%rax,8),%xmm0
  mulsd  %xmm5,%xmm2
  addsd  %xmm2,%xmm0
  divsd  %xmm4,%xmm0
  movsd  %xmm0,(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %r11,%r9
  jne    L115
L72:
  mulsd  (%r10),%xmm1
  add    %rsi,%rcx
  movsd  (%r8),%xmm0
  mulsd  (%rcx),%xmm3
  subsd  %xmm3,%xmm1
  addsd  %xmm0,%xmm1
  mulsd  %xmm5,%xmm0
  addsd  %xmm0,%xmm1
  divsd  %xmm4,%xmm1
  movsd  %xmm1,(%rcx)
  pop    %rbx
  ret
  nopl   0x0(%rax)
L100:
  movslq %ecx,%rax
  movsd  %xmm0,(%rsi)
  lea    -0x2(%r9),%r11d
  lea    -0x8(,%rax,8),%rcx
  lea    (%rdi,%rcx,1),%r8
  movsd  (%r8),%xmm1
  jmp    L60
  nopl   0x0(%rax)
L14:
  movsd  D22(%rip),%xmm1
  movsd  D92(%rip),%xmm3
  movsd  D88(%rip),%xmm5
  movsd  D80(%rip),%xmm4
  jmp    L72
L35:
  movsd  (%rdi),%xmm1
  mov    %rdi,%r8
  xor    %ecx,%ecx
  mov    $0x8,%eax
  mov    $0xffffffff,%r11d
  jmp    L128
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D26
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D78
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D22
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D92
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D88
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D80
type_guessed: f64

tool: {"D26": [0x55,0x55,0x55,0x55,0x55,0x55,0xd5,0xbf]}

tool: D78: 0.3333333333333333

tool: {"D22": [0x00,0x00,0x00,0x00,0x00,0x00,0x20,0xc0]}

tool: {"D92": [0x1c,0xc7,0x71,0x1c,0xc7,0x71,0xec,0x3f]}

tool: {"D88": [0x55,0x55,0x55,0x55,0x55,0x55,0x27,0x40]}

tool: D80: 6

assistant: ```c
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;
}
```