int p_compress(elen, e, h)
int elen;
double *e;
double *h;
{
  double Q, q;
  double Qnew;
  int eindex, hindex;
  double bvirt;
  double enow, hnow;
  int top, bottom;
  bottom = elen - 1;
  Q = e[bottom];
  for (eindex = elen - 2; eindex >= 0; eindex--) {
    enow = e[eindex];
    Qnew = (double)(Q + enow);
    bvirt = Qnew - Q;
    q = enow - bvirt;
    if (q != 0) {
      h[bottom--] = Qnew;
      Q = q;
    } else {
      Q = Qnew;
    }
  }
  top = 0;
  for (hindex = bottom + 1; hindex < elen; hindex++) {
    hnow = h[hindex];
    Qnew = (double)(hnow + Q);
    bvirt = Qnew - hnow;
    q = Q - bvirt;
    if (q != 0) {
      h[top++] = q;
    }
    Q = Qnew;
  }
  h[top] = Q;
  return top + 1;
}


x86 assembly (clang, -O0):
<p_compress>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x58(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jl     11e1 <p_compress+0xe1>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x20(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1199 <p_compress+0x99>
jp     1199 <p_compress+0x99>
jmp    11c4 <p_compress+0xc4>
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x58(%rbp),%ecx
mov    %ecx,%edx
add    $0xffffffff,%edx
mov    %edx,-0x58(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11ce <p_compress+0xce>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11d3 <p_compress+0xd3>
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    1133 <p_compress+0x33>
movl   $0x0,-0x54(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    128d <p_compress+0x18d>
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x20(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1259 <p_compress+0x159>
jp     1259 <p_compress+0x159>
jmp    1275 <p_compress+0x175>
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x54(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x54(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11f1 <p_compress+0xf1>
movsd  -0x20(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<p_compress>:
lea    -0x1(%rdi),%r9d
movslq %edi,%rcx
movsd  -0x8(%rsi,%rcx,8),%xmm2
cmp    $0x2,%ecx
jl     114f <p_compress+0x4f>
lea    -0x2(%rdi),%r8d
test   $0x1,%r8b
jne    1158 <p_compress+0x58>
movsd  (%rsi,%r8,8),%xmm0
movapd %xmm2,%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm3
subsd  %xmm2,%xmm3
subsd  %xmm3,%xmm0
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
jne    1142 <p_compress+0x42>
jnp    1165 <p_compress+0x65>
movslq %r9d,%rax
movsd  %xmm1,(%rdx,%rax,8)
mov    %r8d,%r9d
jmp    1169 <p_compress+0x69>
movapd %xmm2,%xmm0
jmp    120c <p_compress+0x10c>
mov    %r8,%rcx
test   %r8,%r8
jne    117a <p_compress+0x7a>
jmp    120c <p_compress+0x10c>
movapd %xmm1,%xmm0
lea    -0x1(%r8),%rcx
movapd %xmm0,%xmm2
test   %r8,%r8
je     120c <p_compress+0x10c>
add    $0x2,%rcx
xorpd  %xmm1,%xmm1
movapd %xmm2,%xmm0
jmp    11a6 <p_compress+0xa6>
nopl   0x0(%rax,%rax,1)
movslq %r9d,%rax
add    $0xffffffff,%r9d
movsd  %xmm3,(%rdx,%rax,8)
add    $0xfffffffffffffffe,%rcx
cmp    $0x1,%rcx
jle    120c <p_compress+0x10c>
movsd  -0x10(%rsi,%rcx,8),%xmm2
movapd %xmm0,%xmm3
addsd  %xmm2,%xmm3
movapd %xmm3,%xmm4
subsd  %xmm0,%xmm4
subsd  %xmm4,%xmm2
ucomisd %xmm1,%xmm2
jne    11c8 <p_compress+0xc8>
jnp    11e0 <p_compress+0xe0>
movslq %r9d,%rax
add    $0xffffffff,%r9d
movsd  %xmm3,(%rdx,%rax,8)
jmp    11e4 <p_compress+0xe4>
cs nopw 0x0(%rax,%rax,1)
movapd %xmm3,%xmm2
movsd  -0x18(%rsi,%rcx,8),%xmm0
movapd %xmm2,%xmm3
addsd  %xmm0,%xmm3
movapd %xmm3,%xmm4
subsd  %xmm2,%xmm4
subsd  %xmm4,%xmm0
ucomisd %xmm1,%xmm0
jne    1190 <p_compress+0x90>
jp     1190 <p_compress+0x90>
movapd %xmm3,%xmm0
jmp    119c <p_compress+0x9c>
lea    0x1(%r9),%esi
xor    %ecx,%ecx
cmp    %edi,%esi
jge    1268 <p_compress+0x168>
movslq %esi,%rsi
mov    %r9d,%r8d
not    %r8d
add    %edi,%r8d
add    $0xfffffffe,%edi
cmp    %r9d,%edi
jne    1278 <p_compress+0x178>
xor    %ecx,%ecx
movapd %xmm0,%xmm1
test   $0x1,%r8b
je     126c <p_compress+0x16c>
movsd  (%rdx,%rsi,8),%xmm2
movapd %xmm0,%xmm1
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm3
subsd  %xmm2,%xmm3
subsd  %xmm3,%xmm0
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
jne    125b <p_compress+0x15b>
jnp    126c <p_compress+0x16c>
movslq %ecx,%rax
add    $0x1,%ecx
movsd  %xmm0,(%rdx,%rax,8)
jmp    126c <p_compress+0x16c>
movapd %xmm0,%xmm1
movslq %ecx,%rax
movsd  %xmm1,(%rdx,%rax,8)
add    $0x1,%eax
ret
mov    %r8d,%eax
and    $0xfffffffe,%eax
xor    %ecx,%ecx
xorpd  %xmm1,%xmm1
jmp    1299 <p_compress+0x199>
cs nopw 0x0(%rax,%rax,1)
add    $0x2,%rsi
add    $0xfffffffe,%eax
je     122c <p_compress+0x12c>
movsd  (%rdx,%rsi,8),%xmm3
movapd %xmm0,%xmm2
addsd  %xmm3,%xmm2
movapd %xmm2,%xmm4
subsd  %xmm3,%xmm4
subsd  %xmm4,%xmm0
ucomisd %xmm1,%xmm0
jne    12ba <p_compress+0x1ba>
jnp    12c5 <p_compress+0x1c5>
movslq %ecx,%rdi
add    $0x1,%ecx
movsd  %xmm0,(%rdx,%rdi,8)
movsd  0x8(%rdx,%rsi,8),%xmm3
movapd %xmm2,%xmm0
addsd  %xmm3,%xmm0
movapd %xmm0,%xmm4
subsd  %xmm3,%xmm4
subsd  %xmm4,%xmm2
ucomisd %xmm1,%xmm2
jne    12e7 <p_compress+0x1e7>
jnp    1290 <p_compress+0x190>
movslq %ecx,%rdi
add    $0x1,%ecx
movsd  %xmm2,(%rdx,%rdi,8)
jmp    1290 <p_compress+0x190>

LLVM IR (-O0):
define i32 @p_compress(i32 %0, double* %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca double
%8 = alloca double
%9 = alloca double
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca i32
%16 = alloca i32
store i32 %0, i32* %4
store double* %1, double** %5
store double* %2, double** %6
%17 = load i32, i32* %4
%18 = sub nsw i32 %17, 1
store i32 %18, i32* %16
%19 = load double*, double** %5
%20 = load i32, i32* %16
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double, double* %19, i64 %21
%23 = load double, double* %22
store double %23, double* %7
%24 = load i32, i32* %4
%25 = sub nsw i32 %24, 2
store i32 %25, i32* %10
br label %26
26: ; preds = %57, %3
%27 = load i32, i32* %10
%28 = icmp sge i32 %27, 0
br i1 %28, label %29, label %60
29: ; preds = %26
%30 = load double*, double** %5
%31 = load i32, i32* %10
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
%34 = load double, double* %33
store double %34, double* %13
%35 = load double, double* %7
%36 = load double, double* %13
%37 = fadd double %35, %36
store double %37, double* %9
%38 = load double, double* %9
%39 = load double, double* %7
%40 = fsub double %38, %39
store double %40, double* %12
%41 = load double, double* %13
%42 = load double, double* %12
%43 = fsub double %41, %42
store double %43, double* %8
%44 = load double, double* %8
%45 = fcmp une double %44, 0.000000e+00
br i1 %45, label %46, label %54
46: ; preds = %29
%47 = load double, double* %9
%48 = load double*, double** %6
%49 = load i32, i32* %16
%50 = add nsw i32 %49, -1
store i32 %50, i32* %16
%51 = sext i32 %49 to i64
%52 = getelementptr inbounds double, double* %48, i64 %51
store double %47, double* %52
%53 = load double, double* %8
store double %53, double* %7
br label %56
54: ; preds = %29
%55 = load double, double* %9
store double %55, double* %7
br label %56
56: ; preds = %54, %46
br label %57
57: ; preds = %56
%58 = load i32, i32* %10
%59 = add nsw i32 %58, -1
store i32 %59, i32* %10
br label %26
60: ; preds = %26
store i32 0, i32* %15
%61 = load i32, i32* %16
%62 = add nsw i32 %61, 1
store i32 %62, i32* %11
br label %63
63: ; preds = %93, %60
%64 = load i32, i32* %11
%65 = load i32, i32* %4
%66 = icmp slt i32 %64, %65
br i1 %66, label %67, label %96
67: ; preds = %63
%68 = load double*, double** %6
%69 = load i32, i32* %11
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double, double* %68, i64 %70
%72 = load double, double* %71
store double %72, double* %14
%73 = load double, double* %14
%74 = load double, double* %7
%75 = fadd double %73, %74
store double %75, double* %9
%76 = load double, double* %9
%77 = load double, double* %14
%78 = fsub double %76, %77
store double %78, double* %12
%79 = load double, double* %7
%80 = load double, double* %12
%81 = fsub double %79, %80
store double %81, double* %8
%82 = load double, double* %8
%83 = fcmp une double %82, 0.000000e+00
br i1 %83, label %84, label %91
84: ; preds = %67
%85 = load double, double* %8
%86 = load double*, double** %6
%87 = load i32, i32* %15
%88 = add nsw i32 %87, 1
store i32 %88, i32* %15
%89 = sext i32 %87 to i64
%90 = getelementptr inbounds double, double* %86, i64 %89
store double %85, double* %90
br label %91
91: ; preds = %84, %67
%92 = load double, double* %9
store double %92, double* %7
br label %93
93: ; preds = %91
%94 = load i32, i32* %11
%95 = add nsw i32 %94, 1
store i32 %95, i32* %11
br label %63
96: ; preds = %63
%97 = load double, double* %7
%98 = load double*, double** %6
%99 = load i32, i32* %15
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %98, i64 %100
store double %97, double* %101
%102 = load i32, i32* %15
%103 = add nsw i32 %102, 1
ret i32 %103
}

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