#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 M_PI;
double atan2(float, float);
float exp(double);
int /*<<< orphan*/ fft_(float *, float *, int *, int *, int *, int *);
int /*<<< orphan*/ printf(char *, ...);
double sqrt(float);
void transform_gaussian() {
  int n = 64;
  int one = 1, negone = -1;
  int i, ii, j, k;
  float a[64], b[64], c[64], e[64], o[64];
  float sigma = 5.;
  double q, r, t, z;
  for (i = -n / 2; i < n / 2; i++) {
    ii = (i < 0) ? i + n : i;
    z = i - 0.5;
    c[ii] = exp(-0.5 * z * z / (sigma * sigma));
    ;
  }
  if (0)
    for (k = n, i = 1; i < n / 2; i++) {
      c[--k] = -c[i];
      if (1) {
        c[i] *= i;
        c[k] *= -i;
      }
    }
  for (i = 0; i < n; i++) {
    a[i] = c[i];
    b[i] = 0.0;
  }
  for (i = 0; i < n; i++)
    printf("%10d%10.4f%10.4f\n", i, a[i], b[i]);
  printf("OK\n");
  fft_(a, b, &one, &n, &one, &negone);
  for (i = 0; i <= n / 2; i++) {
    r = sqrt(a[i] * a[i] + b[i] * b[i]);
    t = atan2(-b[i], a[i]);
    printf("%10d%10.4f%10.4f%10.4f%10.4f\n", i, a[i], b[i], r,
           t * n / (2. * M_PI));
  }
}


x86 assembly (clang, -O0):
<transform_gaussian>:
push   %rbp
mov    %rsp,%rbp
sub    $0x560,%rsp
movl   $0x40,-0x4(%rbp)
movl   $0x1,-0x8(%rbp)
movl   $0xffffffff,-0xc(%rbp)
movss  0xe88(%rip),%xmm0
movss  %xmm0,-0x524(%rbp)
ff
xor    %eax,%eax
sub    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x54c(%rbp)
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x54c(%rbp),%eax
cmp    %ecx,%eax
jge    124a <transform_gaussian+0xfa>
cmpl   $0x0,-0x10(%rbp)
jge    11cf <transform_gaussian+0x7f>
mov    -0x10(%rbp),%eax
add    -0x4(%rbp),%eax
mov    %eax,-0x550(%rbp)
jmp    11d8 <transform_gaussian+0x88>
mov    -0x10(%rbp),%eax
mov    %eax,-0x550(%rbp)
mov    -0x550(%rbp),%eax
mov    %eax,-0x14(%rbp)
cvtsi2sdl -0x10(%rbp),%xmm0
movsd  0xe2a(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x548(%rbp)
ff
movsd  0xe0e(%rip),%xmm0
mulsd  -0x548(%rbp),%xmm0
ff
mulsd  -0x548(%rbp),%xmm0
ff
movss  -0x524(%rbp),%xmm1
ff
mulss  -0x524(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
call   1050 <exp@plt>
movslq -0x14(%rbp),%rax
movss  %xmm0,-0x320(%rbp,%rax,4)
ff ff
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1190 <transform_gaussian+0x40>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1295 <transform_gaussian+0x145>
movslq -0x10(%rbp),%rax
movss  -0x320(%rbp,%rax,4),%xmm0
ff ff
movslq -0x10(%rbp),%rax
movss  %xmm0,-0x120(%rbp,%rax,4)
ff ff
movslq -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x220(%rbp,%rax,4)
ff ff
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1251 <transform_gaussian+0x101>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12e9 <transform_gaussian+0x199>
mov    -0x10(%rbp),%esi
movslq -0x10(%rbp),%rax
movss  -0x120(%rbp,%rax,4),%xmm0
ff ff
cvtss2sd %xmm0,%xmm0
movslq -0x10(%rbp),%rax
movss  -0x220(%rbp,%rax,4),%xmm1
ff ff
cvtss2sd %xmm1,%xmm1
lea    0xd4c(%rip),%rdi
mov    $0x2,%al
call   1040 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    129c <transform_gaussian+0x14c>
lea    0xd42(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    -0x120(%rbp),%rdi
lea    -0x220(%rbp),%rsi
lea    -0x8(%rbp),%r8
lea    -0x4(%rbp),%rcx
lea    -0xc(%rbp),%r9
mov    %r8,%rdx
call   1060 <fft_@plt>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x554(%rbp)
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x554(%rbp),%eax
cmp    %ecx,%eax
jg     1436 <transform_gaussian+0x2e6>
movslq -0x10(%rbp),%rax
movss  -0x120(%rbp,%rax,4),%xmm0
ff ff
movslq -0x10(%rbp),%rax
movss  -0x120(%rbp,%rax,4),%xmm2
ff ff
movslq -0x10(%rbp),%rax
movss  -0x220(%rbp,%rax,4),%xmm1
ff ff
movslq -0x10(%rbp),%rax
mulss  -0x220(%rbp,%rax,4),%xmm1
ff ff
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
call   1070 <sqrt@plt>
movsd  %xmm0,-0x538(%rbp)
ff
movslq -0x10(%rbp),%rax
movss  -0x220(%rbp,%rax,4),%xmm0
ff ff
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movslq -0x10(%rbp),%rax
movss  -0x120(%rbp,%rax,4),%xmm1
ff ff
call   1030 <atan2@plt>
movsd  %xmm0,-0x540(%rbp)
ff
mov    -0x10(%rbp),%esi
movslq -0x10(%rbp),%rax
movss  -0x120(%rbp,%rax,4),%xmm0
ff ff
cvtss2sd %xmm0,%xmm0
movslq -0x10(%rbp),%rax
movss  -0x220(%rbp,%rax,4),%xmm1
ff ff
cvtss2sd %xmm1,%xmm1
movsd  -0x538(%rbp),%xmm2
ff
movsd  -0x540(%rbp),%xmm3
ff
cvtsi2sdl -0x4(%rbp),%xmm4
mulsd  %xmm4,%xmm3
mov    0x2be2(%rip),%rax
cvtsi2sdl (%rax),%xmm5
movsd  0xbf6(%rip),%xmm4
mulsd  %xmm5,%xmm4
divsd  %xmm4,%xmm3
lea    0xc15(%rip),%rdi
mov    $0x4,%al
call   1040 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1320 <transform_gaussian+0x1d0>
add    $0x560,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<transform_gaussian>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x328,%rsp
movl   $0x40,0xc(%rsp)
movl   $0x1,0x1c(%rsp)
movl   $0xffffffff,0x18(%rsp)
ff
mov    $0xffffffe0,%ebx
movsd  0xe5e(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  %xmm1,0x10(%rsp)
mov    %ebx,%ebp
shr    $0x19,%ebp
and    $0xffffffc0,%ebp
add    %ebx,%ebp
movsd  0x10(%rsp),%xmm1
movsd  0xe3a(%rip),%xmm2
addsd  %xmm2,%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm1,%xmm0
divsd  0xe2a(%rip),%xmm0
call   1070 <exp@plt>
movsd  0x10(%rsp),%xmm1
movslq %ebp,%rax
movss  %xmm0,0x220(%rsp,%rax,4)
addsd  0xe13(%rip),%xmm1
add    $0x1,%ebx
cmp    $0x20,%ebx
jne    11b0 <transform_gaussian+0x40>
xorpd  %xmm0,%xmm0
movapd %xmm0,0x110(%rsp)
movapd %xmm0,0x100(%rsp)
movapd %xmm0,0xf0(%rsp)
movapd %xmm0,0xe0(%rsp)
movapd %xmm0,0xd0(%rsp)
movapd %xmm0,0xc0(%rsp)
movapd %xmm0,0xb0(%rsp)
movapd %xmm0,0xa0(%rsp)
movapd %xmm0,0x90(%rsp)
movapd %xmm0,0x80(%rsp)
movapd %xmm0,0x70(%rsp)
movapd %xmm0,0x60(%rsp)
movapd %xmm0,0x50(%rsp)
movapd %xmm0,0x40(%rsp)
movapd %xmm0,0x30(%rsp)
movapd %xmm0,0x20(%rsp)
lea    0x120(%rsp),%rdi
lea    0x220(%rsp),%rsi
mov    $0x100,%edx
call   1060 <memcpy@plt>
movss  0x120(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xd73(%rip),%r14
xorpd  %xmm1,%xmm1
mov    %r14,%rdi
xor    %esi,%esi
mov    $0x2,%al
call   1050 <printf@plt>
mov    $0x1,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  0x20(%rsp,%rbx,4),%xmm1
movss  0x120(%rsp,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
cvtss2sd %xmm1,%xmm1
mov    %r14,%rdi
mov    %ebx,%esi
mov    $0x2,%al
call   1050 <printf@plt>
add    $0x1,%rbx
cmp    $0x40,%rbx
jne    12e0 <transform_gaussian+0x170>
lea    0xd4c(%rip),%rdi
call   1030 <puts@plt>
lea    0x120(%rsp),%rdi
lea    0x20(%rsp),%rsi
lea    0xc(%rsp),%rcx
lea    0x1c(%rsp),%rdx
lea    0x18(%rsp),%r9
mov    %rdx,%r8
call   1080 <fft_@plt>
cmpl   $0xffffffff,0xc(%rsp)
jl     140e <transform_gaussian+0x29e>
xor    %ebx,%ebx
mov    0x2c97(%rip),%r15
lea    0xcea(%rip),%r14
nopl   0x0(%rax,%rax,1)
movss  0x120(%rsp,%rbx,4),%xmm0
movss  0x20(%rsp,%rbx,4),%xmm1
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
call   1090 <sqrt@plt>
movsd  %xmm0,0x10(%rsp)
movss  0x20(%rsp,%rbx,4),%xmm0
xorpd  0xc8c(%rip),%xmm0
movss  0x120(%rsp,%rbx,4),%xmm1
call   1040 <atan2@plt>
movss  0x120(%rsp,%rbx,4),%xmm1
xorps  %xmm2,%xmm2
cvtss2sd %xmm1,%xmm2
movss  0x20(%rsp,%rbx,4),%xmm1
cvtss2sd %xmm1,%xmm1
xorps  %xmm3,%xmm3
cvtsi2sdl 0xc(%rsp),%xmm3
mulsd  %xmm0,%xmm3
xorps  %xmm0,%xmm0
cvtsi2sdl (%r15),%xmm0
addsd  %xmm0,%xmm0
divsd  %xmm0,%xmm3
mov    %r14,%rdi
mov    %ebx,%esi
movaps %xmm2,%xmm0
movsd  0x10(%rsp),%xmm2
mov    $0x4,%al
call   1050 <printf@plt>
lea    0x1(%rbx),%rax
mov    0xc(%rsp),%ecx
mov    %ecx,%edx
shr    $0x1f,%edx
add    %ecx,%edx
sar    %edx
movslq %edx,%rcx
cmp    %rcx,%rbx
mov    %rax,%rbx
jl     1360 <transform_gaussian+0x1f0>
add    $0x328,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [18 x i8] c"%10d%10.4f%10.4f\0A\00"
@.str.1 = [4 x i8] c"OK\0A\00"
@.str.2 = [30 x i8] c"%10d%10.4f%10.4f%10.4f%10.4f\0A\00"
@M_PI = global i32 0
define void @transform_gaussian() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca [64 x float]
%9 = alloca [64 x float]
%10 = alloca [64 x float]
%11 = alloca [64 x float]
%12 = alloca [64 x float]
%13 = alloca float
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
store i32 64, i32* %1
store i32 1, i32* %2
store i32 -1, i32* %3
store float 5.000000e+00, float* %13
%18 = load i32, i32* %1
%19 = sub nsw i32 0, %18
%20 = sdiv i32 %19, 2
store i32 %20, i32* %4
br label %21
21: ; preds = %53, %0
%22 = load i32, i32* %4
%23 = load i32, i32* %1
%24 = sdiv i32 %23, 2
%25 = icmp slt i32 %22, %24
br i1 %25, label %26, label %56
26: ; preds = %21
%27 = load i32, i32* %4
%28 = icmp slt i32 %27, 0
br i1 %28, label %29, label %33
29: ; preds = %26
%30 = load i32, i32* %4
%31 = load i32, i32* %1
%32 = add nsw i32 %30, %31
br label %35
33: ; preds = %26
%34 = load i32, i32* %4
br label %35
35: ; preds = %33, %29
%36 = phi i32 [ %32, %29 ], [ %34, %33 ]
store i32 %36, i32* %5
%37 = load i32, i32* %4
%38 = sitofp i32 %37 to double
%39 = fsub double %38, 5.000000e-01
store double %39, double* %17
%40 = load double, double* %17
%41 = fmul double -5.000000e-01, %40
%42 = load double, double* %17
%43 = fmul double %41, %42
%44 = load float, float* %13
%45 = load float, float* %13
%46 = fmul float %44, %45
%47 = fpext float %46 to double
%48 = fdiv double %43, %47
%49 = call float @exp(double %48)
%50 = load i32, i32* %5
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds [64 x float], [64 x float]* %10, i64 0, i64 %51
store float %49, float* %52
br label %53
53: ; preds = %35
%54 = load i32, i32* %4
%55 = add nsw i32 %54, 1
store i32 %55, i32* %4
br label %21
56: ; preds = %21
store i32 0, i32* %4
br label %57
57: ; preds = %72, %56
%58 = load i32, i32* %4
%59 = load i32, i32* %1
%60 = icmp slt i32 %58, %59
br i1 %60, label %61, label %75
61: ; preds = %57
%62 = load i32, i32* %4
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds [64 x float], [64 x float]* %10, i64 0, i64 %63
%65 = load float, float* %64
%66 = load i32, i32* %4
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %67
store float %65, float* %68
%69 = load i32, i32* %4
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %70
store float 0.000000e+00, float* %71
br label %72
72: ; preds = %61
%73 = load i32, i32* %4
%74 = add nsw i32 %73, 1
store i32 %74, i32* %4
br label %57
75: ; preds = %57
store i32 0, i32* %4
br label %76
76: ; preds = %93, %75
%77 = load i32, i32* %4
%78 = load i32, i32* %1
%79 = icmp slt i32 %77, %78
br i1 %79, label %80, label %96
80: ; preds = %76
%81 = load i32, i32* %4
%82 = load i32, i32* %4
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %83
%85 = load float, float* %84
%86 = fpext float %85 to double
%87 = load i32, i32* %4
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %88
%90 = load float, float* %89
%91 = fpext float %90 to double
%92 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str, i64 0, i64 0), i32 %81, double %86, double %91)
br label %93
93: ; preds = %80
%94 = load i32, i32* %4
%95 = add nsw i32 %94, 1
store i32 %95, i32* %4
br label %76
96: ; preds = %76
%97 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.1, i64 0, i64 0))
%98 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 0
%99 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 0
%100 = call i32 @fft_(float* %98, float* %99, i32* %2, i32* %1, i32* %2, i32* %3)
store i32 0, i32* %4
br label %101
101: ; preds = %157, %96
%102 = load i32, i32* %4
%103 = load i32, i32* %1
%104 = sdiv i32 %103, 2
%105 = icmp sle i32 %102, %104
br i1 %105, label %106, label %160
106: ; preds = %101
%107 = load i32, i32* %4
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %108
%110 = load float, float* %109
%111 = load i32, i32* %4
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %112
%114 = load float, float* %113
%115 = load i32, i32* %4
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %116
%118 = load float, float* %117
%119 = load i32, i32* %4
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %120
%122 = load float, float* %121
%123 = fmul float %118, %122
%124 = call float @llvm.fmuladd.f32(float %110, float %114, float %123)
%125 = call double @sqrt(float %124)
store double %125, double* %15
%126 = load i32, i32* %4
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %127
%129 = load float, float* %128
%130 = fneg float %129
%131 = load i32, i32* %4
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %132
%134 = load float, float* %133
%135 = call double @atan2(float %130, float %134)
store double %135, double* %16
%136 = load i32, i32* %4
%137 = load i32, i32* %4
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %138
%140 = load float, float* %139
%141 = fpext float %140 to double
%142 = load i32, i32* %4
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %143
%145 = load float, float* %144
%146 = fpext float %145 to double
%147 = load double, double* %15
%148 = load double, double* %16
%149 = load i32, i32* %1
%150 = sitofp i32 %149 to double
%151 = fmul double %148, %150
%152 = load i32, i32* @M_PI
%153 = sitofp i32 %152 to double
%154 = fmul double 2.000000e+00, %153
%155 = fdiv double %151, %154
%156 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.2, i64 0, i64 0), i32 %136, double %141, double %146, double %147, double %155)
br label %157
157: ; preds = %106
%158 = load i32, i32* %4
%159 = add nsw i32 %158, 1
store i32 %159, i32* %4
br label %101
160: ; preds = %101
ret void
}

LLVM IR (-O2):
@.str = [18 x i8] c"%10d%10.4f%10.4f\0A\00"
@.str.2 = [30 x i8] c"%10d%10.4f%10.4f%10.4f%10.4f\0A\00"
@M_PI = global i32 0
@str = [3 x i8] c"OK\00"
define void @transform_gaussian() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca i32
%4 = alloca [64 x float]
%5 = bitcast [64 x float]* %4 to i8*
%6 = alloca [64 x float]
%7 = bitcast [64 x float]* %6 to i8*
%8 = alloca [64 x float]
%9 = bitcast [64 x float]* %8 to i8*
%10 = bitcast i32* %1 to i8*
store i32 64, i32* %1
%11 = bitcast i32* %2 to i8*
store i32 1, i32* %2
%12 = bitcast i32* %3 to i8*
store i32 -1, i32* %3
br label %13
13: ; preds = %0, %13
%14 = phi i32 [ %26, %13 ], [ -32, %0 ]
%15 = lshr i32 %14, 25
%16 = and i32 %15, 64
%17 = add nsw i32 %16, %14
%18 = sitofp i32 %14 to double
%19 = fadd double %18, -5.000000e-01
%20 = fmul double %19, -5.000000e-01
%21 = fmul double %19, %20
%22 = fdiv double %21, 2.500000e+01
%23 = tail call float @exp(double %22)
%24 = sext i32 %17 to i64
%25 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %24
store float %23, float* %25
%26 = add nsw i32 %14, 1
%27 = icmp eq i32 %26, 32
br i1 %27, label %28, label %13
28: ; preds = %13
%29 = getelementptr inbounds [64 x float], [64 x float]* %4, i64 0, i64 0
%30 = load float, float* %29
%31 = fpext float %30 to double
%32 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @.str, i64 0, i64 0), i32 0, double %31, double 0.000000e+00)
br label %33
33: ; preds = %28, %33
%34 = phi i64 [ 1, %28 ], [ %43, %33 ]
%35 = getelementptr inbounds [64 x float], [64 x float]* %6, i64 0, i64 %34
%36 = load float, float* %35
%37 = getelementptr inbounds [64 x float], [64 x float]* %4, i64 0, i64 %34
%38 = load float, float* %37
%39 = fpext float %38 to double
%40 = fpext float %36 to double
%41 = trunc i64 %34 to i32
%42 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @.str, i64 0, i64 0), i32 %41, double %39, double %40)
%43 = add nuw nsw i64 %34, 1
%44 = icmp eq i64 %43, 64
br i1 %44, label %45, label %33
45: ; preds = %33
%46 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([3 x i8], [3 x i8]* @str, i64 0, i64 0))
%47 = getelementptr inbounds [64 x float], [64 x float]* %4, i64 0, i64 0
%48 = getelementptr inbounds [64 x float], [64 x float]* %6, i64 0, i64 0
%49 = call i32 @fft_(float* %47, float* %48, i32* %2, i32* %1, i32* %2, i32* %3)
%50 = load i32, i32* %1
%51 = icmp slt i32 %50, -1
br i1 %51, label %83, label %52
52: ; preds = %45, %52
%53 = phi i64 [ %78, %52 ], [ 0, %45 ]
%54 = getelementptr inbounds [64 x float], [64 x float]* %4, i64 0, i64 %53
%55 = load float, float* %54
%56 = getelementptr inbounds [64 x float], [64 x float]* %6, i64 0, i64 %53
%57 = load float, float* %56
%58 = fmul float %57, %57
%59 = call float @llvm.fmuladd.f32(float %55, float %55, float %58)
%60 = call double @sqrt(float %59)
%61 = load float, float* %56
%62 = fneg float %61
%63 = load float, float* %54
%64 = call double @atan2(float %62, float %63)
%65 = load float, float* %54
%66 = fpext float %65 to double
%67 = load float, float* %56
%68 = fpext float %67 to double
%69 = load i32, i32* %1
%70 = sitofp i32 %69 to double
%71 = fmul double %64, %70
%72 = load i32, i32* @M_PI
%73 = sitofp i32 %72 to double
%74 = fmul double %73, 2.000000e+00
%75 = fdiv double %71, %74
%76 = trunc i64 %53 to i32
%77 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.2, i64 0, i64 0), i32 %76, double %66, double %68, double %60, double %75)
%78 = add nuw nsw i64 %53, 1
%79 = load i32, i32* %1
%80 = sdiv i32 %79, 2
%81 = sext i32 %80 to i64
%82 = icmp slt i64 %53, %81
br i1 %82, label %52, label %83
83: ; preds = %52, %45
ret void
}