user: Translate the below assembly code into C programming language:

```asm
<fast_expansion_sum_zeroelim>:
  endbr64
  mov    %rsi,%r11
  push   %r14
  movq   D5(%rip),%xmm7
  mov    %edx,%esi
  push   %r13
  movsd  (%r11),%xmm3
  mov    %rcx,%r10
  push   %r12
  push   %rbp
  movapd %xmm3,%xmm1
  push   %rbx
  movsd  (%rcx),%xmm0
  xorpd  %xmm7,%xmm1
  comisd %xmm3,%xmm0
  seta   %dl
  comisd %xmm1,%xmm0
  seta   %al
  cmp    %al,%dl
  je     L133
  movsd  0x8(%rcx),%xmm1
  mov    $0x10,%ebp
  xor    %ecx,%ecx
  mov    $0x2,%ebx
  mov    $0x8,%r13d
  mov    $0x1,%r12d
  mov    $0x1,%edx
L146:
  cmp    %edi,%ecx
  setl   %r9b
  cmp    %esi,%edx
  setl   %al
  and    %r9b,%al
  je     L154
  comisd %xmm3,%xmm1
  movapd %xmm3,%xmm2
  movapd %xmm0,%xmm4
  xorpd  %xmm7,%xmm2
  seta   %r14b
  comisd %xmm2,%xmm1
  seta   %r9b
  cmp    %r9b,%r14b
  je     L110
  addsd  %xmm1,%xmm4
  cmp    %ebx,%esi
  mov    %ebx,%edx
  movapd %xmm4,%xmm5
  subsd  %xmm1,%xmm5
  movsd  (%r10,%rbp,1),%xmm1
  setg   %bpl
  subsd  %xmm5,%xmm0
  movapd %xmm0,%xmm2
  movapd %xmm4,%xmm0
L183:
  pxor   %xmm5,%xmm5
  ucomisd %xmm5,%xmm2
  jp     L135
  mov    $0x0,%r9d
  je     L206
L135:
  movsd  %xmm2,(%r8)
  mov    $0x1,%r9d
L206:
  test   %al,%al
  je     L92
  test   %bpl,%bpl
  jne    L111
  jmp    L92
  nopl   0x0(%rax,%rax,1)
L147:
  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
L229:
  ucomisd %xmm5,%xmm2
  jp     L100
  je     L142
L100:
  movslq %r9d,%rax
  add    $0x1,%r9d
  movsd  %xmm2,(%r8,%rax,8)
L142:
  cmp    %edi,%ecx
  jge    L92
  cmp    %esi,%edx
  jge    L92
L111:
  comisd %xmm3,%xmm1
  movapd %xmm3,%xmm2
  movapd %xmm0,%xmm4
  xorpd  %xmm7,%xmm2
  seta   %bl
  comisd %xmm2,%xmm1
  seta   %al
  cmp    %al,%bl
  jne    L147
  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    L229
  nopl   0x0(%rax,%rax,1)
L154:
  xor    %r9d,%r9d
  pxor   %xmm5,%xmm5
L92:
  cmp    %ecx,%edi
  jle    L47
  sub    $0x1,%edi
  movslq %ecx,%rbx
  sub    %ecx,%edi
  lea    (%r11,%rbx,8),%rax
  add    %rbx,%rdi
  lea    0x8(%r11,%rdi,8),%rcx
  nopw   0x0(%rax,%rax,1)
L138:
  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     L37
  je     L213
L37:
  movslq %r9d,%rdi
  add    $0x1,%r9d
  movsd  %xmm2,(%r8,%rdi,8)
L213:
  add    $0x8,%rax
  cmp    %rax,%rcx
  jne    L138
L47:
  cmp    %edx,%esi
  jle    L63
  sub    $0x1,%esi
  movslq %edx,%rcx
  sub    %edx,%esi
  lea    (%r10,%rcx,8),%rax
  add    %rcx,%rsi
  lea    0x8(%r10,%rsi,8),%rdx
  nop
L26:
  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     L254
  je     L197
L254:
  movslq %r9d,%rcx
  add    $0x1,%r9d
  movsd  %xmm2,(%r8,%rcx,8)
L197:
  add    $0x8,%rax
  cmp    %rax,%rdx
  jne    L26
L63:
  ucomisd %xmm5,%xmm0
  mov    $0x1,%edx
  setp   %al
  cmovne %edx,%eax
  test   %al,%al
  jne    L61
  test   %r9d,%r9d
  je     L61
  pop    %rbx
  mov    %r9d,%eax
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  nopl   0x0(%rax,%rax,1)
L61:
  movslq %r9d,%rax
  add    $0x1,%r9d
  pop    %rbx
  pop    %rbp
  movsd  %xmm0,(%r8,%rax,8)
  pop    %r12
  mov    %r9d,%eax
  pop    %r13
  pop    %r14
  ret
  nopl   0x0(%rax)
L133:
  movapd %xmm0,%xmm1
  mov    $0x8,%ebp
  movapd %xmm3,%xmm0
  xor    %edx,%edx
  movsd  0x8(%r11),%xmm3
  mov    $0x1,%ebx
  mov    $0x1,%ecx
  mov    $0x10,%r13d
  mov    $0x2,%r12d
  jmp    L146
L110:
  addsd  %xmm3,%xmm4
  cmp    %r12d,%edi
  mov    %eax,%ebp
  mov    %r12d,%ecx
  setg   %al
  movapd %xmm4,%xmm5
  subsd  %xmm3,%xmm5
  movsd  (%r11,%r13,1),%xmm3
  subsd  %xmm5,%xmm0
  movapd %xmm0,%xmm2
  movapd %xmm4,%xmm0
  jmp    L183
```

tool_calls: name: disassemble_data
arguments: data_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;
}
```