void fft_pass12_1(double *ri, int n) {
  int j;
  double *q = ri + n;
  for (j = 0; j < n; j += 4) {
    double x00 = q[j], x01 = q[j + 1], x02 = q[j + 2], x03 = q[j + 3],
           x0 = x00 + x01, x1 = x00 - x01, x2 = x02 + x03, x3 = x02 - x03;
    ri[2 * (j)] = x0 + x2;
    ri[2 * (j) + 1] = 0.0;
    ri[2 * (j + 2)] = x0 - x2;
    ri[2 * (j + 2) + 1] = 0.0;
    ri[2 * (j + 1)] = ri[2 * (j + 3)] = x1;
    ri[2 * (j + 1) + 1] = x3;
    ri[2 * (j + 3) + 1] = -x3;
  }
}


x86 assembly (clang, -O0):
<fft_pass12_1>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12c0 <fft_pass12_1+0x1c0>
mov    -0x18(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x20(%rbp),%xmm0
subsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x2,%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x2,%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x3,%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    $0x3,%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x10(%rbp)
jmp    1125 <fft_pass12_1+0x25>
pop    %rbp
ret

x86 assembly (clang, -O2):
<fft_pass12_1>:
test   %esi,%esi
jle    11a9 <fft_pass12_1+0xa9>
movslq %esi,%rax
lea    (%rdi,%rax,8),%rcx
add    $0x18,%rcx
add    $0x38,%rdi
xor    %edx,%edx
movapd 0xedf(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x18(%rcx,%rdx,8),%xmm1
movsd  -0x10(%rcx,%rdx,8),%xmm2
movsd  -0x8(%rcx,%rdx,8),%xmm3
movsd  (%rcx,%rdx,8),%xmm4
movapd %xmm1,%xmm5
addsd  %xmm2,%xmm5
subsd  %xmm2,%xmm1
movapd %xmm3,%xmm2
addsd  %xmm4,%xmm2
subsd  %xmm4,%xmm3
movapd %xmm5,%xmm4
addsd  %xmm2,%xmm4
movsd  %xmm4,-0x38(%rdi)
movq   $0x0,-0x30(%rdi)
subsd  %xmm2,%xmm5
movsd  %xmm5,-0x18(%rdi)
movq   $0x0,-0x10(%rdi)
movsd  %xmm1,-0x8(%rdi)
movsd  %xmm1,-0x28(%rdi)
movsd  %xmm3,-0x20(%rdi)
xorpd  %xmm0,%xmm3
movlpd %xmm3,(%rdi)
add    $0x4,%rdx
add    $0x40,%rdi
cmp    %rax,%rdx
jl     1130 <fft_pass12_1+0x30>
ret

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

LLVM IR (-O2):
define void @fft_pass12_1(double* nocapture %0, i32 %1) {
%3 = sext i32 %1 to i64
%4 = getelementptr inbounds double, double* %0, i64 %3
%5 = icmp sgt i32 %1, 0
br i1 %5, label %6, label %48
6: ; preds = %2, %6
%7 = phi i64 [ %46, %6 ], [ 0, %2 ]
%8 = getelementptr inbounds double, double* %4, i64 %7
%9 = load double, double* %8
%10 = or i64 %7, 1
%11 = getelementptr inbounds double, double* %4, i64 %10
%12 = load double, double* %11
%13 = or i64 %7, 2
%14 = getelementptr inbounds double, double* %4, i64 %13
%15 = load double, double* %14
%16 = or i64 %7, 3
%17 = getelementptr inbounds double, double* %4, i64 %16
%18 = load double, double* %17
%19 = fadd double %9, %12
%20 = fsub double %9, %12
%21 = fadd double %15, %18
%22 = fsub double %15, %18
%23 = fadd double %19, %21
%24 = shl nuw nsw i64 %7, 1
%25 = getelementptr inbounds double, double* %0, i64 %24
store double %23, double* %25
%26 = or i64 %24, 1
%27 = getelementptr inbounds double, double* %0, i64 %26
store double 0.000000e+00, double* %27
%28 = fsub double %19, %21
%29 = shl nuw nsw i64 %13, 1
%30 = getelementptr inbounds double, double* %0, i64 %29
store double %28, double* %30
%31 = or i64 %29, 1
%32 = getelementptr inbounds double, double* %0, i64 %31
store double 0.000000e+00, double* %32
%33 = shl nuw nsw i64 %16, 1
%34 = and i64 %33, 4294967294
%35 = getelementptr inbounds double, double* %0, i64 %34
store double %20, double* %35
%36 = shl nuw nsw i64 %10, 1
%37 = and i64 %36, 4294967290
%38 = getelementptr inbounds double, double* %0, i64 %37
store double %20, double* %38
%39 = and i64 %36, 4294967290
%40 = or i64 %39, 1
%41 = getelementptr inbounds double, double* %0, i64 %40
store double %22, double* %41
%42 = fneg double %22
%43 = and i64 %33, 4294967294
%44 = or i64 %43, 1
%45 = getelementptr inbounds double, double* %0, i64 %44
store double %42, double* %45
%46 = add nuw nsw i64 %7, 4
%47 = icmp slt i64 %46, %3
br i1 %47, label %6, label %48
48: ; preds = %6, %2
ret void
}