user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<fast_expansion_sum>:
  endbr64
  mov    %rsi,%r9
  push   %rbx
  movsd  (%rcx),%xmm1
  mov    %rcx,%rsi
  movsd  (%r9),%xmm0
  movq   D8(%rip),%xmm4
  comisd %xmm0,%xmm1
  movapd %xmm0,%xmm2
  xorpd  %xmm4,%xmm2
  seta   %cl
  comisd %xmm2,%xmm1
  seta   %al
  cmp    %al,%cl
  je     L141
  movapd %xmm1,%xmm0
  mov    $0x1,%ecx
  xor    %r10d,%r10d
L151:
  cmp    %edi,%r10d
  setl   %r11b
  cmp    %edx,%ecx
  setl   %al
  and    %al,%r11b
  je     L92
  movslq %ecx,%rax
  movsd  (%rsi,%rax,8),%xmm1
  movslq %r10d,%rax
  movsd  (%r9,%rax,8),%xmm2
  comisd %xmm2,%xmm1
  movapd %xmm2,%xmm3
  xorpd  %xmm4,%xmm3
  seta   %bl
  comisd %xmm3,%xmm1
  seta   %al
  cmp    %al,%bl
  je     L164
  movapd %xmm1,%xmm2
  add    $0x1,%ecx
  addsd  %xmm0,%xmm2
  cmp    %ecx,%edx
  setg   %bl
  movapd %xmm2,%xmm3
  subsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm0
  movapd %xmm0,%xmm1
  movapd %xmm2,%xmm0
L137:
  movsd  %xmm1,(%r8)
  test   %r11b,%r11b
  je     L102
  mov    $0x2,%eax
  test   %bl,%bl
  jne    L165
  jmp    L102
  nopw   0x0(%rax,%rax,1)
L174:
  movapd %xmm1,%xmm2
  add    $0x1,%ecx
  addsd  %xmm0,%xmm2
  movapd %xmm2,%xmm3
  movapd %xmm2,%xmm5
  subsd  %xmm0,%xmm3
  subsd  %xmm3,%xmm5
  subsd  %xmm3,%xmm1
  subsd  %xmm5,%xmm0
  addsd  %xmm0,%xmm1
  movapd %xmm2,%xmm0
L139:
  movsd  %xmm1,-0x8(%r8,%rax,8)
  mov    %eax,%r11d
  add    $0x1,%rax
  cmp    %edi,%r10d
  jge    L150
  cmp    %edx,%ecx
  jge    L150
L165:
  movslq %ecx,%r11
  movsd  (%rsi,%r11,8),%xmm1
  movslq %r10d,%r11
  movsd  (%r9,%r11,8),%xmm2
  comisd %xmm2,%xmm1
  movapd %xmm2,%xmm3
  xorpd  %xmm4,%xmm3
  seta   %bl
  comisd %xmm3,%xmm1
  seta   %r11b
  cmp    %r11b,%bl
  jne    L174
  movapd %xmm2,%xmm3
  add    $0x1,%r10d
  addsd  %xmm0,%xmm3
  movapd %xmm3,%xmm5
  movapd %xmm3,%xmm6
  subsd  %xmm0,%xmm5
  subsd  %xmm5,%xmm6
  subsd  %xmm5,%xmm2
  subsd  %xmm6,%xmm0
  movapd %xmm0,%xmm1
  movapd %xmm3,%xmm0
  addsd  %xmm2,%xmm1
  jmp    L139
  nopw   0x0(%rax,%rax,1)
L92:
  mov    %r8,%rbx
  mov    $0x1,%eax
  xor    %r11d,%r11d
L14:
  cmp    %r10d,%edi
  jle    L87
  movslq %r10d,%rax
  movslq %r11d,%rbx
  sub    %rax,%rbx
  lea    (%r8,%rbx,8),%rbx
  nopl   (%rax)
L121:
  movsd  (%r9,%rax,8),%xmm2
  movapd %xmm0,%xmm1
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm3
  movapd %xmm0,%xmm4
  subsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm4
  subsd  %xmm3,%xmm2
  subsd  %xmm4,%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,(%rbx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%edi
  jg     L121
  add    %edi,%r11d
  sub    %r10d,%r11d
  movslq %r11d,%rax
  lea    (%r8,%rax,8),%rbx
  lea    0x1(%r11),%eax
L87:
  cmp    %ecx,%edx
  jle    L22
  movslq %ecx,%rax
  movslq %r11d,%rdi
  sub    %rax,%rdi
  lea    (%r8,%rdi,8),%rdi
  nopl   0x0(%rax)
L114:
  movsd  (%rsi,%rax,8),%xmm2
  movapd %xmm0,%xmm1
  addsd  %xmm2,%xmm0
  movapd %xmm0,%xmm3
  movapd %xmm0,%xmm4
  subsd  %xmm1,%xmm3
  subsd  %xmm3,%xmm4
  subsd  %xmm3,%xmm2
  subsd  %xmm4,%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,(%rdi,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%edx
  jg     L114
  sub    %ecx,%edx
  add    %r11d,%edx
  movslq %edx,%rax
  lea    (%r8,%rax,8),%rbx
  lea    0x1(%rdx),%eax
L22:
  movsd  %xmm0,(%rbx)
  pop    %rbx
  ret
  nopl   0x0(%rax,%rax,1)
L150:
  movslq %r11d,%rax
  lea    (%r8,%rax,8),%rbx
  lea    0x1(%r11),%eax
  jmp    L14
L141:
  xor    %ecx,%ecx
  mov    $0x1,%r10d
  jmp    L151
  nopl   (%rax)
L164:
  movapd %xmm2,%xmm3
  add    $0x1,%r10d
  mov    %r11d,%ebx
  addsd  %xmm0,%xmm3
  cmp    %r10d,%edi
  setg   %r11b
  movapd %xmm3,%xmm5
  subsd  %xmm2,%xmm5
  subsd  %xmm5,%xmm0
  movapd %xmm0,%xmm1
  movapd %xmm3,%xmm0
  jmp    L137
L102:
  lea    0x8(%r8),%rbx
  mov    $0x2,%eax
  mov    $0x1,%r11d
  jmp    L14
```

tool_calls: name: unpack_data
arguments: label: D8
type_infered: byte[8]

tool: {"D8": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

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