void FFT(double *x, double *y, int n, double w[]) {
  int m, i, j, k, nbits;
  double tempx, tempy;
  m = 0;
  i = n;
  while (i > 0) {
    i /= 2;
    m++;
  }
  m -= 1;
  nbits = n >> 1;
  j = 0;
  for (i = 0; i < n - 1; i++) {
    if (i < j) {
      tempx = x[i];
      tempy = y[i];
      x[i] = x[j];
      y[i] = y[j];
      x[j] = tempx;
      y[j] = tempy;
    }
    k = nbits;
    while (k <= j) {
      j -= k;
      k >>= 1;
    }
    j += k;
  }
  int mj, term_i, stride_i, mi, j2, count2;
  double u1, u2, t1, t2;
  mj = 1;
  stride_i = n / (4 * mj);
  for (k = 0; k < m; k++) {
    mi = 2 * mj;
    term_i = n / mi;
    for (i = 0; i < term_i; i++) {
      count2 = 0;
      for (j = i * mi; count2 < mj; j++, count2++) {
        j %= (n - 1);
        j2 = (j + mj);
        int temp = count2 * n / mi;
        u1 = w[temp * 2 + 0];
        u2 = -w[temp * 2 + 1];
        t1 = u1 * x[j2] - u2 * y[j2];
        t2 = u1 * y[j2] + u2 * x[j2];
        x[j2] = x[j] - t1;
        y[j2] = y[j] - t2;
        x[j] += t1;
        y[j] += t2;
      }
    }
    stride_i /= 2;
    mj *= 2;
  }
  return;
}


x86 assembly (clang, -O0):
<FFT>:
push   %rbp
mov    %rsp,%rbp
push   %rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x28(%rbp)
cmpl   $0x0,-0x28(%rbp)
jle    1147 <FFT+0x47>
mov    -0x28(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x28(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1121 <FFT+0x21>
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x14(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1242 <FFT+0x142>
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1200 <FFT+0x100>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     122b <FFT+0x12b>
mov    -0x30(%rbp),%ecx
mov    -0x2c(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x30(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1206 <FFT+0x106>
mov    -0x30(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1167 <FFT+0x67>
movl   $0x1,-0x4c(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x4c(%rbp),%ecx
shl    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x54(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    145c <FFT+0x35c>
mov    -0x4c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x58(%rbp)
mov    -0x14(%rbp),%eax
cltd
idivl  -0x58(%rbp)
mov    %eax,-0x50(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x50(%rbp),%eax
jge    1437 <FFT+0x337>
movl   $0x0,-0x60(%rbp)
mov    -0x28(%rbp),%eax
imul   -0x58(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x60(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jge    1424 <FFT+0x324>
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x2c(%rbp),%eax
cltd
idiv   %ecx
mov    %edx,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x60(%rbp),%eax
imul   -0x14(%rbp),%eax
cltd
idivl  -0x58(%rbp)
mov    %eax,-0x84(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x84(%rbp),%ecx
shl    $0x1,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x84(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x70(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x68(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x70(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x78(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x80(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    12a2 <FFT+0x1a2>
jmp    1429 <FFT+0x329>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1285 <FFT+0x185>
mov    -0x54(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x54(%rbp)
mov    -0x4c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    125f <FFT+0x15f>
add    $0x8,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<FFT>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    129c <FFT+0x19c>
mov    %edx,%r8d
mov    $0x1,%r10d
mov    %r8d,%eax
xchg   %ax,%ax
shr    %eax
add    $0xffffffff,%r10d
cmp    $0x1,%edx
mov    %eax,%edx
ja     1120 <FFT+0x20>
lea    -0x1(%r8),%r12d
cmp    $0x2,%r8d
jl     1195 <FFT+0x95>
mov    %r12d,%r9d
xor    %edx,%edx
xor    %ebp,%ebp
xchg   %ax,%ax
movslq %ebp,%rax
cmp    %rax,%rdx
jge    1170 <FFT+0x70>
movsd  (%rdi,%rdx,8),%xmm0
movsd  (%rsi,%rdx,8),%xmm1
movsd  (%rdi,%rax,8),%xmm2
movsd  %xmm2,(%rdi,%rdx,8)
movsd  (%rsi,%rax,8),%xmm2
movsd  %xmm2,(%rsi,%rdx,8)
movsd  %xmm0,(%rdi,%rax,8)
movsd  %xmm1,(%rsi,%rax,8)
mov    %r8d,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %ebp,%ebx
sar    %eax
sub    %eax,%ebp
jge    1180 <FFT+0x80>
add    %eax,%ebx
add    $0x1,%rdx
mov    %ebx,%ebp
cmp    %r9,%rdx
jne    1140 <FFT+0x40>
test   %r10d,%r10d
je     129c <FFT+0x19c>
neg    %r10d
mov    $0x1,%ebx
xor    %eax,%eax
mov    %r10d,-0x8(%rsp)
jmp    11c7 <FFT+0xc7>
nop
mov    -0x4(%rsp),%eax
add    $0x1,%eax
mov    %ebp,%ebx
mov    -0x8(%rsp),%r10d
cmp    %eax,%r10d
je     129c <FFT+0x19c>
mov    %eax,-0x4(%rsp)
lea    (%rbx,%rbx,1),%ebp
mov    %r8d,%eax
cltd
idiv   %ebp
test   %eax,%eax
jle    11b0 <FFT+0xb0>
mov    %eax,%r9d
xor    %r13d,%r13d
jmp    11e9 <FFT+0xe9>
add    $0x1,%r13d
cmp    %r9d,%r13d
je     11b0 <FFT+0xb0>
test   %ebx,%ebx
jle    11e0 <FFT+0xe0>
mov    %r13d,%eax
imul   %ebp,%eax
xor    %r10d,%r10d
mov    %ebx,%r14d
nopl   0x0(%rax)
cltd
idiv   %r12d
mov    %edx,%r11d
lea    (%r11,%rbx,1),%r15d
mov    %r10d,%eax
cltd
idiv   %ebp
lea    (%rax,%rax,1),%edx
movslq %edx,%rdx
movsd  (%rcx,%rdx,8),%xmm0
lea    0x1(%rax,%rax,1),%eax
cltq
movsd  (%rcx,%rax,8),%xmm1
movslq %r15d,%rdx
movsd  (%rdi,%rdx,8),%xmm2
movsd  (%rsi,%rdx,8),%xmm3
movapd %xmm1,%xmm4
mulsd  %xmm3,%xmm4
movapd %xmm0,%xmm5
mulsd  %xmm2,%xmm5
addsd  %xmm4,%xmm5
mulsd  %xmm2,%xmm1
mulsd  %xmm3,%xmm0
subsd  %xmm1,%xmm0
movslq %r11d,%rax
movsd  (%rdi,%rax,8),%xmm1
subsd  %xmm5,%xmm1
movsd  %xmm1,(%rdi,%rdx,8)
movsd  (%rsi,%rax,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rdx,8)
addsd  (%rdi,%rax,8),%xmm5
movsd  %xmm5,(%rdi,%rax,8)
addsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rax,8)
add    $0x1,%eax
add    %r8d,%r10d
add    $0xffffffff,%r14d
jne    1200 <FFT+0x100>
jmp    11e0 <FFT+0xe0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @FFT(double* nocapture %0, double* nocapture %1, i32 %2, double* nocapture %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %101
6: ; preds = %4, %6
%7 = phi i32 [ %10, %6 ], [ 0, %4 ]
%8 = phi i32 [ %9, %6 ], [ %2, %4 ]
%9 = lshr i32 %8, 1
%10 = add nuw nsw i32 %7, 1
%11 = icmp ult i32 %8, 2
br i1 %11, label %12, label %6
12: ; preds = %6
%13 = add i32 %2, -1
%14 = icmp sgt i32 %2, 1
br i1 %14, label %15, label %17
15: ; preds = %12
%16 = zext i32 %13 to i64
br label %19
17: ; preds = %40, %12
%18 = icmp eq i32 %7, 0
br i1 %18, label %101, label %44
19: ; preds = %15, %40
%20 = phi i64 [ 0, %15 ], [ %42, %40 ]
%21 = phi i32 [ 0, %15 ], [ %41, %40 ]
%22 = sext i32 %21 to i64
%23 = icmp slt i64 %20, %22
br i1 %23, label %24, label %33
24: ; preds = %19
%25 = getelementptr inbounds double, double* %0, i64 %20
%26 = load double, double* %25
%27 = getelementptr inbounds double, double* %1, i64 %20
%28 = load double, double* %27
%29 = getelementptr inbounds double, double* %0, i64 %22
%30 = load double, double* %29
store double %30, double* %25
%31 = getelementptr inbounds double, double* %1, i64 %22
%32 = load double, double* %31
store double %32, double* %27
store double %26, double* %29
store double %28, double* %31
br label %33
33: ; preds = %24, %19
br label %34
34: ; preds = %33, %34
%35 = phi i32 [ %37, %34 ], [ %2, %33 ]
%36 = phi i32 [ %39, %34 ], [ %21, %33 ]
%37 = ashr i32 %35, 1
%38 = icmp slt i32 %36, %37
%39 = sub nsw i32 %36, %37
br i1 %38, label %40, label %34
40: ; preds = %34
%41 = add nsw i32 %36, %37
%42 = add nuw nsw i64 %20, 1
%43 = icmp eq i64 %42, %16
br i1 %43, label %17, label %19
44: ; preds = %17, %98
%45 = phi i32 [ %99, %98 ], [ 0, %17 ]
%46 = phi i32 [ %47, %98 ], [ 1, %17 ]
%47 = shl nsw i32 %46, 1
%48 = sdiv i32 %2, %47
%49 = icmp sgt i32 %48, 0
br i1 %49, label %50, label %98
50: ; preds = %44
%51 = icmp sgt i32 %46, 0
br label %52
52: ; preds = %50, %95
%53 = phi i32 [ 0, %50 ], [ %96, %95 ]
br i1 %51, label %54, label %95
54: ; preds = %52
%55 = mul nsw i32 %53, %47
br label %56
56: ; preds = %54, %56
%57 = phi i32 [ %92, %56 ], [ %55, %54 ]
%58 = phi i32 [ %93, %56 ], [ 0, %54 ]
%59 = srem i32 %57, %13
%60 = add nsw i32 %59, %46
%61 = mul nsw i32 %58, %2
%62 = sdiv i32 %61, %47
%63 = shl nsw i32 %62, 1
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double, double* %3, i64 %64
%66 = load double, double* %65
%67 = or i32 %63, 1
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %3, i64 %68
%70 = load double, double* %69
%71 = fneg double %70
%72 = sext i32 %60 to i64
%73 = getelementptr inbounds double, double* %0, i64 %72
%74 = load double, double* %73
%75 = getelementptr inbounds double, double* %1, i64 %72
%76 = load double, double* %75
%77 = fmul double %70, %76
%78 = tail call double @llvm.fmuladd.f64(double %66, double %74, double %77)
%79 = fmul double %74, %71
%80 = tail call double @llvm.fmuladd.f64(double %66, double %76, double %79)
%81 = sext i32 %59 to i64
%82 = getelementptr inbounds double, double* %0, i64 %81
%83 = load double, double* %82
%84 = fsub double %83, %78
store double %84, double* %73
%85 = getelementptr inbounds double, double* %1, i64 %81
%86 = load double, double* %85
%87 = fsub double %86, %80
store double %87, double* %75
%88 = load double, double* %82
%89 = fadd double %78, %88
store double %89, double* %82
%90 = load double, double* %85
%91 = fadd double %80, %90
store double %91, double* %85
%92 = add nsw i32 %59, 1
%93 = add nuw nsw i32 %58, 1
%94 = icmp eq i32 %93, %46
br i1 %94, label %95, label %56
95: ; preds = %56, %52
%96 = add nuw nsw i32 %53, 1
%97 = icmp eq i32 %96, %48
br i1 %97, label %98, label %52
98: ; preds = %95, %44
%99 = add nuw nsw i32 %45, 1
%100 = icmp eq i32 %99, %7
br i1 %100, label %101, label %44
101: ; preds = %98, %4, %17
ret void
}