void fftBRInit(long M, short *BRLow) {
  long Mroot_1 = M / 2 - 1;
  long Nroot_1 = ((unsigned long)1 << (Mroot_1));
  long i1;
  long bitsum;
  long bitmask;
  long bit;
  for (i1 = 0; i1 < Nroot_1; i1++) {
    bitsum = 0;
    bitmask = 1;
    for (bit = 1; bit <= Mroot_1; bitmask <<= 1, bit++)
      if (i1 & bitmask)
        bitsum = bitsum + (Nroot_1 >> bit);
    BRLow[i1] = bitsum;
  };
}


x86 assembly (clang, -O0):
<fftBRInit>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x8(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
sub    $0x1,%rax
mov    %rax,-0x18(%rbp)
mov    -0x18(%rbp),%rcx
mov    $0x1,%eax
shl    %cl,%rax
mov    %rax,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
jge    11df <fftBRInit+0xdf>
movq   $0x0,-0x30(%rbp)
movq   $0x1,-0x38(%rbp)
movq   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     11bb <fftBRInit+0xbb>
mov    -0x28(%rbp),%rax
and    -0x38(%rbp),%rax
cmp    $0x0,%rax
je     1199 <fftBRInit+0x99>
mov    -0x30(%rbp),%rax
mov    -0x20(%rbp),%rdx
mov    -0x40(%rbp),%rcx
sar    %cl,%rdx
mov    %rdx,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
jmp    119e <fftBRInit+0x9e>
mov    -0x38(%rbp),%rax
shl    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1160 <fftBRInit+0x60>
mov    -0x30(%rbp),%rax
mov    %ax,%dx
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    113a <fftBRInit+0x3a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<fftBRInit>:
push   %r15
push   %r14
push   %r12
push   %rbx
mov    %rdi,%rax
shr    $0x3f,%rax
add    %rdi,%rax
sar    %rax
lea    -0x1(%rax),%rcx
mov    $0x1,%r8d
mov    $0x1,%r11d
shl    %cl,%r11
cmp    $0x3f,%rcx
jne    1135 <fftBRInit+0x35>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
cmp    $0x3,%rax
mov    $0x2,%r9d
cmovge %rax,%r9
cmp    $0x2,%r11
cmovge %r11,%r8
add    $0xffffffffffffffff,%r9
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %r14d,%r14d
jmp    1170 <fftBRInit+0x70>
nopl   0x0(%rax,%rax,1)
xor    %edx,%edx
mov    %dx,(%rsi,%r14,2)
add    $0x1,%r14
cmp    %r8,%r14
je     112d <fftBRInit+0x2d>
cmp    $0x4,%rdi
jl     1160 <fftBRInit+0x60>
mov    $0x1,%r15d
mov    $0x1,%r12d
xor    %edx,%edx
cmp    $0x6,%rdi
jl     11d1 <fftBRInit+0xd1>
nopw   0x0(%rax,%rax,1)
mov    %r12,%rbx
and    %r14,%rbx
mov    %r11,%rax
mov    %r15d,%ecx
sar    %cl,%rax
test   %rbx,%rbx
cmove  %rbx,%rax
add    %rdx,%rax
lea    (%r12,%r12,1),%rbx
lea    0x1(%r15),%rcx
mov    %r11,%rdx
sar    %cl,%rdx
and    %r14,%rbx
test   %rbx,%rbx
cmove  %rbx,%rdx
add    %rax,%rdx
shl    $0x2,%r12
add    $0x2,%r15
cmp    %r10,%rcx
jne    1190 <fftBRInit+0x90>
test   $0x1,%r9b
je     1162 <fftBRInit+0x62>
and    %r14,%r12
mov    %r11,%rax
mov    %r15d,%ecx
sar    %cl,%rax
test   %r12,%r12
cmove  %r12,%rax
add    %rax,%rdx
jmp    1162 <fftBRInit+0x62>

LLVM IR (-O0):
define void @fftBRInit(i64 %0, i16* %1) {
%3 = alloca i64
%4 = alloca i16*
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
store i64 %0, i64* %3
store i16* %1, i16** %4
%11 = load i64, i64* %3
%12 = sdiv i64 %11, 2
%13 = sub nsw i64 %12, 1
store i64 %13, i64* %5
%14 = load i64, i64* %5
%15 = shl i64 1, %14
store i64 %15, i64* %6
store i64 0, i64* %7
br label %16
16: ; preds = %48, %2
%17 = load i64, i64* %7
%18 = load i64, i64* %6
%19 = icmp slt i64 %17, %18
br i1 %19, label %20, label %51
20: ; preds = %16
store i64 0, i64* %8
store i64 1, i64* %9
store i64 1, i64* %10
br label %21
21: ; preds = %37, %20
%22 = load i64, i64* %10
%23 = load i64, i64* %5
%24 = icmp sle i64 %22, %23
br i1 %24, label %25, label %42
25: ; preds = %21
%26 = load i64, i64* %7
%27 = load i64, i64* %9
%28 = and i64 %26, %27
%29 = icmp ne i64 %28, 0
br i1 %29, label %30, label %36
30: ; preds = %25
%31 = load i64, i64* %8
%32 = load i64, i64* %6
%33 = load i64, i64* %10
%34 = ashr i64 %32, %33
%35 = add nsw i64 %31, %34
store i64 %35, i64* %8
br label %36
36: ; preds = %30, %25
br label %37
37: ; preds = %36
%38 = load i64, i64* %9
%39 = shl i64 %38, 1
store i64 %39, i64* %9
%40 = load i64, i64* %10
%41 = add nsw i64 %40, 1
store i64 %41, i64* %10
br label %21
42: ; preds = %21
%43 = load i64, i64* %8
%44 = trunc i64 %43 to i16
%45 = load i16*, i16** %4
%46 = load i64, i64* %7
%47 = getelementptr inbounds i16, i16* %45, i64 %46
store i16 %44, i16* %47
br label %48
48: ; preds = %42
%49 = load i64, i64* %7
%50 = add nsw i64 %49, 1
store i64 %50, i64* %7
br label %16
51: ; preds = %16
ret void
}

LLVM IR (-O2):
define void @fftBRInit(i64 %0, i16* nocapture writeonly %1) {
%3 = sdiv i64 %0, 2
%4 = add nsw i64 %3, -1
%5 = shl nuw i64 1, %4
%6 = icmp eq i64 %4, 63
br i1 %6, label %57, label %7
7: ; preds = %2
%8 = icmp sgt i64 %0, 3
%9 = call i64 @llvm.smax.i64(i64 %3, i64 2)
%10 = call i64 @llvm.smax.i64(i64 %5, i64 1)
%11 = add nsw i64 %9, -1
%12 = and i64 %11, 1
%13 = icmp slt i64 %0, 6
%14 = and i64 %11, -2
%15 = icmp eq i64 %12, 0
br label %16
16: ; preds = %7, %51
%17 = phi i64 [ 0, %7 ], [ %55, %51 ]
br i1 %8, label %18, label %51
18: ; preds = %16
br i1 %13, label %40, label %19
19: ; preds = %18, %19
%20 = phi i64 [ %37, %19 ], [ 1, %18 ]
%21 = phi i64 [ %36, %19 ], [ 1, %18 ]
%22 = phi i64 [ %35, %19 ], [ 0, %18 ]
%23 = phi i64 [ %38, %19 ], [ 0, %18 ]
%24 = and i64 %21, %17
%25 = icmp eq i64 %24, 0
%26 = ashr i64 %5, %20
%27 = select i1 %25, i64 0, i64 %26
%28 = add nsw i64 %27, %22
%29 = shl i64 %21, 1
%30 = add nuw nsw i64 %20, 1
%31 = and i64 %29, %17
%32 = icmp eq i64 %31, 0
%33 = ashr i64 %5, %30
%34 = select i1 %32, i64 0, i64 %33
%35 = add nsw i64 %34, %28
%36 = shl i64 %21, 2
%37 = add nuw nsw i64 %20, 2
%38 = add i64 %23, 2
%39 = icmp eq i64 %38, %14
br i1 %39, label %40, label %19
40: ; preds = %19, %18
%41 = phi i64 [ undef, %18 ], [ %35, %19 ]
%42 = phi i64 [ 1, %18 ], [ %37, %19 ]
%43 = phi i64 [ 1, %18 ], [ %36, %19 ]
%44 = phi i64 [ 0, %18 ], [ %35, %19 ]
br i1 %15, label %51, label %45
45: ; preds = %40
%46 = and i64 %43, %17
%47 = icmp eq i64 %46, 0
%48 = ashr i64 %5, %42
%49 = select i1 %47, i64 0, i64 %48
%50 = add nsw i64 %49, %44
br label %51
51: ; preds = %45, %40, %16
%52 = phi i64 [ 0, %16 ], [ %41, %40 ], [ %50, %45 ]
%53 = trunc i64 %52 to i16
%54 = getelementptr inbounds i16, i16* %1, i64 %17
store i16 %53, i16* %54
%55 = add nuw nsw i64 %17, 1
%56 = icmp eq i64 %55, %10
br i1 %56, label %57, label %16
57: ; preds = %51, %2
ret void
}