#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */
typedef int uint8_t;
typedef int *p;

/* Variables and functions */
double hllSparseSum(uint8_t *sparse, int sparselen, double *PE, int *ezp,
                    int *invalid) {
  double E = 0;
  int ez = 0, idx = 0, runlen, regval;
  uint8_t *end = sparse + sparselen, *p = sparse;
  while (p < end) {
    if ((((*(p)) & 0xc0) == 0)) {
      runlen = (((*(p)) & 0x3f) + 1);
      idx += runlen;
      ez += runlen;
      p++;
    } else if ((((*(p)) & 0xc0) == 0x40)) {
      runlen = (((((*(p)) & 0x3f) << 8) | (*((p) + 1))) + 1);
      idx += runlen;
      ez += runlen;
      p += 2;
    } else {
      runlen = (((*(p)) & 0x3) + 1);
      regval = ((((*(p)) >> 2) & 0x1f) + 1);
      idx += runlen;
      E += PE[regval] * runlen;
      p++;
    }
  }
  if (idx != (1 << 14) && invalid)
    *invalid = 1;
  E += ez;
  *ezp = ez;
  return E;
}


x86 assembly (clang, -O0):
<hllSparseSum>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x48(%rbp),%rax
jae    1250 <hllSparseSum+0x150>
mov    -0x50(%rbp),%rax
mov    (%rax),%eax
and    $0xc0,%eax
cmp    $0x0,%eax
jne    119c <hllSparseSum+0x9c>
mov    -0x50(%rbp),%rax
mov    (%rax),%eax
and    $0x3f,%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x50(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x50(%rbp)
jmp    124b <hllSparseSum+0x14b>
mov    -0x50(%rbp),%rax
mov    (%rax),%eax
and    $0xc0,%eax
cmp    $0x40,%eax
jne    11ec <hllSparseSum+0xec>
mov    -0x50(%rbp),%rax
mov    (%rax),%eax
and    $0x3f,%eax
shl    $0x8,%eax
mov    -0x50(%rbp),%rcx
or     0x4(%rcx),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x50(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x50(%rbp)
jmp    1246 <hllSparseSum+0x146>
mov    -0x50(%rbp),%rax
mov    (%rax),%eax
and    $0x3,%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x50(%rbp),%rax
mov    (%rax),%eax
sar    $0x2,%eax
and    $0x1f,%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x50(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x50(%rbp)
jmp    124b <hllSparseSum+0x14b>
jmp    1148 <hllSparseSum+0x48>
cmpl   $0x4000,-0x38(%rbp)
je     1272 <hllSparseSum+0x172>
cmpq   $0x0,-0x28(%rbp)
je     1272 <hllSparseSum+0x172>
mov    -0x28(%rbp),%rax
movl   $0x1,(%rax)
cvtsi2sdl -0x34(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x34(%rbp),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
movsd  -0x30(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<hllSparseSum>:
push   %rbx
test   %esi,%esi
jle    119d <hllSparseSum+0x9d>
movslq %esi,%rax
lea    (%rdi,%rax,4),%r9
xor    %r10d,%r10d
xorpd  %xmm0,%xmm0
xor    %r11d,%r11d
jmp    114f <hllSparseSum+0x4f>
nopl   0x0(%rax)
mov    %ebx,%esi
and    $0x3,%esi
add    $0x1,%esi
and    $0x7c,%ebx
add    %ebx,%ebx
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
mulsd  0x8(%rbx,%rdx,1),%xmm1
addsd  %xmm1,%xmm0
mov    $0x1,%eax
add    %esi,%r11d
lea    (%rdi,%rax,4),%rdi
cmp    %r9,%rdi
jae    118c <hllSparseSum+0x8c>
mov    (%rdi),%ebx
mov    %ebx,%eax
and    $0xc0,%al
cmp    $0x40,%al
je     1170 <hllSparseSum+0x70>
test   %al,%al
jne    1120 <hllSparseSum+0x20>
and    $0x3f,%ebx
mov    %ebx,%esi
add    $0x1,%esi
add    %ebx,%r10d
add    $0x1,%r10d
jmp    113e <hllSparseSum+0x3e>
xchg   %ax,%ax
and    $0x3f,%ebx
shl    $0x8,%ebx
or     0x4(%rdi),%ebx
mov    %ebx,%esi
add    $0x1,%esi
add    %ebx,%r10d
add    $0x1,%r10d
mov    $0x2,%eax
jmp    1143 <hllSparseSum+0x43>
cmp    $0x4000,%r11d
setne  %al
test   %r8,%r8
jne    11ab <hllSparseSum+0xab>
jmp    11b6 <hllSparseSum+0xb6>
xorpd  %xmm0,%xmm0
xor    %r10d,%r10d
mov    $0x1,%al
test   %r8,%r8
je     11b6 <hllSparseSum+0xb6>
test   %al,%al
je     11b6 <hllSparseSum+0xb6>
movl   $0x1,(%r8)
xorps  %xmm1,%xmm1
cvtsi2sd %r10d,%xmm1
addsd  %xmm1,%xmm0
mov    %r10d,(%rcx)
pop    %rbx
ret

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

LLVM IR (-O2):
define double @hllSparseSum(i32* %0, i32 %1, double* nocapture %2, i32* nocapture writeonly %3, i32* writeonly %4) {
%6 = sext i32 %1 to i64
%7 = getelementptr inbounds i32, i32* %0, i64 %6
%8 = icmp sgt i32 %1, 0
br i1 %8, label %9, label %50
9: ; preds = %5, %40
%10 = phi i32* [ %45, %40 ], [ %0, %5 ]
%11 = phi double [ %44, %40 ], [ 0.000000e+00, %5 ]
%12 = phi i32 [ %43, %40 ], [ 0, %5 ]
%13 = phi i32 [ %46, %40 ], [ 0, %5 ]
%14 = load i32, i32* %10
%15 = trunc i32 %14 to i8
%16 = and i8 %15, -64
switch i8 %16, label %29 [
i8 0, label %17
i8 64, label %21
]
17: ; preds = %9
%18 = and i32 %14, 63
%19 = add nuw nsw i32 %18, 1
%20 = add nsw i32 %19, %12
br label %40
21: ; preds = %9
%22 = shl i32 %14, 8
%23 = and i32 %22, 16128
%24 = getelementptr inbounds i32, i32* %10, i64 1
%25 = load i32, i32* %24
%26 = or i32 %25, %23
%27 = add nsw i32 %26, 1
%28 = add nsw i32 %27, %12
br label %40
29: ; preds = %9
%30 = and i32 %14, 3
%31 = add nuw nsw i32 %30, 1
%32 = lshr i32 %14, 2
%33 = and i32 %32, 31
%34 = add nuw nsw i32 %33, 1
%35 = zext i32 %34 to i64
%36 = getelementptr inbounds double, double* %2, i64 %35
%37 = load double, double* %36
%38 = sitofp i32 %31 to double
%39 = tail call double @llvm.fmuladd.f64(double %37, double %38, double %11)
br label %40
40: ; preds = %21, %29, %17
%41 = phi i64 [ 2, %21 ], [ 1, %29 ], [ 1, %17 ]
%42 = phi i32 [ %27, %21 ], [ %31, %29 ], [ %19, %17 ]
%43 = phi i32 [ %28, %21 ], [ %12, %29 ], [ %20, %17 ]
%44 = phi double [ %11, %21 ], [ %39, %29 ], [ %11, %17 ]
%45 = getelementptr inbounds i32, i32* %10, i64 %41
%46 = add nsw i32 %42, %13
%47 = icmp ult i32* %45, %7
br i1 %47, label %9, label %48
48: ; preds = %40
%49 = icmp ne i32 %46, 16384
br label %50
50: ; preds = %48, %5
%51 = phi i1 [ true, %5 ], [ %49, %48 ]
%52 = phi i32 [ 0, %5 ], [ %43, %48 ]
%53 = phi double [ 0.000000e+00, %5 ], [ %44, %48 ]
%54 = icmp ne i32* %4, null
%55 = and i1 %54, %51
br i1 %55, label %56, label %57
56: ; preds = %50
store i32 1, i32* %4
br label %57
57: ; preds = %56, %50
%58 = sitofp i32 %52 to double
%59 = fadd double %53, %58
store i32 %52, i32* %3
ret double %59
}