int expansion_sum_zeroelim2(int elen, double *e, int flen, double *f,
                            double *h) {
  double Q, hh;
  double Qnew;
  int eindex, findex, hindex, hlast;
  double enow;
  double bvirt;
  double avirt, bround, around;
  hindex = 0;
  Q = f[0];
  for (eindex = 0; eindex < elen; eindex++) {
    enow = e[eindex];
    Qnew = (double)(Q + enow);
    bvirt = (double)(Qnew - Q);
    avirt = Qnew - bvirt;
    bround = enow - bvirt;
    around = Q - avirt;
    hh = around + bround;
    Q = Qnew;
    if (hh != 0.0) {
      h[hindex++] = hh;
    }
  }
  h[hindex] = Q;
  hlast = hindex;
  for (findex = 1; findex < flen; findex++) {
    hindex = 0;
    Q = f[findex];
    for (eindex = 0; eindex <= hlast; eindex++) {
      enow = h[eindex];
      Qnew = (double)(Q + enow);
      bvirt = (double)(Qnew - Q);
      avirt = Qnew - bvirt;
      bround = enow - bvirt;
      around = Q - avirt;
      hh = around + bround;
      Q = Qnew;
      if (hh != 0) {
        h[hindex++] = hh;
      }
    }
    h[hindex] = Q;
    hlast = hindex;
  }
  return hlast + 1;
}


x86 assembly (clang, -O0):
<expansion_sum_zeroelim2>:
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)
movl   $0x0,-0x4c(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11ff <expansion_sum_zeroelim2+0xff>
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11d0 <expansion_sum_zeroelim2+0xd0>
jp     11d0 <expansion_sum_zeroelim2+0xd0>
jmp    11ec <expansion_sum_zeroelim2+0xec>
movsd  -0x38(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x4c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x4c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    11f1 <expansion_sum_zeroelim2+0xf1>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1131 <expansion_sum_zeroelim2+0x31>
movsd  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x50(%rbp)
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    133e <expansion_sum_zeroelim2+0x23e>
movl   $0x0,-0x4c(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x50(%rbp),%eax
jg     1318 <expansion_sum_zeroelim2+0x218>
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x30(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12e9 <expansion_sum_zeroelim2+0x1e9>
jp     12e9 <expansion_sum_zeroelim2+0x1e9>
jmp    1305 <expansion_sum_zeroelim2+0x205>
movsd  -0x38(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x4c(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x4c(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    130a <expansion_sum_zeroelim2+0x20a>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    124a <expansion_sum_zeroelim2+0x14a>
movsd  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    121e <expansion_sum_zeroelim2+0x11e>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<expansion_sum_zeroelim2>:
movsd  (%rcx),%xmm0
test   %edi,%edi
jle    1168 <expansion_sum_zeroelim2+0x68>
mov    %edi,%r9d
xor    %edi,%edi
xorpd  %xmm1,%xmm1
xor    %eax,%eax
jmp    1129 <expansion_sum_zeroelim2+0x29>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rdi
cmp    %rdi,%r9
je     116a <expansion_sum_zeroelim2+0x6a>
movapd %xmm0,%xmm2
movsd  (%rsi,%rdi,8),%xmm3
addsd  %xmm3,%xmm0
movapd %xmm0,%xmm4
subsd  %xmm2,%xmm4
movapd %xmm0,%xmm5
subsd  %xmm4,%xmm5
subsd  %xmm4,%xmm3
subsd  %xmm5,%xmm2
addsd  %xmm3,%xmm2
ucomisd %xmm1,%xmm2
jne    115a <expansion_sum_zeroelim2+0x5a>
jnp    1120 <expansion_sum_zeroelim2+0x20>
movslq %eax,%r10
add    $0x1,%eax
movsd  %xmm2,(%r8,%r10,8)
jmp    1120 <expansion_sum_zeroelim2+0x20>
xor    %eax,%eax
movslq %eax,%rsi
movsd  %xmm0,(%r8,%rsi,8)
cmp    $0x2,%edx
jl     1209 <expansion_sum_zeroelim2+0x109>
mov    %edx,%r9d
mov    $0x1,%r10d
xorpd  %xmm0,%xmm0
jmp    11a4 <expansion_sum_zeroelim2+0xa4>
nopl   0x0(%rax,%rax,1)
xor    %eax,%eax
movslq %eax,%rdx
movsd  %xmm1,(%r8,%rdx,8)
add    $0x1,%r10
cmp    %r9,%r10
je     1209 <expansion_sum_zeroelim2+0x109>
movsd  (%rcx,%r10,8),%xmm1
test   %eax,%eax
js     1190 <expansion_sum_zeroelim2+0x90>
mov    %eax,%edi
add    $0x1,%edi
xor    %edx,%edx
xor    %eax,%eax
jmp    11c9 <expansion_sum_zeroelim2+0xc9>
nopl   0x0(%rax)
add    $0x1,%rdx
cmp    %rdx,%rdi
je     1192 <expansion_sum_zeroelim2+0x92>
movapd %xmm1,%xmm2
movsd  (%r8,%rdx,8),%xmm3
addsd  %xmm3,%xmm1
movapd %xmm1,%xmm4
subsd  %xmm2,%xmm4
movapd %xmm1,%xmm5
subsd  %xmm4,%xmm5
subsd  %xmm4,%xmm3
subsd  %xmm5,%xmm2
addsd  %xmm3,%xmm2
ucomisd %xmm0,%xmm2
jne    11fb <expansion_sum_zeroelim2+0xfb>
jnp    11c0 <expansion_sum_zeroelim2+0xc0>
movslq %eax,%rsi
add    $0x1,%eax
movsd  %xmm2,(%r8,%rsi,8)
jmp    11c0 <expansion_sum_zeroelim2+0xc0>
add    $0x1,%eax
ret

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

LLVM IR (-O2):
define i32 @expansion_sum_zeroelim2(i32 %0, double* nocapture %1, i32 %2, double* nocapture %3, double* nocapture %4) {
%6 = load double, double* %3
%7 = icmp sgt i32 %0, 0
br i1 %7, label %8, label %31
8: ; preds = %5
%9 = zext i32 %0 to i64
br label %10
10: ; preds = %8, %27
%11 = phi i64 [ 0, %8 ], [ %29, %27 ]
%12 = phi i32 [ 0, %8 ], [ %28, %27 ]
%13 = phi double [ %6, %8 ], [ %16, %27 ]
%14 = getelementptr inbounds double, double* %1, i64 %11
%15 = load double, double* %14
%16 = fadd double %13, %15
%17 = fsub double %16, %13
%18 = fsub double %16, %17
%19 = fsub double %15, %17
%20 = fsub double %13, %18
%21 = fadd double %19, %20
%22 = fcmp une double %21, 0.000000e+00
br i1 %22, label %23, label %27
23: ; preds = %10
%24 = add nsw i32 %12, 1
%25 = sext i32 %12 to i64
%26 = getelementptr inbounds double, double* %4, i64 %25
store double %21, double* %26
br label %27
27: ; preds = %10, %23
%28 = phi i32 [ %24, %23 ], [ %12, %10 ]
%29 = add nuw nsw i64 %11, 1
%30 = icmp eq i64 %29, %9
br i1 %30, label %31, label %10
31: ; preds = %27, %5
%32 = phi double [ %6, %5 ], [ %16, %27 ]
%33 = phi i32 [ 0, %5 ], [ %28, %27 ]
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %4, i64 %34
store double %32, double* %35
%36 = icmp sgt i32 %2, 1
br i1 %36, label %37, label %76
37: ; preds = %31
%38 = zext i32 %2 to i64
br label %39
39: ; preds = %37, %69
%40 = phi i64 [ 1, %37 ], [ %74, %69 ]
%41 = phi i32 [ %33, %37 ], [ %71, %69 ]
%42 = getelementptr inbounds double, double* %3, i64 %40
%43 = load double, double* %42
%44 = icmp slt i32 %41, 0
br i1 %44, label %69, label %45
45: ; preds = %39
%46 = add nuw i32 %41, 1
%47 = zext i32 %46 to i64
br label %48
48: ; preds = %45, %65
%49 = phi i64 [ 0, %45 ], [ %67, %65 ]
%50 = phi i32 [ 0, %45 ], [ %66, %65 ]
%51 = phi double [ %43, %45 ], [ %54, %65 ]
%52 = getelementptr inbounds double, double* %4, i64 %49
%53 = load double, double* %52
%54 = fadd double %51, %53
%55 = fsub double %54, %51
%56 = fsub double %54, %55
%57 = fsub double %53, %55
%58 = fsub double %51, %56
%59 = fadd double %57, %58
%60 = fcmp une double %59, 0.000000e+00
br i1 %60, label %61, label %65
61: ; preds = %48
%62 = add nsw i32 %50, 1
%63 = sext i32 %50 to i64
%64 = getelementptr inbounds double, double* %4, i64 %63
store double %59, double* %64
br label %65
65: ; preds = %48, %61
%66 = phi i32 [ %62, %61 ], [ %50, %48 ]
%67 = add nuw nsw i64 %49, 1
%68 = icmp eq i64 %67, %47
br i1 %68, label %69, label %48
69: ; preds = %65, %39
%70 = phi double [ %43, %39 ], [ %54, %65 ]
%71 = phi i32 [ 0, %39 ], [ %66, %65 ]
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %4, i64 %72
store double %70, double* %73
%74 = add nuw nsw i64 %40, 1
%75 = icmp eq i64 %74, %38
br i1 %75, label %76, label %39
76: ; preds = %69, %31
%77 = phi i32 [ %33, %31 ], [ %71, %69 ]
%78 = add nsw i32 %77, 1
ret i32 %78
}