user: 请从给定的汇编指令推导出C源代码:

```asm
<kernel_align_panel_3_3_lib4>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %edi,0x10(%rsp)
  mov    %rsi,0x8(%rsp)
  test   %edi,%edi
  jle    L38
  mov    %edi,%eax
  mov    %rcx,%r13
  shl    $0x2,%edx
  movslq %edx,%rdx
  lea    (%rsi,%rdx,8),%r14
  cmp    $0x3,%edi
  jle    L35
  lea    0x18(%rsi),%rbp
  sub    $0x4,%eax
  shr    $0x2,%eax
  mov    %eax,0x14(%rsp)
  mov    %eax,%eax
  mov    %rax,0x18(%rsp)
  shl    $0x7,%rax
  lea    0x98(%rsi,%rax,1),%r15
  mov    %r14,%r12
  mov    %rcx,%rbx
L3:
  mov    %rbp,%rdi
  call   <_mm_load_sd@plt>
  mov    %eax,%esi
  mov    %rbx,%rdi
  call   <_mm_store_sd@plt>
  mov    %r12,%rdi
  call   <_mm_loadu_pd@plt>
  mov    %eax,%esi
  lea    0x8(%rbx),%rdi
  call   <_mm_storeu_pd@plt>
  lea    0x20(%rbp),%rdi
  call   <_mm_load_sd@plt>
  mov    %eax,%esi
  lea    0x20(%rbx),%rdi
  call   <_mm_store_sd@plt>
  lea    0x20(%r12),%rdi
  call   <_mm_loadu_pd@plt>
  mov    %eax,%esi
  lea    0x28(%rbx),%rdi
  call   <_mm_storeu_pd@plt>
  lea    0x40(%rbp),%rdi
  call   <_mm_load_sd@plt>
  mov    %eax,%esi
  lea    0x40(%rbx),%rdi
  call   <_mm_store_sd@plt>
  lea    0x40(%r12),%rdi
  call   <_mm_loadu_pd@plt>
  mov    %eax,%esi
  lea    0x48(%rbx),%rdi
  call   <_mm_storeu_pd@plt>
  lea    0x60(%rbp),%rdi
  call   <_mm_load_sd@plt>
  mov    %eax,%esi
  lea    0x60(%rbx),%rdi
  call   <_mm_store_sd@plt>
  lea    0x60(%r12),%rdi
  call   <_mm_loadu_pd@plt>
  mov    %eax,%esi
  lea    0x68(%rbx),%rdi
  call   <_mm_storeu_pd@plt>
  sub    $0xffffffffffffff80,%r12
  sub    $0xffffffffffffff80,%rbx
  sub    $0xffffffffffffff80,%rbp
  cmp    %r15,%rbp
  jne    L3
  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),%ebx
  cmp    %ebx,0x10(%rsp)
  jle    L38
L25:
  mov    0x8(%rsp),%rbp
  add    $0x18,%rbp
L27:
  mov    %rbp,%rdi
  call   <_mm_load_sd@plt>
  mov    %eax,%esi
  mov    %r13,%rdi
  call   <_mm_store_sd@plt>
  mov    %r14,%rdi
  call   <_mm_loadu_pd@plt>
  mov    %eax,%esi
  lea    0x8(%r13),%rdi
  call   <_mm_storeu_pd@plt>
  add    $0x20,%r14
  add    $0x20,%r13
  add    $0x1,%ebx
  add    $0x20,%rbp
  cmp    %ebx,0x10(%rsp)
  jg     L27
L38:
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L35:
  mov    $0x0,%ebx
  jmp    L25
```

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