int fast_expansion_sum(elen, e, flen, f, h)
int elen;
float *e;
int flen;
float *f;
float *h;
{
  float Q;
  float Qnew;
  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;
      h[0] = Q - bvirt;
      enow = e[++eindex];
    } else {
      Qnew = (float)(fnow + Q);
      bvirt = Qnew - fnow;
      h[0] = Q - bvirt;
      fnow = f[++findex];
    }
    Q = Qnew;
    hindex = 1;
    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;
        h[hindex] = around + bround;
        enow = e[++eindex];
      } else {
        Qnew = (float)(Q + fnow);
        bvirt = (float)(Qnew - Q);
        avirt = Qnew - bvirt;
        bround = fnow - bvirt;
        around = Q - avirt;
        h[hindex] = around + bround;
        fnow = f[++findex];
      }
      Q = Qnew;
      hindex++;
    }
  }
  while (eindex < elen) {
    Qnew = (float)(Q + enow);
    bvirt = (float)(Qnew - Q);
    avirt = Qnew - bvirt;
    bround = enow - bvirt;
    around = Q - avirt;
    h[hindex] = around + bround;
    enow = e[++eindex];
    Q = Qnew;
    hindex++;
  }
  while (findex < flen) {
    Qnew = (float)(Q + fnow);
    bvirt = (float)(Qnew - Q);
    avirt = Qnew - bvirt;
    bround = fnow - bvirt;
    around = Q - avirt;
    h[hindex] = around + bround;
    fnow = f[++findex];
    Q = Qnew;
    hindex++;
  }
  h[hindex] = Q;
  return hindex + 1;
}


x86 assembly (clang, -O0):
<fast_expansion_sum>:
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,-0x50(%rbp)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x54(%rbp)
movl   $0x0,-0x48(%rbp)
movl   $0x0,-0x44(%rbp)
movss  -0x54(%rbp),%xmm0
ucomiss -0x50(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movss  -0x54(%rbp),%xmm0
movss  -0x50(%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+0xa4>
movss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x44(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
jmp    11c8 <fast_expansion_sum+0xc8>
movss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x20(%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)
movl   $0x0,-0x4c(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1452 <fast_expansion_sum+0x352>
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1452 <fast_expansion_sum+0x352>
movss  -0x54(%rbp),%xmm0
ucomiss -0x50(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movss  -0x54(%rbp),%xmm0
movss  -0x50(%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    1273 <fast_expansion_sum+0x173>
movss  -0x50(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x44(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
jmp    12bd <fast_expansion_sum+0x1bd>
movss  -0x54(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x20(%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)
movl   $0x1,-0x4c(%rbp)
mov    -0x44(%rbp),%ecx
xor    %eax,%eax
cmp    -0x4(%rbp),%ecx
mov    %al,-0x55(%rbp)
jge    12eb <fast_expansion_sum+0x1eb>
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
setl   %al
mov    %al,-0x55(%rbp)
mov    -0x55(%rbp),%al
test   $0x1,%al
jne    12fb <fast_expansion_sum+0x1fb>
jmp    144d <fast_expansion_sum+0x34d>
movss  -0x54(%rbp),%xmm0
ucomiss -0x50(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movss  -0x54(%rbp),%xmm0
movss  -0x50(%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    13b9 <fast_expansion_sum+0x2b9>
movss  -0x2c(%rbp),%xmm0
addss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x50(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x44(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
jmp    1435 <fast_expansion_sum+0x335>
movss  -0x2c(%rbp),%xmm0
addss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x54(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%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)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12ce <fast_expansion_sum+0x1ce>
jmp    1452 <fast_expansion_sum+0x352>
jmp    1457 <fast_expansion_sum+0x357>
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    14f7 <fast_expansion_sum+0x3f7>
movss  -0x2c(%rbp),%xmm0
addss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x50(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x44(%rbp)
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1457 <fast_expansion_sum+0x357>
jmp    14fc <fast_expansion_sum+0x3fc>
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    159c <fast_expansion_sum+0x49c>
movss  -0x2c(%rbp),%xmm0
addss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x54(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%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)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    14fc <fast_expansion_sum+0x3fc>
movss  -0x2c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<fast_expansion_sum>:
push   %r14
push   %rbx
movss  (%rsi),%xmm2
movss  (%rcx),%xmm0
ucomiss %xmm2,%xmm0
setbe  %al
movaps 0xee8(%rip),%xmm1
xorps  %xmm2,%xmm1
ucomiss %xmm1,%xmm0
seta   %bl
cmp    %bl,%al
je     1138 <fast_expansion_sum+0x38>
movaps %xmm2,%xmm1
movss  0x4(%rsi),%xmm2
xor    %r9d,%r9d
mov    $0x1,%r10d
jmp    1149 <fast_expansion_sum+0x49>
movaps %xmm0,%xmm1
movss  0x4(%rcx),%xmm0
mov    $0x1,%r9d
xor    %r10d,%r10d
xor    %eax,%eax
cmp    %edi,%r10d
jge    1286 <fast_expansion_sum+0x186>
cmp    %edx,%r9d
jge    1286 <fast_expansion_sum+0x186>
ucomiss %xmm2,%xmm0
setbe  %al
movaps 0xe96(%rip),%xmm3
xorps  %xmm2,%xmm3
ucomiss %xmm3,%xmm0
seta   %bl
movaps %xmm1,%xmm3
cmp    %bl,%al
je     119d <fast_expansion_sum+0x9d>
addss  %xmm2,%xmm3
movaps %xmm3,%xmm4
subss  %xmm2,%xmm4
subss  %xmm4,%xmm1
movss  %xmm1,(%r8)
mov    %r10d,%eax
add    $0x1,%r10d
movss  0x4(%rsi,%rax,4),%xmm2
jmp    11be <fast_expansion_sum+0xbe>
addss  %xmm0,%xmm3
movaps %xmm3,%xmm4
subss  %xmm0,%xmm4
subss  %xmm4,%xmm1
movss  %xmm1,(%r8)
mov    %r9d,%eax
add    $0x1,%r9d
movss  0x4(%rcx,%rax,4),%xmm0
movaps %xmm3,%xmm1
mov    $0x1,%eax
cmp    %edi,%r10d
jge    1286 <fast_expansion_sum+0x186>
cmp    %edx,%r9d
jge    1286 <fast_expansion_sum+0x186>
mov    $0x1,%eax
movaps 0xe1c(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
ucomiss %xmm2,%xmm0
setbe  %r11b
movaps %xmm2,%xmm4
xorps  %xmm3,%xmm4
ucomiss %xmm4,%xmm0
seta   %bl
movaps %xmm1,%xmm4
cmp    %bl,%r11b
je     1240 <fast_expansion_sum+0x140>
addss  %xmm2,%xmm4
movaps %xmm4,%xmm5
subss  %xmm1,%xmm5
movaps %xmm4,%xmm6
subss  %xmm5,%xmm6
subss  %xmm5,%xmm2
subss  %xmm6,%xmm1
addss  %xmm2,%xmm1
movss  %xmm1,(%r8,%rax,4)
movslq %r10d,%rbx
add    $0x1,%r10d
movss  0x4(%rsi,%rbx,4),%xmm2
jmp    1271 <fast_expansion_sum+0x171>
xchg   %ax,%ax
addss  %xmm0,%xmm4
movaps %xmm4,%xmm5
subss  %xmm1,%xmm5
movaps %xmm4,%xmm6
subss  %xmm5,%xmm6
subss  %xmm5,%xmm0
subss  %xmm6,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r8,%rax,4)
movslq %r9d,%rbx
add    $0x1,%r9d
movss  0x4(%rcx,%rbx,4),%xmm0
movaps %xmm4,%xmm1
add    $0x1,%rax
cmp    %edi,%r10d
jge    1286 <fast_expansion_sum+0x186>
cmp    %edx,%r9d
jl     11f0 <fast_expansion_sum+0xf0>
cmp    %edi,%r10d
jge    1364 <fast_expansion_sum+0x264>
mov    %eax,%eax
movslq %r10d,%r14
movslq %edi,%r11
sub    %r10d,%edi
mov    %r14,%rbx
not    %rbx
test   $0x1,%dil
je     12df <fast_expansion_sum+0x1df>
movaps %xmm1,%xmm3
addss  %xmm2,%xmm3
movaps %xmm3,%xmm4
subss  %xmm1,%xmm4
movaps %xmm3,%xmm5
subss  %xmm4,%xmm5
subss  %xmm4,%xmm2
subss  %xmm5,%xmm1
addss  %xmm2,%xmm1
movss  %xmm1,(%r8,%rax,4)
movss  0x4(%rsi,%r14,4),%xmm2
add    $0x1,%r14
add    $0x1,%rax
movaps %xmm3,%xmm1
add    %r11,%rbx
je     1364 <fast_expansion_sum+0x264>
sub    %r14,%r11
lea    (%rsi,%r14,4),%rdi
add    $0x8,%rdi
lea    (%r8,%rax,4),%rbx
add    $0x4,%rbx
xor    %esi,%esi
nopl   (%rax)
movaps %xmm1,%xmm3
addss  %xmm2,%xmm3
movaps %xmm3,%xmm4
subss  %xmm1,%xmm4
movaps %xmm3,%xmm5
subss  %xmm4,%xmm5
subss  %xmm4,%xmm2
subss  %xmm5,%xmm1
addss  %xmm2,%xmm1
movss  %xmm1,-0x4(%rbx,%rsi,4)
movss  -0x4(%rdi,%rsi,4),%xmm2
movaps %xmm3,%xmm1
addss  %xmm2,%xmm1
movaps %xmm1,%xmm4
subss  %xmm3,%xmm4
movaps %xmm1,%xmm5
subss  %xmm4,%xmm5
subss  %xmm4,%xmm2
subss  %xmm5,%xmm3
addss  %xmm2,%xmm3
movss  %xmm3,(%rbx,%rsi,4)
movss  (%rdi,%rsi,4),%xmm2
add    $0x2,%rsi
cmp    %rsi,%r11
jne    1300 <fast_expansion_sum+0x200>
add    %rsi,%rax
cmp    %edx,%r9d
jge    1444 <fast_expansion_sum+0x344>
mov    %eax,%eax
movslq %r9d,%rbx
movslq %edx,%rsi
sub    %r9d,%edx
mov    %rbx,%rdi
not    %rdi
test   $0x1,%dl
je     13bb <fast_expansion_sum+0x2bb>
movaps %xmm1,%xmm2
addss  %xmm0,%xmm2
movaps %xmm2,%xmm3
subss  %xmm1,%xmm3
movaps %xmm2,%xmm4
subss  %xmm3,%xmm4
subss  %xmm3,%xmm0
subss  %xmm4,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r8,%rax,4)
movss  0x4(%rcx,%rbx,4),%xmm0
add    $0x1,%rbx
add    $0x1,%rax
movaps %xmm2,%xmm1
add    %rsi,%rdi
je     1444 <fast_expansion_sum+0x344>
sub    %rbx,%rsi
lea    (%rcx,%rbx,4),%rdx
add    $0x8,%rdx
lea    (%r8,%rax,4),%rdi
add    $0x4,%rdi
xor    %ecx,%ecx
nopl   0x0(%rax)
movaps %xmm1,%xmm2
addss  %xmm0,%xmm2
movaps %xmm2,%xmm3
subss  %xmm1,%xmm3
movaps %xmm2,%xmm4
subss  %xmm3,%xmm4
subss  %xmm3,%xmm0
subss  %xmm4,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,-0x4(%rdi,%rcx,4)
movss  -0x4(%rdx,%rcx,4),%xmm0
movaps %xmm2,%xmm1
addss  %xmm0,%xmm1
movaps %xmm1,%xmm3
subss  %xmm2,%xmm3
movaps %xmm1,%xmm4
subss  %xmm3,%xmm4
subss  %xmm3,%xmm0
subss  %xmm4,%xmm2
addss  %xmm0,%xmm2
movss  %xmm2,(%rdi,%rcx,4)
movss  (%rdx,%rcx,4),%xmm0
add    $0x2,%rcx
cmp    %rcx,%rsi
jne    13e0 <fast_expansion_sum+0x2e0>
add    %rcx,%rax
mov    %eax,%ecx
movss  %xmm1,(%r8,%rcx,4)
add    $0x1,%eax
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define i32 @fast_expansion_sum(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 %104
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 %40
32: ; preds = %27
%33 = fadd float %19, %22
%34 = fsub float %33, %22
%35 = fsub float %19, %34
store float %35, float* %4
%36 = add nuw nsw i32 %20, 1
%37 = zext i32 %36 to i64
%38 = getelementptr inbounds float, float* %1, i64 %37
%39 = load float, float* %38
br label %48
40: ; preds = %27
%41 = fadd float %19, %23
%42 = fsub float %41, %23
%43 = fsub float %19, %42
store float %43, float* %4
%44 = add nuw nsw i32 %21, 1
%45 = zext i32 %44 to i64
%46 = getelementptr inbounds float, float* %3, i64 %45
%47 = load float, float* %46
br label %48
48: ; preds = %40, %32
%49 = phi float [ %33, %32 ], [ %41, %40 ]
%50 = phi i32 [ %36, %32 ], [ %20, %40 ]
%51 = phi i32 [ %21, %32 ], [ %44, %40 ]
%52 = phi float [ %39, %32 ], [ %22, %40 ]
%53 = phi float [ %23, %32 ], [ %47, %40 ]
%54 = icmp slt i32 %50, %0
%55 = icmp slt i32 %51, %2
%56 = select i1 %54, i1 %55, i1 false
br i1 %56, label %57, label %104
57: ; preds = %48, %92
%58 = phi i64 [ %98, %92 ], [ 1, %48 ]
%59 = phi float [ %97, %92 ], [ %53, %48 ]
%60 = phi float [ %96, %92 ], [ %52, %48 ]
%61 = phi i32 [ %95, %92 ], [ %51, %48 ]
%62 = phi i32 [ %94, %92 ], [ %50, %48 ]
%63 = phi float [ %93, %92 ], [ %49, %48 ]
%64 = fcmp ule float %59, %60
%65 = fneg float %60
%66 = fcmp ogt float %59, %65
%67 = xor i1 %64, %66
br i1 %67, label %68, label %80
68: ; preds = %57
%69 = fadd float %63, %60
%70 = fsub float %69, %63
%71 = fsub float %69, %70
%72 = fsub float %60, %70
%73 = fsub float %63, %71
%74 = fadd float %72, %73
%75 = getelementptr inbounds float, float* %4, i64 %58
store float %74, float* %75
%76 = add nsw i32 %62, 1
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds float, float* %1, i64 %77
%79 = load float, float* %78
br label %92
80: ; preds = %57
%81 = fadd float %63, %59
%82 = fsub float %81, %63
%83 = fsub float %81, %82
%84 = fsub float %59, %82
%85 = fsub float %63, %83
%86 = fadd float %84, %85
%87 = getelementptr inbounds float, float* %4, i64 %58
store float %86, float* %87
%88 = add nsw i32 %61, 1
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds float, float* %3, i64 %89
%91 = load float, float* %90
br label %92
92: ; preds = %80, %68
%93 = phi float [ %69, %68 ], [ %81, %80 ]
%94 = phi i32 [ %76, %68 ], [ %62, %80 ]
%95 = phi i32 [ %61, %68 ], [ %88, %80 ]
%96 = phi float [ %79, %68 ], [ %60, %80 ]
%97 = phi float [ %59, %68 ], [ %91, %80 ]
%98 = add nuw i64 %58, 1
%99 = icmp slt i32 %94, %0
%100 = icmp slt i32 %95, %2
%101 = select i1 %99, i1 %100, i1 false
br i1 %101, label %57, label %102
102: ; preds = %92
%103 = trunc i64 %98 to i32
br label %104
104: ; preds = %102, %48, %18
%105 = phi float [ %19, %18 ], [ %49, %48 ], [ %93, %102 ]
%106 = phi i32 [ %20, %18 ], [ %50, %48 ], [ %94, %102 ]
%107 = phi i32 [ %21, %18 ], [ %51, %48 ], [ %95, %102 ]
%108 = phi i32 [ 0, %18 ], [ 1, %48 ], [ %103, %102 ]
%109 = phi float [ %22, %18 ], [ %52, %48 ], [ %96, %102 ]
%110 = phi float [ %23, %18 ], [ %53, %48 ], [ %97, %102 ]
%111 = icmp slt i32 %106, %0
br i1 %111, label %112, label %145
112: ; preds = %104
%113 = zext i32 %108 to i64
%114 = sext i32 %106 to i64
%115 = sext i32 %0 to i64
%116 = sub nsw i64 %115, %114
%117 = xor i64 %114, -1
%118 = and i64 %116, 1
%119 = icmp eq i64 %118, 0
br i1 %119, label %132, label %120
120: ; preds = %112
%121 = fadd float %105, %109
%122 = fsub float %121, %105
%123 = fsub float %121, %122
%124 = fsub float %109, %122
%125 = fsub float %105, %123
%126 = fadd float %124, %125
%127 = getelementptr inbounds float, float* %4, i64 %113
store float %126, float* %127
%128 = add nsw i64 %114, 1
%129 = getelementptr inbounds float, float* %1, i64 %128
%130 = load float, float* %129
%131 = add nuw nsw i64 %113, 1
br label %132
132: ; preds = %120, %112
%133 = phi float [ undef, %112 ], [ %121, %120 ]
%134 = phi i64 [ undef, %112 ], [ %131, %120 ]
%135 = phi i64 [ %114, %112 ], [ %128, %120 ]
%136 = phi i64 [ %113, %112 ], [ %131, %120 ]
%137 = phi float [ %109, %112 ], [ %130, %120 ]
%138 = phi float [ %105, %112 ], [ %121, %120 ]
%139 = sub nsw i64 0, %115
%140 = icmp eq i64 %117, %139
br i1 %140, label %141, label %178
141: ; preds = %178, %132
%142 = phi float [ %133, %132 ], [ %194, %178 ]
%143 = phi i64 [ %134, %132 ], [ %204, %178 ]
%144 = trunc i64 %143 to i32
br label %145
145: ; preds = %141, %104
%146 = phi float [ %105, %104 ], [ %142, %141 ]
%147 = phi i32 [ %108, %104 ], [ %144, %141 ]
%148 = icmp slt i32 %107, %2
br i1 %148, label %149, label %238
149: ; preds = %145
%150 = zext i32 %147 to i64
%151 = sext i32 %107 to i64
%152 = sext i32 %2 to i64
%153 = sub nsw i64 %152, %151
%154 = xor i64 %151, -1
%155 = and i64 %153, 1
%156 = icmp eq i64 %155, 0
br i1 %156, label %169, label %157
157: ; preds = %149
%158 = fadd float %146, %110
%159 = fsub float %158, %146
%160 = fsub float %158, %159
%161 = fsub float %110, %159
%162 = fsub float %146, %160
%163 = fadd float %161, %162
%164 = getelementptr inbounds float, float* %4, i64 %150
store float %163, float* %164
%165 = add nsw i64 %151, 1
%166 = getelementptr inbounds float, float* %3, i64 %165
%167 = load float, float* %166
%168 = add nuw nsw i64 %150, 1
br label %169
169: ; preds = %157, %149
%170 = phi float [ undef, %149 ], [ %158, %157 ]
%171 = phi i64 [ undef, %149 ], [ %168, %157 ]
%172 = phi i64 [ %151, %149 ], [ %165, %157 ]
%173 = phi i64 [ %150, %149 ], [ %168, %157 ]
%174 = phi float [ %110, %149 ], [ %167, %157 ]
%175 = phi float [ %146, %149 ], [ %158, %157 ]
%176 = sub nsw i64 0, %152
%177 = icmp eq i64 %154, %176
br i1 %177, label %234, label %206
178: ; preds = %132, %178
%179 = phi i64 [ %201, %178 ], [ %135, %132 ]
%180 = phi i64 [ %204, %178 ], [ %136, %132 ]
%181 = phi float [ %203, %178 ], [ %137, %132 ]
%182 = phi float [ %194, %178 ], [ %138, %132 ]
%183 = fadd float %182, %181
%184 = fsub float %183, %182
%185 = fsub float %183, %184
%186 = fsub float %181, %184
%187 = fsub float %182, %185
%188 = fadd float %186, %187
%189 = getelementptr inbounds float, float* %4, i64 %180
store float %188, float* %189
%190 = add nsw i64 %179, 1
%191 = getelementptr inbounds float, float* %1, i64 %190
%192 = load float, float* %191
%193 = add nuw nsw i64 %180, 1
%194 = fadd float %183, %192
%195 = fsub float %194, %183
%196 = fsub float %194, %195
%197 = fsub float %192, %195
%198 = fsub float %183, %196
%199 = fadd float %197, %198
%200 = getelementptr inbounds float, float* %4, i64 %193
store float %199, float* %200
%201 = add nsw i64 %179, 2
%202 = getelementptr inbounds float, float* %1, i64 %201
%203 = load float, float* %202
%204 = add nuw nsw i64 %180, 2
%205 = icmp eq i64 %201, %115
br i1 %205, label %141, label %178
206: ; preds = %169, %206
%207 = phi i64 [ %229, %206 ], [ %172, %169 ]
%208 = phi i64 [ %232, %206 ], [ %173, %169 ]
%209 = phi float [ %231, %206 ], [ %174, %169 ]
%210 = phi float [ %222, %206 ], [ %175, %169 ]
%211 = fadd float %210, %209
%212 = fsub float %211, %210
%213 = fsub float %211, %212
%214 = fsub float %209, %212
%215 = fsub float %210, %213
%216 = fadd float %214, %215
%217 = getelementptr inbounds float, float* %4, i64 %208
store float %216, float* %217
%218 = add nsw i64 %207, 1
%219 = getelementptr inbounds float, float* %3, i64 %218
%220 = load float, float* %219
%221 = add nuw nsw i64 %208, 1
%222 = fadd float %211, %220
%223 = fsub float %222, %211
%224 = fsub float %222, %223
%225 = fsub float %220, %223
%226 = fsub float %211, %224
%227 = fadd float %225, %226
%228 = getelementptr inbounds float, float* %4, i64 %221
store float %227, float* %228
%229 = add nsw i64 %207, 2
%230 = getelementptr inbounds float, float* %3, i64 %229
%231 = load float, float* %230
%232 = add nuw nsw i64 %208, 2
%233 = icmp eq i64 %229, %152
br i1 %233, label %234, label %206
234: ; preds = %206, %169
%235 = phi float [ %170, %169 ], [ %222, %206 ]
%236 = phi i64 [ %171, %169 ], [ %232, %206 ]
%237 = trunc i64 %236 to i32
br label %238
238: ; preds = %234, %145
%239 = phi float [ %146, %145 ], [ %235, %234 ]
%240 = phi i32 [ %147, %145 ], [ %237, %234 ]
%241 = zext i32 %240 to i64
%242 = getelementptr inbounds float, float* %4, i64 %241
store float %239, float* %242
%243 = add nuw nsw i32 %240, 1
ret i32 %243
}