int expansion_sum_zeroelim2(elen, e, flen, f, h)
int elen;
float *e;
int flen;
float *f;
float *h;
{
  float Q, hh;
  float Qnew;
  int eindex, findex, hindex, hlast;
  float enow;
  float bvirt;
  float avirt, bround, around;
  hindex = 0;
  Q = f[0];
  for (eindex = 0; eindex < elen; eindex++) {
    enow = e[eindex];
    Qnew = (float)(Q + enow);
    bvirt = (float)(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 = (float)(Q + enow);
      bvirt = (float)(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,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1203 <expansion_sum_zeroelim2+0x103>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x2c(%rbp),%xmm0
addss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x34(%rbp),%xmm0
subss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x48(%rbp),%xmm0
subss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x58(%rbp),%xmm0
addss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11d4 <expansion_sum_zeroelim2+0xd4>
jp     11d4 <expansion_sum_zeroelim2+0xd4>
jmp    11f0 <expansion_sum_zeroelim2+0xf0>
movss  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x40(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    11f5 <expansion_sum_zeroelim2+0xf5>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1131 <expansion_sum_zeroelim2+0x31>
movss  -0x2c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
mov    %eax,-0x44(%rbp)
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1341 <expansion_sum_zeroelim2+0x241>
movl   $0x0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     131b <expansion_sum_zeroelim2+0x21b>
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x2c(%rbp),%xmm0
addss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x34(%rbp),%xmm0
subss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x48(%rbp),%xmm0
subss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x58(%rbp),%xmm0
addss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    12ec <expansion_sum_zeroelim2+0x1ec>
jp     12ec <expansion_sum_zeroelim2+0x1ec>
jmp    1308 <expansion_sum_zeroelim2+0x208>
movss  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x40(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    130d <expansion_sum_zeroelim2+0x20d>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    124e <expansion_sum_zeroelim2+0x14e>
movss  -0x2c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1222 <expansion_sum_zeroelim2+0x122>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<expansion_sum_zeroelim2>:
movss  (%rcx),%xmm0
test   %edi,%edi
jle    1164 <expansion_sum_zeroelim2+0x64>
mov    %edi,%r9d
xor    %edi,%edi
xorps  %xmm1,%xmm1
xor    %eax,%eax
jmp    1129 <expansion_sum_zeroelim2+0x29>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rdi
cmp    %rdi,%r9
je     1166 <expansion_sum_zeroelim2+0x66>
movaps %xmm0,%xmm2
movss  (%rsi,%rdi,4),%xmm3
addss  %xmm3,%xmm0
movaps %xmm0,%xmm4
subss  %xmm2,%xmm4
movaps %xmm0,%xmm5
subss  %xmm4,%xmm5
subss  %xmm4,%xmm3
subss  %xmm5,%xmm2
addss  %xmm3,%xmm2
ucomiss %xmm1,%xmm2
jne    1156 <expansion_sum_zeroelim2+0x56>
jnp    1120 <expansion_sum_zeroelim2+0x20>
movslq %eax,%r10
add    $0x1,%eax
movss  %xmm2,(%r8,%r10,4)
jmp    1120 <expansion_sum_zeroelim2+0x20>
xor    %eax,%eax
movslq %eax,%rsi
movss  %xmm0,(%r8,%rsi,4)
cmp    $0x2,%edx
jl     1205 <expansion_sum_zeroelim2+0x105>
mov    %edx,%r9d
mov    $0x1,%r10d
xorps  %xmm0,%xmm0
jmp    11a4 <expansion_sum_zeroelim2+0xa4>
cs nopw 0x0(%rax,%rax,1)
xor    %eax,%eax
movslq %eax,%rdx
movss  %xmm1,(%r8,%rdx,4)
add    $0x1,%r10
cmp    %r9,%r10
je     1205 <expansion_sum_zeroelim2+0x105>
movss  (%rcx,%r10,4),%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>
movaps %xmm1,%xmm2
movss  (%r8,%rdx,4),%xmm3
addss  %xmm3,%xmm1
movaps %xmm1,%xmm4
subss  %xmm2,%xmm4
movaps %xmm1,%xmm5
subss  %xmm4,%xmm5
subss  %xmm4,%xmm3
subss  %xmm5,%xmm2
addss  %xmm3,%xmm2
ucomiss %xmm0,%xmm2
jne    11f7 <expansion_sum_zeroelim2+0xf7>
jnp    11c0 <expansion_sum_zeroelim2+0xc0>
movslq %eax,%rsi
add    $0x1,%eax
movss  %xmm2,(%r8,%rsi,4)
jmp    11c0 <expansion_sum_zeroelim2+0xc0>
add    $0x1,%eax
ret

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

LLVM IR (-O2):
define i32 @expansion_sum_zeroelim2(i32 %0, float* nocapture %1, i32 %2, float* nocapture %3, float* nocapture %4) {
%6 = load float, float* %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 float [ %6, %8 ], [ %16, %27 ]
%14 = getelementptr inbounds float, float* %1, i64 %11
%15 = load float, float* %14
%16 = fadd float %13, %15
%17 = fsub float %16, %13
%18 = fsub float %16, %17
%19 = fsub float %15, %17
%20 = fsub float %13, %18
%21 = fadd float %19, %20
%22 = fcmp une float %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 float, float* %4, i64 %25
store float %21, float* %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 float [ %6, %5 ], [ %16, %27 ]
%33 = phi i32 [ 0, %5 ], [ %28, %27 ]
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds float, float* %4, i64 %34
store float %32, float* %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 float, float* %3, i64 %40
%43 = load float, float* %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 float [ %43, %45 ], [ %54, %65 ]
%52 = getelementptr inbounds float, float* %4, i64 %49
%53 = load float, float* %52
%54 = fadd float %51, %53
%55 = fsub float %54, %51
%56 = fsub float %54, %55
%57 = fsub float %53, %55
%58 = fsub float %51, %56
%59 = fadd float %57, %58
%60 = fcmp une float %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 float, float* %4, i64 %63
store float %59, float* %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 float [ %43, %39 ], [ %54, %65 ]
%71 = phi i32 [ 0, %39 ], [ %66, %65 ]
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds float, float* %4, i64 %72
store float %70, float* %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
}