#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 */

/* Variables and functions */
int /*<<< orphan*/ bitreverse(float *, int);
double pow(float, double);
float sin(float);
void JGcfft(float x[], int NC, int forward) {
  float wr, wi, wpr, wpi, theta, scale;
  int mmax, ND, m, i, j, delta;
  ND = NC << 1;
  bitreverse(x, ND);
  for (mmax = 2; mmax < ND; mmax = delta) {
    delta = mmax << 1;
    theta = (2.0 * 3.14159265358979323846264338327) / (forward ? mmax : -mmax);
    wpr = -2.0 * pow(sin(0.5 * theta), 2.0);
    wpi = sin(theta);
    wr = 1.0;
    wi = 0.0;
    for (m = 0; m < mmax; m += 2) {
      float rtemp, itemp;
      for (i = m; i < ND; i += delta) {
        j = i + mmax;
        rtemp = wr * x[j] - wi * x[j + 1];
        itemp = wr * x[j + 1] + wi * x[j];
        x[j] = x[i] - rtemp;
        x[j + 1] = x[i + 1] - itemp;
        x[i] += rtemp;
        x[i + 1] += itemp;
      }
      rtemp = wr;
      wr = wr * wpr - wi * wpi + wr;
      wi = wi * wpr + rtemp * wpi + wi;
    }
  }
  scale = forward ? 1.0 / ND : 2.0;
  {
    float *xi = x, *xe = x + ND;
    while (xi < xe)
      *xi++ *= scale;
  }
}


x86 assembly (clang, -O0):
<JGcfft>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x30(%rbp),%esi
call   1040 <bitreverse@plt>
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    13c4 <JGcfft+0x294>
mov    -0x2c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x0,-0x10(%rbp)
je     1188 <JGcfft+0x58>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x5c(%rbp)
jmp    1190 <JGcfft+0x60>
xor    %eax,%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xe81(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe5b(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1030 <sin@plt>
movsd  0xe2e(%rip),%xmm1
call   1050 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xe2e(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x24(%rbp),%xmm0
call   1030 <sin@plt>
movss  %xmm0,-0x20(%rbp)
movss  0xe22(%rip),%xmm0
movss  %xmm0,-0x14(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    13b4 <JGcfft+0x284>
mov    -0x34(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    1343 <JGcfft+0x213>
mov    -0x38(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x14(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x44(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x44(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    122c <JGcfft+0xfc>
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x14(%rbp),%xmm0
movss  -0x1c(%rbp),%xmm2
movss  -0x18(%rbp),%xmm1
mulss  -0x20(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
addss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x18(%rbp),%xmm0
movss  -0x1c(%rbp),%xmm2
movss  -0x44(%rbp),%xmm1
mulss  -0x20(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
addss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x34(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x34(%rbp)
jmp    121a <JGcfft+0xea>
jmp    13b9 <JGcfft+0x289>
mov    -0x40(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    115e <JGcfft+0x2e>
cmpl   $0x0,-0x10(%rbp)
je     13e9 <JGcfft+0x2b9>
cvtsi2sdl -0x30(%rbp),%xmm1
movsd  0xc2d(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    13fb <JGcfft+0x2cb>
movsd  0xc0f(%rip),%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    13fb <JGcfft+0x2cb>
movsd  -0x68(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x58(%rbp),%rax
jae    1453 <JGcfft+0x323>
movss  -0x28(%rbp),%xmm0
mov    -0x50(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x50(%rbp)
mulss  (%rax),%xmm0
movss  %xmm0,(%rax)
jmp    1424 <JGcfft+0x2f4>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<JGcfft>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %edx,%r14d
mov    %esi,%r15d
mov    %rdi,%rbx
add    %esi,%r15d
mov    %r15d,%esi
call   1040 <bitreverse@plt>
cmp    $0x3,%r15d
jl     12ce <JGcfft+0x19e>
mov    %r15d,%r12d
mov    $0x2,%ebp
jmp    1179 <JGcfft+0x49>
cs nopw 0x0(%rax,%rax,1)
cmp    %r15d,%ebp
jge    12ce <JGcfft+0x19e>
mov    %ebp,%r13d
lea    0x0(%r13,%r13,1),%ebp
mov    %r13d,%eax
neg    %eax
test   %r14d,%r14d
cmovne %r13d,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0xe64(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0xc(%rsp)
mulss  0xe6b(%rip),%xmm0
call   1030 <sin@plt>
movsd  0xe46(%rip),%xmm1
call   1050 <pow@plt>
movsd  %xmm0,0x10(%rsp)
movss  0xc(%rsp),%xmm0
call   1030 <sin@plt>
movaps 0xe50(%rip),%xmm9
test   %r13d,%r13d
jle    1170 <JGcfft+0x40>
movsd  0x10(%rsp),%xmm1
mulsd  0xe1d(%rip),%xmm1
xorps  %xmm8,%xmm8
cvtsd2ss %xmm1,%xmm8
movslq %ebp,%rax
mov    %r13d,%ecx
shufps $0x0,%xmm8,%xmm8
shufps $0x0,%xmm0,%xmm0
lea    0x1(%r13),%r8d
lea    (%rbx,%rcx,4),%rdx
xor    %ecx,%ecx
movss  0xe07(%rip),%xmm4
jmp    1244 <JGcfft+0x114>
nop
unpcklps %xmm4,%xmm3
mulps  %xmm0,%xmm3
movaps %xmm4,%xmm1
mulps  %xmm8,%xmm1
addps  %xmm3,%xmm1
addps  %xmm1,%xmm4
add    $0x2,%rcx
add    $0x2,%r8d
cmp    %ecx,%r13d
jle    1170 <JGcfft+0x40>
movaps %xmm4,%xmm1
shufps $0x55,%xmm4,%xmm1
movaps %xmm1,%xmm3
xorps  %xmm9,%xmm3
cmp    %ecx,%r15d
jle    1220 <JGcfft+0xf0>
mov    %r8d,%edi
mov    %rcx,%rsi
nopl   (%rax)
movss  (%rdx,%rsi,4),%xmm5
movslq %edi,%rdi
movss  (%rbx,%rdi,4),%xmm6
movaps %xmm6,%xmm7
mulss  %xmm3,%xmm7
movaps %xmm4,%xmm2
mulss  %xmm5,%xmm2
addss  %xmm7,%xmm2
mulss  %xmm1,%xmm5
mulss  %xmm4,%xmm6
addss  %xmm5,%xmm6
movss  (%rbx,%rsi,4),%xmm5
subss  %xmm2,%xmm5
movss  %xmm5,(%rdx,%rsi,4)
movss  0x4(%rbx,%rsi,4),%xmm5
movaps %xmm5,%xmm7
subss  %xmm6,%xmm7
movss  %xmm7,(%rbx,%rdi,4)
addss  (%rbx,%rsi,4),%xmm2
movss  %xmm2,(%rbx,%rsi,4)
addss  %xmm5,%xmm6
movss  %xmm6,0x4(%rbx,%rsi,4)
add    %rax,%rsi
add    %ebp,%edi
cmp    %r12,%rsi
jl     1260 <JGcfft+0x130>
jmp    1220 <JGcfft+0xf0>
test   %r14d,%r14d
je     12f8 <JGcfft+0x1c8>
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
movsd  0xd35(%rip),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
test   %r15d,%r15d
jg     1309 <JGcfft+0x1d9>
jmp    13f8 <JGcfft+0x2c8>
movss  0xd28(%rip),%xmm0
test   %r15d,%r15d
jle    13f8 <JGcfft+0x2c8>
movslq %r15d,%rax
lea    (%rbx,%rax,4),%rax
lea    0x4(%rbx),%rdx
cmp    %rdx,%rax
cmova  %rax,%rdx
mov    %rbx,%rcx
not    %rcx
add    %rdx,%rcx
cmp    $0x1c,%rcx
jb     13e0 <JGcfft+0x2b0>
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rcx,%rdx
and    $0xfffffffffffffff8,%rdx
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
lea    -0x8(%rdx),%rsi
mov    %rsi,%rbp
shr    $0x3,%rbp
add    $0x1,%rbp
test   %rsi,%rsi
je     1407 <JGcfft+0x2d7>
mov    %rbp,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movups (%rbx,%rsi,4),%xmm2
movups 0x10(%rbx,%rsi,4),%xmm3
movups 0x20(%rbx,%rsi,4),%xmm4
movups 0x30(%rbx,%rsi,4),%xmm5
mulps  %xmm1,%xmm2
mulps  %xmm1,%xmm3
movups %xmm2,(%rbx,%rsi,4)
movups %xmm3,0x10(%rbx,%rsi,4)
mulps  %xmm1,%xmm4
mulps  %xmm1,%xmm5
movups %xmm4,0x20(%rbx,%rsi,4)
movups %xmm5,0x30(%rbx,%rsi,4)
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1370 <JGcfft+0x240>
test   $0x1,%bpl
je     13ca <JGcfft+0x29a>
movups (%rbx,%rsi,4),%xmm2
movups 0x10(%rbx,%rsi,4),%xmm3
mulps  %xmm1,%xmm2
mulps  %xmm1,%xmm3
movups %xmm2,(%rbx,%rsi,4)
movups %xmm3,0x10(%rbx,%rsi,4)
cmp    %rdx,%rcx
je     13f8 <JGcfft+0x2c8>
lea    (%rbx,%rdx,4),%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  (%rbx),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rbx)
lea    0x4(%rbx),%rcx
mov    %rcx,%rbx
cmp    %rax,%rcx
jb     13e0 <JGcfft+0x2b0>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %esi,%esi
test   $0x1,%bpl
jne    13b2 <JGcfft+0x282>
jmp    13ca <JGcfft+0x29a>

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

LLVM IR (-O2):
define void @JGcfft(float* %0, i32 %1, i32 %2) {
%4 = ptrtoint float* %0 to i64
%5 = shl i32 %1, 1
%6 = tail call i32 @bitreverse(float* %0, i32 %5)
%7 = icmp sgt i32 %5, 2
br i1 %7, label %8, label %80
8: ; preds = %3
%9 = icmp eq i32 %2, 0
%10 = zext i32 %5 to i64
br label %13
11: ; preds = %71, %13
%12 = icmp slt i32 %15, %5
br i1 %12, label %13, label %80
13: ; preds = %8, %11
%14 = phi i32 [ 2, %8 ], [ %15, %11 ]
%15 = shl i32 %14, 1
%16 = sub nsw i32 0, %14
%17 = select i1 %9, i32 %16, i32 %14
%18 = sitofp i32 %17 to double
%19 = fdiv double 0x401921FB54442D18, %18
%20 = fptrunc double %19 to float
%21 = fmul float %20, 5.000000e-01
%22 = tail call float @sin(float %21)
%23 = tail call double @pow(float %22, double 2.000000e+00)
%24 = tail call float @sin(float %20)
%25 = icmp sgt i32 %14, 0
br i1 %25, label %26, label %11
26: ; preds = %13
%27 = fmul double %23, -2.000000e+00
%28 = fptrunc double %27 to float
%29 = sext i32 %15 to i64
%30 = zext i32 %14 to i64
%31 = insertelement <2 x float> poison, float %28, i64 0
%32 = shufflevector <2 x float> %31, <2 x float> poison, <2 x i32> zeroinitializer
%33 = insertelement <2 x float> poison, float %24, i64 0
%34 = shufflevector <2 x float> %33, <2 x float> poison, <2 x i32> zeroinitializer
br label %35
35: ; preds = %26, %71
%36 = phi i64 [ 0, %26 ], [ %77, %71 ]
%37 = phi <2 x float> [ <float 1.000000e+00, float 0.000000e+00>, %26 ], [ %76, %71 ]
%38 = trunc i64 %36 to i32
%39 = icmp sgt i32 %5, %38
%40 = extractelement <2 x float> %37, i64 1
%41 = fneg float %40
br i1 %39, label %42, label %71
42: ; preds = %35
%43 = extractelement <2 x float> %37, i64 0
br label %44
44: ; preds = %42, %44
%45 = phi i64 [ %69, %44 ], [ %36, %42 ]
%46 = trunc i64 %45 to i32
%47 = add nsw i64 %45, %30
%48 = add nsw i32 %14, %46
%49 = getelementptr inbounds float, float* %0, i64 %47
%50 = load float, float* %49
%51 = or i32 %48, 1
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds float, float* %0, i64 %52
%54 = load float, float* %53
%55 = fmul float %54, %41
%56 = tail call float @llvm.fmuladd.f32(float %43, float %50, float %55)
%57 = fmul float %40, %50
%58 = tail call float @llvm.fmuladd.f32(float %43, float %54, float %57)
%59 = getelementptr inbounds float, float* %0, i64 %45
%60 = load float, float* %59
%61 = fsub float %60, %56
store float %61, float* %49
%62 = or i64 %45, 1
%63 = getelementptr inbounds float, float* %0, i64 %62
%64 = load float, float* %63
%65 = fsub float %64, %58
store float %65, float* %53
%66 = load float, float* %59
%67 = fadd float %56, %66
store float %67, float* %59
%68 = fadd float %58, %64
store float %68, float* %63
%69 = add nsw i64 %45, %29
%70 = icmp slt i64 %69, %10
br i1 %70, label %44, label %71
71: ; preds = %44, %35
%72 = insertelement <2 x float> poison, float %41, i64 0
%73 = shufflevector <2 x float> %72, <2 x float> %37, <2 x i32> <i32 0, i32 2>
%74 = fmul <2 x float> %34, %73
%75 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %37, <2 x float> %32, <2 x float> %74)
%76 = fadd <2 x float> %37, %75
%77 = add nuw nsw i64 %36, 2
%78 = trunc i64 %77 to i32
%79 = icmp sgt i32 %14, %78
br i1 %79, label %35, label %11
80: ; preds = %11, %3
%81 = icmp eq i32 %2, 0
br i1 %81, label %86, label %82
82: ; preds = %80
%83 = sitofp i32 %5 to double
%84 = fdiv double 1.000000e+00, %83
%85 = fptrunc double %84 to float
br label %86
86: ; preds = %80, %82
%87 = phi float [ %85, %82 ], [ 2.000000e+00, %80 ]
%88 = sext i32 %5 to i64
%89 = getelementptr inbounds float, float* %0, i64 %88
%90 = icmp sgt i32 %5, 0
br i1 %90, label %91, label %166
91: ; preds = %86
%92 = shl nsw i64 %88, 2
%93 = add i64 %92, %4
%94 = add i64 %4, 4
%95 = call i64 @llvm.umax.i64(i64 %93, i64 %94)
%96 = xor i64 %4, -1
%97 = add i64 %95, %96
%98 = lshr i64 %97, 2
%99 = add nuw nsw i64 %98, 1
%100 = icmp ult i64 %97, 28
br i1 %100, label %158, label %101
101: ; preds = %91
%102 = and i64 %99, 9223372036854775800
%103 = getelementptr float, float* %0, i64 %102
%104 = insertelement <4 x float> poison, float %87, i64 0
%105 = shufflevector <4 x float> %104, <4 x float> poison, <4 x i32> zeroinitializer
%106 = insertelement <4 x float> poison, float %87, i64 0
%107 = shufflevector <4 x float> %106, <4 x float> poison, <4 x i32> zeroinitializer
%108 = add nsw i64 %102, -8
%109 = lshr exact i64 %108, 3
%110 = add nuw nsw i64 %109, 1
%111 = and i64 %110, 1
%112 = icmp eq i64 %108, 0
br i1 %112, label %142, label %113
113: ; preds = %101
%114 = and i64 %110, 4611686018427387902
br label %115
115: ; preds = %115, %113
%116 = phi i64 [ 0, %113 ], [ %139, %115 ]
%117 = phi i64 [ 0, %113 ], [ %140, %115 ]
%118 = getelementptr float, float* %0, i64 %116
%119 = bitcast float* %118 to <4 x float>*
%120 = load <4 x float>, <4 x float>* %119
%121 = getelementptr float, float* %118, i64 4
%122 = bitcast float* %121 to <4 x float>*
%123 = load <4 x float>, <4 x float>* %122
%124 = fmul <4 x float> %105, %120
%125 = fmul <4 x float> %107, %123
%126 = bitcast float* %118 to <4 x float>*
store <4 x float> %124, <4 x float>* %126
%127 = bitcast float* %121 to <4 x float>*
store <4 x float> %125, <4 x float>* %127
%128 = or i64 %116, 8
%129 = getelementptr float, float* %0, i64 %128
%130 = bitcast float* %129 to <4 x float>*
%131 = load <4 x float>, <4 x float>* %130
%132 = getelementptr float, float* %129, i64 4
%133 = bitcast float* %132 to <4 x float>*
%134 = load <4 x float>, <4 x float>* %133
%135 = fmul <4 x float> %105, %131
%136 = fmul <4 x float> %107, %134
%137 = bitcast float* %129 to <4 x float>*
store <4 x float> %135, <4 x float>* %137
%138 = bitcast float* %132 to <4 x float>*
store <4 x float> %136, <4 x float>* %138
%139 = add nuw i64 %116, 16
%140 = add i64 %117, 2
%141 = icmp eq i64 %140, %114
br i1 %141, label %142, label %115
142: ; preds = %115, %101
%143 = phi i64 [ 0, %101 ], [ %139, %115 ]
%144 = icmp eq i64 %111, 0
br i1 %144, label %156, label %145
145: ; preds = %142
%146 = getelementptr float, float* %0, i64 %143
%147 = bitcast float* %146 to <4 x float>*
%148 = load <4 x float>, <4 x float>* %147
%149 = getelementptr float, float* %146, i64 4
%150 = bitcast float* %149 to <4 x float>*
%151 = load <4 x float>, <4 x float>* %150
%152 = fmul <4 x float> %105, %148
%153 = fmul <4 x float> %107, %151
%154 = bitcast float* %146 to <4 x float>*
store <4 x float> %152, <4 x float>* %154
%155 = bitcast float* %149 to <4 x float>*
store <4 x float> %153, <4 x float>* %155
br label %156
156: ; preds = %142, %145
%157 = icmp eq i64 %99, %102
br i1 %157, label %166, label %158
158: ; preds = %91, %156
%159 = phi float* [ %0, %91 ], [ %103, %156 ]
br label %160
160: ; preds = %158, %160
%161 = phi float* [ %162, %160 ], [ %159, %158 ]
%162 = getelementptr inbounds float, float* %161, i64 1
%163 = load float, float* %161
%164 = fmul float %87, %163
store float %164, float* %161
%165 = icmp ult float* %162, %89
br i1 %165, label %160, label %166
166: ; preds = %160, %156, %86
ret void
}