int expansion_sum(elen, e, flen, f, h)
int elen;
float *e;
int flen;
float *f;
float *h;
{
  float Q;
  float Qnew;
  int findex, hindex, hlast;
  float hnow;
  float bvirt;
  float avirt, bround, around;
  Q = f[0];
  for (hindex = 0; hindex < elen; hindex++) {
    hnow = e[hindex];
    Qnew = (float)(Q + hnow);
    bvirt = (float)(Qnew - Q);
    avirt = Qnew - bvirt;
    bround = hnow - bvirt;
    around = Q - avirt;
    h[hindex] = around + bround;
    Q = Qnew;
  }
  h[hindex] = Q;
  hlast = hindex;
  for (findex = 1; findex < flen; findex++) {
    Q = f[findex];
    for (hindex = findex; hindex <= hlast; hindex++) {
      hnow = h[hindex];
      Qnew = (float)(Q + hnow);
      bvirt = (float)(Qnew - Q);
      avirt = Qnew - bvirt;
      bround = hnow - bvirt;
      around = Q - avirt;
      h[hindex] = around + bround;
      Q = Qnew;
    }
    h[++hlast] = Q;
  }
  return hlast + 1;
}


x86 assembly (clang, -O0):
<expansion_sum>:
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)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11c2 <expansion_sum+0xc2>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x2c(%rbp),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x40(%rbp),%xmm0
subss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
addss  -0x4c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    112a <expansion_sum+0x2a>
movss  -0x2c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12c5 <expansion_sum+0x1c5>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     129d <expansion_sum+0x19d>
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x2c(%rbp),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x30(%rbp),%xmm0
subss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x40(%rbp),%xmm0
subss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x2c(%rbp),%xmm0
subss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
addss  -0x4c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1205 <expansion_sum+0x105>
movss  -0x2c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x3c(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11e1 <expansion_sum+0xe1>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<expansion_sum>:
movss  (%rcx),%xmm0
test   %edi,%edi
jle    1116 <expansion_sum+0x16>
mov    %edi,%eax
mov    %edi,%r9d
cmp    $0x1,%edi
jne    111d <expansion_sum+0x1d>
xor    %edi,%edi
jmp    1193 <expansion_sum+0x93>
xor    %eax,%eax
jmp    11c8 <expansion_sum+0xc8>
mov    %r9d,%r10d
and    $0xfffffffe,%r10d
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movss  (%rsi,%rdi,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,(%r8,%rdi,4)
movss  0x4(%rsi,%rdi,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(%r8,%rdi,4)
add    $0x2,%rdi
cmp    %rdi,%r10
jne    1130 <expansion_sum+0x30>
test   $0x1,%r9b
je     11c8 <expansion_sum+0xc8>
movss  (%rsi,%rdi,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,(%r8,%rdi,4)
movaps %xmm2,%xmm0
mov    %eax,%esi
movss  %xmm0,(%r8,%rsi,4)
cmp    $0x2,%edx
jl     125c <expansion_sum+0x15c>
mov    %edx,%r9d
mov    $0x1,%edx
xor    %eax,%eax
jmp    120b <expansion_sum+0x10b>
cs nopw 0x0(%rax,%rax,1)
nop
movaps %xmm1,%xmm0
movss  %xmm0,0x4(%r8,%rsi,4)
add    $0x1,%rsi
add    $0x1,%rdx
add    $0x1,%rax
cmp    %r9,%rdx
je     125a <expansion_sum+0x15a>
movss  (%rcx,%rdx,4),%xmm1
cmp    %rsi,%rdx
ja     11f0 <expansion_sum+0xf0>
mov    %rax,%rdi
movaps %xmm1,%xmm0
nopl   0x0(%rax,%rax,1)
movss  0x4(%r8,%rdi,4),%xmm2
addss  %xmm2,%xmm0
movaps %xmm0,%xmm3
subss  %xmm1,%xmm3
movaps %xmm0,%xmm4
subss  %xmm3,%xmm4
subss  %xmm3,%xmm2
subss  %xmm4,%xmm1
addss  %xmm2,%xmm1
movss  %xmm1,0x4(%r8,%rdi,4)
add    $0x1,%rdi
movaps %xmm0,%xmm1
cmp    %rsi,%rdi
jb     1220 <expansion_sum+0x120>
jmp    11f3 <expansion_sum+0xf3>
mov    %esi,%eax
add    $0x1,%eax
ret

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

LLVM IR (-O2):
define i32 @expansion_sum(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 %55
8: ; preds = %5
%9 = zext i32 %0 to i64
%10 = and i64 %9, 1
%11 = icmp eq i32 %0, 1
br i1 %11, label %40, label %12
12: ; preds = %8
%13 = and i64 %9, 4294967294
br label %14
14: ; preds = %14, %12
%15 = phi i64 [ 0, %12 ], [ %37, %14 ]
%16 = phi float [ %6, %12 ], [ %30, %14 ]
%17 = phi i64 [ 0, %12 ], [ %38, %14 ]
%18 = getelementptr inbounds float, float* %1, i64 %15
%19 = load float, float* %18
%20 = fadd float %16, %19
%21 = fsub float %20, %16
%22 = fsub float %20, %21
%23 = fsub float %19, %21
%24 = fsub float %16, %22
%25 = fadd float %23, %24
%26 = getelementptr inbounds float, float* %4, i64 %15
store float %25, float* %26
%27 = or i64 %15, 1
%28 = getelementptr inbounds float, float* %1, i64 %27
%29 = load float, float* %28
%30 = fadd float %20, %29
%31 = fsub float %30, %20
%32 = fsub float %30, %31
%33 = fsub float %29, %31
%34 = fsub float %20, %32
%35 = fadd float %33, %34
%36 = getelementptr inbounds float, float* %4, i64 %27
store float %35, float* %36
%37 = add nuw nsw i64 %15, 2
%38 = add i64 %17, 2
%39 = icmp eq i64 %38, %13
br i1 %39, label %40, label %14
40: ; preds = %14, %8
%41 = phi float [ undef, %8 ], [ %30, %14 ]
%42 = phi i64 [ 0, %8 ], [ %37, %14 ]
%43 = phi float [ %6, %8 ], [ %30, %14 ]
%44 = icmp eq i64 %10, 0
br i1 %44, label %55, label %45
45: ; preds = %40
%46 = getelementptr inbounds float, float* %1, i64 %42
%47 = load float, float* %46
%48 = fadd float %43, %47
%49 = fsub float %48, %43
%50 = fsub float %48, %49
%51 = fsub float %47, %49
%52 = fsub float %43, %50
%53 = fadd float %51, %52
%54 = getelementptr inbounds float, float* %4, i64 %42
store float %53, float* %54
br label %55
55: ; preds = %45, %40, %5
%56 = phi float [ %6, %5 ], [ %41, %40 ], [ %48, %45 ]
%57 = phi i32 [ 0, %5 ], [ %0, %40 ], [ %0, %45 ]
%58 = zext i32 %57 to i64
%59 = getelementptr inbounds float, float* %4, i64 %58
store float %56, float* %59
%60 = icmp sgt i32 %2, 1
br i1 %60, label %61, label %90
61: ; preds = %55
%62 = zext i32 %2 to i64
br label %63
63: ; preds = %61, %82
%64 = phi i64 [ %58, %61 ], [ %84, %82 ]
%65 = phi i64 [ 1, %61 ], [ %86, %82 ]
%66 = getelementptr inbounds float, float* %3, i64 %65
%67 = load float, float* %66
%68 = icmp ugt i64 %65, %64
br i1 %68, label %82, label %69
69: ; preds = %63, %69
%70 = phi i64 [ %80, %69 ], [ %65, %63 ]
%71 = phi float [ %74, %69 ], [ %67, %63 ]
%72 = getelementptr inbounds float, float* %4, i64 %70
%73 = load float, float* %72
%74 = fadd float %71, %73
%75 = fsub float %74, %71
%76 = fsub float %74, %75
%77 = fsub float %73, %75
%78 = fsub float %71, %76
%79 = fadd float %77, %78
store float %79, float* %72
%80 = add nuw nsw i64 %70, 1
%81 = icmp ult i64 %70, %64
br i1 %81, label %69, label %82
82: ; preds = %69, %63
%83 = phi float [ %67, %63 ], [ %74, %69 ]
%84 = add nuw nsw i64 %64, 1
%85 = getelementptr inbounds float, float* %4, i64 %84
store float %83, float* %85
%86 = add nuw nsw i64 %65, 1
%87 = icmp eq i64 %86, %62
br i1 %87, label %88, label %63
88: ; preds = %82
%89 = trunc i64 %84 to i32
br label %90
90: ; preds = %88, %55
%91 = phi i32 [ %57, %55 ], [ %89, %88 ]
%92 = add nuw nsw i32 %91, 1
ret i32 %92
}