int linear_expansion_sum_zeroelim(int elen, double *e, int flen, double *f,
                                  double *h) {
  double Q, q, hh;
  double Qnew;
  double R;
  double bvirt;
  double avirt, bround, around;
  int eindex, findex, hindex;
  int count;
  double g0;
  eindex = findex = 0;
  hindex = 0;
  if (((f[findex]) > (e[eindex])) == ((f[findex]) > -(e[eindex]))) {
    g0 = (e[eindex]);
    ++eindex;
  } else {
    g0 = (f[findex]);
    ++findex;
  }
  if ((eindex < elen) && ((findex >= flen) || (((f[findex]) > (e[eindex])) ==
                                               ((f[findex]) > -(e[eindex]))))) {
    Qnew = (double)((e[eindex]) + g0);
    bvirt = Qnew - (e[eindex]);
    q = g0 - bvirt;
    ++eindex;
  } else {
    Qnew = (double)((f[findex]) + g0);
    bvirt = Qnew - (f[findex]);
    q = g0 - bvirt;
    ++findex;
  }
  Q = Qnew;
  for (count = 2; count < elen + flen; count++) {
    if ((eindex < elen) &&
        ((findex >= flen) ||
         (((f[findex]) > (e[eindex])) == ((f[findex]) > -(e[eindex]))))) {
      R = (double)((e[eindex]) + q);
      bvirt = R - (e[eindex]);
      hh = q - bvirt;
      ++eindex;
    } else {
      R = (double)((f[findex]) + q);
      bvirt = R - (f[findex]);
      hh = q - bvirt;
      ++findex;
    }
    Qnew = (double)(Q + R);
    bvirt = (double)(Qnew - Q);
    avirt = Qnew - bvirt;
    bround = R - bvirt;
    around = Q - avirt;
    q = around + bround;
    Q = Qnew;
    if (hh != 0) {
      h[hindex++] = hh;
    }
  }
  if (q != 0) {
    h[hindex++] = q;
  }
  if ((Q != 0.0) || (hindex == 0)) {
    h[hindex++] = Q;
  }
  return hindex;
}


x86 assembly (clang, -O0):
<linear_expansion_sum_zeroelim>:
push   %rbp
mov    %rsp,%rbp
push   %rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x0,-0x78(%rbp)
movl   $0x0,-0x74(%rbp)
movl   $0x0,-0x7c(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    -0x20(%rbp),%rcx
movslq -0x78(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
mov    -0x10(%rbp),%rcx
movslq -0x74(%rbp),%rdx
movsd  (%rcx,%rdx,8),%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    11b7 <linear_expansion_sum_zeroelim+0xb7>
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    11d5 <linear_expansion_sum_zeroelim+0xd5>
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x78(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x78(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12a6 <linear_expansion_sum_zeroelim+0x1a6>
mov    -0x78(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1255 <linear_expansion_sum_zeroelim+0x155>
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    -0x20(%rbp),%rcx
movslq -0x78(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
mov    -0x10(%rbp),%rcx
movslq -0x74(%rbp),%rdx
movsd  (%rcx,%rdx,8),%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    12a6 <linear_expansion_sum_zeroelim+0x1a6>
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    12f2 <linear_expansion_sum_zeroelim+0x1f2>
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x78(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x78(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x2,-0x80(%rbp)
mov    -0x80(%rbp),%eax
mov    -0x4(%rbp),%ecx
add    -0x14(%rbp),%ecx
cmp    %ecx,%eax
jge    14d5 <linear_expansion_sum_zeroelim+0x3d5>
mov    -0x74(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13df <linear_expansion_sum_zeroelim+0x2df>
mov    -0x78(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1394 <linear_expansion_sum_zeroelim+0x294>
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    -0x20(%rbp),%rcx
movslq -0x78(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
mov    -0x10(%rbp),%rcx
movslq -0x74(%rbp),%rdx
movsd  (%rcx,%rdx,8),%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    13df <linear_expansion_sum_zeroelim+0x2df>
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    1425 <linear_expansion_sum_zeroelim+0x325>
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x78(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x78(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x78(%rbp)
movsd  -0x30(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    14a6 <linear_expansion_sum_zeroelim+0x3a6>
jp     14a6 <linear_expansion_sum_zeroelim+0x3a6>
jmp    14c2 <linear_expansion_sum_zeroelim+0x3c2>
movsd  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x7c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x7c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    14c7 <linear_expansion_sum_zeroelim+0x3c7>
mov    -0x80(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x80(%rbp)
jmp    1303 <linear_expansion_sum_zeroelim+0x203>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    14f2 <linear_expansion_sum_zeroelim+0x3f2>
jp     14f2 <linear_expansion_sum_zeroelim+0x3f2>
jmp    150e <linear_expansion_sum_zeroelim+0x40e>
movsd  -0x38(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x7c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x7c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1530 <linear_expansion_sum_zeroelim+0x430>
jp     1530 <linear_expansion_sum_zeroelim+0x430>
cmpl   $0x0,-0x7c(%rbp)
jne    154c <linear_expansion_sum_zeroelim+0x44c>
movsd  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x7c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x7c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x7c(%rbp),%eax
add    $0x8,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<linear_expansion_sum_zeroelim>:
push   %rbp
push   %rbx
movsd  (%rcx),%xmm1
movsd  (%rsi),%xmm0
ucomisd %xmm0,%xmm1
setbe  %r9b
movapd 0xee6(%rip),%xmm2
xorpd  %xmm0,%xmm2
ucomisd %xmm2,%xmm1
seta   %al
xor    %r9b,%al
movzbl %al,%r9d
mov    %r9d,%eax
xor    $0x1,%al
test   %r9b,%r9b
jne    113a <linear_expansion_sum_zeroelim+0x3a>
movapd %xmm1,%xmm0
movzbl %al,%r10d
cmp    %edi,%r9d
jge    117a <linear_expansion_sum_zeroelim+0x7a>
cmp    %edx,%r10d
jge    119d <linear_expansion_sum_zeroelim+0x9d>
mov    %r10d,%eax
movsd  (%rcx,%rax,8),%xmm2
mov    %r9d,%eax
movsd  (%rsi,%rax,8),%xmm4
ucomisd %xmm4,%xmm2
setbe  %r11b
movapd 0xe98(%rip),%xmm1
xorpd  %xmm4,%xmm1
ucomisd %xmm1,%xmm2
seta   %al
cmp    %al,%r11b
je     1182 <linear_expansion_sum_zeroelim+0x82>
jmp    11a5 <linear_expansion_sum_zeroelim+0xa5>
mov    %r10d,%eax
movsd  (%rcx,%rax,8),%xmm2
movapd %xmm0,%xmm1
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm3
subsd  %xmm2,%xmm3
mov    $0x2,%r10d
sub    %r9d,%r10d
jmp    11b9 <linear_expansion_sum_zeroelim+0xb9>
mov    %r9d,%eax
movsd  (%rsi,%rax,8),%xmm4
movapd %xmm0,%xmm1
addsd  %xmm4,%xmm1
movapd %xmm1,%xmm3
subsd  %xmm4,%xmm3
add    $0x1,%r9d
subsd  %xmm3,%xmm0
lea    (%rdx,%rdi,1),%r11d
xor    %eax,%eax
cmp    $0x3,%r11d
jl     12bd <linear_expansion_sum_zeroelim+0x1bd>
lea    (%rdx,%rdi,1),%r11d
add    $0xfffffffe,%r11d
xor    %eax,%eax
movapd 0xe21(%rip),%xmm2
xorpd  %xmm3,%xmm3
jmp    1226 <linear_expansion_sum_zeroelim+0x126>
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm1,%xmm4
addsd  %xmm5,%xmm4
movapd %xmm4,%xmm0
subsd  %xmm1,%xmm0
movapd %xmm4,%xmm6
subsd  %xmm0,%xmm6
subsd  %xmm0,%xmm5
subsd  %xmm6,%xmm1
movapd %xmm1,%xmm0
addsd  %xmm5,%xmm0
movapd %xmm4,%xmm1
add    $0xffffffff,%r11d
je     12c1 <linear_expansion_sum_zeroelim+0x1c1>
cmp    %edi,%r9d
jge    1260 <linear_expansion_sum_zeroelim+0x160>
cmp    %edx,%r10d
jge    1280 <linear_expansion_sum_zeroelim+0x180>
movslq %r10d,%rbx
movsd  (%rcx,%rbx,8),%xmm6
movslq %r9d,%rbx
movsd  (%rsi,%rbx,8),%xmm7
ucomisd %xmm7,%xmm6
setbe  %bpl
movapd %xmm7,%xmm4
xorpd  %xmm2,%xmm4
ucomisd %xmm4,%xmm6
seta   %bl
cmp    %bl,%bpl
je     1268 <linear_expansion_sum_zeroelim+0x168>
jmp    1288 <linear_expansion_sum_zeroelim+0x188>
xchg   %ax,%ax
movslq %r10d,%rbx
movsd  (%rcx,%rbx,8),%xmm6
movapd %xmm0,%xmm5
addsd  %xmm6,%xmm5
movapd %xmm5,%xmm4
subsd  %xmm6,%xmm4
add    $0x1,%r10d
jmp    129c <linear_expansion_sum_zeroelim+0x19c>
xchg   %ax,%ax
movslq %r9d,%rbx
movsd  (%rsi,%rbx,8),%xmm7
movapd %xmm0,%xmm5
addsd  %xmm7,%xmm5
movapd %xmm5,%xmm4
subsd  %xmm7,%xmm4
add    $0x1,%r9d
subsd  %xmm4,%xmm0
ucomisd %xmm3,%xmm0
jne    12ac <linear_expansion_sum_zeroelim+0x1ac>
jnp    11f0 <linear_expansion_sum_zeroelim+0xf0>
movslq %eax,%rbx
add    $0x1,%eax
movsd  %xmm0,(%r8,%rbx,8)
jmp    11f0 <linear_expansion_sum_zeroelim+0xf0>
movapd %xmm1,%xmm4
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12cd <linear_expansion_sum_zeroelim+0x1cd>
jnp    12d9 <linear_expansion_sum_zeroelim+0x1d9>
movslq %eax,%rcx
add    $0x1,%eax
movsd  %xmm0,(%r8,%rcx,8)
ucomisd %xmm1,%xmm4
jne    12e5 <linear_expansion_sum_zeroelim+0x1e5>
jp     12e5 <linear_expansion_sum_zeroelim+0x1e5>
test   %eax,%eax
jne    12f1 <linear_expansion_sum_zeroelim+0x1f1>
movslq %eax,%rcx
add    $0x1,%eax
movsd  %xmm4,(%r8,%rcx,8)
pop    %rbx
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @linear_expansion_sum_zeroelim(i32 %0, double* nocapture %1, i32 %2, double* nocapture %3, double* nocapture writeonly %4) {
%6 = load double, double* %3
%7 = load double, double* %1
%8 = fcmp ule double %6, %7
%9 = fneg double %7
%10 = fcmp ogt double %6, %9
%11 = xor i1 %8, %10
%12 = zext i1 %11 to i32
%13 = xor i1 %11, true
%14 = zext i1 %13 to i32
%15 = select i1 %11, double %7, double %6
%16 = icmp slt i32 %12, %0
br i1 %16, label %21, label %17
17: ; preds = %5
%18 = zext i1 %13 to i64
%19 = getelementptr inbounds double, double* %3, i64 %18
%20 = load double, double* %19
br label %43
21: ; preds = %5
%22 = icmp slt i32 %14, %2
br i1 %22, label %27, label %23
23: ; preds = %21
%24 = zext i1 %11 to i64
%25 = getelementptr inbounds double, double* %1, i64 %24
%26 = load double, double* %25
br label %38
27: ; preds = %21
%28 = zext i1 %13 to i64
%29 = getelementptr inbounds double, double* %3, i64 %28
%30 = load double, double* %29
%31 = zext i1 %11 to i64
%32 = getelementptr inbounds double, double* %1, i64 %31
%33 = load double, double* %32
%34 = fcmp ule double %30, %33
%35 = fneg double %33
%36 = fcmp ogt double %30, %35
%37 = xor i1 %34, %36
br i1 %37, label %38, label %43
38: ; preds = %23, %27
%39 = phi double [ %26, %23 ], [ %33, %27 ]
%40 = fadd double %15, %39
%41 = fsub double %40, %39
%42 = select i1 %11, i32 2, i32 1
br label %48
43: ; preds = %17, %27
%44 = phi double [ %20, %17 ], [ %30, %27 ]
%45 = fadd double %15, %44
%46 = fsub double %45, %44
%47 = select i1 %11, i32 1, i32 2
br label %48
48: ; preds = %43, %38
%49 = phi double [ %40, %38 ], [ %45, %43 ]
%50 = phi double [ %41, %38 ], [ %46, %43 ]
%51 = phi i32 [ %42, %38 ], [ %12, %43 ]
%52 = phi i32 [ %14, %38 ], [ %47, %43 ]
%53 = fsub double %15, %50
%54 = add i32 %2, %0
%55 = icmp sgt i32 %54, 2
br i1 %55, label %56, label %116
56: ; preds = %48, %112
%57 = phi i32 [ %114, %112 ], [ 2, %48 ]
%58 = phi i32 [ %113, %112 ], [ 0, %48 ]
%59 = phi i32 [ %99, %112 ], [ %52, %48 ]
%60 = phi i32 [ %98, %112 ], [ %51, %48 ]
%61 = phi double [ %101, %112 ], [ %49, %48 ]
%62 = phi double [ %106, %112 ], [ %53, %48 ]
%63 = icmp slt i32 %60, %0
br i1 %63, label %68, label %64
64: ; preds = %56
%65 = sext i32 %59 to i64
%66 = getelementptr inbounds double, double* %3, i64 %65
%67 = load double, double* %66
br label %90
68: ; preds = %56
%69 = icmp slt i32 %59, %2
br i1 %69, label %74, label %70
70: ; preds = %68
%71 = sext i32 %60 to i64
%72 = getelementptr inbounds double, double* %1, i64 %71
%73 = load double, double* %72
br label %85
74: ; preds = %68
%75 = sext i32 %59 to i64
%76 = getelementptr inbounds double, double* %3, i64 %75
%77 = load double, double* %76
%78 = sext i32 %60 to i64
%79 = getelementptr inbounds double, double* %1, i64 %78
%80 = load double, double* %79
%81 = fcmp ule double %77, %80
%82 = fneg double %80
%83 = fcmp ogt double %77, %82
%84 = xor i1 %81, %83
br i1 %84, label %85, label %90
85: ; preds = %70, %74
%86 = phi double [ %73, %70 ], [ %80, %74 ]
%87 = fadd double %62, %86
%88 = fsub double %87, %86
%89 = add nsw i32 %60, 1
br label %95
90: ; preds = %64, %74
%91 = phi double [ %67, %64 ], [ %77, %74 ]
%92 = fadd double %62, %91
%93 = fsub double %92, %91
%94 = add nsw i32 %59, 1
br label %95
95: ; preds = %90, %85
%96 = phi double [ %87, %85 ], [ %92, %90 ]
%97 = phi double [ %88, %85 ], [ %93, %90 ]
%98 = phi i32 [ %89, %85 ], [ %60, %90 ]
%99 = phi i32 [ %59, %85 ], [ %94, %90 ]
%100 = fsub double %62, %97
%101 = fadd double %61, %96
%102 = fsub double %101, %61
%103 = fsub double %101, %102
%104 = fsub double %96, %102
%105 = fsub double %61, %103
%106 = fadd double %104, %105
%107 = fcmp une double %100, 0.000000e+00
br i1 %107, label %108, label %112
108: ; preds = %95
%109 = add nsw i32 %58, 1
%110 = sext i32 %58 to i64
%111 = getelementptr inbounds double, double* %4, i64 %110
store double %100, double* %111
br label %112
112: ; preds = %95, %108
%113 = phi i32 [ %109, %108 ], [ %58, %95 ]
%114 = add nuw nsw i32 %57, 1
%115 = icmp eq i32 %114, %54
br i1 %115, label %116, label %56
116: ; preds = %112, %48
%117 = phi double [ %53, %48 ], [ %106, %112 ]
%118 = phi double [ %49, %48 ], [ %101, %112 ]
%119 = phi i32 [ 0, %48 ], [ %113, %112 ]
%120 = fcmp une double %117, 0.000000e+00
br i1 %120, label %121, label %125
121: ; preds = %116
%122 = add nsw i32 %119, 1
%123 = sext i32 %119 to i64
%124 = getelementptr inbounds double, double* %4, i64 %123
store double %117, double* %124
br label %125
125: ; preds = %121, %116
%126 = phi i32 [ %122, %121 ], [ %119, %116 ]
%127 = fcmp une double %118, 0.000000e+00
%128 = icmp eq i32 %126, 0
%129 = select i1 %127, i1 true, i1 %128
br i1 %129, label %130, label %134
130: ; preds = %125
%131 = add nsw i32 %126, 1
%132 = sext i32 %126 to i64
%133 = getelementptr inbounds double, double* %4, i64 %132
store double %118, double* %133
br label %134
134: ; preds = %125, %130
%135 = phi i32 [ %131, %130 ], [ %126, %125 ]
ret i32 %135
}