int grow_expansion_zeroelim(elen, e, b, h)
int elen;
float *e;
float b;
float *h;
{
  float Q, hh;
  float Qnew;
  int eindex, hindex;
  float enow;
  float bvirt;
  float avirt, bround, around;
  hindex = 0;
  Q = b;
  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;
    }
  }
  if ((Q != 0.0) || (hindex == 0)) {
    h[hindex++] = Q;
  }
  return hindex;
}


x86 assembly (clang, -O0):
<grow_expansion_zeroelim>:
push   %rbp
mov    %rsp,%rbp
cvtsd2ss %xmm0,%xmm0
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
movl   $0x0,-0x34(%rbp)
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1202 <grow_expansion_zeroelim+0x102>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x24(%rbp),%xmm0
addss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x38(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x24(%rbp),%xmm0
subss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11d3 <grow_expansion_zeroelim+0xd3>
jp     11d3 <grow_expansion_zeroelim+0xd3>
jmp    11ef <grow_expansion_zeroelim+0xef>
movss  -0x28(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x34(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    11f4 <grow_expansion_zeroelim+0xf4>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1130 <grow_expansion_zeroelim+0x30>
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1228 <grow_expansion_zeroelim+0x128>
jp     1228 <grow_expansion_zeroelim+0x128>
cmpl   $0x0,-0x34(%rbp)
jne    1244 <grow_expansion_zeroelim+0x144>
movss  -0x24(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x34(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<grow_expansion_zeroelim>:
cvtsd2ss %xmm0,%xmm0
test   %edi,%edi
jle    1172 <grow_expansion_zeroelim+0x72>
mov    %edi,%r8d
xor    %edi,%edi
xorps  %xmm1,%xmm1
xor    %eax,%eax
jmp    1129 <grow_expansion_zeroelim+0x29>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rdi
cmp    %rdi,%r8
je     1163 <grow_expansion_zeroelim+0x63>
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 <grow_expansion_zeroelim+0x56>
jnp    1120 <grow_expansion_zeroelim+0x20>
movslq %eax,%rcx
add    $0x1,%eax
movss  %xmm2,(%rdx,%rcx,4)
jmp    1120 <grow_expansion_zeroelim+0x20>
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1174 <grow_expansion_zeroelim+0x74>
jp     1174 <grow_expansion_zeroelim+0x74>
test   %eax,%eax
je     1174 <grow_expansion_zeroelim+0x74>
ret
xor    %eax,%eax
movslq %eax,%rcx
add    $0x1,%eax
movss  %xmm0,(%rdx,%rcx,4)
ret

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

LLVM IR (-O2):
define i32 @grow_expansion_zeroelim(i32 %0, float* nocapture %1, double %2, float* nocapture writeonly %3) {
%5 = fptrunc double %2 to float
%6 = icmp sgt i32 %0, 0
br i1 %6, label %7, label %34
7: ; preds = %4
%8 = zext i32 %0 to i64
br label %9
9: ; preds = %7, %26
%10 = phi i64 [ 0, %7 ], [ %28, %26 ]
%11 = phi float [ %5, %7 ], [ %15, %26 ]
%12 = phi i32 [ 0, %7 ], [ %27, %26 ]
%13 = getelementptr inbounds float, float* %1, i64 %10
%14 = load float, float* %13
%15 = fadd float %11, %14
%16 = fsub float %15, %11
%17 = fsub float %15, %16
%18 = fsub float %14, %16
%19 = fsub float %11, %17
%20 = fadd float %18, %19
%21 = fcmp une float %20, 0.000000e+00
br i1 %21, label %22, label %26
22: ; preds = %9
%23 = add nsw i32 %12, 1
%24 = sext i32 %12 to i64
%25 = getelementptr inbounds float, float* %3, i64 %24
store float %20, float* %25
br label %26
26: ; preds = %9, %22
%27 = phi i32 [ %23, %22 ], [ %12, %9 ]
%28 = add nuw nsw i64 %10, 1
%29 = icmp eq i64 %28, %8
br i1 %29, label %30, label %9
30: ; preds = %26
%31 = fcmp une float %15, 0.000000e+00
%32 = icmp eq i32 %27, 0
%33 = select i1 %31, i1 true, i1 %32
br i1 %33, label %34, label %40
34: ; preds = %4, %30
%35 = phi float [ %15, %30 ], [ %5, %4 ]
%36 = phi i32 [ %27, %30 ], [ 0, %4 ]
%37 = add nsw i32 %36, 1
%38 = sext i32 %36 to i64
%39 = getelementptr inbounds float, float* %3, i64 %38
store float %35, float* %39
br label %40
40: ; preds = %30, %34
%41 = phi i32 [ %37, %34 ], [ %27, %30 ]
ret i32 %41
}