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

```asm
<rpmva>:
  endbr64
  push   %r15
  pxor   %xmm0,%xmm0
  push   %r14
  cvtsi2sd %esi,%xmm0
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x68,%rsp
  subsd  D38(%rip),%xmm0
  mov    %esi,0x2c(%rsp)
  mov    %rdx,0x30(%rsp)
  mov    %ecx,0x58(%rsp)
  mov    %r8,0x38(%rsp)
  call   <log@plt>
  movsd  D0(%rip),%xmm3
  divsd  D6(%rip),%xmm0
  addsd  %xmm3,%xmm0
  cvttsd2si %xmm0,%r13d
  test   %r13d,%r13d
  jle    L14
  mov    $0x1,%edx
  xor    %eax,%eax
  nopl   (%rax)
L110:
  add    $0x1,%eax
  add    %edx,%edx
  cmp    %eax,%r13d
  jne    L110
  mov    %edx,0x40(%rsp)
  sar    %edx
  mov    %edx,%r12d
L112:
  movslq 0x2c(%rsp),%rbp
  shl    $0x3,%rbp
  movsd  -0x8(%rbx,%rbp,1),%xmm2
  movapd %xmm2,%xmm0
  addsd  %xmm3,%xmm0
  ucomisd %xmm3,%xmm0
  jp     L170
  je     L283
L170:
  movslq 0x40(%rsp),%rdi
  movsd  %xmm2,(%rsp)
  shl    $0x3,%rdi
  call   <malloc@plt>
  movslq %r12d,%rdi
  shl    $0x3,%rdi
  mov    %rax,%r14
  call   <malloc@plt>
  mov    0x2c(%rsp),%edx
  movsd  (%rsp),%xmm2
  mov    %rax,0x50(%rsp)
  mov    D0(%rip),%rax
  test   %edx,%edx
  movq   %rax,%xmm3
  jle    L253
  mov    0x2c(%rsp),%eax
  cmp    $0x1,%eax
  je     L156
  shr    %eax
  movapd %xmm2,%xmm1
  mov    %eax,%edx
  unpcklpd %xmm1,%xmm1
  xor    %eax,%eax
  shl    $0x4,%rdx
L187:
  movupd (%rbx,%rax,1),%xmm0
  divpd  %xmm1,%xmm0
  movups %xmm0,(%r14,%rax,1)
  add    $0x10,%rax
  cmp    %rdx,%rax
  jne    L187
  mov    0x2c(%rsp),%edi
  mov    %edi,%eax
  and    $0xfffffffe,%eax
  and    $0x1,%edi
  je     L253
L310:
  cltq
  movsd  (%rbx,%rax,8),%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%r14,%rax,8)
L253:
  mov    0x40(%rsp),%ebx
  cmp    %ebx,0x2c(%rsp)
  jge    L176
  mov    %ebx,%eax
  sub    0x2c(%rsp),%eax
  lea    (%r14,%rbp,1),%rdi
  xor    %esi,%esi
  sub    $0x1,%eax
  movsd  %xmm2,(%rsp)
  lea    0x8(,%rax,8),%rdx
  call   <memset@plt>
  mov    D0(%rip),%rax
  movsd  (%rsp),%xmm2
  movq   %rax,%xmm3
L176:
  lea    -0x8(%r14),%rax
  movl   $0x1,0x28(%rsp)
  mov    %r12d,%ebp
  mov    $0x1,%ebx
  mov    %rax,0x48(%rsp)
  lea    -0x18(%r14),%rax
  mov    %rax,0x18(%rsp)
  cmp    $0x1,%r13d
  jle    L113
  mov    %r13d,0x44(%rsp)
  mov    %ebx,%r11d
  mov    %r12d,0x5c(%rsp)
  nopl   0x0(%rax,%rax,1)
L182:
  lea    -0x1(%rbp),%r13d
  test   %r11d,%r11d
  jle    L246
  movslq %ebp,%r10
  lea    -0x2(%rbp),%ebx
  lea    0x0(%rbp,%rbp,1),%r12d
  movslq %r13d,%rcx
  mov    %r10,%rax
  mov    %ebx,0xc(%rsp)
  movslq %r12d,%r12
  neg    %rax
  shl    $0x3,%rax
  cmp    $0xfffffffffffffff8,%rax
  setne  %sil
  cmp    $0x3,%ebp
  setne  %dl
  xor    %r9d,%r9d
  and    %edx,%esi
  mov    %r14,%rdx
  sub    %rax,%rdx
  mov    %sil,(%rsp)
  mov    %ebx,%eax
  mov    %rbx,%rsi
  and    $0xfffffffe,%esi
  shr    %eax
  not    %rbx
  lea    -0x1(%rax),%edi
  mov    %esi,0x8(%rsp)
  lea    0x2(%rsi),%r8d
  mov    0x48(%rsp),%rsi
  lea    -0x3(%rbp),%eax
  not    %rdi
  shl    $0x3,%rbx
  shl    $0x3,%rax
  shl    $0x4,%rdi
  sub    %rax,%rsi
  lea    -0x18(%rdx),%rax
  mov    %rsi,0x20(%rsp)
  mov    %rax,0x10(%rsp)
  nopl   0x0(%rax)
L53:
  movsd  (%r14,%rcx,8),%xmm0
  mov    %ecx,%eax
  subsd  %xmm3,%xmm0
  movsd  %xmm0,(%r14,%rcx,8)
  cmp    $0x1,%ebp
  jle    L199
  lea    -0x1(%rcx),%edx
  movslq %edx,%rdx
  lea    (%r14,%rdx,8),%rsi
  lea    0x0(%r13,%rcx,1),%edx
  movslq %edx,%rdx
  movsd  (%rsi),%xmm1
  mulsd  (%r14,%rdx,8),%xmm0
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%rsi)
  cmp    $0x2,%ebp
  je     L199
  lea    0x0(,%rcx,8),%rdx
  lea    (%rbx,%rdx,1),%r15
  lea    0x8(%rdx),%rsi
  cmp    %rsi,%r15
  lea    -0x8(%rdx),%rsi
  setge  %r15b
  cmp    %rsi,%rdx
  setge  %sil
  or     %r15b,%sil
  je     L26
  cmpb   $0x0,(%rsp)
  je     L26
  movsd  (%r14,%rcx,8),%xmm4
  mov    0x18(%rsp),%rsi
  mov    0x10(%rsp),%r15
  add    %rdx,%rsi
  unpcklpd %xmm4,%xmm4
  add    %rdx,%r15
  xor    %edx,%edx
  nopl   0x0(%rax)
L237:
  movupd (%r15,%rdx,1),%xmm1
  movupd (%rsi,%rdx,1),%xmm0
  shufpd $0x1,%xmm1,%xmm1
  mulpd  %xmm4,%xmm1
  shufpd $0x1,%xmm0,%xmm0
  subpd  %xmm1,%xmm0
  shufpd $0x1,%xmm0,%xmm0
  movups %xmm0,(%rsi,%rdx,1)
  sub    $0x10,%rdx
  cmp    %rdi,%rdx
  jne    L237
  mov    0xc(%rsp),%edx
  cmp    %edx,0x8(%rsp)
  je     L199
  sub    %r8d,%eax
  movslq %eax,%rdx
  add    %ebp,%eax
  cltq
  lea    (%r14,%rdx,8),%rdx
  movsd  (%r14,%rax,8),%xmm1
  mulsd  (%r14,%rcx,8),%xmm1
  movsd  (%rdx),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx)
L199:
  add    $0x1,%r9d
  add    %r12,%rcx
  cmp    %r11d,%r9d
  jne    L53
L246:
  addl   $0x1,0x28(%rsp)
  sar    %ebp
  mov    0x28(%rsp),%eax
  add    %r11d,%r11d
  cmp    %eax,0x44(%rsp)
  jne    L182
  mov    0x44(%rsp),%r13d
  mov    0x5c(%rsp),%r12d
L113:
  mov    0x58(%rsp),%eax
  test   %eax,%eax
  jle    L259
  movslq 0x58(%rsp),%rax
  lea    -0x1(%r12),%r15d
  xor    %ebp,%ebp
  mov    %r15d,0x8(%rsp)
  mov    %rax,(%rsp)
  mov    %r15d,%eax
  shr    %eax
  lea    -0x1(%rax),%ebx
  mov    %r15d,%eax
  mov    0x50(%rsp),%r15
  and    $0xfffffffe,%eax
  add    $0x1,%rbx
  mov    %eax,0xc(%rsp)
  shl    $0x4,%rbx
  nopl   0x0(%rax)
L254:
  mov    0x30(%rsp),%rax
  movsd  (%rax,%rbp,8),%xmm3
  test   %r12d,%r12d
  jle    L266
  cmpl   $0x1,0x8(%rsp)
  jbe    L214
  movapd %xmm3,%xmm1
  xor    %eax,%eax
  unpcklpd %xmm1,%xmm1
  nopl   0x0(%rax)
L93:
  movupd (%r14,%rax,2),%xmm0
  movhpd 0x10(%r14,%rax,2),%xmm0
  addpd  %xmm1,%xmm0
  movups %xmm0,(%r15,%rax,1)
  add    $0x10,%rax
  cmp    %rbx,%rax
  jne    L93
  mov    0xc(%rsp),%eax
L178:
  lea    (%rax,%rax,1),%edx
  movslq %eax,%rcx
  add    $0x1,%eax
  movslq %edx,%rdx
  lea    0x0(,%rcx,8),%rsi
  movsd  (%r14,%rdx,8),%xmm0
  lea    0x0(,%rdx,8),%rdi
  addsd  %xmm3,%xmm0
  movsd  %xmm0,(%r15,%rcx,8)
  cmp    %r12d,%eax
  jge    L266
  movsd  0x10(%r14,%rdi,1),%xmm0
  addsd  %xmm3,%xmm0
  movsd  %xmm0,0x8(%r15,%rsi,1)
L266:
  movapd %xmm3,%xmm1
  mov    $0x1,%r9d
  mov    $0x2,%esi
  mov    $0x1,%edi
  cmp    $0x1,%r13d
  jle    L289
L101:
  mulsd  %xmm1,%xmm1
  test   %r12d,%r12d
  jle    L118
  lea    (%rsi,%rsi,1),%r8d
  lea    (%r9,%r9,1),%edx
  movslq %esi,%r10
  mov    %r15,%rax
  movslq %edx,%rdx
  movslq %r8d,%r8
  shl    $0x3,%r10
  movslq %r9d,%r11
  lea    -0x8(%r14,%rdx,8),%rcx
  shl    $0x3,%r8
  xor    %edx,%edx
  nop
L275:
  movsd  (%rcx),%xmm0
  add    %esi,%edx
  add    %r8,%rcx
  addsd  %xmm1,%xmm0
  mulsd  (%rax,%r11,8),%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
  add    %r10,%rax
  cmp    %r12d,%edx
  jl     L275
L118:
  add    $0x1,%edi
  add    %r9d,%r9d
  add    %esi,%esi
  cmp    %edi,%r13d
  jne    L101
L289:
  movsd  (%r15),%xmm0
  mov    0x40(%rsp),%edi
  cmp    %edi,0x2c(%rsp)
  je     L309
  mulsd  %xmm1,%xmm1
  divsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%r15)
L309:
  mulsd  %xmm2,%xmm0
  mov    0x38(%rsp),%rax
  movsd  %xmm0,(%rax,%rbp,8)
  add    $0x1,%rbp
  cmp    %rbp,(%rsp)
  jne    L254
L259:
  mov    %r14,%rdi
  call   <free@plt>
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  add    $0x68,%rsp
  mov    $0x1,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L26:
  lea    (%r14,%rdx,1),%rax
  add    0x20(%rsp),%rdx
  nopl   0x0(%rax)
L129:
  movsd  -0x10(%rax,%r10,8),%xmm1
  movsd  -0x10(%rax),%xmm0
  sub    $0x8,%rax
  mulsd  (%r14,%rcx,8),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rdx,%rax
  jne    L129
  jmp    L199
L214:
  xor    %eax,%eax
  jmp    L178
L283:
  lea    D53(%rip),%rdi
  call   <puts@plt>
  add    $0x68,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L14:
  movl   $0x1,0x40(%rsp)
  xor    %r12d,%r12d
  jmp    L112
L156:
  xor    %eax,%eax
  jmp    L310
```

tool_calls: name: parse_data
arguments: tag: D38
type_guessed: f64
name: parse_data
arguments: tag: D0
type_guessed: f64
name: parse_data
arguments: tag: D6
type_guessed: byte[8]
name: parse_data
arguments: tag: D53
type_guessed: string

tool: D38: 0.5

tool: D0: 1

tool: {"D6": [0xef,0x39,0xfa,0xfe,0x42,0x2e,0xe6,0x3f]}

tool: D53: (rpmva)coefficient too small,please check.

assistant: ```c
int rpmva(a, n, x, m, y)
int n, m;
double *a, *x, *y;
{
  int i, j, k, nn, nk, kk, l, count, tail;
  double *ax, *b;
  double a1, t1, t2;
  k = log(n - 0.5) / log(2.0) + 1;
  nn = 1;
  for (i = 0; i < k; i++)
    nn = nn << 1;
  nk = nn >> 1;
  a1 = a[n - 1];
  if (a1 + 1.0 == 1.0) {
    printf("(rpmva)coefficient too small,please check.\n");
    return (0);
  }
  b = (double *)malloc(nn * sizeof(double));
  ax = (double *)malloc(nk * sizeof(double));
  for (i = 0; i < n; i++)
    b[i] = a[i] / a1;
  for (i = n; i < nn; i++)
    b[i] = 0.0;
  l = nn >> 1;
  count = 1;
  for (i = 1; i < k; i++) {
    tail = l - 1;
    for (j = 0; j < count; j++) {
      b[tail] = b[tail] - 1.0;
      for (kk = 1; kk < l; kk++)
        b[tail - kk] = b[tail - kk] - b[tail] * b[tail - kk + l];
      tail = tail + 2 * l;
    }
    l = l >> 1;
    count = count << 1;
  }
  for (kk = 0; kk < m; kk++) {
    t1 = x[kk];
    for (i = 0; i < nk; i++)
      ax[i] = t1 + b[2 * i];
    tail = 1;
    l = 2;
    t2 = t1;
    for (i = 1; i < k; i++) {
      t2 = t2 * t2;
      for (j = 0; j < nk; j = j + l)
        ax[j] = ax[j] + ax[j + tail] * (t2 + b[2 * (j + tail) - 1]);
      tail = tail << 1;
      l = l << 1;
    }
    t2 = t2 * t2 / t1;
    if (nn != n)
      ax[0] -= t2;
    y[kk] = ax[0] * a1;
  }
  free(b);
  free(ax);
  return (1);
}
```