user: What is the c source code of the assembly code below?

```asm
<kernel_dgetr_2_lib4>:
  endbr64
  mov    %esi,%r9d
  mov    %rdx,%rax
  mov    %rcx,%r10
  test   %esi,%esi
  jle    L129
  movslq %r9d,%rcx
  lea    -0x1(%rsi),%esi
  mov    %rcx,%rdx
  lea    0x20(%r10,%rcx,8),%rcx
  shl    $0x5,%rdx
  lea    -0x10(%rax,%rdx,1),%rdx
  cmp    %rdx,%r10
  setae  %dl
  cmp    %rcx,%rax
  setae  %cl
  or     %cl,%dl
  je     L54
  cmp    $0x3,%esi
  jbe    L54
  mov    %esi,%r11d
  mov    %rax,%rdx
  mov    %r10,%rcx
  shr    %r11d
  sub    $0x1,%r11d
  shl    $0x6,%r11
  lea    0x40(%rax,%r11,1),%r11
  nopw   0x0(%rax,%rax,1)
L53:
  movupd (%rdx),%xmm0
  movupd 0x20(%rdx),%xmm1
  add    $0x40,%rdx
  add    $0x10,%rcx
  movupd -0x30(%rdx),%xmm4
  movupd -0x10(%rdx),%xmm3
  movapd %xmm0,%xmm2
  movapd %xmm1,%xmm5
  unpcklpd %xmm4,%xmm2
  unpcklpd %xmm3,%xmm5
  unpckhpd %xmm4,%xmm0
  unpckhpd %xmm3,%xmm1
  unpcklpd %xmm5,%xmm2
  unpcklpd %xmm1,%xmm0
  movups %xmm2,-0x10(%rcx)
  movups %xmm0,0x10(%rcx)
  cmp    %rdx,%r11
  jne    L53
  and    $0xfffffffe,%esi
  mov    %r9d,%r11d
  mov    %esi,%ecx
  add    $0x1,%esi
  mov    %rcx,%rdx
  lea    (%r10,%rcx,8),%rcx
  shl    $0x5,%rdx
  add    %rax,%rdx
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rcx)
  movsd  0x8(%rdx),%xmm0
  movsd  %xmm0,0x20(%rcx)
  cmp    %esi,%r9d
  jle    L119
  movsd  0x20(%rdx),%xmm0
  movsd  %xmm0,0x8(%rcx)
  movsd  0x28(%rdx),%xmm0
  movsd  %xmm0,0x28(%rcx)
L119:
  mov    %r11,%rdx
  shl    $0x5,%rdx
  add    %rdx,%rax
  lea    -0x4(,%r8,4),%edx
  movslq %edx,%rdx
  add    %r11,%rdx
  lea    (%r10,%rdx,8),%r10
  lea    -0x3(%rdi),%edx
  cmp    %edx,%r9d
  jge    L75
L13:
  lea    0x0(,%r8,4),%ecx
  lea    -0x4(%rdi),%r8d
  mov    %r10,%rdx
  sub    %r9d,%r8d
  movslq %ecx,%rcx
  shr    $0x2,%r8d
  shl    $0x3,%rcx
  lea    0x1(%r8),%r11d
  mov    %r11,%rsi
  shl    $0x7,%rsi
  add    %rax,%rsi
  nopw   0x0(%rax,%rax,1)
L99:
  movsd  (%rax),%xmm0
  sub    $0xffffffffffffff80,%rax
  movsd  %xmm0,(%rdx)
  movsd  -0x78(%rax),%xmm0
  movsd  %xmm0,0x20(%rdx)
  movsd  -0x60(%rax),%xmm0
  movsd  %xmm0,0x8(%rdx)
  movsd  -0x58(%rax),%xmm0
  movsd  %xmm0,0x28(%rdx)
  movsd  -0x40(%rax),%xmm0
  movsd  %xmm0,0x10(%rdx)
  movsd  -0x38(%rax),%xmm0
  movsd  %xmm0,0x30(%rdx)
  movsd  -0x20(%rax),%xmm0
  movsd  %xmm0,0x18(%rdx)
  movsd  -0x18(%rax),%xmm0
  movsd  %xmm0,0x38(%rdx)
  add    %rcx,%rdx
  cmp    %rsi,%rax
  jne    L99
  imul   %r11,%rcx
  lea    0x4(%r9,%r8,4),%r9d
  add    %rcx,%r10
L79:
  cmp    %r9d,%edi
  jle    L56
  movsd  (%rsi),%xmm0
  lea    0x1(%r9),%eax
  movsd  %xmm0,(%r10)
  movsd  0x8(%rsi),%xmm0
  movsd  %xmm0,0x20(%r10)
  cmp    %eax,%edi
  jle    L56
  movsd  0x20(%rsi),%xmm0
  lea    0x2(%r9),%eax
  movsd  %xmm0,0x8(%r10)
  movsd  0x28(%rsi),%xmm0
  movsd  %xmm0,0x28(%r10)
  cmp    %eax,%edi
  jle    L56
  movsd  0x40(%rsi),%xmm0
  add    $0x3,%r9d
  movsd  %xmm0,0x10(%r10)
  movsd  0x48(%rsi),%xmm0
  movsd  %xmm0,0x30(%r10)
  cmp    %r9d,%edi
  jle    L56
  movsd  0x60(%rsi),%xmm0
  movsd  %xmm0,0x18(%r10)
  movsd  0x68(%rsi),%xmm0
  movsd  %xmm0,0x38(%r10)
L56:
  ret
  nopl   0x0(%rax,%rax,1)
L129:
  xor    %r9d,%r9d
  lea    -0x3(%rdi),%edx
  cmp    %edx,%r9d
  jl     L13
L75:
  mov    %rax,%rsi
  jmp    L79
  nopw   0x0(%rax,%rax,1)
L54:
  mov    %r9d,%r11d
  mov    %r10,%rcx
  mov    %rax,%rdx
  mov    %r11,%rsi
  shl    $0x5,%rsi
  add    %rax,%rsi
  nopl   0x0(%rax,%rax,1)
L102:
  movsd  (%rdx),%xmm0
  add    $0x20,%rdx
  add    $0x8,%rcx
  movsd  %xmm0,-0x8(%rcx)
  movsd  -0x18(%rdx),%xmm0
  movsd  %xmm0,0x18(%rcx)
  cmp    %rdx,%rsi
  jne    L102
  jmp    L119
```

assistant: ```c
void kernel_dgetr_2_lib4(int kmax, int kna, double *A, double *C, int sdc) {
  const int bs = 4;
  int k;
  k = 0;
  if (kna > 0) {
    for (; k < kna; k++) {
      C[0 + bs * 0] = A[0 + bs * 0];
      C[0 + bs * 1] = A[1 + bs * 0];
      C += 1;
      A += bs;
    }
    C += bs * (sdc - 1);
  }
  for (; k < kmax - 3; k += 4) {
    C[0 + bs * 0] = A[0 + bs * 0];
    C[0 + bs * 1] = A[1 + bs * 0];
    C[1 + bs * 0] = A[0 + bs * 1];
    C[1 + bs * 1] = A[1 + bs * 1];
    C[2 + bs * 0] = A[0 + bs * 2];
    C[2 + bs * 1] = A[1 + bs * 2];
    C[3 + bs * 0] = A[0 + bs * 3];
    C[3 + bs * 1] = A[1 + bs * 3];
    C += bs * sdc;
    A += bs * bs;
  }
  for (; k < kmax; k++) {
    C[0 + bs * 0] = A[0 + bs * 0];
    C[0 + bs * 1] = A[1 + bs * 0];
    C += 1;
    A += bs;
  }
}
```