#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*/ fht(double *);
void psycho_1_fft(double *x_real, double *energy, int N) {
  double a, b;
  int i, j;
  fht(x_real);
  energy[0] = x_real[0] * x_real[0];
  for (i = 1, j = N - 1; i < N / 2; i++, j--) {
    a = x_real[i];
    b = x_real[j];
    energy[i] = (a * a + b * b) / 2.0;
  }
  energy[N / 2] = x_real[N / 2] * x_real[N / 2];
}


x86 assembly (clang, -O0):
<psycho_1_fft>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x8(%rbp),%rdi
call   1030 <fht@plt>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x2c(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x34(%rbp),%eax
cmp    %ecx,%eax
jge    11e2 <psycho_1_fft+0xd2>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xe46(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    1154 <psycho_1_fft+0x44>
mov    -0x8(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x50(%rbp),%rax
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x48(%rbp),%rax
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<psycho_1_fft>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %edx,%ebp
mov    %rsi,%r14
mov    %rdi,%rbx
call   1030 <fht@plt>
movsd  (%rbx),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,(%r14)
mov    %ebp,%eax
shr    $0x1f,%eax
add    %ebp,%eax
sar    %eax
cmp    $0x4,%ebp
jl     12fc <psycho_1_fft+0x1ec>
cmp    $0x3,%eax
mov    $0x2,%ecx
cmovge %eax,%ecx
mov    %ebp,%r9d
lea    -0x1(%rcx),%r8
mov    $0x1,%esi
cmp    $0x4,%r8
jb     1244 <psycho_1_fft+0x134>
lea    0x8(%r14),%rdx
lea    (%r14,%rcx,8),%rdi
mov    %r9,%rbp
sub    %rcx,%rbp
lea    (%rbx,%rbp,8),%r10
add    $0x8,%r10
lea    (%rbx,%r9,8),%rbp
lea    0x8(%rbx),%r11
lea    (%rbx,%rcx,8),%r15
cmp    %rbp,%rdx
setb   %bpl
cmp    %rdi,%r10
setb   %r10b
cmp    %r15,%rdx
setb   %dl
cmp    %rdi,%r11
setb   %dil
test   %r10b,%bpl
jne    1244 <psycho_1_fft+0x134>
and    %dil,%dl
jne    1244 <psycho_1_fft+0x134>
mov    %r8,%rdi
and    $0xfffffffffffffffc,%rdi
lea    0x1(%rdi),%rsi
lea    (%rbx,%r9,8),%rbp
add    $0xfffffffffffffff0,%rbp
sub    %rdi,%r9
xor    %edx,%edx
movapd 0xe2d(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd 0x8(%rbx,%rdx,8),%xmm1
movupd 0x18(%rbx,%rdx,8),%xmm2
movupd -0x10(%rbp),%xmm3
movupd 0x0(%rbp),%xmm4
mulpd  %xmm4,%xmm4
shufpd $0x1,%xmm4,%xmm4
mulpd  %xmm3,%xmm3
shufpd $0x1,%xmm3,%xmm3
mulpd  %xmm1,%xmm1
addpd  %xmm4,%xmm1
mulpd  %xmm2,%xmm2
addpd  %xmm3,%xmm2
mulpd  %xmm0,%xmm1
mulpd  %xmm0,%xmm2
movupd %xmm1,0x8(%r14,%rdx,8)
movupd %xmm2,0x18(%r14,%rdx,8)
add    $0x4,%rdx
add    $0xffffffffffffffe0,%rbp
cmp    %rdx,%rdi
jne    11e0 <psycho_1_fft+0xd0>
cmp    %rdi,%r8
je     12fc <psycho_1_fft+0x1ec>
mov    %rsi,%rdx
not    %rdx
test   $0x1,%cl
jne    127d <psycho_1_fft+0x16d>
movsd  (%rbx,%rsi,8),%xmm0
movsd  -0x8(%rbx,%r9,8),%xmm1
add    $0xffffffffffffffff,%r9
mulsd  %xmm1,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
mulsd  0xd9d(%rip),%xmm0
movsd  %xmm0,(%r14,%rsi,8)
add    $0x1,%rsi
add    %rcx,%rdx
je     12fc <psycho_1_fft+0x1ec>
sub    %rsi,%rcx
lea    (%r14,%rsi,8),%rdx
add    $0x8,%rdx
lea    (%rbx,%rsi,8),%rsi
add    $0x8,%rsi
lea    (%rbx,%r9,8),%rdi
add    $0xfffffffffffffff8,%rdi
xor    %ebp,%ebp
movsd  0xd69(%rip),%xmm0
nopw   0x0(%rax,%rax,1)
movsd  -0x8(%rsi,%rbp,8),%xmm1
movsd  (%rdi),%xmm2
mulsd  %xmm2,%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%rdx,%rbp,8)
movsd  (%rsi,%rbp,8),%xmm1
movsd  -0x8(%rdi),%xmm2
mulsd  %xmm2,%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rbp,8)
add    $0x2,%rbp
add    $0xfffffffffffffff0,%rdi
cmp    %rbp,%rcx
jne    12b0 <psycho_1_fft+0x1a0>
cltq
movsd  (%rbx,%rax,8),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,(%r14,%rax,8)
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @psycho_1_fft(double* %0, double* %1, i32 %2) {
%4 = alloca double*
%5 = alloca double*
%6 = alloca i32
%7 = alloca double
%8 = alloca double
%9 = alloca i32
%10 = alloca i32
store double* %0, double** %4
store double* %1, double** %5
store i32 %2, i32* %6
%11 = load double*, double** %4
%12 = call i32 @fht(double* %11)
%13 = load double*, double** %4
%14 = getelementptr inbounds double, double* %13, i64 0
%15 = load double, double* %14
%16 = load double*, double** %4
%17 = getelementptr inbounds double, double* %16, i64 0
%18 = load double, double* %17
%19 = fmul double %15, %18
%20 = load double*, double** %5
%21 = getelementptr inbounds double, double* %20, i64 0
store double %19, double* %21
store i32 1, i32* %9
%22 = load i32, i32* %6
%23 = sub nsw i32 %22, 1
store i32 %23, i32* %10
br label %24
24: ; preds = %51, %3
%25 = load i32, i32* %9
%26 = load i32, i32* %6
%27 = sdiv i32 %26, 2
%28 = icmp slt i32 %25, %27
br i1 %28, label %29, label %56
29: ; preds = %24
%30 = load double*, double** %4
%31 = load i32, i32* %9
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
%34 = load double, double* %33
store double %34, double* %7
%35 = load double*, double** %4
%36 = load i32, i32* %10
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %35, i64 %37
%39 = load double, double* %38
store double %39, double* %8
%40 = load double, double* %7
%41 = load double, double* %7
%42 = load double, double* %8
%43 = load double, double* %8
%44 = fmul double %42, %43
%45 = call double @llvm.fmuladd.f64(double %40, double %41, double %44)
%46 = fdiv double %45, 2.000000e+00
%47 = load double*, double** %5
%48 = load i32, i32* %9
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
store double %46, double* %50
br label %51
51: ; preds = %29
%52 = load i32, i32* %9
%53 = add nsw i32 %52, 1
store i32 %53, i32* %9
%54 = load i32, i32* %10
%55 = add nsw i32 %54, -1
store i32 %55, i32* %10
br label %24
56: ; preds = %24
%57 = load double*, double** %4
%58 = load i32, i32* %6
%59 = sdiv i32 %58, 2
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %57, i64 %60
%62 = load double, double* %61
%63 = load double*, double** %4
%64 = load i32, i32* %6
%65 = sdiv i32 %64, 2
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %63, i64 %66
%68 = load double, double* %67
%69 = fmul double %62, %68
%70 = load double*, double** %5
%71 = load i32, i32* %6
%72 = sdiv i32 %71, 2
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %70, i64 %73
store double %69, double* %74
ret void
}

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