user: Please show me the C implementation that matches the assembly code listed below:

```asm
<interpolate_product_single>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %rsi,-0x50(%rbp)
  mov    %edx,-0x54(%rbp)
  mov    %ecx,-0x58(%rbp)
  mov    %r8,-0x60(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  mov    $0x0,%eax
  call   <_mm_setzero_ps@plt>
  mov    %eax,-0x3c(%rbp)
  mov    -0x60(%rbp),%rax
  mov    %rax,%rdi
  call   <_mm_loadu_ps@plt>
  mov    %eax,-0x34(%rbp)
  movl   $0x0,-0x38(%rbp)
  jmp    L27
L5:
  mov    -0x38(%rbp),%eax
  imul   -0x58(%rbp),%eax
  mov    %eax,%eax
  lea    0x0(,%rax,4),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,%rdi
  call   <_mm_loadu_ps@plt>
  mov    %eax,%ebx
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,%rdi
  call   <_mm_load1_ps@plt>
  mov    %ebx,%esi
  mov    %eax,%edi
  call   <_mm_mul_ps@plt>
  mov    %eax,%edx
  mov    -0x3c(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_add_ps@plt>
  mov    %eax,-0x3c(%rbp)
  mov    -0x38(%rbp),%eax
  add    $0x1,%eax
  imul   -0x58(%rbp),%eax
  mov    %eax,%eax
  lea    0x0(,%rax,4),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,%rdi
  call   <_mm_loadu_ps@plt>
  mov    %eax,%ebx
  mov    -0x38(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,%rdi
  call   <_mm_load1_ps@plt>
  mov    %ebx,%esi
  mov    %eax,%edi
  call   <_mm_mul_ps@plt>
  mov    %eax,%edx
  mov    -0x3c(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_add_ps@plt>
  mov    %eax,-0x3c(%rbp)
  addl   $0x2,-0x38(%rbp)
L27:
  mov    -0x38(%rbp),%eax
  cmp    %eax,-0x54(%rbp)
  ja     L5
  mov    -0x3c(%rbp),%edx
  mov    -0x34(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_mul_ps@plt>
  mov    %eax,-0x3c(%rbp)
  mov    -0x3c(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_movehl_ps@plt>
  mov    %eax,%edx
  mov    -0x3c(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_add_ps@plt>
  mov    %eax,-0x3c(%rbp)
  mov    -0x3c(%rbp),%ecx
  mov    -0x3c(%rbp),%eax
  mov    $0x55,%edx
  mov    %ecx,%esi
  mov    %eax,%edi
  call   <_mm_shuffle_ps@plt>
  mov    %eax,%edx
  mov    -0x3c(%rbp),%eax
  mov    %edx,%esi
  mov    %eax,%edi
  call   <_mm_add_ss@plt>
  mov    %eax,-0x3c(%rbp)
  mov    -0x3c(%rbp),%edx
  lea    -0x40(%rbp),%rax
  mov    %edx,%esi
  mov    %rax,%rdi
  call   <_mm_store_ss@plt>
  movss  -0x40(%rbp),%xmm0
  mov    -0x18(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L8
  call   <__stack_chk_fail@plt>
L8:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

assistant: ```c
float interpolate_product_single(const float *a, const float *b, unsigned int len, const unsigned int oversample, float *frac) {
  int i;
  float ret;
  if (1) {
    __m128 sum = _mm_setzero_ps();
    __m128 f = _mm_loadu_ps(frac);
    for (i = 0; i < len; i += 2) {
      sum = _mm_add_ps(sum, _mm_mul_ps(_mm_load1_ps(a + i), _mm_loadu_ps(b + i * oversample)));
      sum = _mm_add_ps(sum, _mm_mul_ps(_mm_load1_ps(a + i + 1), _mm_loadu_ps(b + (i + 1) * oversample)));
    }
    sum = _mm_mul_ps(f, sum);
    sum = _mm_add_ps(sum, _mm_movehl_ps(sum, sum));
    sum = _mm_add_ss(sum, _mm_shuffle_ps(sum, sum, 0x55));
    _mm_store_ss(&ret, sum);
  } else {
    float accum[4] = {0, 0, 0, 0};
    for (i = 0; i < len; i++) {
      const float curr_in = a[i];
      accum[0] += curr_in * b[i * oversample + 0];
      accum[1] += curr_in * b[i * oversample + 1];
      accum[2] += curr_in * b[i * oversample + 2];
      accum[3] += curr_in * b[i * oversample + 3];
    }
    ret = accum[0] * frac[0] + accum[1] * frac[1] + accum[2] * frac[2] + accum[3] * frac[3];
  }
  return ret;
}
```