user: What would the C source look like for the assembly code provided here?

```asm
<kernel_align_panel_3_3_lib4>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %rsi,0x8(%rsp)
  test   %edi,%edi
  jle    L16
  shl    $0x2,%edx
  mov    %edi,%r15d
  mov    %rcx,%r13
  movslq %edx,%rdx
  lea    (%rsi,%rdx,8),%r14
  cmp    $0x3,%edi
  jle    L38
  lea    -0x4(%rdi),%eax
  lea    0x18(%rsi),%rbp
  mov    %r14,%r12
  mov    %rcx,%rbx
  shr    $0x2,%eax
  mov    %eax,0x14(%rsp)
  mov    %rax,0x18(%rsp)
  shl    $0x7,%rax
  lea    0x98(%rsi,%rax,1),%rax
  mov    %rax,(%rsp)
  nopl   0x0(%rax)
L5:
  mov    %rbp,%rdi
  call   <_mm_load_sd@plt>
  mov    %rbx,%rdi
  mov    %eax,%esi
  call   <_mm_store_sd@plt>
  mov    %r12,%rdi
  call   <_mm_loadu_pd@plt>
  lea    0x8(%rbx),%rdi
  mov    %eax,%esi
  call   <_mm_storeu_pd@plt>
  lea    0x20(%rbp),%rdi
  call   <_mm_load_sd@plt>
  lea    0x20(%rbx),%rdi
  mov    %eax,%esi
  call   <_mm_store_sd@plt>
  lea    0x20(%r12),%rdi
  call   <_mm_loadu_pd@plt>
  lea    0x28(%rbx),%rdi
  mov    %eax,%esi
  call   <_mm_storeu_pd@plt>
  lea    0x40(%rbp),%rdi
  call   <_mm_load_sd@plt>
  lea    0x40(%rbx),%rdi
  mov    %eax,%esi
  call   <_mm_store_sd@plt>
  lea    0x40(%r12),%rdi
  call   <_mm_loadu_pd@plt>
  lea    0x48(%rbx),%rdi
  mov    %eax,%esi
  call   <_mm_storeu_pd@plt>
  lea    0x60(%rbp),%rdi
  sub    $0xffffffffffffff80,%rbp
  call   <_mm_load_sd@plt>
  lea    0x60(%rbx),%rdi
  mov    %eax,%esi
  call   <_mm_store_sd@plt>
  lea    0x60(%r12),%rdi
  sub    $0xffffffffffffff80,%r12
  call   <_mm_loadu_pd@plt>
  lea    0x68(%rbx),%rdi
  sub    $0xffffffffffffff80,%rbx
  mov    %eax,%esi
  call   <_mm_storeu_pd@plt>
  cmp    (%rsp),%rbp
  jne    L5
  mov    0x18(%rsp),%rax
  add    $0x1,%rax
  shl    $0x7,%rax
  add    %rax,0x8(%rsp)
  add    %rax,%r14
  add    %rax,%r13
  mov    0x14(%rsp),%eax
  lea    0x4(,%rax,4),%ebp
  cmp    %ebp,%r15d
  jle    L16
L27:
  mov    0x8(%rsp),%rbx
  add    $0x18,%rbx
L41:
  mov    %rbx,%rdi
  add    $0x1,%ebp
  add    $0x20,%rbx
  call   <_mm_load_sd@plt>
  mov    %r13,%rdi
  mov    %eax,%esi
  call   <_mm_store_sd@plt>
  mov    %r14,%rdi
  add    $0x20,%r14
  call   <_mm_loadu_pd@plt>
  lea    0x8(%r13),%rdi
  add    $0x20,%r13
  mov    %eax,%esi
  call   <_mm_storeu_pd@plt>
  cmp    %ebp,%r15d
  jg     L41
L16:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L38:
  xor    %ebp,%ebp
  jmp    L27
```

assistant: ```c
void kernel_align_panel_3_3_lib4(int kmax, double *A0, int sda, double *B) {
  if (kmax <= 0)
    return;
  const int bs = 4;
  double *A1 = A0 + bs * sda;
  __m128d a_0, a_1;
  int k;
  for (k = 0; k < kmax - 3; k += 4) {
    a_0 = _mm_load_sd(&A0[3 + bs * 0]);
    _mm_store_sd(&B[0 + bs * 0], a_0);
    a_1 = _mm_loadu_pd(&A1[0 + bs * 0]);
    _mm_storeu_pd(&B[1 + bs * 0], a_1);
    a_0 = _mm_load_sd(&A0[3 + bs * 1]);
    _mm_store_sd(&B[0 + bs * 1], a_0);
    a_1 = _mm_loadu_pd(&A1[0 + bs * 1]);
    _mm_storeu_pd(&B[1 + bs * 1], a_1);
    a_0 = _mm_load_sd(&A0[3 + bs * 2]);
    _mm_store_sd(&B[0 + bs * 2], a_0);
    a_1 = _mm_loadu_pd(&A1[0 + bs * 2]);
    _mm_storeu_pd(&B[1 + bs * 2], a_1);
    a_0 = _mm_load_sd(&A0[3 + bs * 3]);
    _mm_store_sd(&B[0 + bs * 3], a_0);
    a_1 = _mm_loadu_pd(&A1[0 + bs * 3]);
    _mm_storeu_pd(&B[1 + bs * 3], a_1);
    A0 += 16;
    A1 += 16;
    B += 16;
  }
  for (; k < kmax; k++) {
    a_0 = _mm_load_sd(&A0[3 + bs * 0]);
    _mm_store_sd(&B[0 + bs * 0], a_0);
    a_1 = _mm_loadu_pd(&A1[0 + bs * 0]);
    _mm_storeu_pd(&B[1 + bs * 0], a_1);
    A0 += 4;
    A1 += 4;
    B += 4;
  }
}
```