void countup_subsites(int len, int pos, long **countsp) {
  long *src = countsp[0];
  long *dst = countsp[1];
  int depth = 0;
  int max_depth = len - pos - 1;
  unsigned long site;
  unsigned long max = 1ul << len * 2;
  while (depth < max_depth) {
    for (site = 0ul; site < max; site++)
      dst[site >> 2] += src[site];
    max >>= 2;
    depth++;
    src = countsp[depth];
    dst = countsp[depth + 1];
  }
}


x86 assembly (clang, -O0):
<countup_subsites>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x4(%rbp),%eax
sub    -0x8(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
mov    %eax,%eax
mov    %eax,%ecx
mov    $0x1,%eax
shl    %cl,%rax
mov    %rax,-0x38(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    11e0 <countup_subsites+0xe0>
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x38(%rbp),%rax
jae    11a1 <countup_subsites+0xa1>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%rcx
shr    $0x2,%rcx
add    (%rax,%rcx,8),%rdx
mov    %rdx,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1162 <countup_subsites+0x62>
mov    -0x38(%rbp),%rax
shr    $0x2,%rax
mov    %rax,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x20(%rbp)
jmp    114e <countup_subsites+0x4e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<countup_subsites>:
push   %rbx
not    %esi
add    %edi,%esi
test   %esi,%esi
jle    11a9 <countup_subsites+0xa9>
add    %edi,%edi
mov    $0x1,%r9d
mov    %edi,%ecx
shl    %cl,%r9
lea    0x8(%rdx),%rax
mov    %esi,%r8d
xor    %r10d,%r10d
jmp    1145 <countup_subsites+0x45>
cs nopw 0x0(%rax,%rax,1)
shr    $0x2,%r9
lea    (%rdx,%r10,8),%rax
add    $0x10,%rax
add    $0x1,%r10
cmp    %r8,%r10
je     11a9 <countup_subsites+0xa9>
test   %r9,%r9
je     1130 <countup_subsites+0x30>
mov    (%rdx,%r10,8),%rcx
mov    (%rax),%rdi
cmp    $0x1,%r9
jne    1160 <countup_subsites+0x60>
xor    %eax,%eax
jmp    1195 <countup_subsites+0x95>
nopl   0x0(%rax,%rax,1)
mov    %r9,%r11
and    $0xfffffffffffffffe,%r11
xor    %eax,%eax
nopl   0x0(%rax)
mov    %rax,%rsi
and    $0xfffffffffffffffc,%rsi
mov    (%rdi,%rsi,2),%rbx
add    (%rcx,%rax,8),%rbx
mov    %rbx,(%rdi,%rsi,2)
add    0x8(%rcx,%rax,8),%rbx
mov    %rbx,(%rdi,%rsi,2)
add    $0x2,%rax
cmp    %rax,%r11
jne    1170 <countup_subsites+0x70>
test   $0x1,%r9b
je     1130 <countup_subsites+0x30>
mov    (%rcx,%rax,8),%rcx
and    $0xfffffffffffffffc,%rax
add    %rcx,(%rdi,%rax,2)
jmp    1130 <countup_subsites+0x30>
pop    %rbx
ret

LLVM IR (-O0):
define void @countup_subsites(i32 %0, i32 %1, i64** %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca i64**
%7 = alloca i64*
%8 = alloca i64*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i64
%12 = alloca i64
store i32 %0, i32* %4
store i32 %1, i32* %5
store i64** %2, i64*** %6
%13 = load i64**, i64*** %6
%14 = getelementptr inbounds i64*, i64** %13, i64 0
%15 = load i64*, i64** %14
store i64* %15, i64** %7
%16 = load i64**, i64*** %6
%17 = getelementptr inbounds i64*, i64** %16, i64 1
%18 = load i64*, i64** %17
store i64* %18, i64** %8
store i32 0, i32* %9
%19 = load i32, i32* %4
%20 = load i32, i32* %5
%21 = sub nsw i32 %19, %20
%22 = sub nsw i32 %21, 1
store i32 %22, i32* %10
%23 = load i32, i32* %4
%24 = mul nsw i32 %23, 2
%25 = zext i32 %24 to i64
%26 = shl i64 1, %25
store i64 %26, i64* %12
br label %27
27: ; preds = %50, %3
%28 = load i32, i32* %9
%29 = load i32, i32* %10
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %66
31: ; preds = %27
store i64 0, i64* %11
br label %32
32: ; preds = %47, %31
%33 = load i64, i64* %11
%34 = load i64, i64* %12
%35 = icmp ult i64 %33, %34
br i1 %35, label %36, label %50
36: ; preds = %32
%37 = load i64*, i64** %7
%38 = load i64, i64* %11
%39 = getelementptr inbounds i64, i64* %37, i64 %38
%40 = load i64, i64* %39
%41 = load i64*, i64** %8
%42 = load i64, i64* %11
%43 = lshr i64 %42, 2
%44 = getelementptr inbounds i64, i64* %41, i64 %43
%45 = load i64, i64* %44
%46 = add nsw i64 %45, %40
store i64 %46, i64* %44
br label %47
47: ; preds = %36
%48 = load i64, i64* %11
%49 = add i64 %48, 1
store i64 %49, i64* %11
br label %32
50: ; preds = %32
%51 = load i64, i64* %12
%52 = lshr i64 %51, 2
store i64 %52, i64* %12
%53 = load i32, i32* %9
%54 = add nsw i32 %53, 1
store i32 %54, i32* %9
%55 = load i64**, i64*** %6
%56 = load i32, i32* %9
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds i64*, i64** %55, i64 %57
%59 = load i64*, i64** %58
store i64* %59, i64** %7
%60 = load i64**, i64*** %6
%61 = load i32, i32* %9
%62 = add nsw i32 %61, 1
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds i64*, i64** %60, i64 %63
%65 = load i64*, i64** %64
store i64* %65, i64** %8
br label %27
66: ; preds = %27
ret void
}

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