int fast_expansion_sum_zeroelim(elen, e, flen, f, h)
int elen;
float *e;
int flen;
float *f;
float *h;
{
  float Q;
  float Qnew;
  float hh;
  float bvirt;
  float avirt, bround, around;
  int eindex, findex, hindex;
  float 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 = (float)(enow + Q);
      bvirt = Qnew - enow;
      hh = Q - bvirt;
      enow = e[++eindex];
    } else {
      Qnew = (float)(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 = (float)(Q + enow);
        bvirt = (float)(Qnew - Q);
        avirt = Qnew - bvirt;
        bround = enow - bvirt;
        around = Q - avirt;
        hh = around + bround;
        enow = e[++eindex];
      } else {
        Qnew = (float)(Q + fnow);
        bvirt = (float)(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 = (float)(Q + enow);
    bvirt = (float)(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 = (float)(Q + fnow);
    bvirt = (float)(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
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x58(%rbp)
movl   $0x0,-0x4c(%rbp)
movl   $0x0,-0x48(%rbp)
movss  -0x58(%rbp),%xmm0
ucomiss -0x54(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movss  -0x58(%rbp),%xmm0
movss  -0x54(%rbp),%xmm1
movd   %xmm1,%ecx
xor    $0x80000000,%ecx
movd   %ecx,%xmm1
ucomiss %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    11a4 <fast_expansion_sum_zeroelim+0xa4>
movss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
jmp    11c8 <fast_expansion_sum_zeroelim+0xc8>
movss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x4c(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x58(%rbp)
movl   $0x0,-0x50(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    14ab <fast_expansion_sum_zeroelim+0x3ab>
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    14ab <fast_expansion_sum_zeroelim+0x3ab>
movss  -0x58(%rbp),%xmm0
ucomiss -0x54(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movss  -0x58(%rbp),%xmm0
movss  -0x54(%rbp),%xmm1
movd   %xmm1,%ecx
xor    $0x80000000,%ecx
movd   %ecx,%xmm1
ucomiss %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    1270 <fast_expansion_sum_zeroelim+0x170>
movss  -0x54(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
jmp    12b7 <fast_expansion_sum_zeroelim+0x1b7>
movss  -0x58(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x4c(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12e2 <fast_expansion_sum_zeroelim+0x1e2>
jp     12e2 <fast_expansion_sum_zeroelim+0x1e2>
jmp    12fe <fast_expansion_sum_zeroelim+0x1fe>
movss  -0x34(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x50(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1303 <fast_expansion_sum_zeroelim+0x203>
mov    -0x48(%rbp),%ecx
xor    %eax,%eax
cmp    -0x4(%rbp),%ecx
mov    %al,-0x59(%rbp)
jge    1320 <fast_expansion_sum_zeroelim+0x220>
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
setl   %al
mov    %al,-0x59(%rbp)
mov    -0x59(%rbp),%al
test   $0x1,%al
jne    1330 <fast_expansion_sum_zeroelim+0x230>
jmp    14a6 <fast_expansion_sum_zeroelim+0x3a6>
movss  -0x58(%rbp),%xmm0
ucomiss -0x54(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movss  -0x58(%rbp),%xmm0
movss  -0x54(%rbp),%xmm1
movd   %xmm1,%ecx
xor    $0x80000000,%ecx
movd   %ecx,%xmm1
ucomiss %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    13e6 <fast_expansion_sum_zeroelim+0x2e6>
movss  -0x2c(%rbp),%xmm0
addss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x54(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
jmp    145a <fast_expansion_sum_zeroelim+0x35a>
movss  -0x2c(%rbp),%xmm0
addss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x58(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x4c(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1485 <fast_expansion_sum_zeroelim+0x385>
jp     1485 <fast_expansion_sum_zeroelim+0x385>
jmp    14a1 <fast_expansion_sum_zeroelim+0x3a1>
movss  -0x34(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x50(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1303 <fast_expansion_sum_zeroelim+0x203>
jmp    14ab <fast_expansion_sum_zeroelim+0x3ab>
jmp    14b0 <fast_expansion_sum_zeroelim+0x3b0>
mov    -0x48(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    157c <fast_expansion_sum_zeroelim+0x47c>
movss  -0x2c(%rbp),%xmm0
addss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x54(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    155b <fast_expansion_sum_zeroelim+0x45b>
jp     155b <fast_expansion_sum_zeroelim+0x45b>
jmp    1577 <fast_expansion_sum_zeroelim+0x477>
movss  -0x34(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x50(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    14b0 <fast_expansion_sum_zeroelim+0x3b0>
jmp    1581 <fast_expansion_sum_zeroelim+0x481>
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    164d <fast_expansion_sum_zeroelim+0x54d>
movss  -0x2c(%rbp),%xmm0
addss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x58(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x44(%rbp),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x4c(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x34(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    162c <fast_expansion_sum_zeroelim+0x52c>
jp     162c <fast_expansion_sum_zeroelim+0x52c>
jmp    1648 <fast_expansion_sum_zeroelim+0x548>
movss  -0x34(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x50(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1581 <fast_expansion_sum_zeroelim+0x481>
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1673 <fast_expansion_sum_zeroelim+0x573>
jp     1673 <fast_expansion_sum_zeroelim+0x573>
cmpl   $0x0,-0x50(%rbp)
jne    168f <fast_expansion_sum_zeroelim+0x58f>
movss  -0x2c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x50(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x50(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<fast_expansion_sum_zeroelim>:
push   %rbx
movss  (%rsi),%xmm2
movss  (%rcx),%xmm0
ucomiss %xmm2,%xmm0
setbe  %r9b
movaps 0xee9(%rip),%xmm1
xorps  %xmm2,%xmm1
ucomiss %xmm1,%xmm0
seta   %al
cmp    %al,%r9b
je     1138 <fast_expansion_sum_zeroelim+0x38>
movaps %xmm2,%xmm3
movss  0x4(%rsi),%xmm2
xor    %r9d,%r9d
mov    $0x1,%r10d
jmp    1149 <fast_expansion_sum_zeroelim+0x49>
movaps %xmm0,%xmm3
movss  0x4(%rcx),%xmm0
mov    $0x1,%r9d
xor    %r10d,%r10d
xor    %eax,%eax
cmp    %edi,%r10d
jge    1195 <fast_expansion_sum_zeroelim+0x95>
cmp    %edx,%r9d
jge    1195 <fast_expansion_sum_zeroelim+0x95>
ucomiss %xmm2,%xmm0
setbe  %r11b
movaps 0xe9d(%rip),%xmm1
xorps  %xmm2,%xmm1
ucomiss %xmm1,%xmm0
seta   %al
movaps %xmm3,%xmm1
cmp    %al,%r11b
je     1273 <fast_expansion_sum_zeroelim+0x173>
addss  %xmm2,%xmm1
movaps %xmm1,%xmm4
subss  %xmm2,%xmm4
mov    %r10d,%eax
add    $0x1,%r10d
movss  0x4(%rsi,%rax,4),%xmm2
jmp    128b <fast_expansion_sum_zeroelim+0x18b>
movaps %xmm3,%xmm1
cmp    %edi,%r10d
jge    11f6 <fast_expansion_sum_zeroelim+0xf6>
movslq %r10d,%r10
movslq %edi,%rdi
xorps  %xmm3,%xmm3
jmp    11b5 <fast_expansion_sum_zeroelim+0xb5>
nopl   0x0(%rax,%rax,1)
cmp    %r10,%rdi
je     11f6 <fast_expansion_sum_zeroelim+0xf6>
movaps %xmm1,%xmm4
addss  %xmm2,%xmm1
movaps %xmm1,%xmm5
subss  %xmm4,%xmm5
movaps %xmm1,%xmm6
subss  %xmm5,%xmm6
subss  %xmm5,%xmm2
subss  %xmm6,%xmm4
addss  %xmm2,%xmm4
movss  0x4(%rsi,%r10,4),%xmm2
add    $0x1,%r10
ucomiss %xmm3,%xmm4
jne    11e8 <fast_expansion_sum_zeroelim+0xe8>
jnp    11b0 <fast_expansion_sum_zeroelim+0xb0>
movslq %eax,%rbx
add    $0x1,%eax
movss  %xmm4,(%r8,%rbx,4)
jmp    11b0 <fast_expansion_sum_zeroelim+0xb0>
cmp    %edx,%r9d
jge    1255 <fast_expansion_sum_zeroelim+0x155>
movslq %r9d,%rsi
movslq %edx,%rdx
xorps  %xmm2,%xmm2
jmp    1215 <fast_expansion_sum_zeroelim+0x115>
cs nopw 0x0(%rax,%rax,1)
cmp    %rsi,%rdx
je     1255 <fast_expansion_sum_zeroelim+0x155>
movaps %xmm1,%xmm3
addss  %xmm0,%xmm1
movaps %xmm1,%xmm4
subss  %xmm3,%xmm4
movaps %xmm1,%xmm5
subss  %xmm4,%xmm5
subss  %xmm4,%xmm0
subss  %xmm5,%xmm3
addss  %xmm0,%xmm3
movss  0x4(%rcx,%rsi,4),%xmm0
add    $0x1,%rsi
ucomiss %xmm2,%xmm3
jne    1247 <fast_expansion_sum_zeroelim+0x147>
jnp    1210 <fast_expansion_sum_zeroelim+0x110>
movslq %eax,%rdi
add    $0x1,%eax
movss  %xmm3,(%r8,%rdi,4)
jmp    1210 <fast_expansion_sum_zeroelim+0x110>
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm1
jne    1265 <fast_expansion_sum_zeroelim+0x165>
jp     1265 <fast_expansion_sum_zeroelim+0x165>
test   %eax,%eax
je     1265 <fast_expansion_sum_zeroelim+0x165>
pop    %rbx
ret
movslq %eax,%rcx
add    $0x1,%eax
movss  %xmm1,(%r8,%rcx,4)
pop    %rbx
ret
addss  %xmm0,%xmm1
movaps %xmm1,%xmm4
subss  %xmm0,%xmm4
mov    %r9d,%eax
add    $0x1,%r9d
movss  0x4(%rcx,%rax,4),%xmm0
subss  %xmm4,%xmm3
xorps  %xmm4,%xmm4
ucomiss %xmm4,%xmm3
jne    1299 <fast_expansion_sum_zeroelim+0x199>
jnp    12ad <fast_expansion_sum_zeroelim+0x1ad>
movss  %xmm3,(%r8)
mov    $0x1,%eax
cmp    %edi,%r10d
jl     12b8 <fast_expansion_sum_zeroelim+0x1b8>
jmp    1198 <fast_expansion_sum_zeroelim+0x98>
xor    %eax,%eax
cmp    %edi,%r10d
jge    1198 <fast_expansion_sum_zeroelim+0x98>
cmp    %edx,%r9d
jge    1198 <fast_expansion_sum_zeroelim+0x98>
movaps 0xd38(%rip),%xmm3
xorps  %xmm4,%xmm4
nopl   0x0(%rax,%rax,1)
ucomiss %xmm2,%xmm0
setbe  %r11b
movaps %xmm2,%xmm5
xorps  %xmm3,%xmm5
ucomiss %xmm5,%xmm0
seta   %bl
movaps %xmm1,%xmm6
cmp    %bl,%r11b
je     1320 <fast_expansion_sum_zeroelim+0x220>
addss  %xmm2,%xmm6
movaps %xmm6,%xmm5
subss  %xmm1,%xmm5
movaps %xmm6,%xmm7
subss  %xmm5,%xmm7
subss  %xmm5,%xmm2
subss  %xmm7,%xmm1
movaps %xmm1,%xmm5
addss  %xmm2,%xmm5
movslq %r10d,%rbx
add    $0x1,%r10d
movss  0x4(%rsi,%rbx,4),%xmm2
jmp    134e <fast_expansion_sum_zeroelim+0x24e>
nopl   0x0(%rax,%rax,1)
addss  %xmm0,%xmm6
movaps %xmm6,%xmm5
subss  %xmm1,%xmm5
movaps %xmm6,%xmm7
subss  %xmm5,%xmm7
subss  %xmm5,%xmm0
subss  %xmm7,%xmm1
movaps %xmm1,%xmm5
addss  %xmm0,%xmm5
movslq %r9d,%rbx
add    $0x1,%r9d
movss  0x4(%rcx,%rbx,4),%xmm0
movaps %xmm6,%xmm1
ucomiss %xmm4,%xmm5
jne    1358 <fast_expansion_sum_zeroelim+0x258>
jnp    1364 <fast_expansion_sum_zeroelim+0x264>
movslq %eax,%rbx
add    $0x1,%eax
movss  %xmm5,(%r8,%rbx,4)
cmp    %edi,%r10d
jge    1198 <fast_expansion_sum_zeroelim+0x98>
cmp    %edx,%r9d
jl     12d0 <fast_expansion_sum_zeroelim+0x1d0>
jmp    1198 <fast_expansion_sum_zeroelim+0x98>

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

LLVM IR (-O2):
define i32 @fast_expansion_sum_zeroelim(i32 %0, float* nocapture %1, i32 %2, float* nocapture %3, float* nocapture writeonly %4) {
%6 = load float, float* %1
%7 = load float, float* %3
%8 = fcmp ule float %7, %6
%9 = fneg float %6
%10 = fcmp ogt float %7, %9
%11 = xor i1 %8, %10
br i1 %11, label %12, label %15
12: ; preds = %5
%13 = getelementptr inbounds float, float* %1, i64 1
%14 = load float, float* %13
br label %18
15: ; preds = %5
%16 = getelementptr inbounds float, float* %3, i64 1
%17 = load float, float* %16
br label %18
18: ; preds = %15, %12
%19 = phi float [ %6, %12 ], [ %7, %15 ]
%20 = phi i32 [ 1, %12 ], [ 0, %15 ]
%21 = phi i32 [ 0, %12 ], [ 1, %15 ]
%22 = phi float [ %14, %12 ], [ %6, %15 ]
%23 = phi float [ %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 float %23, %22
%29 = fneg float %22
%30 = fcmp ogt float %23, %29
%31 = xor i1 %28, %30
br i1 %31, label %32, label %39
32: ; preds = %27
%33 = fadd float %19, %22
%34 = fsub float %33, %22
%35 = add nuw nsw i32 %20, 1
%36 = zext i32 %35 to i64
%37 = getelementptr inbounds float, float* %1, i64 %36
%38 = load float, float* %37
br label %46
39: ; preds = %27
%40 = fadd float %19, %23
%41 = fsub float %40, %23
%42 = add nuw nsw i32 %21, 1
%43 = zext i32 %42 to i64
%44 = getelementptr inbounds float, float* %3, i64 %43
%45 = load float, float* %44
br label %46
46: ; preds = %39, %32
%47 = phi float [ %33, %32 ], [ %40, %39 ]
%48 = phi float [ %34, %32 ], [ %41, %39 ]
%49 = phi i32 [ %35, %32 ], [ %20, %39 ]
%50 = phi i32 [ %21, %32 ], [ %42, %39 ]
%51 = phi float [ %38, %32 ], [ %22, %39 ]
%52 = phi float [ %23, %32 ], [ %45, %39 ]
%53 = fsub float %19, %48
%54 = fcmp une float %53, 0.000000e+00
br i1 %54, label %55, label %56
55: ; preds = %46
store float %53, float* %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 float [ %100, %106 ], [ %52, %56 ]
%63 = phi float [ %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 float [ %95, %106 ], [ %47, %56 ]
%68 = fcmp ule float %62, %63
%69 = fneg float %63
%70 = fcmp ogt float %62, %69
%71 = xor i1 %68, %70
br i1 %71, label %72, label %83
72: ; preds = %61
%73 = fadd float %67, %63
%74 = fsub float %73, %67
%75 = fsub float %73, %74
%76 = fsub float %63, %74
%77 = fsub float %67, %75
%78 = fadd float %76, %77
%79 = add nsw i32 %66, 1
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds float, float* %1, i64 %80
%82 = load float, float* %81
br label %94
83: ; preds = %61
%84 = fadd float %67, %62
%85 = fsub float %84, %67
%86 = fsub float %84, %85
%87 = fsub float %62, %85
%88 = fsub float %67, %86
%89 = fadd float %87, %88
%90 = add nsw i32 %65, 1
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds float, float* %3, i64 %91
%93 = load float, float* %92
br label %94
94: ; preds = %83, %72
%95 = phi float [ %73, %72 ], [ %84, %83 ]
%96 = phi float [ %78, %72 ], [ %89, %83 ]
%97 = phi i32 [ %79, %72 ], [ %66, %83 ]
%98 = phi i32 [ %65, %72 ], [ %90, %83 ]
%99 = phi float [ %82, %72 ], [ %63, %83 ]
%100 = phi float [ %62, %72 ], [ %93, %83 ]
%101 = fcmp une float %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 float, float* %4, i64 %104
store float %96, float* %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 float [ %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 float [ %22, %18 ], [ %51, %56 ], [ %99, %106 ]
%117 = phi float [ %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 float [ %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 float [ %116, %119 ], [ %142, %148 ]
%132 = phi i32 [ %115, %119 ], [ %149, %148 ]
%133 = phi float [ %112, %119 ], [ %134, %148 ]
%134 = fadd float %133, %131
%135 = fsub float %134, %133
%136 = fsub float %134, %135
%137 = fsub float %131, %135
%138 = fsub float %133, %136
%139 = fadd float %137, %138
%140 = add nsw i64 %130, 1
%141 = getelementptr inbounds float, float* %1, i64 %140
%142 = load float, float* %141
%143 = fcmp une float %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 float, float* %4, i64 %146
store float %139, float* %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 float [ %117, %126 ], [ %164, %170 ]
%154 = phi i32 [ %124, %126 ], [ %171, %170 ]
%155 = phi float [ %123, %126 ], [ %156, %170 ]
%156 = fadd float %155, %153
%157 = fsub float %156, %155
%158 = fsub float %156, %157
%159 = fsub float %153, %157
%160 = fsub float %155, %158
%161 = fadd float %159, %160
%162 = add nsw i64 %152, 1
%163 = getelementptr inbounds float, float* %3, i64 %162
%164 = load float, float* %163
%165 = fcmp une float %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 float, float* %4, i64 %168
store float %161, float* %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 float [ %123, %122 ], [ %156, %170 ]
%175 = phi i32 [ %124, %122 ], [ %171, %170 ]
%176 = fcmp une float %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 float, float* %4, i64 %181
store float %174, float* %182
br label %183
183: ; preds = %173, %179
%184 = phi i32 [ %180, %179 ], [ %175, %173 ]
ret i32 %184
}