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

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

tool_calls: name: read_data
arguments: label: D18
data_type: byte[16]
name: read_data
arguments: label: D6
data_type: f64

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

tool: D6: 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);
  }
}
```