int fast_expansion_sum(elen, e, flen, f, h)
int elen;
double *e;
int flen;
double *f;
double *h;
{
  double Q;
  double Qnew;
  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;
      h[0] = Q - bvirt;
      enow = e[++eindex];
    } else {
      Qnew = (double)(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 = (double)(Q + enow);
        bvirt = (double)(Qnew - Q);
        avirt = Qnew - bvirt;
        bround = enow - bvirt;
        around = Q - avirt;
        h[hindex] = around + bround;
        enow = e[++eindex];
      } else {
        Qnew = (double)(Q + fnow);
        bvirt = (double)(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 = (double)(Q + enow);
    bvirt = (double)(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 = (double)(Q + fnow);
    bvirt = (double)(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
movsd  (%rax),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x60(%rbp)
movl   $0x0,-0x5c(%rbp)
movsd  -0x78(%rbp),%xmm0
ucomisd -0x70(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x78(%rbp),%xmm0
movsd  -0x70(%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    11af <fast_expansion_sum+0xaf>
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x5c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    11d3 <fast_expansion_sum+0xd3>
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x60(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x60(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1473 <fast_expansion_sum+0x373>
mov    -0x60(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1473 <fast_expansion_sum+0x373>
movsd  -0x78(%rbp),%xmm0
ucomisd -0x70(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x78(%rbp),%xmm0
movsd  -0x70(%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    1289 <fast_expansion_sum+0x189>
movsd  -0x70(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x5c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    12d3 <fast_expansion_sum+0x1d3>
movsd  -0x78(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
mov    -0x60(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x60(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x64(%rbp)
mov    -0x5c(%rbp),%ecx
xor    %eax,%eax
cmp    -0x4(%rbp),%ecx
mov    %al,-0x79(%rbp)
jge    1301 <fast_expansion_sum+0x201>
mov    -0x60(%rbp),%eax
cmp    -0x14(%rbp),%eax
setl   %al
mov    %al,-0x79(%rbp)
mov    -0x79(%rbp),%al
test   $0x1,%al
jne    1311 <fast_expansion_sum+0x211>
jmp    146e <fast_expansion_sum+0x36e>
movsd  -0x78(%rbp),%xmm0
ucomisd -0x70(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x78(%rbp),%xmm0
movsd  -0x70(%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    13da <fast_expansion_sum+0x2da>
movsd  -0x30(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x70(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x5c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    1456 <fast_expansion_sum+0x356>
movsd  -0x30(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x78(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x60(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x60(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    12e4 <fast_expansion_sum+0x1e4>
jmp    1473 <fast_expansion_sum+0x373>
jmp    1478 <fast_expansion_sum+0x378>
mov    -0x5c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1518 <fast_expansion_sum+0x418>
movsd  -0x30(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x70(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x5c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1478 <fast_expansion_sum+0x378>
jmp    151d <fast_expansion_sum+0x41d>
mov    -0x60(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    15bd <fast_expansion_sum+0x4bd>
movsd  -0x30(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x78(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x60(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x60(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    151d <fast_expansion_sum+0x41d>
movsd  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<fast_expansion_sum>:
push   %r14
push   %rbx
movsd  (%rsi),%xmm2
movsd  (%rcx),%xmm0
ucomisd %xmm2,%xmm0
setbe  %al
movapd 0xee6(%rip),%xmm1
xorpd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
seta   %bl
cmp    %bl,%al
je     113d <fast_expansion_sum+0x3d>
movapd %xmm2,%xmm1
movsd  0x8(%rsi),%xmm2
xor    %r9d,%r9d
mov    $0x1,%r10d
jmp    114f <fast_expansion_sum+0x4f>
movapd %xmm0,%xmm1
movsd  0x8(%rcx),%xmm0
mov    $0x1,%r9d
xor    %r10d,%r10d
xor    %eax,%eax
cmp    %edi,%r10d
jge    12a9 <fast_expansion_sum+0x1a9>
cmp    %edx,%r9d
jge    12a9 <fast_expansion_sum+0x1a9>
ucomisd %xmm2,%xmm0
setbe  %al
movapd 0xe8e(%rip),%xmm3
xorpd  %xmm2,%xmm3
ucomisd %xmm3,%xmm0
seta   %bl
movapd %xmm1,%xmm3
cmp    %bl,%al
je     11a9 <fast_expansion_sum+0xa9>
addsd  %xmm2,%xmm3
movapd %xmm3,%xmm4
subsd  %xmm2,%xmm4
subsd  %xmm4,%xmm1
movsd  %xmm1,(%r8)
mov    %r10d,%eax
add    $0x1,%r10d
movsd  0x8(%rsi,%rax,8),%xmm2
jmp    11cb <fast_expansion_sum+0xcb>
addsd  %xmm0,%xmm3
movapd %xmm3,%xmm4
subsd  %xmm0,%xmm4
subsd  %xmm4,%xmm1
movsd  %xmm1,(%r8)
mov    %r9d,%eax
add    $0x1,%r9d
movsd  0x8(%rcx,%rax,8),%xmm0
movapd %xmm3,%xmm1
mov    $0x1,%eax
cmp    %edi,%r10d
jge    12a9 <fast_expansion_sum+0x1a9>
cmp    %edx,%r9d
jge    12a9 <fast_expansion_sum+0x1a9>
mov    $0x1,%eax
movapd 0xe0d(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
ucomisd %xmm2,%xmm0
setbe  %r11b
movapd %xmm2,%xmm4
xorpd  %xmm3,%xmm4
ucomisd %xmm4,%xmm0
seta   %bl
movapd %xmm1,%xmm4
cmp    %bl,%r11b
je     1260 <fast_expansion_sum+0x160>
addsd  %xmm2,%xmm4
movapd %xmm4,%xmm5
subsd  %xmm1,%xmm5
movapd %xmm4,%xmm6
subsd  %xmm5,%xmm6
subsd  %xmm5,%xmm2
subsd  %xmm6,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%r8,%rax,8)
movslq %r10d,%rbx
add    $0x1,%r10d
movsd  0x8(%rsi,%rbx,8),%xmm2
jmp    1293 <fast_expansion_sum+0x193>
cs nopw 0x0(%rax,%rax,1)
nop
addsd  %xmm0,%xmm4
movapd %xmm4,%xmm5
subsd  %xmm1,%xmm5
movapd %xmm4,%xmm6
subsd  %xmm5,%xmm6
subsd  %xmm5,%xmm0
subsd  %xmm6,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r8,%rax,8)
movslq %r9d,%rbx
add    $0x1,%r9d
movsd  0x8(%rcx,%rbx,8),%xmm0
movapd %xmm4,%xmm1
add    $0x1,%rax
cmp    %edi,%r10d
jge    12a9 <fast_expansion_sum+0x1a9>
cmp    %edx,%r9d
jl     1200 <fast_expansion_sum+0x100>
cmp    %edi,%r10d
jge    139a <fast_expansion_sum+0x29a>
mov    %eax,%eax
movslq %r10d,%r14
movslq %edi,%r11
sub    %r10d,%edi
mov    %r14,%rbx
not    %rbx
test   $0x1,%dil
je     1306 <fast_expansion_sum+0x206>
movapd %xmm1,%xmm3
addsd  %xmm2,%xmm3
movapd %xmm3,%xmm4
subsd  %xmm1,%xmm4
movapd %xmm3,%xmm5
subsd  %xmm4,%xmm5
subsd  %xmm4,%xmm2
subsd  %xmm5,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,(%r8,%rax,8)
movsd  0x8(%rsi,%r14,8),%xmm2
add    $0x1,%r14
add    $0x1,%rax
movapd %xmm3,%xmm1
add    %r11,%rbx
je     139a <fast_expansion_sum+0x29a>
sub    %r14,%r11
lea    (%rsi,%r14,8),%rdi
add    $0x10,%rdi
lea    (%r8,%rax,8),%rbx
add    $0x8,%rbx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm1,%xmm3
addsd  %xmm2,%xmm3
movapd %xmm3,%xmm4
subsd  %xmm1,%xmm4
movapd %xmm3,%xmm5
subsd  %xmm4,%xmm5
subsd  %xmm4,%xmm2
subsd  %xmm5,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,-0x8(%rbx,%rsi,8)
movsd  -0x8(%rdi,%rsi,8),%xmm2
movapd %xmm3,%xmm1
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm4
subsd  %xmm3,%xmm4
movapd %xmm1,%xmm5
subsd  %xmm4,%xmm5
subsd  %xmm4,%xmm2
subsd  %xmm5,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rbx,%rsi,8)
movsd  (%rdi,%rsi,8),%xmm2
add    $0x2,%rsi
cmp    %rsi,%r11
jne    1330 <fast_expansion_sum+0x230>
add    %rsi,%rax
cmp    %edx,%r9d
jge    148a <fast_expansion_sum+0x38a>
mov    %eax,%eax
movslq %r9d,%rbx
movslq %edx,%rsi
sub    %r9d,%edx
mov    %rbx,%rdi
not    %rdi
test   $0x1,%dl
je     13f5 <fast_expansion_sum+0x2f5>
movapd %xmm1,%xmm2
addsd  %xmm0,%xmm2
movapd %xmm2,%xmm3
subsd  %xmm1,%xmm3
movapd %xmm2,%xmm4
subsd  %xmm3,%xmm4
subsd  %xmm3,%xmm0
subsd  %xmm4,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r8,%rax,8)
movsd  0x8(%rcx,%rbx,8),%xmm0
add    $0x1,%rbx
add    $0x1,%rax
movapd %xmm2,%xmm1
add    %rsi,%rdi
je     148a <fast_expansion_sum+0x38a>
sub    %rbx,%rsi
lea    (%rcx,%rbx,8),%rdx
add    $0x10,%rdx
lea    (%r8,%rax,8),%rdi
add    $0x8,%rdi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd %xmm1,%xmm2
addsd  %xmm0,%xmm2
movapd %xmm2,%xmm3
subsd  %xmm1,%xmm3
movapd %xmm2,%xmm4
subsd  %xmm3,%xmm4
subsd  %xmm3,%xmm0
subsd  %xmm4,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%rdi,%rcx,8)
movsd  -0x8(%rdx,%rcx,8),%xmm0
movapd %xmm2,%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm3
subsd  %xmm2,%xmm3
movapd %xmm1,%xmm4
subsd  %xmm3,%xmm4
subsd  %xmm3,%xmm0
subsd  %xmm4,%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,(%rdi,%rcx,8)
movsd  (%rdx,%rcx,8),%xmm0
add    $0x2,%rcx
cmp    %rcx,%rsi
jne    1420 <fast_expansion_sum+0x320>
add    %rcx,%rax
mov    %eax,%ecx
movsd  %xmm1,(%r8,%rcx,8)
add    $0x1,%eax
pop    %rbx
pop    %r14
ret

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