user: Converte this assembly snippet into its equivalent C code:

```asm
<filter_mem16_8>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,%r14
  mov    %rsi,%r13
  mov    %rdx,%rbp
  mov    %rcx,%r15
  mov    %r8d,%r12d
  mov    0x70(%rsp),%rdi
  call   <_mm_loadu_ps@plt>
  mov    %eax,%ebx
  mov    %r13,%rdi
  call   <_mm_loadu_ps@plt>
  mov    %eax,0x18(%rsp)
  mov    %rbp,%rdi
  call   <_mm_loadu_ps@plt>
  mov    %eax,0x1c(%rsp)
  mov    0x70(%rsp),%rax
  add    $0x10,%rax
  mov    %rax,0x28(%rsp)
  mov    %rax,%rdi
  call   <_mm_loadu_ps@plt>
  shl    $0x20,%rax
  or     %rax,%rbx
  lea    0x10(%r13),%rdi
  call   <_mm_loadu_ps@plt>
  mov    %eax,0x10(%rsp)
  lea    0x10(%rbp),%rdi
  call   <_mm_loadu_ps@plt>
  mov    %eax,0x14(%rsp)
  test   %r12d,%r12d
  jle    L14
  mov    %r14,%r13
  lea    -0x1(%r12),%eax
  lea    0x4(%r14,%rax,4),%rax
  mov    %rax,0x20(%rsp)
L4:
  mov    %r13,%rdi
  call   <_mm_load_ps1@plt>
  mov    %eax,%ebp
  mov    %ebx,%esi
  mov    %eax,%edi
  call   <_mm_add_ss@plt>
  mov    %eax,%r12d
  mov    %eax,%esi
  mov    %r15,%rdi
  call   <_mm_store_ss@plt>
  mov    $0x0,%edx
  mov    %r12d,%esi
  mov    %r12d,%edi
  call   <_mm_shuffle_ps@plt>
  mov    %eax,0xc(%rsp)
  mov    %rbx,%rcx
  shr    $0x20,%rcx
  mov    %rcx,%r12
  mov    %ecx,%esi
  mov    %ebx,%edi
  call   <_mm_move_ss@plt>
  mov    %eax,%edi
  mov    $0x39,%edx
  mov    %eax,%esi
  call   <_mm_shuffle_ps@plt>
  mov    %eax,%r14d
  mov    0x18(%rsp),%esi
  mov    %ebp,%edi
  call   <_mm_mul_ps@plt>
  mov    %eax,%esi
  mov    %r14d,%edi
  call   <_mm_add_ps@plt>
  mov    %eax,%r14d
  mov    0x1c(%rsp),%esi
  mov    0xc(%rsp),%edi
  call   <_mm_mul_ps@plt>
  mov    %eax,%esi
  mov    %r14d,%edi
  call   <_mm_sub_ps@plt>
  mov    %eax,%eax
  movabs $0xffffffff00000000,%rdx
  and    %rdx,%rbx
  or     %rax,%rbx
  mov    %r12d,%esi
  mov    %r12d,%edi
  call   <_mm_sub_ss@plt>
  mov    %eax,%edi
  mov    $0x39,%edx
  mov    %eax,%esi
  call   <_mm_shuffle_ps@plt>
  mov    %eax,%r12d
  mov    0x10(%rsp),%esi
  mov    %ebp,%edi
  call   <_mm_mul_ps@plt>
  mov    %eax,%esi
  mov    %r12d,%edi
  call   <_mm_add_ps@plt>
  mov    %eax,%ebp
  mov    0x14(%rsp),%esi
  mov    0xc(%rsp),%edi
  call   <_mm_mul_ps@plt>
  mov    %eax,%esi
  mov    %ebp,%edi
  call   <_mm_sub_ps@plt>
  shl    $0x20,%rax
  mov    %ebx,%ebx
  or     %rax,%rbx
  add    $0x4,%r13
  add    $0x4,%r15
  cmp    0x20(%rsp),%r13
  jne    L4
L14:
  mov    %ebx,%esi
  mov    0x70(%rsp),%rdi
  call   <_mm_storeu_ps@plt>
  shr    $0x20,%rbx
  mov    %rbx,%rsi
  mov    0x28(%rsp),%rdi
  call   <_mm_storeu_ps@plt>
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

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]);
}
```