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

```asm
<filter_mem16_8>:
  endbr64
  push   %r15
  mov    %rdx,%r15
  push   %r14
  push   %r13
  mov    %rdi,%r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x38,%rsp
  mov    0x70(%rsp),%r14
  mov    %r8d,0x10(%rsp)
  mov    %rcx,0x28(%rsp)
  mov    %r14,%rdi
  mov    %r14,0x20(%rsp)
  call   <_mm_loadu_ps@plt>
  mov    %rbx,%rdi
  mov    %eax,%ebp
  call   <_mm_loadu_ps@plt>
  mov    %r15,%rdi
  mov    %eax,(%rsp)
  call   <_mm_loadu_ps@plt>
  mov    %eax,0x4(%rsp)
  lea    0x10(%r14),%rax
  mov    %rax,%rdi
  mov    %rax,0x18(%rsp)
  call   <_mm_loadu_ps@plt>
  lea    0x10(%rbx),%rdi
  mov    %eax,%r12d
  call   <_mm_loadu_ps@plt>
  lea    0x10(%r15),%rdi
  mov    %eax,0x8(%rsp)
  call   <_mm_loadu_ps@plt>
  mov    0x10(%rsp),%r8d
  mov    %eax,0xc(%rsp)
  test   %r8d,%r8d
  jle    L17
  lea    -0x1(%r8),%eax
  mov    0x28(%rsp),%rbx
  lea    0x4(%r13,%rax,4),%rax
  mov    %rax,0x10(%rsp)
  nop
L18:
  mov    %r13,%rdi
  add    $0x4,%r13
  call   <_mm_load_ps1@plt>
  mov    %ebp,%esi
  mov    %eax,%edi
  mov    %eax,%r15d
  call   <_mm_add_ss@plt>
  mov    %rbx,%rdi
  add    $0x4,%rbx
  mov    %eax,%r14d
  mov    %eax,%esi
  call   <_mm_store_ss@plt>
  xor    %edx,%edx
  mov    %r14d,%esi
  mov    %r14d,%edi
  call   <_mm_shuffle_ps@plt>
  mov    %r12d,%esi
  mov    %ebp,%edi
  mov    %eax,%r14d
  call   <_mm_move_ss@plt>
  mov    $0x39,%edx
  mov    %eax,%edi
  mov    %eax,%esi
  call   <_mm_shuffle_ps@plt>
  mov    (%rsp),%esi
  mov    %r15d,%edi
  mov    %eax,%ebp
  call   <_mm_mul_ps@plt>
  mov    %ebp,%edi
  mov    %eax,%esi
  call   <_mm_add_ps@plt>
  mov    0x4(%rsp),%esi
  mov    %r14d,%edi
  mov    %eax,%ebp
  call   <_mm_mul_ps@plt>
  mov    %ebp,%edi
  mov    %eax,%esi
  call   <_mm_sub_ps@plt>
  mov    %r12d,%esi
  mov    %r12d,%edi
  mov    %eax,%ebp
  call   <_mm_sub_ss@plt>
  mov    $0x39,%edx
  mov    %eax,%edi
  mov    %eax,%esi
  call   <_mm_shuffle_ps@plt>
  mov    0x8(%rsp),%esi
  mov    %r15d,%edi
  mov    %eax,%r12d
  call   <_mm_mul_ps@plt>
  mov    %r12d,%edi
  mov    %eax,%esi
  call   <_mm_add_ps@plt>
  mov    0xc(%rsp),%esi
  mov    %r14d,%edi
  mov    %eax,%r12d
  call   <_mm_mul_ps@plt>
  mov    %r12d,%edi
  mov    %eax,%esi
  call   <_mm_sub_ps@plt>
  mov    %eax,%r12d
  cmp    0x10(%rsp),%r13
  jne    L18
L17:
  mov    0x20(%rsp),%rdi
  mov    %ebp,%esi
  call   <_mm_storeu_ps@plt>
  mov    0x18(%rsp),%rdi
  add    $0x38,%rsp
  mov    %r12d,%esi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1190 <_mm_storeu_ps@plt>
```

assistant: ```c
void filter_mem16_8(const float *x, const float *_num, const float *_den, float *y, int N, int ord, float *_mem) {
  __m128 num[2], den[2], mem[2];
  int i;
  for (i = 0; i < 2; i++) {
    mem[i] = _mm_loadu_ps(_mem + 4 * i);
    num[i] = _mm_loadu_ps(_num + 4 * i);
    den[i] = _mm_loadu_ps(_den + 4 * i);
  }
  for (i = 0; i < N; i++) {
    __m128 xx;
    __m128 yy;
    xx = _mm_load_ps1(x + i);
    yy = _mm_add_ss(xx, mem[0]);
    _mm_store_ss(y + i, yy);
    yy = _mm_shuffle_ps(yy, yy, 0);
    mem[0] = _mm_move_ss(mem[0], mem[1]);
    mem[0] = _mm_shuffle_ps(mem[0], mem[0], 0x39);
    mem[0] = _mm_add_ps(mem[0], _mm_mul_ps(xx, num[0]));
    mem[0] = _mm_sub_ps(mem[0], _mm_mul_ps(yy, den[0]));
    mem[1] = _mm_sub_ss(mem[1], mem[1]);
    mem[1] = _mm_shuffle_ps(mem[1], mem[1], 0x39);
    mem[1] = _mm_add_ps(mem[1], _mm_mul_ps(xx, num[1]));
    mem[1] = _mm_sub_ps(mem[1], _mm_mul_ps(yy, den[1]));
  }
  _mm_storeu_ps(_mem, mem[0]);
  _mm_storeu_ps(_mem + 4, mem[1]);
}
```