int grow_expansion(elen, e, b, h)
int elen;
float *e;
float b;
float *h;
{
  float Q;
  float Qnew;
  int eindex;
  float enow;
  float bvirt;
  float avirt, bround, around;
  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;
    h[eindex] = around + bround;
    Q = Qnew;
  }
  h[eindex] = Q;
  return eindex + 1;
}


x86 assembly (clang, -O0):
<grow_expansion>:
push   %rbp
mov    %rsp,%rbp
cvtsd2ss %xmm0,%xmm0
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11c1 <grow_expansion+0xc1>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x24(%rbp),%xmm0
addss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
subss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x28(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x34(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x24(%rbp),%xmm0
subss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1129 <grow_expansion+0x29>
movss  -0x24(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<grow_expansion>:
cvtsd2ss %xmm0,%xmm0
test   %edi,%edi
jle    1116 <grow_expansion+0x16>
mov    %edi,%eax
mov    %edi,%r8d
cmp    $0x1,%edi
jne    111d <grow_expansion+0x1d>
xor    %ecx,%ecx
jmp    1191 <grow_expansion+0x91>
xor    %eax,%eax
jmp    11c5 <grow_expansion+0xc5>
mov    %r8d,%edi
and    $0xfffffffe,%edi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rsi,%rcx,4),%xmm1
movaps %xmm0,%xmm2
addss  %xmm1,%xmm2
movaps %xmm2,%xmm3
subss  %xmm0,%xmm3
movaps %xmm2,%xmm4
subss  %xmm3,%xmm4
subss  %xmm3,%xmm1
subss  %xmm4,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rdx,%rcx,4)
movss  0x4(%rsi,%rcx,4),%xmm1
movaps %xmm2,%xmm0
addss  %xmm1,%xmm0
movaps %xmm0,%xmm3
subss  %xmm2,%xmm3
movaps %xmm0,%xmm4
subss  %xmm3,%xmm4
subss  %xmm3,%xmm1
subss  %xmm4,%xmm2
addss  %xmm1,%xmm2
movss  %xmm2,0x4(%rdx,%rcx,4)
add    $0x2,%rcx
cmp    %rcx,%rdi
jne    1130 <grow_expansion+0x30>
test   $0x1,%r8b
je     11c5 <grow_expansion+0xc5>
movss  (%rsi,%rcx,4),%xmm1
movaps %xmm0,%xmm2
addss  %xmm1,%xmm2
movaps %xmm2,%xmm3
subss  %xmm0,%xmm3
movaps %xmm2,%xmm4
subss  %xmm3,%xmm4
subss  %xmm3,%xmm1
subss  %xmm4,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rdx,%rcx,4)
movaps %xmm2,%xmm0
mov    %eax,%ecx
movss  %xmm0,(%rdx,%rcx,4)
add    $0x1,%eax
ret

LLVM IR (-O0):
define i32 @grow_expansion(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 i32
%12 = alloca float
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = fptrunc double %2 to float
store i32 %0, i32* %5
store float* %1, float** %6
store float %17, float* %7
store float* %3, float** %8
%18 = load float, float* %7
store float %18, float* %9
store i32 0, i32* %11
br label %19
19: ; preds = %52, %4
%20 = load i32, i32* %11
%21 = load i32, i32* %5
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %55
23: ; preds = %19
%24 = load float*, float** %6
%25 = load i32, i32* %11
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds float, float* %24, i64 %26
%28 = load float, float* %27
store float %28, float* %12
%29 = load float, float* %9
%30 = load float, float* %12
%31 = fadd float %29, %30
store float %31, float* %10
%32 = load float, float* %10
%33 = load float, float* %9
%34 = fsub float %32, %33
store float %34, float* %13
%35 = load float, float* %10
%36 = load float, float* %13
%37 = fsub float %35, %36
store float %37, float* %14
%38 = load float, float* %12
%39 = load float, float* %13
%40 = fsub float %38, %39
store float %40, float* %15
%41 = load float, float* %9
%42 = load float, float* %14
%43 = fsub float %41, %42
store float %43, float* %16
%44 = load float, float* %16
%45 = load float, float* %15
%46 = fadd float %44, %45
%47 = load float*, float** %8
%48 = load i32, i32* %11
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds float, float* %47, i64 %49
store float %46, float* %50
%51 = load float, float* %10
store float %51, float* %9
br label %52
52: ; preds = %23
%53 = load i32, i32* %11
%54 = add nsw i32 %53, 1
store i32 %54, i32* %11
br label %19
55: ; preds = %19
%56 = load float, float* %9
%57 = load float*, float** %8
%58 = load i32, i32* %11
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %57, i64 %59
store float %56, float* %60
%61 = load i32, i32* %11
%62 = add nsw i32 %61, 1
ret i32 %62
}

LLVM IR (-O2):
define i32 @grow_expansion(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 %54
7: ; preds = %4
%8 = zext i32 %0 to i64
%9 = and i64 %8, 1
%10 = icmp eq i32 %0, 1
br i1 %10, label %39, label %11
11: ; preds = %7
%12 = and i64 %8, 4294967294
br label %13
13: ; preds = %13, %11
%14 = phi i64 [ 0, %11 ], [ %36, %13 ]
%15 = phi float [ %5, %11 ], [ %29, %13 ]
%16 = phi i64 [ 0, %11 ], [ %37, %13 ]
%17 = getelementptr inbounds float, float* %1, i64 %14
%18 = load float, float* %17
%19 = fadd float %15, %18
%20 = fsub float %19, %15
%21 = fsub float %19, %20
%22 = fsub float %18, %20
%23 = fsub float %15, %21
%24 = fadd float %22, %23
%25 = getelementptr inbounds float, float* %3, i64 %14
store float %24, float* %25
%26 = or i64 %14, 1
%27 = getelementptr inbounds float, float* %1, i64 %26
%28 = load float, float* %27
%29 = fadd float %19, %28
%30 = fsub float %29, %19
%31 = fsub float %29, %30
%32 = fsub float %28, %30
%33 = fsub float %19, %31
%34 = fadd float %32, %33
%35 = getelementptr inbounds float, float* %3, i64 %26
store float %34, float* %35
%36 = add nuw nsw i64 %14, 2
%37 = add i64 %16, 2
%38 = icmp eq i64 %37, %12
br i1 %38, label %39, label %13
39: ; preds = %13, %7
%40 = phi float [ undef, %7 ], [ %29, %13 ]
%41 = phi i64 [ 0, %7 ], [ %36, %13 ]
%42 = phi float [ %5, %7 ], [ %29, %13 ]
%43 = icmp eq i64 %9, 0
br i1 %43, label %54, label %44
44: ; preds = %39
%45 = getelementptr inbounds float, float* %1, i64 %41
%46 = load float, float* %45
%47 = fadd float %42, %46
%48 = fsub float %47, %42
%49 = fsub float %47, %48
%50 = fsub float %46, %48
%51 = fsub float %42, %49
%52 = fadd float %50, %51
%53 = getelementptr inbounds float, float* %3, i64 %41
store float %52, float* %53
br label %54
54: ; preds = %44, %39, %4
%55 = phi i32 [ 0, %4 ], [ %0, %39 ], [ %0, %44 ]
%56 = phi float [ %5, %4 ], [ %40, %39 ], [ %47, %44 ]
%57 = zext i32 %55 to i64
%58 = getelementptr inbounds float, float* %3, i64 %57
store float %56, float* %58
%59 = add nuw nsw i32 %55, 1
ret i32 %59
}