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

```asm
<dtrset_lib>:
  endbr64
  mov    %rdx,%r8
  mov    %ecx,%r10d
  mov    %esi,%eax
  sar    $0x1f,%eax
  shr    $0x1e,%eax
  add    %eax,%esi
  and    $0x3,%esi
  sub    %esi,%eax
  lea    0x4(%rax),%eax
  cltd
  shr    $0x1e,%edx
  lea    (%rax,%rdx,1),%r9d
  and    $0x3,%r9d
  sub    %edx,%r9d
  cmp    %r9d,%edi
  cmovle %edi,%r9d
  test   %r9d,%r9d
  jle    L61
  lea    0x20(%r8),%rdx
  mov    %r8,%rsi
  mov    $0x0,%ecx
  jmp    L185
L121:
  movsd  %xmm0,(%rax)
  add    $0x20,%rax
  cmp    %rdx,%rax
  jne    L121
L178:
  add    $0x8,%rsi
  add    $0x1,%ecx
  add    $0x28,%rdx
  cmp    %ecx,%r9d
  je     L6
L185:
  mov    %rsi,%rax
  test   %ecx,%ecx
  jns    L121
  jmp    L178
L6:
  lea    -0x4(,%r10,4),%eax
  cltq
  mov    %r9d,%edx
  add    %rdx,%rax
  lea    (%r8,%rax,8),%r8
L79:
  lea    -0x3(%rdi),%r11d
  cmp    %r11d,%r9d
  jge    L151
  push   %rbp
  push   %rbx
  shl    $0x2,%r10d
  movslq %r10d,%r10
  shl    $0x3,%r10
  mov    %r9d,%esi
  mov    %r8,%rdx
  mov    $0x0,%ebx
  jmp    L175
L61:
  mov    $0x0,%r9d
  jmp    L79
L166:
  mov    %rdx,%rax
  lea    -0x1(%rsi),%ecx
  shl    $0x5,%rcx
  lea    0x20(%rdx,%rcx,1),%rcx
L62:
  movsd  %xmm0,(%rax)
  movsd  %xmm0,0x8(%rax)
  movsd  %xmm0,0x10(%rax)
  movsd  %xmm0,0x18(%rax)
  add    $0x20,%rax
  cmp    %rcx,%rax
  jne    L62
  mov    %esi,%eax
L144:
  shl    $0x2,%eax
  movslq %eax,%rcx
  movsd  %xmm0,(%rdx,%rcx,8)
  lea    0x1(%rax),%ecx
  movslq %ecx,%rcx
  lea    0x0(,%rcx,8),%rbp
  movsd  %xmm0,(%rdx,%rcx,8)
  movsd  %xmm0,0x8(%rdx,%rbp,1)
  movsd  %xmm0,0x10(%rdx,%rbp,1)
  lea    0x5(%rax),%ecx
  movslq %ecx,%rcx
  lea    0x0(,%rcx,8),%rbp
  movsd  %xmm0,(%rdx,%rcx,8)
  movsd  %xmm0,0x8(%rdx,%rbp,1)
  movsd  %xmm0,0x10(%rdx,%rbp,1)
  lea    0xa(%rax),%ecx
  movslq %ecx,%rcx
  movsd  %xmm0,(%rdx,%rcx,8)
  movsd  %xmm0,0x8(%rdx,%rcx,8)
  add    $0xf,%eax
  cltq
  movsd  %xmm0,(%rdx,%rax,8)
  add    %r10,%rdx
  add    $0x4,%esi
  cmp    %r11d,%esi
  jge    L46
L175:
  mov    %ebx,%eax
  test   %esi,%esi
  jg     L166
  jmp    L144
L46:
  lea    -0x4(%rdi),%eax
  sub    %r9d,%eax
  shr    $0x2,%eax
  lea    0x1(%rax),%edx
  mov    %edx,%edx
  imul   %r10,%rdx
  add    %rdx,%r8
  lea    0x4(%r9,%rax,4),%r9d
  cmp    %r9d,%edi
  jle    L171
L150:
  test   %r9d,%r9d
  js     L162
  mov    %r8,%rax
  mov    %r9d,%edx
  shl    $0x5,%rdx
  lea    0x20(%r8,%rdx,1),%rdx
L195:
  movsd  %xmm0,(%rax)
  add    $0x20,%rax
  cmp    %rdx,%rax
  jne    L195
L162:
  add    $0x8,%r8
  add    $0x1,%r9d
  cmp    %r9d,%edi
  jne    L150
L171:
  pop    %rbx
  pop    %rbp
  ret
L151:
  cmp    %r9d,%edi
  jle    L31
L88:
  test   %r9d,%r9d
  js     L169
  mov    %r8,%rax
  mov    %r9d,%edx
  shl    $0x5,%rdx
  lea    0x20(%r8,%rdx,1),%rdx
L114:
  movsd  %xmm0,(%rax)
  add    $0x20,%rax
  cmp    %rdx,%rax
  jne    L114
L169:
  add    $0x8,%r8
  add    $0x1,%r9d
  cmp    %r9d,%edi
  jne    L88
  ret
L31:
  ret
```

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