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

```asm
<dedisperse_kernel_lagged>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %rdi,%rbp
  mov    %rsi,%r15
  mov    %rsi,0x18(%rsp)
  mov    %edx,%ebx
  mov    %ecx,%r12d
  mov    %ecx,0xc(%rsp)
  mov    %r8d,%r14d
  mov    %r8d,0x14(%rsp)
  mov    %edx,%edi
  xor    $0x1,%edi
  and    $0x1,%edi
  call   <assert@plt>
  mov    %ebx,%eax
  shr    $0x1f,%eax
  add    %ebx,%eax
  sar    %eax
  mov    %eax,0x28(%rsp)
  cmp    $0x1,%ebx
  jle    L43
  mov    %eax,%edi
  mov    %rbp,%rax
  mov    %r15,%rcx
  mov    %edi,0x10(%rsp)
  mov    %r14d,%r11d
  lea    -0x1(%r14),%r8d
  lea    0x4(,%r8,4),%r15
  lea    -0x1(%r12),%edx
  sub    %r14d,%edx
  lea    0x4(,%rdx,4),%r14
  movslq %r11d,%rsi
  lea    0x0(,%rsi,4),%rdi
  mov    %r12d,%edx
  add    %r11d,%edx
  mov    %edx,0x24(%rsp)
  movslq %r12d,%rdx
  lea    0x0(,%rdx,4),%r10
  mov    %r10,0x38(%rsp)
  lea    0x1(%rdx,%r8,1),%r8
  lea    0x0(,%r8,4),%r10
  mov    %rsi,%rbp
  neg    %rbp
  shl    $0x2,%rbp
  mov    %r11d,%r9d
  lea    0x4(,%r9,4),%r13
  lea    -0x2(%r12),%r8d
  sub    %r11d,%r8d
  lea    0x4(,%r8,4),%r12
  sub    %rsi,%rdx
  lea    -0x4(,%rdx,4),%rbx
  mov    %rbx,0x30(%rsp)
  add    %rdx,%r9
  lea    0x0(,%r9,4),%rbx
  movl   $0x0,0x20(%rsp)
  lea    0x1(%r11),%edx
  mov    %edx,0x2c(%rsp)
  lea    0x4(%rdi),%r9
  mov    %r9,%r11
  mov    %rbx,%r9
  jmp    L54
L52:
  mov    0x8(%rax),%rsi
  movss  (%rsi,%rdx,1),%xmm0
  mov    (%rcx),%rsi
  movss  %xmm0,(%rsi,%rdx,1)
  add    $0x4,%rdx
  cmp    %r15,%rdx
  jne    L52
L58:
  mov    0xc(%rsp),%edx
  cmp    %edx,0x14(%rsp)
  jge    L148
  mov    $0x0,%edx
  mov    %r9,%r8
L28:
  lea    (%rdx,%rdi,1),%rsi
  mov    (%rax),%rbx
  mov    0x8(%rax),%r9
  movss  (%rbx,%rdx,1),%xmm0
  addss  (%r9,%rsi,1),%xmm0
  mov    (%rcx),%rbx
  movss  %xmm0,(%rbx,%rsi,1)
  add    $0x4,%rdx
  cmp    %r14,%rdx
  jne    L28
  mov    %r8,%r9
L148:
  mov    0x24(%rsp),%edx
  cmp    %edx,0xc(%rsp)
  jge    L87
  mov    0x38(%rsp),%rdx
L18:
  mov    %rdx,%rsi
  add    (%rax),%rsi
  movss  (%rsi,%rbp,1),%xmm0
  mov    (%rcx),%rsi
  movss  %xmm0,(%rsi,%rdx,1)
  add    $0x4,%rdx
  cmp    %r10,%rdx
  jne    L18
L87:
  cmpl   $0x0,0x14(%rsp)
  js     L99
  movslq 0x10(%rsp),%rdx
  mov    0x18(%rsp),%rbx
  lea    (%rbx,%rdx,8),%r8
  mov    $0x0,%edx
L7:
  mov    0x8(%rax),%rsi
  movss  (%rsi,%rdx,1),%xmm0
  mov    (%r8),%rsi
  movss  %xmm0,(%rsi,%rdx,1)
  add    $0x4,%rdx
  cmp    %r13,%rdx
  jne    L7
L99:
  mov    0x2c(%rsp),%edx
  cmp    %edx,0xc(%rsp)
  jle    L123
  movslq 0x10(%rsp),%rdx
  mov    0x18(%rsp),%rbx
  lea    (%rbx,%rdx,8),%rbx
  mov    $0x0,%edx
  mov    %r9,%r8
  mov    %r13,%r9
L98:
  lea    (%r11,%rdx,1),%rsi
  mov    (%rax),%r13
  movss  0x0(%r13,%rdx,1),%xmm0
  mov    0x8(%rax),%r13
  addss  0x0(%r13,%rsi,1),%xmm0
  mov    (%rbx),%r13
  movss  %xmm0,0x0(%r13,%rsi,1)
  add    $0x4,%rdx
  cmp    %r12,%rdx
  jne    L98
  mov    %r9,%r13
  mov    %r8,%r9
L123:
  mov    0x24(%rsp),%edx
  cmp    %edx,0xc(%rsp)
  jg     L62
  movslq 0x10(%rsp),%rdx
  mov    0x18(%rsp),%rbx
  lea    (%rbx,%rdx,8),%r8
  mov    0x30(%rsp),%rdx
L46:
  mov    (%rax),%rsi
  movss  (%rsi,%rdx,1),%xmm0
  mov    %rdx,%rsi
  add    (%r8),%rsi
  movss  %xmm0,0x4(%rdi,%rsi,1)
  add    $0x4,%rdx
  cmp    %rdx,%r9
  jne    L46
L62:
  addl   $0x1,0x20(%rsp)
  mov    0x20(%rsp),%ebx
  add    $0x10,%rax
  add    $0x8,%rcx
  addl   $0x1,0x10(%rsp)
  cmp    %ebx,0x28(%rsp)
  jle    L43
L54:
  mov    $0x0,%edx
  cmpl   $0x0,0x14(%rsp)
  jg     L52
  jmp    L58
L43:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void dedisperse_kernel_lagged(float **in, float **out, int n, int m, int s) {
  assert(n % 2 == 0);
  int npair = n / 2;
  for (int jj = 0; jj < npair; jj++) {
    for (int i = 0; i < s; i++)
      out[jj][i] = in[2 * jj + 1][i];
    for (int i = s; i < m; i++)
      out[jj][i] = in[2 * jj][i - s] + in[2 * jj + 1][i];
    for (int i = m; i < m + s; i++)
      out[jj][i] = in[2 * jj][i - s];
    for (int i = 0; i < s + 1; i++)
      out[jj + npair][i] = in[2 * jj + 1][i];
    for (int i = s + 1; i < m; i++)
      out[jj + npair][i] = in[2 * jj][i - s - 1] + in[2 * jj + 1][i];
    for (int i = m; i < m + s + 1; i++)
      out[jj + npair][i] = in[2 * jj][i - s - 1];
  }
}
```