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;
}


x86 assembly (clang, -O0):
<fast_expansion_sum_zeroelim>:
push   %rbp
mov    %rsp,%rbp
push   %rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
movl   $0x0,-0x68(%rbp)
movl   $0x0,-0x64(%rbp)
movsd  -0x80(%rbp),%xmm0
ucomisd -0x78(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x80(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    11b0 <fast_expansion_sum_zeroelim+0xb0>
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x64(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    11d4 <fast_expansion_sum_zeroelim+0xd4>
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x68(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    14ce <fast_expansion_sum_zeroelim+0x3ce>
mov    -0x68(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    14ce <fast_expansion_sum_zeroelim+0x3ce>
movsd  -0x80(%rbp),%xmm0
ucomisd -0x78(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x80(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    1287 <fast_expansion_sum_zeroelim+0x187>
movsd  -0x78(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x64(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    12ce <fast_expansion_sum_zeroelim+0x1ce>
movsd  -0x80(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x68(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12f5 <fast_expansion_sum_zeroelim+0x1f5>
jp     12f5 <fast_expansion_sum_zeroelim+0x1f5>
jmp    1311 <fast_expansion_sum_zeroelim+0x211>
movsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x6c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1316 <fast_expansion_sum_zeroelim+0x216>
mov    -0x64(%rbp),%ecx
xor    %eax,%eax
cmp    -0x4(%rbp),%ecx
mov    %al,-0x81(%rbp)
jge    1339 <fast_expansion_sum_zeroelim+0x239>
mov    -0x68(%rbp),%eax
cmp    -0x14(%rbp),%eax
setl   %al
mov    %al,-0x81(%rbp)
mov    -0x81(%rbp),%al
test   $0x1,%al
jne    134c <fast_expansion_sum_zeroelim+0x24c>
jmp    14c9 <fast_expansion_sum_zeroelim+0x3c9>
movsd  -0x80(%rbp),%xmm0
ucomisd -0x78(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x80(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    140d <fast_expansion_sum_zeroelim+0x30d>
movsd  -0x30(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x78(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x64(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    1481 <fast_expansion_sum_zeroelim+0x381>
movsd  -0x30(%rbp),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x80(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x68(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    14a8 <fast_expansion_sum_zeroelim+0x3a8>
jp     14a8 <fast_expansion_sum_zeroelim+0x3a8>
jmp    14c4 <fast_expansion_sum_zeroelim+0x3c4>
movsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x6c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1316 <fast_expansion_sum_zeroelim+0x216>
jmp    14ce <fast_expansion_sum_zeroelim+0x3ce>
jmp    14d3 <fast_expansion_sum_zeroelim+0x3d3>
mov    -0x64(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    159b <fast_expansion_sum_zeroelim+0x49b>
movsd  -0x30(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x78(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x64(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    157a <fast_expansion_sum_zeroelim+0x47a>
jp     157a <fast_expansion_sum_zeroelim+0x47a>
jmp    1596 <fast_expansion_sum_zeroelim+0x496>
movsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x6c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    14d3 <fast_expansion_sum_zeroelim+0x3d3>
jmp    15a0 <fast_expansion_sum_zeroelim+0x4a0>
mov    -0x68(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1668 <fast_expansion_sum_zeroelim+0x568>
movsd  -0x30(%rbp),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x80(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x68(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1647 <fast_expansion_sum_zeroelim+0x547>
jp     1647 <fast_expansion_sum_zeroelim+0x547>
jmp    1663 <fast_expansion_sum_zeroelim+0x563>
movsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x6c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    15a0 <fast_expansion_sum_zeroelim+0x4a0>
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    168a <fast_expansion_sum_zeroelim+0x58a>
jp     168a <fast_expansion_sum_zeroelim+0x58a>
cmpl   $0x0,-0x6c(%rbp)
jne    16a6 <fast_expansion_sum_zeroelim+0x5a6>
movsd  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x6c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x8,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fast_expansion_sum_zeroelim>:
push   %rbx
movsd  (%rsi),%xmm2
movsd  (%rcx),%xmm0
ucomisd %xmm2,%xmm0
setbe  %r9b
movapd 0xee7(%rip),%xmm1
xorpd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
cmp    %al,%r9b
je     113d <fast_expansion_sum_zeroelim+0x3d>
movapd %xmm2,%xmm3
movsd  0x8(%rsi),%xmm2
xor    %r9d,%r9d
mov    $0x1,%r10d
jmp    114f <fast_expansion_sum_zeroelim+0x4f>
movapd %xmm0,%xmm3
movsd  0x8(%rcx),%xmm0
mov    $0x1,%r9d
xor    %r10d,%r10d
xor    %eax,%eax
cmp    %edi,%r10d
jge    11a1 <fast_expansion_sum_zeroelim+0xa1>
cmp    %edx,%r9d
jge    11a1 <fast_expansion_sum_zeroelim+0xa1>
ucomisd %xmm2,%xmm0
setbe  %r11b
movapd 0xe95(%rip),%xmm1
xorpd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
movapd %xmm3,%xmm1
cmp    %al,%r11b
je     1289 <fast_expansion_sum_zeroelim+0x189>
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm4
subsd  %xmm2,%xmm4
mov    %r10d,%eax
add    $0x1,%r10d
movsd  0x8(%rsi,%rax,8),%xmm2
jmp    12a2 <fast_expansion_sum_zeroelim+0x1a2>
movapd %xmm3,%xmm1
cmp    %edi,%r10d
jge    120a <fast_expansion_sum_zeroelim+0x10a>
movslq %r10d,%r10
movslq %edi,%rdi
xorpd  %xmm3,%xmm3
jmp    11c5 <fast_expansion_sum_zeroelim+0xc5>
cs nopw 0x0(%rax,%rax,1)
cmp    %r10,%rdi
je     120a <fast_expansion_sum_zeroelim+0x10a>
movapd %xmm1,%xmm4
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm5
subsd  %xmm4,%xmm5
movapd %xmm1,%xmm6
subsd  %xmm5,%xmm6
subsd  %xmm5,%xmm2
subsd  %xmm6,%xmm4
addsd  %xmm2,%xmm4
movsd  0x8(%rsi,%r10,8),%xmm2
add    $0x1,%r10
ucomisd %xmm3,%xmm4
jne    11fc <fast_expansion_sum_zeroelim+0xfc>
jnp    11c0 <fast_expansion_sum_zeroelim+0xc0>
movslq %eax,%rbx
add    $0x1,%eax
movsd  %xmm4,(%r8,%rbx,8)
jmp    11c0 <fast_expansion_sum_zeroelim+0xc0>
cmp    %edx,%r9d
jge    1269 <fast_expansion_sum_zeroelim+0x169>
movslq %r9d,%rsi
movslq %edx,%rdx
xorpd  %xmm2,%xmm2
jmp    1225 <fast_expansion_sum_zeroelim+0x125>
nopl   0x0(%rax,%rax,1)
cmp    %rsi,%rdx
je     1269 <fast_expansion_sum_zeroelim+0x169>
movapd %xmm1,%xmm3
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm4
subsd  %xmm3,%xmm4
movapd %xmm1,%xmm5
subsd  %xmm4,%xmm5
subsd  %xmm4,%xmm0
subsd  %xmm5,%xmm3
addsd  %xmm0,%xmm3
movsd  0x8(%rcx,%rsi,8),%xmm0
add    $0x1,%rsi
ucomisd %xmm2,%xmm3
jne    125b <fast_expansion_sum_zeroelim+0x15b>
jnp    1220 <fast_expansion_sum_zeroelim+0x120>
movslq %eax,%rdi
add    $0x1,%eax
movsd  %xmm3,(%r8,%rdi,8)
jmp    1220 <fast_expansion_sum_zeroelim+0x120>
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    127b <fast_expansion_sum_zeroelim+0x17b>
jp     127b <fast_expansion_sum_zeroelim+0x17b>
test   %eax,%eax
je     127b <fast_expansion_sum_zeroelim+0x17b>
pop    %rbx
ret
movslq %eax,%rcx
add    $0x1,%eax
movsd  %xmm1,(%r8,%rcx,8)
pop    %rbx
ret
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm4
subsd  %xmm0,%xmm4
mov    %r9d,%eax
add    $0x1,%r9d
movsd  0x8(%rcx,%rax,8),%xmm0
subsd  %xmm4,%xmm3
xorpd  %xmm4,%xmm4
ucomisd %xmm4,%xmm3
jne    12b2 <fast_expansion_sum_zeroelim+0x1b2>
jnp    12c6 <fast_expansion_sum_zeroelim+0x1c6>
movsd  %xmm3,(%r8)
mov    $0x1,%eax
cmp    %edi,%r10d
jl     12d1 <fast_expansion_sum_zeroelim+0x1d1>
jmp    11a5 <fast_expansion_sum_zeroelim+0xa5>
xor    %eax,%eax
cmp    %edi,%r10d
jge    11a5 <fast_expansion_sum_zeroelim+0xa5>
cmp    %edx,%r9d
jge    11a5 <fast_expansion_sum_zeroelim+0xa5>
movapd 0xd1e(%rip),%xmm3
xorpd  %xmm4,%xmm4
cs nopw 0x0(%rax,%rax,1)
ucomisd %xmm2,%xmm0
setbe  %r11b
movapd %xmm2,%xmm5
xorpd  %xmm3,%xmm5
ucomisd %xmm5,%xmm0
seta   %bl
movapd %xmm1,%xmm6
cmp    %bl,%r11b
je     1350 <fast_expansion_sum_zeroelim+0x250>
addsd  %xmm2,%xmm6
movapd %xmm6,%xmm5
subsd  %xmm1,%xmm5
movapd %xmm6,%xmm7
subsd  %xmm5,%xmm7
subsd  %xmm5,%xmm2
subsd  %xmm7,%xmm1
movapd %xmm1,%xmm5
addsd  %xmm2,%xmm5
movslq %r10d,%rbx
add    $0x1,%r10d
movsd  0x8(%rsi,%rbx,8),%xmm2
jmp    1381 <fast_expansion_sum_zeroelim+0x281>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  %xmm0,%xmm6
movapd %xmm6,%xmm5
subsd  %xmm1,%xmm5
movapd %xmm6,%xmm7
subsd  %xmm5,%xmm7
subsd  %xmm5,%xmm0
subsd  %xmm7,%xmm1
movapd %xmm1,%xmm5
addsd  %xmm0,%xmm5
movslq %r9d,%rbx
add    $0x1,%r9d
movsd  0x8(%rcx,%rbx,8),%xmm0
movapd %xmm6,%xmm1
ucomisd %xmm4,%xmm5
jne    138d <fast_expansion_sum_zeroelim+0x28d>
jnp    1399 <fast_expansion_sum_zeroelim+0x299>
movslq %eax,%rbx
add    $0x1,%eax
movsd  %xmm5,(%r8,%rbx,8)
cmp    %edi,%r10d
jge    11a5 <fast_expansion_sum_zeroelim+0xa5>
cmp    %edx,%r9d
jl     12f0 <fast_expansion_sum_zeroelim+0x1f0>
jmp    11a5 <fast_expansion_sum_zeroelim+0xa5>

LLVM IR (-O0):
define i32 @fast_expansion_sum_zeroelim(i32 %0, double* %1, i32 %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca double*
%8 = alloca i32
%9 = alloca double*
%10 = alloca double*
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca double
%22 = alloca double
store i32 %0, i32* %6
store double* %1, double** %7
store i32 %2, i32* %8
store double* %3, double** %9
store double* %4, double** %10
%23 = load double*, double** %7
%24 = getelementptr inbounds double, double* %23, i64 0
%25 = load double, double* %24
store double %25, double* %21
%26 = load double*, double** %9
%27 = getelementptr inbounds double, double* %26, i64 0
%28 = load double, double* %27
store double %28, double* %22
store i32 0, i32* %19
store i32 0, i32* %18
%29 = load double, double* %22
%30 = load double, double* %21
%31 = fcmp ogt double %29, %30
%32 = zext i1 %31 to i32
%33 = load double, double* %22
%34 = load double, double* %21
%35 = fneg double %34
%36 = fcmp ogt double %33, %35
%37 = zext i1 %36 to i32
%38 = icmp eq i32 %32, %37
br i1 %38, label %39, label %47
39: ; preds = %5
%40 = load double, double* %21
store double %40, double* %11
%41 = load double*, double** %7
%42 = load i32, i32* %18
%43 = add nsw i32 %42, 1
store i32 %43, i32* %18
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %41, i64 %44
%46 = load double, double* %45
store double %46, double* %21
br label %55
47: ; preds = %5
%48 = load double, double* %22
store double %48, double* %11
%49 = load double*, double** %9
%50 = load i32, i32* %19
%51 = add nsw i32 %50, 1
store i32 %51, i32* %19
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %49, i64 %52
%54 = load double, double* %53
store double %54, double* %22
br label %55
55: ; preds = %47, %39
store i32 0, i32* %20
%56 = load i32, i32* %18
%57 = load i32, i32* %6
%58 = icmp slt i32 %56, %57
br i1 %58, label %59, label %202
59: ; preds = %55
%60 = load i32, i32* %19
%61 = load i32, i32* %8
%62 = icmp slt i32 %60, %61
br i1 %62, label %63, label %202
63: ; preds = %59
%64 = load double, double* %22
%65 = load double, double* %21
%66 = fcmp ogt double %64, %65
%67 = zext i1 %66 to i32
%68 = load double, double* %22
%69 = load double, double* %21
%70 = fneg double %69
%71 = fcmp ogt double %68, %70
%72 = zext i1 %71 to i32
%73 = icmp eq i32 %67, %72
br i1 %73, label %74, label %90
74: ; preds = %63
%75 = load double, double* %21
%76 = load double, double* %11
%77 = fadd double %75, %76
store double %77, double* %12
%78 = load double, double* %12
%79 = load double, double* %21
%80 = fsub double %78, %79
store double %80, double* %14
%81 = load double, double* %11
%82 = load double, double* %14
%83 = fsub double %81, %82
store double %83, double* %13
%84 = load double*, double** %7
%85 = load i32, i32* %18
%86 = add nsw i32 %85, 1
store i32 %86, i32* %18
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %84, i64 %87
%89 = load double, double* %88
store double %89, double* %21
br label %106
90: ; preds = %63
%91 = load double, double* %22
%92 = load double, double* %11
%93 = fadd double %91, %92
store double %93, double* %12
%94 = load double, double* %12
%95 = load double, double* %22
%96 = fsub double %94, %95
store double %96, double* %14
%97 = load double, double* %11
%98 = load double, double* %14
%99 = fsub double %97, %98
store double %99, double* %13
%100 = load double*, double** %9
%101 = load i32, i32* %19
%102 = add nsw i32 %101, 1
store i32 %102, i32* %19
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %100, i64 %103
%105 = load double, double* %104
store double %105, double* %22
br label %106
106: ; preds = %90, %74
%107 = load double, double* %12
store double %107, double* %11
%108 = load double, double* %13
%109 = fcmp une double %108, 0.000000e+00
br i1 %109, label %110, label %117
110: ; preds = %106
%111 = load double, double* %13
%112 = load double*, double** %10
%113 = load i32, i32* %20
%114 = add nsw i32 %113, 1
store i32 %114, i32* %20
%115 = sext i32 %113 to i64
%116 = getelementptr inbounds double, double* %112, i64 %115
store double %111, double* %116
br label %117
117: ; preds = %110, %106
br label %118
118: ; preds = %200, %117
%119 = load i32, i32* %18
%120 = load i32, i32* %6
%121 = icmp slt i32 %119, %120
br i1 %121, label %122, label %126
122: ; preds = %118
%123 = load i32, i32* %19
%124 = load i32, i32* %8
%125 = icmp slt i32 %123, %124
br label %126
126: ; preds = %122, %118
%127 = phi i1 [ false, %118 ], [ %125, %122 ]
br i1 %127, label %128, label %201
128: ; preds = %126
%129 = load double, double* %22
%130 = load double, double* %21
%131 = fcmp ogt double %129, %130
%132 = zext i1 %131 to i32
%133 = load double, double* %22
%134 = load double, double* %21
%135 = fneg double %134
%136 = fcmp ogt double %133, %135
%137 = zext i1 %136 to i32
%138 = icmp eq i32 %132, %137
br i1 %138, label %139, label %164
139: ; preds = %128
%140 = load double, double* %11
%141 = load double, double* %21
%142 = fadd double %140, %141
store double %142, double* %12
%143 = load double, double* %12
%144 = load double, double* %11
%145 = fsub double %143, %144
store double %145, double* %14
%146 = load double, double* %12
%147 = load double, double* %14
%148 = fsub double %146, %147
store double %148, double* %15
%149 = load double, double* %21
%150 = load double, double* %14
%151 = fsub double %149, %150
store double %151, double* %16
%152 = load double, double* %11
%153 = load double, double* %15
%154 = fsub double %152, %153
store double %154, double* %17
%155 = load double, double* %17
%156 = load double, double* %16
%157 = fadd double %155, %156
store double %157, double* %13
%158 = load double*, double** %7
%159 = load i32, i32* %18
%160 = add nsw i32 %159, 1
store i32 %160, i32* %18
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds double, double* %158, i64 %161
%163 = load double, double* %162
store double %163, double* %21
br label %189
164: ; preds = %128
%165 = load double, double* %11
%166 = load double, double* %22
%167 = fadd double %165, %166
store double %167, double* %12
%168 = load double, double* %12
%169 = load double, double* %11
%170 = fsub double %168, %169
store double %170, double* %14
%171 = load double, double* %12
%172 = load double, double* %14
%173 = fsub double %171, %172
store double %173, double* %15
%174 = load double, double* %22
%175 = load double, double* %14
%176 = fsub double %174, %175
store double %176, double* %16
%177 = load double, double* %11
%178 = load double, double* %15
%179 = fsub double %177, %178
store double %179, double* %17
%180 = load double, double* %17
%181 = load double, double* %16
%182 = fadd double %180, %181
store double %182, double* %13
%183 = load double*, double** %9
%184 = load i32, i32* %19
%185 = add nsw i32 %184, 1
store i32 %185, i32* %19
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds double, double* %183, i64 %186
%188 = load double, double* %187
store double %188, double* %22
br label %189
189: ; preds = %164, %139
%190 = load double, double* %12
store double %190, double* %11
%191 = load double, double* %13
%192 = fcmp une double %191, 0.000000e+00
br i1 %192, label %193, label %200
193: ; preds = %189
%194 = load double, double* %13
%195 = load double*, double** %10
%196 = load i32, i32* %20
%197 = add nsw i32 %196, 1
store i32 %197, i32* %20
%198 = sext i32 %196 to i64
%199 = getelementptr inbounds double, double* %195, i64 %198
store double %194, double* %199
br label %200
200: ; preds = %193, %189
br label %118
201: ; preds = %126
br label %202
202: ; preds = %201, %59, %55
br label %203
203: ; preds = %242, %202
%204 = load i32, i32* %18
%205 = load i32, i32* %6
%206 = icmp slt i32 %204, %205
br i1 %206, label %207, label %243
207: ; preds = %203
%208 = load double, double* %11
%209 = load double, double* %21
%210 = fadd double %208, %209
store double %210, double* %12
%211 = load double, double* %12
%212 = load double, double* %11
%213 = fsub double %211, %212
store double %213, double* %14
%214 = load double, double* %12
%215 = load double, double* %14
%216 = fsub double %214, %215
store double %216, double* %15
%217 = load double, double* %21
%218 = load double, double* %14
%219 = fsub double %217, %218
store double %219, double* %16
%220 = load double, double* %11
%221 = load double, double* %15
%222 = fsub double %220, %221
store double %222, double* %17
%223 = load double, double* %17
%224 = load double, double* %16
%225 = fadd double %223, %224
store double %225, double* %13
%226 = load double*, double** %7
%227 = load i32, i32* %18
%228 = add nsw i32 %227, 1
store i32 %228, i32* %18
%229 = sext i32 %228 to i64
%230 = getelementptr inbounds double, double* %226, i64 %229
%231 = load double, double* %230
store double %231, double* %21
%232 = load double, double* %12
store double %232, double* %11
%233 = load double, double* %13
%234 = fcmp une double %233, 0.000000e+00
br i1 %234, label %235, label %242
235: ; preds = %207
%236 = load double, double* %13
%237 = load double*, double** %10
%238 = load i32, i32* %20
%239 = add nsw i32 %238, 1
store i32 %239, i32* %20
%240 = sext i32 %238 to i64
%241 = getelementptr inbounds double, double* %237, i64 %240
store double %236, double* %241
br label %242
242: ; preds = %235, %207
br label %203
243: ; preds = %203
br label %244
244: ; preds = %283, %243
%245 = load i32, i32* %19
%246 = load i32, i32* %8
%247 = icmp slt i32 %245, %246
br i1 %247, label %248, label %284
248: ; preds = %244
%249 = load double, double* %11
%250 = load double, double* %22
%251 = fadd double %249, %250
store double %251, double* %12
%252 = load double, double* %12
%253 = load double, double* %11
%254 = fsub double %252, %253
store double %254, double* %14
%255 = load double, double* %12
%256 = load double, double* %14
%257 = fsub double %255, %256
store double %257, double* %15
%258 = load double, double* %22
%259 = load double, double* %14
%260 = fsub double %258, %259
store double %260, double* %16
%261 = load double, double* %11
%262 = load double, double* %15
%263 = fsub double %261, %262
store double %263, double* %17
%264 = load double, double* %17
%265 = load double, double* %16
%266 = fadd double %264, %265
store double %266, double* %13
%267 = load double*, double** %9
%268 = load i32, i32* %19
%269 = add nsw i32 %268, 1
store i32 %269, i32* %19
%270 = sext i32 %269 to i64
%271 = getelementptr inbounds double, double* %267, i64 %270
%272 = load double, double* %271
store double %272, double* %22
%273 = load double, double* %12
store double %273, double* %11
%274 = load double, double* %13
%275 = fcmp une double %274, 0.000000e+00
br i1 %275, label %276, label %283
276: ; preds = %248
%277 = load double, double* %13
%278 = load double*, double** %10
%279 = load i32, i32* %20
%280 = add nsw i32 %279, 1
store i32 %280, i32* %20
%281 = sext i32 %279 to i64
%282 = getelementptr inbounds double, double* %278, i64 %281
store double %277, double* %282
br label %283
283: ; preds = %276, %248
br label %244
284: ; preds = %244
%285 = load double, double* %11
%286 = fcmp une double %285, 0.000000e+00
br i1 %286, label %290, label %287
287: ; preds = %284
%288 = load i32, i32* %20
%289 = icmp eq i32 %288, 0
br i1 %289, label %290, label %297
290: ; preds = %287, %284
%291 = load double, double* %11
%292 = load double*, double** %10
%293 = load i32, i32* %20
%294 = add nsw i32 %293, 1
store i32 %294, i32* %20
%295 = sext i32 %293 to i64
%296 = getelementptr inbounds double, double* %292, i64 %295
store double %291, double* %296
br label %297
297: ; preds = %290, %287
%298 = load i32, i32* %20
ret i32 %298
}

LLVM IR (-O2):
define i32 @fast_expansion_sum_zeroelim(i32 %0, double* nocapture %1, i32 %2, double* nocapture %3, double* nocapture writeonly %4) {
%6 = load double, double* %1
%7 = load double, double* %3
%8 = fcmp ule double %7, %6
%9 = fneg double %6
%10 = fcmp ogt double %7, %9
%11 = xor i1 %8, %10
br i1 %11, label %12, label %15
12: ; preds = %5
%13 = getelementptr inbounds double, double* %1, i64 1
%14 = load double, double* %13
br label %18
15: ; preds = %5
%16 = getelementptr inbounds double, double* %3, i64 1
%17 = load double, double* %16
br label %18
18: ; preds = %15, %12
%19 = phi double [ %6, %12 ], [ %7, %15 ]
%20 = phi i32 [ 1, %12 ], [ 0, %15 ]
%21 = phi i32 [ 0, %12 ], [ 1, %15 ]
%22 = phi double [ %14, %12 ], [ %6, %15 ]
%23 = phi double [ %7, %12 ], [ %17, %15 ]
%24 = icmp slt i32 %20, %0
%25 = icmp slt i32 %21, %2
%26 = and i1 %24, %25
br i1 %26, label %27, label %111
27: ; preds = %18
%28 = fcmp ule double %23, %22
%29 = fneg double %22
%30 = fcmp ogt double %23, %29
%31 = xor i1 %28, %30
br i1 %31, label %32, label %39
32: ; preds = %27
%33 = fadd double %19, %22
%34 = fsub double %33, %22
%35 = add nuw nsw i32 %20, 1
%36 = zext i32 %35 to i64
%37 = getelementptr inbounds double, double* %1, i64 %36
%38 = load double, double* %37
br label %46
39: ; preds = %27
%40 = fadd double %19, %23
%41 = fsub double %40, %23
%42 = add nuw nsw i32 %21, 1
%43 = zext i32 %42 to i64
%44 = getelementptr inbounds double, double* %3, i64 %43
%45 = load double, double* %44
br label %46
46: ; preds = %39, %32
%47 = phi double [ %33, %32 ], [ %40, %39 ]
%48 = phi double [ %34, %32 ], [ %41, %39 ]
%49 = phi i32 [ %35, %32 ], [ %20, %39 ]
%50 = phi i32 [ %21, %32 ], [ %42, %39 ]
%51 = phi double [ %38, %32 ], [ %22, %39 ]
%52 = phi double [ %23, %32 ], [ %45, %39 ]
%53 = fsub double %19, %48
%54 = fcmp une double %53, 0.000000e+00
br i1 %54, label %55, label %56
55: ; preds = %46
store double %53, double* %4
br label %56
56: ; preds = %55, %46
%57 = phi i32 [ 1, %55 ], [ 0, %46 ]
%58 = icmp slt i32 %49, %0
%59 = icmp slt i32 %50, %2
%60 = select i1 %58, i1 %59, i1 false
br i1 %60, label %61, label %111
61: ; preds = %56, %106
%62 = phi double [ %100, %106 ], [ %52, %56 ]
%63 = phi double [ %99, %106 ], [ %51, %56 ]
%64 = phi i32 [ %107, %106 ], [ %57, %56 ]
%65 = phi i32 [ %98, %106 ], [ %50, %56 ]
%66 = phi i32 [ %97, %106 ], [ %49, %56 ]
%67 = phi double [ %95, %106 ], [ %47, %56 ]
%68 = fcmp ule double %62, %63
%69 = fneg double %63
%70 = fcmp ogt double %62, %69
%71 = xor i1 %68, %70
br i1 %71, label %72, label %83
72: ; preds = %61
%73 = fadd double %67, %63
%74 = fsub double %73, %67
%75 = fsub double %73, %74
%76 = fsub double %63, %74
%77 = fsub double %67, %75
%78 = fadd double %76, %77
%79 = add nsw i32 %66, 1
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %1, i64 %80
%82 = load double, double* %81
br label %94
83: ; preds = %61
%84 = fadd double %67, %62
%85 = fsub double %84, %67
%86 = fsub double %84, %85
%87 = fsub double %62, %85
%88 = fsub double %67, %86
%89 = fadd double %87, %88
%90 = add nsw i32 %65, 1
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %3, i64 %91
%93 = load double, double* %92
br label %94
94: ; preds = %83, %72
%95 = phi double [ %73, %72 ], [ %84, %83 ]
%96 = phi double [ %78, %72 ], [ %89, %83 ]
%97 = phi i32 [ %79, %72 ], [ %66, %83 ]
%98 = phi i32 [ %65, %72 ], [ %90, %83 ]
%99 = phi double [ %82, %72 ], [ %63, %83 ]
%100 = phi double [ %62, %72 ], [ %93, %83 ]
%101 = fcmp une double %96, 0.000000e+00
br i1 %101, label %102, label %106
102: ; preds = %94
%103 = add nsw i32 %64, 1
%104 = sext i32 %64 to i64
%105 = getelementptr inbounds double, double* %4, i64 %104
store double %96, double* %105
br label %106
106: ; preds = %102, %94
%107 = phi i32 [ %103, %102 ], [ %64, %94 ]
%108 = icmp slt i32 %97, %0
%109 = icmp slt i32 %98, %2
%110 = select i1 %108, i1 %109, i1 false
br i1 %110, label %61, label %111
111: ; preds = %106, %56, %18
%112 = phi double [ %19, %18 ], [ %47, %56 ], [ %95, %106 ]
%113 = phi i32 [ %20, %18 ], [ %49, %56 ], [ %97, %106 ]
%114 = phi i32 [ %21, %18 ], [ %50, %56 ], [ %98, %106 ]
%115 = phi i32 [ 0, %18 ], [ %57, %56 ], [ %107, %106 ]
%116 = phi double [ %22, %18 ], [ %51, %56 ], [ %99, %106 ]
%117 = phi double [ %23, %18 ], [ %52, %56 ], [ %100, %106 ]
%118 = icmp slt i32 %113, %0
br i1 %118, label %119, label %122
119: ; preds = %111
%120 = sext i32 %113 to i64
%121 = sext i32 %0 to i64
br label %129
122: ; preds = %148, %111
%123 = phi double [ %112, %111 ], [ %134, %148 ]
%124 = phi i32 [ %115, %111 ], [ %149, %148 ]
%125 = icmp slt i32 %114, %2
br i1 %125, label %126, label %173
126: ; preds = %122
%127 = sext i32 %114 to i64
%128 = sext i32 %2 to i64
br label %151
129: ; preds = %119, %148
%130 = phi i64 [ %120, %119 ], [ %140, %148 ]
%131 = phi double [ %116, %119 ], [ %142, %148 ]
%132 = phi i32 [ %115, %119 ], [ %149, %148 ]
%133 = phi double [ %112, %119 ], [ %134, %148 ]
%134 = fadd double %133, %131
%135 = fsub double %134, %133
%136 = fsub double %134, %135
%137 = fsub double %131, %135
%138 = fsub double %133, %136
%139 = fadd double %137, %138
%140 = add nsw i64 %130, 1
%141 = getelementptr inbounds double, double* %1, i64 %140
%142 = load double, double* %141
%143 = fcmp une double %139, 0.000000e+00
br i1 %143, label %144, label %148
144: ; preds = %129
%145 = add nsw i32 %132, 1
%146 = sext i32 %132 to i64
%147 = getelementptr inbounds double, double* %4, i64 %146
store double %139, double* %147
br label %148
148: ; preds = %144, %129
%149 = phi i32 [ %145, %144 ], [ %132, %129 ]
%150 = icmp eq i64 %140, %121
br i1 %150, label %122, label %129
151: ; preds = %126, %170
%152 = phi i64 [ %127, %126 ], [ %162, %170 ]
%153 = phi double [ %117, %126 ], [ %164, %170 ]
%154 = phi i32 [ %124, %126 ], [ %171, %170 ]
%155 = phi double [ %123, %126 ], [ %156, %170 ]
%156 = fadd double %155, %153
%157 = fsub double %156, %155
%158 = fsub double %156, %157
%159 = fsub double %153, %157
%160 = fsub double %155, %158
%161 = fadd double %159, %160
%162 = add nsw i64 %152, 1
%163 = getelementptr inbounds double, double* %3, i64 %162
%164 = load double, double* %163
%165 = fcmp une double %161, 0.000000e+00
br i1 %165, label %166, label %170
166: ; preds = %151
%167 = add nsw i32 %154, 1
%168 = sext i32 %154 to i64
%169 = getelementptr inbounds double, double* %4, i64 %168
store double %161, double* %169
br label %170
170: ; preds = %166, %151
%171 = phi i32 [ %167, %166 ], [ %154, %151 ]
%172 = icmp eq i64 %162, %128
br i1 %172, label %173, label %151
173: ; preds = %170, %122
%174 = phi double [ %123, %122 ], [ %156, %170 ]
%175 = phi i32 [ %124, %122 ], [ %171, %170 ]
%176 = fcmp une double %174, 0.000000e+00
%177 = icmp eq i32 %175, 0
%178 = select i1 %176, i1 true, i1 %177
br i1 %178, label %179, label %183
179: ; preds = %173
%180 = add nsw i32 %175, 1
%181 = sext i32 %175 to i64
%182 = getelementptr inbounds double, double* %4, i64 %181
store double %174, double* %182
br label %183
183: ; preds = %173, %179
%184 = phi i32 [ %180, %179 ], [ %175, %173 ]
ret i32 %184
}