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

```asm
<_lnplm_array>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %edi,0x1c(%rsp)
  mov    %esi,%r14d
  mov    %esi,0x18(%rsp)
  movaps %xmm0,(%rsp)
  mov    %rdx,%rbx
  mov    %rcx,%rbp
  mov    $0x1,%eax
  call   <logq@plt>
  mov    %eax,%r12d
  test   %r14d,%r14d
  jne    L22
  movl   $0x1,0x0(%rbp)
  pxor   %xmm0,%xmm0
L36:
  movaps %xmm0,(%rbx)
  mov    0x18(%rsp),%r15d
  mov    0x1c(%rsp),%r14d
  cmp    %r14d,%r15d
  je     L26
  mov    %r12d,%edi
  call   <__floatsitf>
  movaps %xmm0,0x20(%rsp)
  mov    0x0(%rbp),%eax
  mov    %eax,0x4(%rbp)
  movdqa (%rbx),%xmm1
  call   <__addtf3>
  movaps %xmm0,(%rsp)
  lea    0x1(%r15,%r15,1),%edi
  mov    $0x0,%eax
  call   <logq@plt>
  mov    %eax,%edi
  call   <__floatsitf>
  movdqa (%rsp),%xmm1
  call   <__addtf3>
  movaps %xmm0,0x10(%rbx)
  lea    0x2(%r15),%eax
  cmp    %eax,%r14d
  jl     L26
  add    $0x8,%rbp
  mov    $0x2,%r12d
  lea    -0x1(%r15,%r15,1),%r14d
L40:
  mov    -0x8(%rbp),%r13d
  neg    %r13d
  lea    (%r14,%r12,1),%edi
  mov    $0x0,%eax
  call   <logq@plt>
  mov    %eax,%edi
  call   <__floatsitf>
  movdqa (%rbx),%xmm1
  call   <__addtf3>
  movaps %xmm0,(%rsp)
  mov    -0x4(%rbp),%r15d
  lea    (%r14,%r12,2),%edi
  mov    $0x0,%eax
  call   <logq@plt>
  mov    %eax,%edi
  call   <__floatsitf>
  movdqa 0x20(%rsp),%xmm1
  call   <__addtf3>
  movdqa 0x10(%rbx),%xmm1
  call   <__addtf3>
  mov    %rbp,%rdx
  mov    %r13d,%esi
  movdqa (%rsp),%xmm1
  mov    %r15d,%edi
  mov    $0x2,%eax
  call   <logadd_s@plt>
  mov    %eax,%edi
  call   <__floatsitf>
  movaps %xmm0,0x20(%rbx)
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <logq@plt>
  mov    %eax,%edi
  call   <__floatsitf>
  movdqa %xmm0,%xmm1
  movdqa 0x20(%rbx),%xmm0
  call   <__subtf3>
  movaps %xmm0,0x20(%rbx)
  add    $0x10,%rbx
  add    $0x4,%rbp
  add    $0x1,%r12d
  mov    0x18(%rsp),%eax
  add    %r12d,%eax
  cmp    %eax,0x1c(%rsp)
  jge    L40
L26:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L22:
  mov    0x18(%rsp),%r14d
  mov    %r14d,%eax
  shr    $0x1f,%eax
  add    %r14d,%eax
  sar    %eax
  mov    %r14d,%ecx
  shr    $0x1f,%ecx
  lea    (%r14,%rcx,1),%edx
  and    $0x1,%edx
  sub    %ecx,%edx
  add    %edx,%eax
  and    $0x1,%eax
  add    %eax,%eax
  mov    $0x1,%edx
  sub    %eax,%edx
  mov    %edx,0x0(%rbp)
  lea    -0x1(%r14,%r14,1),%edi
  mov    $0x0,%eax
  call   <ln_doublefact@plt>
  mov    %eax,%r13d
  movdqa (%rsp),%xmm1
  movdqa %xmm1,%xmm0
  call   <__multf3>
  movdqa D13(%rip),%xmm1
  call   <__subtf3>
  mov    $0x1,%eax
  call   <logq@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %r14d,%xmm0
  mulsd  D11(%rip),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  mulsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sd %r13d,%xmm1
  addsd  %xmm1,%xmm0
  call   <__extenddftf2>
  jmp    L36
```

tool_calls: name: read_data
arguments: data_tag: D13
type_guessed: byte[16]
name: read_data
arguments: data_tag: D11
type_guessed: f64

tool: {"D13": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xff,0x3f]}

tool: D11: 0.5

assistant: ```c
void _lnplm_array(int lmax, int m, __float128 x, __float128 lnplm[], int sign[]) {
  int l;
  __float128 logx = logq(x);
  if (m == 0) {
    sign[0] = +1;
    lnplm[0] = 0;
  } else {
    sign[0] = (1 - 2 * (((int)(m / 2) + m % 2) & 1));
    lnplm[0] = ln_doublefact(2 * m - 1) + m * 0.5 * logq(((x) * (x)) - 1);
  }
  if (lmax == m)
    return;
  sign[1] = sign[0];
  lnplm[1] = lnplm[0] + logx + logq(2 * m + 1);
  for (l = m + 2; l <= lmax; l++) {
    lnplm[l - m] = logadd_s(logq(2 * l - 1) + logx + lnplm[l - m - 1], sign[l - m - 1], logq(l + m - 1) + lnplm[l - m - 2], -sign[l - m - 2], &sign[l - m]);
    lnplm[l - m] -= logq(l - m);
  }
}
```