void hc_fftw_to_mx(register double *Xr, register double *Xi, register double *X,
                   long N) {
  register long k;
  register double re, im;
  if (N % 2) {
    Xr[0] = X[0];
    Xi[0] = 0.0;
    for (k = 1; k < ((N + 1) / 2); ++k) {
      re = X[k];
      im = X[N - k];
      Xr[k] = re;
      Xi[k] = im;
      Xr[N - k] = re;
      Xi[N - k] = -im;
    }
  } else {
    Xr[0] = X[0];
    Xi[0] = 0.0;
    for (k = 1; k < (N / 2); ++k) {
      re = X[k];
      im = X[N - k];
      Xr[k] = re;
      Xi[k] = im;
      Xr[N - k] = re;
      Xi[N - k] = -im;
    }
    Xr[N / 2] = X[N / 2];
    Xi[N / 2] = 0.0;
  }
}


x86 assembly (clang, -O0):
<hc_fftw_to_mx>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x20(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
cmp    $0x0,%rdx
je     121e <hc_fftw_to_mx+0x11e>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movq   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x40(%rbp),%rax
cmp    %rcx,%rax
jge    1219 <hc_fftw_to_mx+0x119>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    114f <hc_fftw_to_mx+0x4f>
jmp    1370 <hc_fftw_to_mx+0x270>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movq   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x48(%rbp),%rax
cmp    %rcx,%rax
jge    1307 <hc_fftw_to_mx+0x207>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
sub    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1241 <hc_fftw_to_mx+0x141>
mov    -0x18(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x20(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x60(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x20(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x58(%rbp),%rax
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x20(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x50(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<hc_fftw_to_mx>:
push   %rbx
movsd  (%rdx),%xmm0
movsd  %xmm0,(%rdi)
movq   $0x0,(%rsi)
test   $0x1,%cl
jne    121b <hc_fftw_to_mx+0x11b>
mov    %rcx,%r9
shr    $0x3f,%r9
add    %rcx,%r9
sar    %r9
cmp    $0x4,%rcx
jl     1202 <hc_fftw_to_mx+0x102>
cmp    $0x3,%r9
mov    $0x2,%r8d
cmovge %r9,%r8
add    $0xffffffffffffffff,%r8
mov    $0x1,%r10d
cmp    $0x6,%rcx
jl     11d2 <hc_fftw_to_mx+0xd2>
lea    -0x1(%rcx),%rax
mov    %r8,%r11
and    $0xfffffffffffffffe,%r11
neg    %r11
mov    $0x1,%r10d
movaps 0xe93(%rip),%xmm0
nopl   (%rax)
movsd  (%rdx,%r10,8),%xmm1
movsd  (%rdx,%rax,8),%xmm2
movsd  %xmm1,(%rdi,%r10,8)
movsd  %xmm2,(%rsi,%r10,8)
movsd  %xmm1,(%rdi,%rax,8)
xorps  %xmm0,%xmm2
movlps %xmm2,(%rsi,%rax,8)
movsd  0x8(%rdx,%r10,8),%xmm1
movsd  -0x8(%rdx,%rax,8),%xmm2
movsd  %xmm1,0x8(%rdi,%r10,8)
movsd  %xmm2,0x8(%rsi,%r10,8)
movsd  %xmm1,-0x8(%rdi,%rax,8)
xorps  %xmm0,%xmm2
movlps %xmm2,-0x8(%rsi,%rax,8)
add    $0xfffffffffffffffe,%rax
lea    (%r11,%r10,1),%rbx
add    $0x2,%rbx
add    $0x2,%r10
cmp    $0x1,%rbx
jne    1170 <hc_fftw_to_mx+0x70>
test   $0x1,%r8b
je     1202 <hc_fftw_to_mx+0x102>
movsd  (%rdx,%r10,8),%xmm0
sub    %r10,%rcx
movsd  (%rdx,%rcx,8),%xmm1
movsd  %xmm0,(%rdi,%r10,8)
movsd  %xmm1,(%rsi,%r10,8)
movsd  %xmm0,(%rdi,%rcx,8)
xorps  0xe02(%rip),%xmm1
movlps %xmm1,(%rsi,%rcx,8)
movsd  (%rdx,%r9,8),%xmm0
movsd  %xmm0,(%rdi,%r9,8)
movq   $0x0,(%rsi,%r9,8)
jmp    12c8 <hc_fftw_to_mx+0x1c8>
cmp    $0x3,%rcx
jl     12c8 <hc_fftw_to_mx+0x1c8>
lea    0x1(%rcx),%rax
shr    %rax
lea    -0x1(%rax),%r8
mov    $0x1,%r9d
cmp    $0x2,%rax
je     12c2 <hc_fftw_to_mx+0x1c2>
lea    -0x1(%rcx),%rax
mov    %r8,%r10
and    $0xfffffffffffffffe,%r10
neg    %r10
mov    $0x1,%r9d
movaps 0xda5(%rip),%xmm0
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%r9,8),%xmm1
movsd  (%rdx,%rax,8),%xmm2
movsd  %xmm1,(%rdi,%r9,8)
movsd  %xmm2,(%rsi,%r9,8)
movsd  %xmm1,(%rdi,%rax,8)
xorps  %xmm0,%xmm2
movlps %xmm2,(%rsi,%rax,8)
movsd  0x8(%rdx,%r9,8),%xmm1
movsd  -0x8(%rdx,%rax,8),%xmm2
movsd  %xmm1,0x8(%rdi,%r9,8)
movsd  %xmm2,0x8(%rsi,%r9,8)
movsd  %xmm1,-0x8(%rdi,%rax,8)
xorps  %xmm0,%xmm2
movlps %xmm2,-0x8(%rsi,%rax,8)
add    $0xfffffffffffffffe,%rax
lea    (%r10,%r9,1),%r11
add    $0x2,%r11
add    $0x2,%r9
cmp    $0x1,%r11
jne    1260 <hc_fftw_to_mx+0x160>
test   $0x1,%r8b
jne    12ca <hc_fftw_to_mx+0x1ca>
pop    %rbx
ret
movsd  (%rdx,%r9,8),%xmm0
sub    %r9,%rcx
movsd  (%rdx,%rcx,8),%xmm1
movsd  %xmm0,(%rdi,%r9,8)
movsd  %xmm1,(%rsi,%r9,8)
movsd  %xmm0,(%rdi,%rcx,8)
xorps  0xd10(%rip),%xmm1
movlps %xmm1,(%rsi,%rcx,8)
pop    %rbx
ret

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

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