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

```asm
<fast_expansion_sum_zeroelim>:
  endbr64
  push   %rbx
  mov    %rsi,%r10
  mov    %edx,%esi
  mov    %rcx,%r9
  movsd  (%r10),%xmm3
  movsd  (%rcx),%xmm0
  comisd %xmm3,%xmm0
  seta   %dl
  movapd %xmm3,%xmm1
  xorpd  D20(%rip),%xmm1
  comisd %xmm1,%xmm0
  seta   %al
  cmp    %al,%dl
  je     L158
  movsd  0x8(%rcx),%xmm2
  mov    $0x1,%ecx
  mov    $0x0,%r11d
L219:
  cmp    %edi,%r11d
  jge    L41
  cmp    %esi,%ecx
  jge    L41
  comisd %xmm3,%xmm2
  seta   %dl
  movapd %xmm3,%xmm1
  xorpd  D20(%rip),%xmm1
  comisd %xmm1,%xmm2
  seta   %al
  cmp    %al,%dl
  je     L333
  movapd %xmm0,%xmm4
  addsd  %xmm2,%xmm4
  movapd %xmm4,%xmm5
  subsd  %xmm2,%xmm5
  subsd  %xmm5,%xmm0
  movapd %xmm0,%xmm1
  add    $0x1,%ecx
  movslq %ecx,%rax
  movsd  (%r9,%rax,8),%xmm2
  movapd %xmm4,%xmm0
L82:
  ucomisd D14(%rip),%xmm1
  jp     L10
  mov    $0x0,%edx
  jne    L10
L54:
  cmp    %edi,%r11d
  jge    L165
  cmp    %esi,%ecx
  jge    L165
  movq   D20(%rip),%xmm8
  pxor   %xmm7,%xmm7
  movapd %xmm7,%xmm9
  jmp    L252
L158:
  movapd %xmm0,%xmm2
  movapd %xmm3,%xmm0
  movsd  0x8(%r10),%xmm3
  mov    $0x0,%ecx
  mov    $0x1,%r11d
  jmp    L219
L333:
  movapd %xmm0,%xmm4
  addsd  %xmm3,%xmm4
  movapd %xmm4,%xmm5
  subsd  %xmm3,%xmm5
  subsd  %xmm5,%xmm0
  movapd %xmm0,%xmm1
  add    $0x1,%r11d
  movslq %r11d,%rax
  movsd  (%r10,%rax,8),%xmm3
  movapd %xmm4,%xmm0
  jmp    L82
L10:
  movsd  %xmm1,(%r8)
  mov    $0x1,%edx
  jmp    L54
L162:
  movapd %xmm0,%xmm4
  addsd  %xmm3,%xmm4
  movapd %xmm4,%xmm5
  subsd  %xmm0,%xmm5
  movapd %xmm4,%xmm6
  subsd  %xmm5,%xmm6
  subsd  %xmm6,%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm5,%xmm3
  addsd  %xmm3,%xmm1
  add    $0x1,%r11d
  movslq %r11d,%rax
  movsd  (%r10,%rax,8),%xmm3
  movapd %xmm4,%xmm0
  jmp    L277
L16:
  movslq %edx,%rax
  movsd  %xmm1,(%r8,%rax,8)
  lea    0x1(%rdx),%edx
L66:
  cmp    %edi,%r11d
  jge    L165
  cmp    %esi,%ecx
  jge    L165
L252:
  comisd %xmm3,%xmm2
  seta   %bl
  movapd %xmm3,%xmm1
  xorpd  %xmm8,%xmm1
  comisd %xmm1,%xmm2
  seta   %al
  cmp    %al,%bl
  je     L162
  movapd %xmm0,%xmm4
  addsd  %xmm2,%xmm4
  movapd %xmm4,%xmm5
  subsd  %xmm0,%xmm5
  movapd %xmm4,%xmm6
  subsd  %xmm5,%xmm6
  subsd  %xmm6,%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm5,%xmm2
  addsd  %xmm2,%xmm1
  add    $0x1,%ecx
  movslq %ecx,%rax
  movsd  (%r9,%rax,8),%xmm2
  movapd %xmm4,%xmm0
L277:
  ucomisd %xmm7,%xmm1
  jp     L16
  ucomisd %xmm9,%xmm1
  je     L66
  jmp    L16
L11:
  movslq %edx,%r10
  movsd  %xmm1,(%r8,%r10,8)
  lea    0x1(%rdx),%edx
L283:
  add    $0x8,%rax
  cmp    %rax,%rdi
  je     L152
L160:
  movapd %xmm0,%xmm1
  addsd  %xmm3,%xmm0
  movapd %xmm0,%xmm4
  subsd  %xmm1,%xmm4
  movapd %xmm0,%xmm5
  subsd  %xmm4,%xmm5
  subsd  %xmm5,%xmm1
  subsd  %xmm4,%xmm3
  addsd  %xmm3,%xmm1
  movsd  0x8(%rax),%xmm3
  ucomisd %xmm6,%xmm1
  jp     L11
  ucomisd %xmm7,%xmm1
  je     L283
  jmp    L11
L41:
  mov    $0x0,%edx
L165:
  cmp    %r11d,%edi
  jle    L152
  movslq %r11d,%rbx
  lea    (%r10,%rbx,8),%rax
  sub    $0x1,%edi
  sub    %r11d,%edi
  add    %rbx,%rdi
  lea    0x8(%r10,%rdi,8),%rdi
  pxor   %xmm6,%xmm6
  movapd %xmm6,%xmm7
  jmp    L160
L152:
  cmp    %ecx,%esi
  jle    L279
  movslq %ecx,%rdi
  lea    (%r9,%rdi,8),%rax
  sub    $0x1,%esi
  sub    %ecx,%esi
  add    %rdi,%rsi
  lea    0x8(%r9,%rsi,8),%rcx
  pxor   %xmm5,%xmm5
  movapd %xmm5,%xmm6
  jmp    L327
L262:
  movslq %edx,%rsi
  movsd  %xmm1,(%r8,%rsi,8)
  lea    0x1(%rdx),%edx
L331:
  add    $0x8,%rax
  cmp    %rcx,%rax
  je     L279
L327:
  movapd %xmm0,%xmm1
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm3
  subsd  %xmm1,%xmm3
  movapd %xmm0,%xmm4
  subsd  %xmm3,%xmm4
  subsd  %xmm4,%xmm1
  subsd  %xmm3,%xmm2
  addsd  %xmm2,%xmm1
  movsd  0x8(%rax),%xmm2
  ucomisd %xmm5,%xmm1
  jp     L262
  ucomisd %xmm6,%xmm1
  je     L331
  jmp    L262
L279:
  ucomisd D14(%rip),%xmm0
  setp   %al
  mov    $0x1,%ecx
  cmovne %ecx,%eax
  test   %al,%al
  jne    L4
  test   %edx,%edx
  je     L4
L306:
  mov    %edx,%eax
  pop    %rbx
  ret
L4:
  movslq %edx,%rax
  movsd  %xmm0,(%r8,%rax,8)
  lea    0x1(%rdx),%edx
  jmp    L306
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D20
data_type: f32[4]
name: extract_symbolic_data
arguments: data_label: D14
data_type: f32[2]

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

tool: {"D14": [0.0,0.0]}

assistant: ```c
int fast_expansion_sum_zeroelim(int elen, double *e, int flen, double *f, double *h) {
  double Q;
  double Qnew;
  double hh;
  double bvirt;
  double avirt, bround, around;
  int eindex, findex, hindex;
  double enow, fnow;
  enow = e[0];
  fnow = f[0];
  eindex = findex = 0;
  if ((fnow > enow) == (fnow > -enow)) {
    Q = enow;
    enow = e[++eindex];
  } else {
    Q = fnow;
    fnow = f[++findex];
  }
  hindex = 0;
  if ((eindex < elen) && (findex < flen)) {
    if ((fnow > enow) == (fnow > -enow)) {
      Qnew = (double)(enow + Q);
      bvirt = Qnew - enow;
      hh = Q - bvirt;
      enow = e[++eindex];
    } else {
      Qnew = (double)(fnow + Q);
      bvirt = Qnew - fnow;
      hh = Q - bvirt;
      fnow = f[++findex];
    }
    Q = Qnew;
    if (hh != 0.0) {
      h[hindex++] = hh;
    }
    while ((eindex < elen) && (findex < flen)) {
      if ((fnow > enow) == (fnow > -enow)) {
        Qnew = (double)(Q + enow);
        bvirt = (double)(Qnew - Q);
        avirt = Qnew - bvirt;
        bround = enow - bvirt;
        around = Q - avirt;
        hh = around + bround;
        enow = e[++eindex];
      } else {
        Qnew = (double)(Q + fnow);
        bvirt = (double)(Qnew - Q);
        avirt = Qnew - bvirt;
        bround = fnow - bvirt;
        around = Q - avirt;
        hh = around + bround;
        fnow = f[++findex];
      }
      Q = Qnew;
      if (hh != 0.0) {
        h[hindex++] = hh;
      }
    }
  }
  while (eindex < elen) {
    Qnew = (double)(Q + enow);
    bvirt = (double)(Qnew - Q);
    avirt = Qnew - bvirt;
    bround = enow - bvirt;
    around = Q - avirt;
    hh = around + bround;
    enow = e[++eindex];
    Q = Qnew;
    if (hh != 0.0) {
      h[hindex++] = hh;
    }
  }
  while (findex < flen) {
    Qnew = (double)(Q + fnow);
    bvirt = (double)(Qnew - Q);
    avirt = Qnew - bvirt;
    bround = fnow - bvirt;
    around = Q - avirt;
    hh = around + bround;
    fnow = f[++findex];
    Q = Qnew;
    if (hh != 0.0) {
      h[hindex++] = hh;
    }
  }
  if ((Q != 0.0) || (hindex == 0)) {
    h[hindex++] = Q;
  }
  return hindex;
}
```