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

```asm
<fast_expansion_sum_zeroelim>:
  endbr64
  mov    %rsi,%r11
  push   %rbx
  movsd  (%rcx),%xmm0
  mov    %edx,%esi
  movsd  (%r11),%xmm3
  movq   D5(%rip),%xmm7
  mov    %rcx,%r10
  comisd %xmm3,%xmm0
  movapd %xmm3,%xmm1
  xorpd  %xmm7,%xmm1
  seta   %dl
  comisd %xmm1,%xmm0
  seta   %al
  cmp    %al,%dl
  je     L192
  movsd  0x8(%rcx),%xmm1
  mov    $0x1,%edx
  xor    %ecx,%ecx
L46:
  cmp    %edi,%ecx
  setl   %r9b
  cmp    %esi,%edx
  setl   %al
  and    %r9b,%al
  je     L114
  comisd %xmm3,%xmm1
  movapd %xmm3,%xmm2
  movapd %xmm0,%xmm4
  xorpd  %xmm7,%xmm2
  seta   %bl
  comisd %xmm2,%xmm1
  seta   %r9b
  cmp    %r9b,%bl
  je     L187
  addsd  %xmm1,%xmm4
  add    $0x1,%edx
  movslq %edx,%r9
  cmp    %edx,%esi
  setg   %bl
  movapd %xmm4,%xmm5
  subsd  %xmm1,%xmm5
  movsd  (%r10,%r9,8),%xmm1
  subsd  %xmm5,%xmm0
  movapd %xmm0,%xmm2
  movapd %xmm4,%xmm0
L203:
  pxor   %xmm5,%xmm5
  ucomisd %xmm5,%xmm2
  jp     L32
  mov    $0x0,%r9d
  je     L54
L32:
  movsd  %xmm2,(%r8)
  mov    $0x1,%r9d
L54:
  test   %al,%al
  je     L181
  test   %bl,%bl
  jne    L90
  jmp    L181
  nopw   0x0(%rax,%rax,1)
L79:
  addsd  %xmm1,%xmm4
  add    $0x1,%edx
  movslq %edx,%rax
  movapd %xmm4,%xmm6
  movapd %xmm4,%xmm8
  subsd  %xmm0,%xmm6
  subsd  %xmm6,%xmm8
  subsd  %xmm6,%xmm1
  subsd  %xmm8,%xmm0
  movapd %xmm0,%xmm2
  movapd %xmm4,%xmm0
  addsd  %xmm1,%xmm2
  movsd  (%r10,%rax,8),%xmm1
L144:
  ucomisd %xmm5,%xmm2
  jp     L8
  je     L17
L8:
  movslq %r9d,%rax
  add    $0x1,%r9d
  movsd  %xmm2,(%r8,%rax,8)
L17:
  cmp    %edi,%ecx
  jge    L181
  cmp    %esi,%edx
  jge    L181
L90:
  comisd %xmm3,%xmm1
  movapd %xmm3,%xmm2
  movapd %xmm0,%xmm4
  xorpd  %xmm7,%xmm2
  seta   %bl
  comisd %xmm2,%xmm1
  seta   %al
  cmp    %al,%bl
  jne    L79
  addsd  %xmm3,%xmm4
  add    $0x1,%ecx
  movslq %ecx,%rax
  movapd %xmm4,%xmm6
  movapd %xmm4,%xmm8
  subsd  %xmm0,%xmm6
  subsd  %xmm6,%xmm8
  subsd  %xmm6,%xmm3
  subsd  %xmm8,%xmm0
  movapd %xmm0,%xmm2
  movapd %xmm4,%xmm0
  addsd  %xmm3,%xmm2
  movsd  (%r11,%rax,8),%xmm3
  jmp    L144
  nopl   0x0(%rax,%rax,1)
L114:
  xor    %r9d,%r9d
  pxor   %xmm5,%xmm5
L181:
  cmp    %ecx,%edi
  jle    L94
  sub    $0x1,%edi
  movslq %ecx,%rbx
  sub    %ecx,%edi
  lea    (%r11,%rbx,8),%rax
  add    %rbx,%rdi
  lea    0x8(%r11,%rdi,8),%rcx
  nop
L206:
  movapd %xmm0,%xmm2
  addsd  %xmm3,%xmm0
  movapd %xmm0,%xmm4
  movapd %xmm0,%xmm6
  subsd  %xmm2,%xmm4
  subsd  %xmm4,%xmm6
  subsd  %xmm4,%xmm3
  subsd  %xmm6,%xmm2
  addsd  %xmm3,%xmm2
  movsd  0x8(%rax),%xmm3
  ucomisd %xmm5,%xmm2
  jp     L170
  je     L124
L170:
  movslq %r9d,%rdi
  add    $0x1,%r9d
  movsd  %xmm2,(%r8,%rdi,8)
L124:
  add    $0x8,%rax
  cmp    %rax,%rcx
  jne    L206
L94:
  cmp    %edx,%esi
  jle    L133
  sub    $0x1,%esi
  movslq %edx,%rcx
  sub    %edx,%esi
  lea    (%r10,%rcx,8),%rax
  add    %rcx,%rsi
  lea    0x8(%r10,%rsi,8),%rdx
  nop
L225:
  movapd %xmm0,%xmm2
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm3
  movapd %xmm0,%xmm4
  subsd  %xmm2,%xmm3
  subsd  %xmm3,%xmm4
  subsd  %xmm3,%xmm1
  subsd  %xmm4,%xmm2
  addsd  %xmm1,%xmm2
  movsd  0x8(%rax),%xmm1
  ucomisd %xmm5,%xmm2
  jp     L112
  je     L23
L112:
  movslq %r9d,%rcx
  add    $0x1,%r9d
  movsd  %xmm2,(%r8,%rcx,8)
L23:
  add    $0x8,%rax
  cmp    %rax,%rdx
  jne    L225
L133:
  ucomisd %xmm5,%xmm0
  mov    $0x1,%edx
  setp   %al
  cmovne %edx,%eax
  test   %al,%al
  jne    L250
  test   %r9d,%r9d
  je     L250
  mov    %r9d,%eax
  pop    %rbx
  ret
  nopl   0x0(%rax)
L250:
  movslq %r9d,%rax
  add    $0x1,%r9d
  pop    %rbx
  movsd  %xmm0,(%r8,%rax,8)
  mov    %r9d,%eax
  ret
  nopw   0x0(%rax,%rax,1)
L192:
  movapd %xmm0,%xmm1
  xor    %edx,%edx
  movapd %xmm3,%xmm0
  mov    $0x1,%ecx
  movsd  0x8(%r11),%xmm3
  jmp    L46
  nopw   0x0(%rax,%rax,1)
L187:
  addsd  %xmm3,%xmm4
  add    $0x1,%ecx
  mov    %eax,%ebx
  movslq %ecx,%r9
  cmp    %ecx,%edi
  setg   %al
  movapd %xmm4,%xmm5
  subsd  %xmm3,%xmm5
  movsd  (%r11,%r9,8),%xmm3
  subsd  %xmm5,%xmm0
  movapd %xmm0,%xmm2
  movapd %xmm4,%xmm0
  jmp    L203
```

tool_calls: name: retrieve_data
arguments: tag: D5
datatype: f32[2]

tool: {"D5": [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;
}
```