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

```asm
<dtrset_lib>:
  endbr64
  neg    %esi
  push   %r15
  push   %r14
  mov    %rdx,%r14
  push   %r13
  mov    %esi,%r13d
  and    $0x3,%r13d
  push   %r12
  cmp    %r13d,%edi
  push   %rbp
  push   %rbx
  cmovle %edi,%r13d
  xor    %eax,%eax
  mov    %edi,-0x10(%rsp)
  test   %r13d,%r13d
  jle    L130
L107:
  movsd  %xmm0,(%r14,%rax,8)
  test   %rax,%rax
  je     L137
  movsd  %xmm0,0x20(%r14,%rax,8)
  cmp    $0x2,%eax
  jne    L137
  movsd  %xmm0,0x40(%r14,%rax,8)
L137:
  add    $0x1,%rax
  cmp    %eax,%r13d
  jg     L107
  lea    -0x4(,%rcx,4),%eax
  movslq %r13d,%rdx
  cltq
  add    %rdx,%rax
  lea    (%r14,%rax,8),%r14
L11:
  mov    -0x10(%rsp),%eax
  lea    -0x3(%rax),%r15d
  cmp    %r15d,%r13d
  jge    L51
  lea    0x0(,%rcx,4),%eax
  mov    %r13d,-0xc(%rsp)
  mov    %r13d,%r12d
  lea    0x0(,%r13,4),%edx
  cltq
  mov    %r14,-0x8(%rsp)
  lea    0x0(,%rax,8),%rbx
  mov    %r14,%rax
  mov    %rbx,-0x18(%rsp)
  mov    %rbx,%r14
  nopl   0x0(%rax)
L70:
  test   %r12d,%r12d
  je     L68
  lea    -0x1(%r12),%esi
  mov    %rax,%rcx
  shl    $0x5,%rsi
  lea    0x20(%rax,%rsi,1),%rsi
  nopw   0x0(%rax,%rax,1)
L61:
  movsd  %xmm0,(%rcx)
  add    $0x20,%rcx
  movsd  %xmm0,-0x18(%rcx)
  movsd  %xmm0,-0x10(%rcx)
  movsd  %xmm0,-0x8(%rcx)
  cmp    %rsi,%rcx
  jne    L61
  movslq %edx,%rcx
  lea    0x2(%rdx),%esi
  lea    0x3(%rdx),%edi
  lea    (%rax,%rcx,8),%r13
  lea    0x5(%rdx),%r8d
  movslq %esi,%rsi
  movslq %edi,%rdi
  lea    0x1(%rdx),%ecx
  lea    0x6(%rdx),%r9d
  movslq %r8d,%r8
  shl    $0x3,%rsi
  lea    0x7(%rdx),%r10d
  lea    0xa(%rdx),%r11d
  movslq %ecx,%rcx
  movslq %r9d,%r9
  lea    0xb(%rdx),%ebx
  lea    0xf(%rdx),%ebp
  movslq %r10d,%r10
  movslq %r11d,%r11
  shl    $0x3,%rcx
  movslq %ebx,%rbx
  movslq %ebp,%rbp
  shl    $0x3,%rdi
  shl    $0x3,%r8
  shl    $0x3,%r9
  shl    $0x3,%r10
  shl    $0x3,%r11
  shl    $0x3,%rbx
  shl    $0x3,%rbp
L134:
  add    $0x4,%r12d
  movsd  %xmm0,0x0(%r13)
  add    $0x10,%edx
  movsd  %xmm0,(%rax,%rcx,1)
  movsd  %xmm0,(%rax,%rsi,1)
  movsd  %xmm0,(%rax,%rdi,1)
  movsd  %xmm0,(%rax,%r8,1)
  movsd  %xmm0,(%rax,%r9,1)
  movsd  %xmm0,(%rax,%r10,1)
  movsd  %xmm0,(%rax,%r11,1)
  movsd  %xmm0,(%rax,%rbx,1)
  movsd  %xmm0,(%rax,%rbp,1)
  add    %r14,%rax
  cmp    %r15d,%r12d
  jl     L70
  mov    -0x10(%rsp),%eax
  mov    -0xc(%rsp),%r13d
  mov    -0x8(%rsp),%r14
  sub    $0x4,%eax
  sub    %r13d,%eax
  shr    $0x2,%eax
  lea    0x1(%rax),%edx
  lea    0x4(%r13,%rax,4),%r13d
  imul   -0x18(%rsp),%rdx
  add    %rdx,%r14
L51:
  mov    -0x10(%rsp),%ecx
  cmp    -0x10(%rsp),%r13d
  jge    L16
L127:
  test   %r13d,%r13d
  js     L92
  movslq %r13d,%rdx
  mov    %r14,%rax
  shl    $0x5,%rdx
  lea    0x20(%r14,%rdx,1),%rdx
  nopw   0x0(%rax,%rax,1)
L27:
  movsd  %xmm0,(%rax)
  add    $0x20,%rax
  cmp    %rdx,%rax
  jne    L27
L92:
  add    $0x1,%r13d
  add    $0x8,%r14
  cmp    %r13d,%ecx
  jne    L127
L16:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L68:
  mov    %rax,%r13
  mov    $0x78,%ebp
  mov    $0x58,%ebx
  mov    $0x50,%r11d
  mov    $0x38,%r10d
  mov    $0x30,%r9d
  mov    $0x28,%r8d
  mov    $0x18,%edi
  mov    $0x10,%esi
  mov    $0x8,%ecx
  jmp    L134
L130:
  xor    %r13d,%r13d
  jmp    L11
```

assistant: ```c
void dtrset_lib(int row, double alpha, int offset, double *pA, int sda) {
  const int bs = 4;
  int ii, jj;
  int col = row;
  int na = (bs - offset % bs) % bs;
  na = row < na ? row : na;
  ii = 0;
  if (na > 0) {
    for (; ii < na; ii++) {
      for (jj = 0; jj <= ii; jj++) {
        pA[jj * bs] = alpha;
      }
      pA += 1;
    }
    pA += bs * (sda - 1);
  }
  for (; ii < row - 3; ii += 4) {
    for (jj = 0; jj < ii; jj++) {
      pA[0 + jj * bs] = alpha;
      pA[1 + jj * bs] = alpha;
      pA[2 + jj * bs] = alpha;
      pA[3 + jj * bs] = alpha;
    }
    pA[0 + (jj + 0) * bs] = alpha;
    pA[1 + (jj + 0) * bs] = alpha;
    pA[2 + (jj + 0) * bs] = alpha;
    pA[3 + (jj + 0) * bs] = alpha;
    pA[1 + (jj + 1) * bs] = alpha;
    pA[2 + (jj + 1) * bs] = alpha;
    pA[3 + (jj + 1) * bs] = alpha;
    pA[2 + (jj + 2) * bs] = alpha;
    pA[3 + (jj + 2) * bs] = alpha;
    pA[3 + (jj + 3) * bs] = alpha;
    pA += bs * sda;
  }
  for (; ii < row; ii++) {
    for (jj = 0; jj < ii + 1; jj++) {
      pA[jj * bs] = alpha;
    }
    pA += 1;
  }
}
```