void quickfwt(double p[], int n, int k, double value[]) {
  int m, l, it, i, j, is;
  double t;
  m = 1;
  l = n;
  it = 2;
  value[0] = 1;
  value[n >> 1] = 2;
  for (i = 1; i < k; i++) {
    m <<= 1;
    l >> 1;
    it <<= 1;
    for (j = 0; j < m; j++)
      value[j * l + l / 2] = it + 1 - value[j * l];
  }
  for (i = 0; i < n; i++)
    value[i] = p[(int)(value)[i] - 1];
  for (i = l = 1; i <= k; i++, l <<= 1) {
    m = n / (2 * l) - 1;
    for (j = 0; j <= m; j++) {
      it = 2 * l * j;
      for (is = 0; is < l; is++) {
        t = value[it + is] - value[it + is + l];
        value[it + is] = value[it + is] + value[it + is + l];
        value[it + is + l] = t;
      }
    }
  }
}


x86 assembly (clang, -O0):
<quickfwt>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x20(%rbp)
movl   $0x2,-0x24(%rbp)
mov    -0x18(%rbp),%rax
movsd  0xed6(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
sar    $0x1,%ecx
movslq %ecx,%rcx
movsd  0xeb5(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11f0 <quickfwt+0xf0>
mov    -0x1c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x24(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11dd <quickfwt+0xdd>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x2c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
mov    -0x20(%rbp),%eax
mov    $0x2,%esi
cltd
idiv   %esi
mov    %eax,%edx
mov    -0x40(%rbp),%rax
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    117c <quickfwt+0x7c>
jmp    11e2 <quickfwt+0xe2>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1157 <quickfwt+0x57>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    123a <quickfwt+0x13a>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x28(%rbp),%rdx
cvttsd2si (%rcx,%rdx,8),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11f7 <quickfwt+0xf7>
movl   $0x1,-0x20(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1355 <quickfwt+0x255>
mov    -0xc(%rbp),%eax
mov    -0x20(%rbp),%ecx
shl    $0x1,%ecx
cltd
idiv   %ecx
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     1339 <quickfwt+0x239>
mov    -0x20(%rbp),%eax
shl    $0x1,%eax
imul   -0x2c(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1326 <quickfwt+0x226>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    128d <quickfwt+0x18d>
jmp    132b <quickfwt+0x22b>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    126d <quickfwt+0x16d>
jmp    133e <quickfwt+0x23e>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x20(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1248 <quickfwt+0x148>
pop    %rbp
ret

x86 assembly (clang, -O2):
<quickfwt>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,%r15
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rcx)
mov    %esi,-0x54(%rsp)
mov    %esi,%eax
sar    %eax
cltq
movabs $0x4000000000000000,%rcx
mov    %rcx,(%r15,%rax,8)
mov    %edx,-0x50(%rsp)
cmp    $0x2,%edx
jl     11ca <quickfwt+0xca>
mov    -0x54(%rsp),%eax
mov    %eax,%ecx
shr    $0x1f,%ecx
add    %eax,%ecx
sar    %ecx
cltq
movslq %ecx,%rcx
lea    (%rcx,%rax,1),%rdx
mov    %rax,%rsi
shl    $0x4,%rsi
mov    $0x2,%r8d
mov    $0x1,%r9d
mov    $0x1,%r10d
jmp    117b <quickfwt+0x7b>
add    $0x1,%r10d
cmp    -0x50(%rsp),%r10d
je     11ca <quickfwt+0xca>
add    %r9d,%r9d
add    %r8d,%r8d
test   %r9d,%r9d
jle    1170 <quickfwt+0x70>
mov    %r8d,%ebp
or     $0x1,%ebp
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
mov    %r9d,%ebx
mov    %r15,%rbp
nopl   0x0(%rax)
movapd %xmm0,%xmm1
subsd  0x0(%rbp),%xmm1
movsd  %xmm1,0x0(%rbp,%rcx,8)
movapd %xmm0,%xmm1
subsd  0x0(%rbp,%rax,8),%xmm1
movsd  %xmm1,0x0(%rbp,%rdx,8)
add    %rsi,%rbp
add    $0xfffffffffffffffe,%rbx
jne    11a0 <quickfwt+0xa0>
jmp    1170 <quickfwt+0x70>
cmpl   $0x0,-0x54(%rsp)
jle    123d <quickfwt+0x13d>
mov    -0x54(%rsp),%ecx
mov    %ecx,%eax
cmp    $0x1,%ecx
jne    11e0 <quickfwt+0xe0>
xor    %ecx,%ecx
jmp    1225 <quickfwt+0x125>
mov    %eax,%edx
and    $0xfffffffe,%edx
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
cvttsd2si (%r15,%rcx,8),%esi
movslq %esi,%rsi
movsd  -0x8(%rdi,%rsi,8),%xmm0
movsd  %xmm0,(%r15,%rcx,8)
cvttsd2si 0x8(%r15,%rcx,8),%esi
movslq %esi,%rsi
movsd  -0x8(%rdi,%rsi,8),%xmm0
movsd  %xmm0,0x8(%r15,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rdx
jne    11f0 <quickfwt+0xf0>
test   $0x1,%al
je     123d <quickfwt+0x13d>
cvttsd2si (%r15,%rcx,8),%eax
cltq
movsd  -0x8(%rdi,%rax,8),%xmm0
movsd  %xmm0,(%r15,%rcx,8)
cmpl   $0x0,-0x50(%rsp)
jle    14e9 <quickfwt+0x3e9>
lea    0x10(%r15),%rax
mov    %rax,-0x40(%rsp)
lea    0x8(%r15),%rax
mov    %rax,-0x48(%rsp)
mov    $0x1,%eax
mov    $0x1,%edi
jmp    1288 <quickfwt+0x188>
cs nopw 0x0(%rax,%rax,1)
mov    -0x38(%rsp),%rcx
lea    0x1(%rcx),%eax
mov    -0x4c(%rsp),%edx
mov    %edx,%edi
cmp    -0x50(%rsp),%ecx
je     14e9 <quickfwt+0x3e9>
mov    %rax,-0x38(%rsp)
lea    (%rdi,%rdi,1),%ecx
mov    -0x54(%rsp),%eax
cltd
mov    %ecx,-0x4c(%rsp)
idiv   %ecx
test   %eax,%eax
jle    1270 <quickfwt+0x170>
movslq %edi,%rdx
movslq -0x4c(%rsp),%rcx
mov    %eax,%eax
mov    %edi,%ebp
lea    (%rdx,%rbp,1),%rsi
mov    %rsi,-0x10(%rsp)
mov    %ebp,%esi
and    $0xfffffffe,%esi
mov    %rsi,-0x18(%rsp)
add    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x20(%rsp)
mov    %rsi,%rbx
shr    %rbx
add    $0x1,%rbx
mov    %rbx,%rsi
mov    %rbx,-0x28(%rsp)
and    $0xfffffffffffffffe,%rbx
mov    %rbx,-0x30(%rsp)
mov    %rbp,%rsi
neg    %rsi
mov    %rsi,-0x8(%rsp)
mov    -0x40(%rsp),%r9
lea    (%r9,%rdx,8),%r8
mov    %rcx,%r11
lea    0x0(,%rcx,8),%r14
mov    -0x48(%rsp),%rcx
lea    (%rcx,%rdx,8),%r10
mov    %r15,%rcx
xor    %r12d,%r12d
jmp    1339 <quickfwt+0x239>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r12
add    %r14,%r8
add    %r14,%r9
add    %r14,%r10
add    %r14,%rcx
cmp    %rax,%r12
je     1270 <quickfwt+0x170>
test   %edi,%edi
jle    1320 <quickfwt+0x220>
mov    %r12,%r13
imul   %r11,%r13
cmp    $0x1,%edi
je     137d <quickfwt+0x27d>
mov    -0x10(%rsp),%rsi
add    %r13,%rsi
lea    (%r15,%r13,8),%rbx
lea    (%r15,%rsi,8),%rsi
cmp    %rsi,%rbx
jae    1426 <quickfwt+0x326>
mov    %r13,%rsi
add    %rbp,%rsi
lea    (%rdx,%r13,1),%rbx
lea    (%r15,%rsi,8),%rsi
lea    (%r15,%rbx,8),%rbx
cmp    %rsi,%rbx
jae    1426 <quickfwt+0x326>
xor    %ebx,%ebx
mov    %rbx,%rsi
test   $0x1,%bpl
je     13ba <quickfwt+0x2ba>
add    %rbx,%r13
movsd  (%r15,%r13,8),%xmm0
lea    (%rdx,%r13,1),%rsi
movsd  (%r15,%rsi,8),%xmm1
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r15,%r13,8)
movsd  %xmm2,(%r15,%rsi,8)
mov    %rbx,%rsi
or     $0x1,%rsi
not    %rbx
cmp    -0x8(%rsp),%rbx
je     1320 <quickfwt+0x220>
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rsi,8),%xmm0
movsd  -0x8(%r10,%rsi,8),%xmm1
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rcx,%rsi,8)
movsd  %xmm2,-0x8(%r10,%rsi,8)
movsd  0x8(%rcx,%rsi,8),%xmm0
movsd  (%r10,%rsi,8),%xmm1
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
addsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rcx,%rsi,8)
movsd  %xmm2,(%r10,%rsi,8)
add    $0x2,%rsi
cmp    %rsi,%rbp
jne    13d0 <quickfwt+0x2d0>
jmp    1320 <quickfwt+0x220>
cmpq   $0x0,-0x20(%rsp)
je     14de <quickfwt+0x3de>
mov    -0x30(%rsp),%rbx
xor    %esi,%esi
nopl   0x0(%rax)
movupd -0x10(%r9,%rsi,8),%xmm0
movupd -0x10(%r8,%rsi,8),%xmm1
movapd %xmm0,%xmm2
subpd  %xmm1,%xmm2
addpd  %xmm1,%xmm0
movupd %xmm0,-0x10(%r9,%rsi,8)
movupd %xmm2,-0x10(%r8,%rsi,8)
movupd (%r9,%rsi,8),%xmm0
movupd (%r8,%rsi,8),%xmm1
movapd %xmm0,%xmm2
subpd  %xmm1,%xmm2
addpd  %xmm1,%xmm0
movupd %xmm0,(%r9,%rsi,8)
movupd %xmm2,(%r8,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rbx
jne    1440 <quickfwt+0x340>
testb  $0x1,-0x28(%rsp)
je     14c8 <quickfwt+0x3c8>
add    %r13,%rsi
movupd (%r15,%rsi,8),%xmm0
lea    (%rsi,%rdx,1),%rbx
movupd (%r15,%rbx,8),%xmm1
movapd %xmm0,%xmm2
subpd  %xmm1,%xmm2
addpd  %xmm1,%xmm0
movupd %xmm0,(%r15,%rsi,8)
movupd %xmm2,(%r15,%rbx,8)
mov    -0x18(%rsp),%rsi
mov    %rsi,%rbx
cmp    %rbp,%rsi
je     1320 <quickfwt+0x220>
jmp    137f <quickfwt+0x27f>
xor    %esi,%esi
testb  $0x1,-0x28(%rsp)
jne    149d <quickfwt+0x39d>
jmp    14c8 <quickfwt+0x3c8>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @quickfwt(double* %0, i32 %1, i32 %2, double* %3) {
%5 = alloca double*
%6 = alloca i32
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca double
store double* %0, double** %5
store i32 %1, i32* %6
store i32 %2, i32* %7
store double* %3, double** %8
store i32 1, i32* %9
%16 = load i32, i32* %6
store i32 %16, i32* %10
store i32 2, i32* %11
%17 = load double*, double** %8
%18 = getelementptr inbounds double, double* %17, i64 0
store double 1.000000e+00, double* %18
%19 = load double*, double** %8
%20 = load i32, i32* %6
%21 = ashr i32 %20, 1
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %19, i64 %22
store double 2.000000e+00, double* %23
store i32 1, i32* %12
br label %24
24: ; preds = %64, %4
%25 = load i32, i32* %12
%26 = load i32, i32* %7
%27 = icmp slt i32 %25, %26
br i1 %27, label %28, label %67
28: ; preds = %24
%29 = load i32, i32* %9
%30 = shl i32 %29, 1
store i32 %30, i32* %9
%31 = load i32, i32* %10
%32 = ashr i32 %31, 1
%33 = load i32, i32* %11
%34 = shl i32 %33, 1
store i32 %34, i32* %11
store i32 0, i32* %13
br label %35
35: ; preds = %60, %28
%36 = load i32, i32* %13
%37 = load i32, i32* %9
%38 = icmp slt i32 %36, %37
br i1 %38, label %39, label %63
39: ; preds = %35
%40 = load i32, i32* %11
%41 = add nsw i32 %40, 1
%42 = sitofp i32 %41 to double
%43 = load double*, double** %8
%44 = load i32, i32* %13
%45 = load i32, i32* %10
%46 = mul nsw i32 %44, %45
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %43, i64 %47
%49 = load double, double* %48
%50 = fsub double %42, %49
%51 = load double*, double** %8
%52 = load i32, i32* %13
%53 = load i32, i32* %10
%54 = mul nsw i32 %52, %53
%55 = load i32, i32* %10
%56 = sdiv i32 %55, 2
%57 = add nsw i32 %54, %56
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %51, i64 %58
store double %50, double* %59
br label %60
60: ; preds = %39
%61 = load i32, i32* %13
%62 = add nsw i32 %61, 1
store i32 %62, i32* %13
br label %35
63: ; preds = %35
br label %64
64: ; preds = %63
%65 = load i32, i32* %12
%66 = add nsw i32 %65, 1
store i32 %66, i32* %12
br label %24
67: ; preds = %24
store i32 0, i32* %12
br label %68
68: ; preds = %88, %67
%69 = load i32, i32* %12
%70 = load i32, i32* %6
%71 = icmp slt i32 %69, %70
br i1 %71, label %72, label %91
72: ; preds = %68
%73 = load double*, double** %5
%74 = load double*, double** %8
%75 = load i32, i32* %12
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %74, i64 %76
%78 = load double, double* %77
%79 = fptosi double %78 to i32
%80 = sub nsw i32 %79, 1
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %73, i64 %81
%83 = load double, double* %82
%84 = load double*, double** %8
%85 = load i32, i32* %12
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
store double %83, double* %87
br label %88
88: ; preds = %72
%89 = load i32, i32* %12
%90 = add nsw i32 %89, 1
store i32 %90, i32* %12
br label %68
91: ; preds = %68
store i32 1, i32* %10
store i32 1, i32* %12
br label %92
92: ; preds = %173, %91
%93 = load i32, i32* %12
%94 = load i32, i32* %7
%95 = icmp sle i32 %93, %94
br i1 %95, label %96, label %178
96: ; preds = %92
%97 = load i32, i32* %6
%98 = load i32, i32* %10
%99 = mul nsw i32 2, %98
%100 = sdiv i32 %97, %99
%101 = sub nsw i32 %100, 1
store i32 %101, i32* %9
store i32 0, i32* %13
br label %102
102: ; preds = %169, %96
%103 = load i32, i32* %13
%104 = load i32, i32* %9
%105 = icmp sle i32 %103, %104
br i1 %105, label %106, label %172
106: ; preds = %102
%107 = load i32, i32* %10
%108 = mul nsw i32 2, %107
%109 = load i32, i32* %13
%110 = mul nsw i32 %108, %109
store i32 %110, i32* %11
store i32 0, i32* %14
br label %111
111: ; preds = %165, %106
%112 = load i32, i32* %14
%113 = load i32, i32* %10
%114 = icmp slt i32 %112, %113
br i1 %114, label %115, label %168
115: ; preds = %111
%116 = load double*, double** %8
%117 = load i32, i32* %11
%118 = load i32, i32* %14
%119 = add nsw i32 %117, %118
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds double, double* %116, i64 %120
%122 = load double, double* %121
%123 = load double*, double** %8
%124 = load i32, i32* %11
%125 = load i32, i32* %14
%126 = add nsw i32 %124, %125
%127 = load i32, i32* %10
%128 = add nsw i32 %126, %127
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %123, i64 %129
%131 = load double, double* %130
%132 = fsub double %122, %131
store double %132, double* %15
%133 = load double*, double** %8
%134 = load i32, i32* %11
%135 = load i32, i32* %14
%136 = add nsw i32 %134, %135
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %133, i64 %137
%139 = load double, double* %138
%140 = load double*, double** %8
%141 = load i32, i32* %11
%142 = load i32, i32* %14
%143 = add nsw i32 %141, %142
%144 = load i32, i32* %10
%145 = add nsw i32 %143, %144
%146 = sext i32 %145 to i64
%147 = getelementptr inbounds double, double* %140, i64 %146
%148 = load double, double* %147
%149 = fadd double %139, %148
%150 = load double*, double** %8
%151 = load i32, i32* %11
%152 = load i32, i32* %14
%153 = add nsw i32 %151, %152
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds double, double* %150, i64 %154
store double %149, double* %155
%156 = load double, double* %15
%157 = load double*, double** %8
%158 = load i32, i32* %11
%159 = load i32, i32* %14
%160 = add nsw i32 %158, %159
%161 = load i32, i32* %10
%162 = add nsw i32 %160, %161
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %157, i64 %163
store double %156, double* %164
br label %165
165: ; preds = %115
%166 = load i32, i32* %14
%167 = add nsw i32 %166, 1
store i32 %167, i32* %14
br label %111
168: ; preds = %111
br label %169
169: ; preds = %168
%170 = load i32, i32* %13
%171 = add nsw i32 %170, 1
store i32 %171, i32* %13
br label %102
172: ; preds = %102
br label %173
173: ; preds = %172
%174 = load i32, i32* %12
%175 = add nsw i32 %174, 1
store i32 %175, i32* %12
%176 = load i32, i32* %10
%177 = shl i32 %176, 1
store i32 %177, i32* %10
br label %92
178: ; preds = %92
ret void
}

LLVM IR (-O2):
define void @quickfwt(double* nocapture %0, i32 %1, i32 %2, double* nocapture %3) {
store double 1.000000e+00, double* %3
%5 = ashr i32 %1, 1
%6 = sext i32 %5 to i64
%7 = getelementptr inbounds double, double* %3, i64 %6
store double 2.000000e+00, double* %7
%8 = icmp sgt i32 %2, 1
br i1 %8, label %9, label %13
9: ; preds = %4
%10 = sdiv i32 %1, 2
%11 = sext i32 %1 to i64
%12 = sext i32 %10 to i64
br label %21
13: ; preds = %49, %4
%14 = icmp sgt i32 %1, 0
br i1 %14, label %15, label %63
15: ; preds = %13
%16 = zext i32 %1 to i64
%17 = and i64 %16, 1
%18 = icmp eq i32 %1, 1
br i1 %18, label %52, label %19
19: ; preds = %15
%20 = and i64 %16, 4294967294
br label %65
21: ; preds = %9, %49
%22 = phi i32 [ 1, %9 ], [ %25, %49 ]
%23 = phi i32 [ 1, %9 ], [ %50, %49 ]
%24 = phi i32 [ 2, %9 ], [ %26, %49 ]
%25 = shl i32 %22, 1
%26 = shl i32 %24, 1
%27 = icmp sgt i32 %25, 0
br i1 %27, label %28, label %49
28: ; preds = %21
%29 = or i32 %26, 1
%30 = sitofp i32 %29 to double
%31 = zext i32 %25 to i64
br label %32
32: ; preds = %32, %28
%33 = phi i64 [ 0, %28 ], [ %47, %32 ]
%34 = mul nsw i64 %33, %11
%35 = getelementptr inbounds double, double* %3, i64 %34
%36 = load double, double* %35
%37 = fsub double %30, %36
%38 = add nsw i64 %34, %12
%39 = getelementptr inbounds double, double* %3, i64 %38
store double %37, double* %39
%40 = or i64 %33, 1
%41 = mul nsw i64 %40, %11
%42 = getelementptr inbounds double, double* %3, i64 %41
%43 = load double, double* %42
%44 = fsub double %30, %43
%45 = add nsw i64 %41, %12
%46 = getelementptr inbounds double, double* %3, i64 %45
store double %44, double* %46
%47 = add nuw nsw i64 %33, 2
%48 = icmp eq i64 %47, %31
br i1 %48, label %49, label %32
49: ; preds = %32, %21
%50 = add nuw nsw i32 %23, 1
%51 = icmp eq i32 %50, %2
br i1 %51, label %13, label %21
52: ; preds = %65, %15
%53 = phi i64 [ 0, %15 ], [ %83, %65 ]
%54 = icmp eq i64 %17, 0
br i1 %54, label %63, label %55
55: ; preds = %52
%56 = getelementptr inbounds double, double* %3, i64 %53
%57 = load double, double* %56
%58 = fptosi double %57 to i32
%59 = add nsw i32 %58, -1
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %0, i64 %60
%62 = load double, double* %61
store double %62, double* %56
br label %63
63: ; preds = %55, %52, %13
%64 = icmp slt i32 %2, 1
br i1 %64, label %220, label %86
65: ; preds = %65, %19
%66 = phi i64 [ 0, %19 ], [ %83, %65 ]
%67 = phi i64 [ 0, %19 ], [ %84, %65 ]
%68 = getelementptr inbounds double, double* %3, i64 %66
%69 = load double, double* %68
%70 = fptosi double %69 to i32
%71 = add nsw i32 %70, -1
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %0, i64 %72
%74 = load double, double* %73
store double %74, double* %68
%75 = or i64 %66, 1
%76 = getelementptr inbounds double, double* %3, i64 %75
%77 = load double, double* %76
%78 = fptosi double %77 to i32
%79 = add nsw i32 %78, -1
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %0, i64 %80
%82 = load double, double* %81
store double %82, double* %76
%83 = add nuw nsw i64 %66, 2
%84 = add i64 %67, 2
%85 = icmp eq i64 %84, %20
br i1 %85, label %52, label %65
86: ; preds = %63, %217
%87 = phi i32 [ %218, %217 ], [ 1, %63 ]
%88 = phi i32 [ %89, %217 ], [ 1, %63 ]
%89 = shl nsw i32 %88, 1
%90 = sdiv i32 %1, %89
%91 = icmp sgt i32 %90, 0
br i1 %91, label %92, label %217
92: ; preds = %86
%93 = icmp sgt i32 %88, 0
%94 = sext i32 %88 to i64
%95 = sext i32 %89 to i64
%96 = zext i32 %90 to i64
%97 = zext i32 %88 to i64
%98 = add nsw i64 %94, %97
%99 = and i64 %97, 4294967294
%100 = add nsw i64 %99, -2
%101 = lshr exact i64 %100, 1
%102 = add nuw i64 %101, 1
%103 = icmp eq i32 %88, 1
%104 = and i64 %97, 4294967294
%105 = and i64 %102, 1
%106 = icmp eq i64 %100, 0
%107 = and i64 %102, -2
%108 = icmp eq i64 %105, 0
%109 = icmp eq i64 %104, %97
%110 = and i64 %97, 1
%111 = icmp eq i64 %110, 0
%112 = sub nsw i64 0, %97
br label %113
113: ; preds = %92, %214
%114 = phi i64 [ 0, %92 ], [ %215, %214 ]
%115 = mul i64 %114, %95
%116 = getelementptr double, double* %3, i64 %115
%117 = add i64 %115, %97
%118 = getelementptr double, double* %3, i64 %117
%119 = add i64 %115, %94
%120 = getelementptr double, double* %3, i64 %119
%121 = add i64 %98, %115
%122 = getelementptr double, double* %3, i64 %121
%123 = mul nsw i64 %114, %95
br i1 %93, label %124, label %214
124: ; preds = %113
br i1 %103, label %177, label %125
125: ; preds = %124
%126 = icmp ult double* %116, %122
%127 = icmp ult double* %120, %118
%128 = and i1 %126, %127
br i1 %128, label %177, label %129
129: ; preds = %125
br i1 %106, label %161, label %130
130: ; preds = %129, %130
%131 = phi i64 [ %158, %130 ], [ 0, %129 ]
%132 = phi i64 [ %159, %130 ], [ 0, %129 ]
%133 = add nsw i64 %131, %123
%134 = getelementptr inbounds double, double* %3, i64 %133
%135 = bitcast double* %134 to <2 x double>*
%136 = load <2 x double>, <2 x double>* %135
%137 = add nsw i64 %133, %94
%138 = getelementptr inbounds double, double* %3, i64 %137
%139 = bitcast double* %138 to <2 x double>*
%140 = load <2 x double>, <2 x double>* %139
%141 = fsub <2 x double> %136, %140
%142 = fadd <2 x double> %136, %140
%143 = bitcast double* %134 to <2 x double>*
store <2 x double> %142, <2 x double>* %143
%144 = bitcast double* %138 to <2 x double>*
store <2 x double> %141, <2 x double>* %144
%145 = or i64 %131, 2
%146 = add nsw i64 %145, %123
%147 = getelementptr inbounds double, double* %3, i64 %146
%148 = bitcast double* %147 to <2 x double>*
%149 = load <2 x double>, <2 x double>* %148
%150 = add nsw i64 %146, %94
%151 = getelementptr inbounds double, double* %3, i64 %150
%152 = bitcast double* %151 to <2 x double>*
%153 = load <2 x double>, <2 x double>* %152
%154 = fsub <2 x double> %149, %153
%155 = fadd <2 x double> %149, %153
%156 = bitcast double* %147 to <2 x double>*
store <2 x double> %155, <2 x double>* %156
%157 = bitcast double* %151 to <2 x double>*
store <2 x double> %154, <2 x double>* %157
%158 = add nuw i64 %131, 4
%159 = add nuw i64 %132, 2
%160 = icmp eq i64 %159, %107
br i1 %160, label %161, label %130
161: ; preds = %130, %129
%162 = phi i64 [ 0, %129 ], [ %158, %130 ]
br i1 %108, label %176, label %163
163: ; preds = %161
%164 = add nsw i64 %162, %123
%165 = getelementptr inbounds double, double* %3, i64 %164
%166 = bitcast double* %165 to <2 x double>*
%167 = load <2 x double>, <2 x double>* %166
%168 = add nsw i64 %164, %94
%169 = getelementptr inbounds double, double* %3, i64 %168
%170 = bitcast double* %169 to <2 x double>*
%171 = load <2 x double>, <2 x double>* %170
%172 = fsub <2 x double> %167, %171
%173 = fadd <2 x double> %167, %171
%174 = bitcast double* %165 to <2 x double>*
store <2 x double> %173, <2 x double>* %174
%175 = bitcast double* %169 to <2 x double>*
store <2 x double> %172, <2 x double>* %175
br label %176
176: ; preds = %161, %163
br i1 %109, label %214, label %177
177: ; preds = %125, %124, %176
%178 = phi i64 [ 0, %125 ], [ 0, %124 ], [ %104, %176 ]
%179 = xor i64 %178, -1
br i1 %111, label %190, label %180
180: ; preds = %177
%181 = add nsw i64 %178, %123
%182 = getelementptr inbounds double, double* %3, i64 %181
%183 = load double, double* %182
%184 = add nsw i64 %181, %94
%185 = getelementptr inbounds double, double* %3, i64 %184
%186 = load double, double* %185
%187 = fsub double %183, %186
%188 = fadd double %183, %186
store double %188, double* %182
store double %187, double* %185
%189 = or i64 %178, 1
br label %190
190: ; preds = %180, %177
%191 = phi i64 [ %178, %177 ], [ %189, %180 ]
%192 = icmp eq i64 %179, %112
br i1 %192, label %214, label %193
193: ; preds = %190, %193
%194 = phi i64 [ %212, %193 ], [ %191, %190 ]
%195 = add nsw i64 %194, %123
%196 = getelementptr inbounds double, double* %3, i64 %195
%197 = load double, double* %196
%198 = add nsw i64 %195, %94
%199 = getelementptr inbounds double, double* %3, i64 %198
%200 = load double, double* %199
%201 = fsub double %197, %200
%202 = fadd double %197, %200
store double %202, double* %196
store double %201, double* %199
%203 = add nuw nsw i64 %194, 1
%204 = add nsw i64 %203, %123
%205 = getelementptr inbounds double, double* %3, i64 %204
%206 = load double, double* %205
%207 = add nsw i64 %204, %94
%208 = getelementptr inbounds double, double* %3, i64 %207
%209 = load double, double* %208
%210 = fsub double %206, %209
%211 = fadd double %206, %209
store double %211, double* %205
store double %210, double* %208
%212 = add nuw nsw i64 %194, 2
%213 = icmp eq i64 %212, %97
br i1 %213, label %214, label %193
214: ; preds = %190, %193, %176, %113
%215 = add nuw nsw i64 %114, 1
%216 = icmp eq i64 %215, %96
br i1 %216, label %217, label %113
217: ; preds = %214, %86
%218 = add nuw i32 %87, 1
%219 = icmp eq i32 %87, %2
br i1 %219, label %220, label %86
220: ; preds = %217, %63
ret void
}