user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<cpoly_cauchy>:
  endbr64
  push   %r14
  push   %r13
  lea    -0x1(%rdi),%r13d
  push   %r12
  mov    %edi,%r12d
  movslq %r13d,%rdi
  push   %rbp
  lea    0x0(,%rdi,8),%r14
  mov    %rsi,%rbp
  push   %rbx
  lea    (%rsi,%r14,1),%rax
  mov    %rdx,%rbx
  sub    $0x20,%rsp
  movsd  (%rax),%xmm5
  movapd %xmm5,%xmm0
  xorpd  D11(%rip),%xmm0
  movsd  %xmm5,0x18(%rsp)
  movsd  %xmm0,(%rax)
  movapd %xmm5,%xmm0
  call   <log@plt>
  movsd  0x0(%rbp),%xmm1
  movsd  %xmm0,0x10(%rsp)
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x8(%rsp)
  call   <log@plt>
  movsd  0x10(%rsp),%xmm2
  subsd  %xmm0,%xmm2
  pxor   %xmm0,%xmm0
  cvtsi2sd %r13d,%xmm0
  divsd  %xmm0,%xmm2
  movapd %xmm2,%xmm0
  call   <exp@plt>
  pxor   %xmm4,%xmm4
  movsd  0x8(%rsp),%xmm1
  movsd  0x18(%rsp),%xmm5
  movapd %xmm0,%xmm3
  movsd  -0x8(%rbp,%r14,1),%xmm0
  ucomisd %xmm4,%xmm0
  jp     L63
  jne    L63
L28:
  lea    -0x2(%r12),%ecx
  movsd  D13(%rip),%xmm5
  lea    0x10(%rbp,%rcx,8),%rdx
  cs nopw 0x0(%rax,%rax,1)
L93:
  movapd %xmm3,%xmm0
  mulsd  %xmm5,%xmm3
  cmp    $0x1,%r12d
  jle    L69
  lea    0x8(%rbp),%rax
  movapd %xmm1,%xmm2
  nopw   0x0(%rax,%rax,1)
L137:
  mulsd  %xmm3,%xmm2
  add    $0x8,%rax
  addsd  -0x8(%rax),%xmm2
  cmp    %rax,%rdx
  jne    L137
  comisd %xmm2,%xmm4
  jb     L93
L50:
  movapd %xmm0,%xmm2
  movq   D24(%rip),%xmm3
  movsd  D23(%rip),%xmm4
  divsd  %xmm0,%xmm2
  andpd  %xmm3,%xmm2
  comisd %xmm4,%xmm2
  jbe    L108
  lea    -0x3(%r12),%eax
  lea    (%rbx,%r14,1),%rdi
  lea    0x10(%rbx,%rax,8),%rsi
  lea    0x1(%rcx),%rdx
L156:
  movsd  %xmm1,(%rbx)
  cmp    $0x1,%r12d
  jle    L7
  mov    $0x1,%eax
  jmp    L139
  nopl   0x0(%rax)
L64:
  add    $0x1,%rax
L139:
  mulsd  %xmm0,%xmm1
  addsd  0x0(%rbp,%rax,8),%xmm1
  movsd  %xmm1,(%rbx,%rax,8)
  cmp    %rax,%rdx
  jne    L64
  movsd  (%rbx),%xmm1
L7:
  movsd  (%rdi),%xmm2
  cmp    $0x1,%r13d
  jle    L22
  lea    0x8(%rbx),%rax
  nopw   0x0(%rax,%rax,1)
L98:
  mulsd  %xmm0,%xmm1
  add    $0x8,%rax
  addsd  -0x8(%rax),%xmm1
  cmp    %rax,%rsi
  jne    L98
L22:
  divsd  %xmm1,%xmm2
  subsd  %xmm2,%xmm0
  divsd  %xmm0,%xmm2
  andpd  %xmm3,%xmm2
  comisd %xmm4,%xmm2
  jbe    L108
  movsd  0x0(%rbp),%xmm1
  jmp    L156
  nopw   0x0(%rax,%rax,1)
L108:
  add    $0x20,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   (%rax)
L69:
  movapd %xmm1,%xmm2
  comisd %xmm2,%xmm4
  jb     L93
  jmp    L50
L63:
  divsd  %xmm0,%xmm5
  minsd  %xmm3,%xmm5
  movapd %xmm5,%xmm3
  jmp    L28
```

tool_calls: name: parse_data
arguments: label: D11
type_infered: f32[4]
name: parse_data
arguments: label: D13
type_infered: f64
name: parse_data
arguments: label: D24
type_infered: qword
name: parse_data
arguments: label: D23
type_infered: f64

tool: {"D11": [0.0,-0.0,0.0,0.0]}

tool: D13: 0.1

tool: D24: 9223372036854775807

tool: D23: 0.005

assistant: ```c
double cpoly_cauchy(int n, double *pot, double *q) {
  double f, x, delf, dx, xm;
  int i, n1 = n - 1;
  pot[n1] = -pot[n1];
  x = exp((log(-pot[n1]) - log(pot[0])) / (double)n1);
  if (pot[n1 - 1] != 0.) {
    xm = -pot[n1] / pot[n1 - 1];
    if (xm < x)
      x = xm;
  }
  for (;;) {
    xm = x * 0.1;
    f = pot[0];
    for (i = 1; i < n; i++)
      f = f * xm + pot[i];
    if (f <= 0.0) {
      break;
    }
    x = xm;
  }
  dx = x;
  while (fabs(dx / x) > 0.005) {
    q[0] = pot[0];
    for (i = 1; i < n; i++)
      q[i] = q[i - 1] * x + pot[i];
    f = q[n1];
    delf = q[0];
    for (i = 1; i < n1; i++)
      delf = delf * x + q[i];
    dx = f / delf;
    x -= dx;
  }
  return x;
}
```