void FFT_2(int sofarRadix, int radix, int remainRadix, int groupOffset,
           double yRe[], double yIm[]) {
  int yIndx[2];
  int groupNo, blockNo, adr;
  double ar0, ar1;
  double ai0, ai1;
  adr = groupOffset;
  for (groupNo = 0; groupNo < remainRadix; groupNo++) {
    for (blockNo = 0; blockNo < radix; blockNo++) {
      yIndx[blockNo] = adr;
      adr = adr + sofarRadix;
    }
    ar0 = yRe[yIndx[0]];
    ai0 = yIm[yIndx[0]];
    ar1 = yRe[yIndx[1]];
    ai1 = yIm[yIndx[1]];
    yRe[yIndx[0]] = ar0 + ar1;
    yRe[yIndx[1]] = ar0 - ar1;
    yIm[yIndx[0]] = ai0 + ai1;
    yIm[yIndx[1]] = ai0 - ai1;
    groupOffset = groupOffset + sofarRadix * radix;
    adr = groupOffset;
  }
}


x86 assembly (clang, -O0):
<FFT_2>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    122d <FFT_2+0x12d>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1166 <FFT_2+0x66>
mov    -0x34(%rbp),%ecx
movslq -0x30(%rbp),%rax
mov    %ecx,-0x28(%rbp,%rax,4)
mov    -0x34(%rbp),%eax
add    -0x4(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1138 <FFT_2+0x38>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
subsd  -0x58(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1125 <FFT_2+0x25>
pop    %rbp
ret

x86 assembly (clang, -O2):
<FFT_2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    1304 <FFT_2+0x204>
mov    %esi,%r10d
imul   %edi,%r10d
mov    %esi,%eax
mov    %eax,%r11d
and    $0xfffffff8,%r11d
lea    -0x8(%r11),%rbx
mov    %rbx,-0x8(%rsp)
shr    $0x3,%rbx
add    $0x1,%rbx
mov    %r11d,%ebp
imul   %edi,%ebp
mov    %rbp,-0x20(%rsp)
movd   %edi,%xmm0
pshufd $0x0,%xmm0,%xmm0
movdqa 0xeb1(%rip),%xmm1
pmuludq %xmm0,%xmm1
pshufd $0xe8,%xmm1,%xmm2
movdqa 0xeb0(%rip),%xmm1
pmuludq %xmm0,%xmm1
pshufd $0xe8,%xmm1,%xmm1
punpckldq %xmm2,%xmm1
pslld  $0x2,%xmm0
mov    %ebx,%r12d
and    $0x3,%r12d
and    $0xfffffffffffffffc,%rbx
mov    %rbx,-0x18(%rsp)
mov    %r12,-0x10(%rsp)
shl    $0x5,%r12
xor    %r13d,%r13d
jmp    11f7 <FFT_2+0xf7>
mov    -0x28(%rsp),%r15d
mov    -0x24(%rsp),%ebx
movslq %r15d,%rbp
movsd  (%r8,%rbp,8),%xmm2
movsd  (%r9,%rbp,8),%xmm3
movslq %ebx,%r14
movsd  (%r8,%r14,8),%xmm4
movsd  (%r9,%r14,8),%xmm5
movapd %xmm2,%xmm6
addsd  %xmm4,%xmm6
movsd  %xmm6,(%r8,%rbp,8)
subsd  %xmm4,%xmm2
movsd  %xmm2,(%r8,%r14,8)
movapd %xmm3,%xmm2
addsd  %xmm5,%xmm2
movsd  %xmm2,(%r9,%rbp,8)
subsd  %xmm5,%xmm3
movsd  %xmm3,(%r9,%r14,8)
add    %r10d,%ecx
add    $0x1,%r13d
cmp    %edx,%r13d
je     1304 <FFT_2+0x204>
test   %esi,%esi
jle    1199 <FFT_2+0x99>
cmp    $0x8,%esi
jae    1210 <FFT_2+0x110>
xor    %ebx,%ebx
mov    %ecx,%ebp
jmp    12f0 <FFT_2+0x1f0>
nopl   0x0(%rax)
movd   %ecx,%xmm2
pshufd $0x0,%xmm2,%xmm2
paddd  %xmm1,%xmm2
cmpq   $0x18,-0x8(%rsp)
jae    122a <FFT_2+0x12a>
xor    %r15d,%r15d
jmp    12a2 <FFT_2+0x1a2>
mov    -0x18(%rsp),%rbx
xor    %r15d,%r15d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movdqu %xmm2,-0x28(%rsp,%r15,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,-0x18(%rsp,%r15,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,-0x8(%rsp,%r15,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,0x8(%rsp,%r15,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,0x18(%rsp,%r15,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,0x28(%rsp,%r15,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,0x38(%rsp,%r15,4)
paddd  %xmm0,%xmm2
movdqu %xmm2,0x48(%rsp,%r15,4)
add    $0x20,%r15
paddd  %xmm0,%xmm2
add    $0xfffffffffffffffc,%rbx
jne    1240 <FFT_2+0x140>
cmpq   $0x0,-0x10(%rsp)
je     12dc <FFT_2+0x1dc>
lea    -0x18(%rsp),%rbx
lea    (%rbx,%r15,4),%rbx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
movdqu %xmm2,-0x10(%rbx,%rbp,1)
paddd  %xmm0,%xmm2
movdqu %xmm2,(%rbx,%rbp,1)
paddd  %xmm0,%xmm2
add    $0x20,%rbp
cmp    %rbp,%r12
jne    12c0 <FFT_2+0x1c0>
cmp    %rax,%r11
je     1190 <FFT_2+0x90>
mov    -0x20(%rsp),%rbx
lea    (%rcx,%rbx,1),%ebp
mov    %r11,%rbx
mov    %ebp,-0x28(%rsp,%rbx,4)
add    %edi,%ebp
add    $0x1,%rbx
cmp    %rbx,%rax
jne    12f0 <FFT_2+0x1f0>
jmp    1190 <FFT_2+0x90>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @FFT_2(i32 %0, i32 %1, i32 %2, i32 %3, double* nocapture %4, double* nocapture %5) {
%7 = alloca [2 x i32]
%8 = bitcast [2 x i32]* %7 to i8*
%9 = icmp sgt i32 %2, 0
br i1 %9, label %10, label %130
10: ; preds = %6
%11 = icmp sgt i32 %1, 0
%12 = getelementptr inbounds [2 x i32], [2 x i32]* %7, i64 0, i64 0
%13 = getelementptr inbounds [2 x i32], [2 x i32]* %7, i64 0, i64 1
%14 = mul nsw i32 %1, %0
%15 = zext i32 %1 to i64
%16 = and i64 %15, 4294967288
%17 = add nsw i64 %16, -8
%18 = lshr exact i64 %17, 3
%19 = add nuw nsw i64 %18, 1
%20 = icmp ult i32 %1, 8
%21 = and i64 %15, 4294967288
%22 = trunc i64 %21 to i32
%23 = mul i32 %22, %0
%24 = insertelement <4 x i32> poison, i32 %0, i64 0
%25 = shufflevector <4 x i32> %24, <4 x i32> poison, <4 x i32> zeroinitializer
%26 = mul <4 x i32> %25, <i32 0, i32 1, i32 2, i32 3>
%27 = shl i32 %0, 2
%28 = insertelement <4 x i32> poison, i32 %27, i64 0
%29 = shufflevector <4 x i32> %28, <4 x i32> poison, <4 x i32> zeroinitializer
%30 = and i64 %19, 3
%31 = icmp ult i64 %17, 24
%32 = and i64 %19, 4611686018427387900
%33 = icmp eq i64 %30, 0
%34 = icmp eq i64 %21, %15
br label %35
35: ; preds = %10, %110
%36 = phi i32 [ undef, %10 ], [ %111, %110 ]
%37 = phi i32 [ undef, %10 ], [ %112, %110 ]
%38 = phi i32 [ %3, %10 ], [ %127, %110 ]
%39 = phi i32 [ 0, %10 ], [ %128, %110 ]
br i1 %11, label %40, label %110
40: ; preds = %35
br i1 %20, label %97, label %41
41: ; preds = %40
%42 = add i32 %38, %23
%43 = insertelement <4 x i32> poison, i32 %38, i64 0
%44 = shufflevector <4 x i32> %43, <4 x i32> poison, <4 x i32> zeroinitializer
%45 = add <4 x i32> %44, %26
br i1 %31, label %80, label %46
46: ; preds = %41, %46
%47 = phi i64 [ %76, %46 ], [ 0, %41 ]
%48 = phi <4 x i32> [ %77, %46 ], [ %45, %41 ]
%49 = phi i64 [ %78, %46 ], [ 0, %41 ]
%50 = add <4 x i32> %48, %29
%51 = getelementptr inbounds [2 x i32], [2 x i32]* %7, i64 0, i64 %47
%52 = bitcast i32* %51 to <4 x i32>*
store <4 x i32> %48, <4 x i32>* %52
%53 = getelementptr inbounds i32, i32* %51, i64 4
%54 = bitcast i32* %53 to <4 x i32>*
store <4 x i32> %50, <4 x i32>* %54
%55 = or i64 %47, 8
%56 = add <4 x i32> %50, %29
%57 = add <4 x i32> %56, %29
%58 = getelementptr inbounds [2 x i32], [2 x i32]* %7, i64 0, i64 %55
%59 = bitcast i32* %58 to <4 x i32>*
store <4 x i32> %56, <4 x i32>* %59
%60 = getelementptr inbounds i32, i32* %58, i64 4
%61 = bitcast i32* %60 to <4 x i32>*
store <4 x i32> %57, <4 x i32>* %61
%62 = or i64 %47, 16
%63 = add <4 x i32> %57, %29
%64 = add <4 x i32> %63, %29
%65 = getelementptr inbounds [2 x i32], [2 x i32]* %7, i64 0, i64 %62
%66 = bitcast i32* %65 to <4 x i32>*
store <4 x i32> %63, <4 x i32>* %66
%67 = getelementptr inbounds i32, i32* %65, i64 4
%68 = bitcast i32* %67 to <4 x i32>*
store <4 x i32> %64, <4 x i32>* %68
%69 = or i64 %47, 24
%70 = add <4 x i32> %64, %29
%71 = add <4 x i32> %70, %29
%72 = getelementptr inbounds [2 x i32], [2 x i32]* %7, i64 0, i64 %69
%73 = bitcast i32* %72 to <4 x i32>*
store <4 x i32> %70, <4 x i32>* %73
%74 = getelementptr inbounds i32, i32* %72, i64 4
%75 = bitcast i32* %74 to <4 x i32>*
store <4 x i32> %71, <4 x i32>* %75
%76 = add nuw i64 %47, 32
%77 = add <4 x i32> %71, %29
%78 = add i64 %49, 4
%79 = icmp eq i64 %78, %32
br i1 %79, label %80, label %46
80: ; preds = %46, %41
%81 = phi i64 [ 0, %41 ], [ %76, %46 ]
%82 = phi <4 x i32> [ %45, %41 ], [ %77, %46 ]
br i1 %33, label %96, label %83
83: ; preds = %80, %83
%84 = phi i64 [ %92, %83 ], [ %81, %80 ]
%85 = phi <4 x i32> [ %93, %83 ], [ %82, %80 ]
%86 = phi i64 [ %94, %83 ], [ 0, %80 ]
%87 = add <4 x i32> %85, %29
%88 = getelementptr inbounds [2 x i32], [2 x i32]* %7, i64 0, i64 %84
%89 = bitcast i32* %88 to <4 x i32>*
store <4 x i32> %85, <4 x i32>* %89
%90 = getelementptr inbounds i32, i32* %88, i64 4
%91 = bitcast i32* %90 to <4 x i32>*
store <4 x i32> %87, <4 x i32>* %91
%92 = add nuw i64 %84, 8
%93 = add <4 x i32> %87, %29
%94 = add i64 %86, 1
%95 = icmp eq i64 %94, %30
br i1 %95, label %96, label %83
96: ; preds = %83, %80
br i1 %34, label %107, label %97
97: ; preds = %40, %96
%98 = phi i64 [ 0, %40 ], [ %21, %96 ]
%99 = phi i32 [ %38, %40 ], [ %42, %96 ]
br label %100
100: ; preds = %97, %100
%101 = phi i64 [ %105, %100 ], [ %98, %97 ]
%102 = phi i32 [ %104, %100 ], [ %99, %97 ]
%103 = getelementptr inbounds [2 x i32], [2 x i32]* %7, i64 0, i64 %101
store i32 %102, i32* %103
%104 = add nsw i32 %102, %0
%105 = add nuw nsw i64 %101, 1
%106 = icmp eq i64 %105, %15
br i1 %106, label %107, label %100
107: ; preds = %100, %96
%108 = load i32, i32* %12
%109 = load i32, i32* %13
br label %110
110: ; preds = %107, %35
%111 = phi i32 [ %109, %107 ], [ %36, %35 ]
%112 = phi i32 [ %108, %107 ], [ %37, %35 ]
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %4, i64 %113
%115 = load double, double* %114
%116 = getelementptr inbounds double, double* %5, i64 %113
%117 = load double, double* %116
%118 = sext i32 %111 to i64
%119 = getelementptr inbounds double, double* %4, i64 %118
%120 = load double, double* %119
%121 = getelementptr inbounds double, double* %5, i64 %118
%122 = load double, double* %121
%123 = fadd double %115, %120
store double %123, double* %114
%124 = fsub double %115, %120
store double %124, double* %119
%125 = fadd double %117, %122
store double %125, double* %116
%126 = fsub double %117, %122
store double %126, double* %121
%127 = add nsw i32 %38, %14
%128 = add nuw nsw i32 %39, 1
%129 = icmp eq i32 %128, %2
br i1 %129, label %130, label %35
130: ; preds = %110, %6
ret void
}