#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*/ exit(int /*<<< orphan*/);
int /*<<< orphan*/ fft_(float *, float *, int *, int *, int *, int *);
int /*<<< orphan*/ printf(char *, int, float, float, float, float, float,
                          float);
void precision_test() {
  extern double dnormal();
  int n = 64;
  int one = 1, negone = -1;
  int i, j, k;
  float a[64], b[64], c[64], d[64];
  for (i = 0; i < n; i++) {
    a[i] = c[i] = dnormal();
    b[i] = d[i] = dnormal();
  }
  fft_(a, b, &one, &n, &one, &negone);
  fft_(a, b, &one, &n, &one, &one);
  for (i = 0; i < n; i++) {
    printf("%10d%10.6f%10.6f%10.6f%10.6f%10.6f%10.6f\n", i, a[i], b[i], c[i],
           d[i], c[i] - a[i], d[i] - b[i]);
  }
  exit(0);
}


x86 assembly (clang, -O0):
<precision_test>:
push   %rbp
mov    %rsp,%rbp
sub    $0x420,%rsp
movl   $0x40,-0x4(%rbp)
movl   $0x1,-0x8(%rbp)
movl   $0xffffffff,-0xc(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11cb <precision_test+0x8b>
mov    $0x0,%al
call   1030 <dnormal@plt>
cvtsd2ss %xmm0,%xmm0
movslq -0x10(%rbp),%rax
movss  %xmm0,-0x320(%rbp,%rax,4)
ff ff
movslq -0x10(%rbp),%rax
movss  %xmm0,-0x120(%rbp,%rax,4)
ff ff
mov    $0x0,%al
call   1030 <dnormal@plt>
cvtsd2ss %xmm0,%xmm0
movslq -0x10(%rbp),%rax
movss  %xmm0,-0x420(%rbp,%rax,4)
ff ff
movslq -0x10(%rbp),%rax
movss  %xmm0,-0x220(%rbp,%rax,4)
ff ff
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1167 <precision_test+0x27>
lea    -0x120(%rbp),%rdi
lea    -0x220(%rbp),%rsi
lea    -0x8(%rbp),%r8
lea    -0x4(%rbp),%rcx
lea    -0xc(%rbp),%r9
mov    %r8,%rdx
call   1050 <fft_@plt>
lea    -0x120(%rbp),%rdi
lea    -0x220(%rbp),%rsi
lea    -0x8(%rbp),%r9
lea    -0x4(%rbp),%rcx
mov    %r9,%rdx
mov    %r9,%r8
call   1050 <fft_@plt>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12a6 <precision_test+0x166>
mov    -0x10(%rbp),%esi
movslq -0x10(%rbp),%rax
movss  -0x120(%rbp,%rax,4),%xmm0
ff ff
movslq -0x10(%rbp),%rax
movss  -0x220(%rbp,%rax,4),%xmm1
ff ff
movslq -0x10(%rbp),%rax
movss  -0x320(%rbp,%rax,4),%xmm2
ff ff
movslq -0x10(%rbp),%rax
movss  -0x420(%rbp,%rax,4),%xmm3
ff ff
movslq -0x10(%rbp),%rax
movss  -0x320(%rbp,%rax,4),%xmm4
ff ff
movslq -0x10(%rbp),%rax
subss  -0x120(%rbp,%rax,4),%xmm4
ff ff
movslq -0x10(%rbp),%rax
movss  -0x420(%rbp,%rax,4),%xmm5
ff ff
movslq -0x10(%rbp),%rax
subss  -0x220(%rbp,%rax,4),%xmm5
ff ff
lea    0xd6d(%rip),%rdi
call   1040 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1215 <precision_test+0xd5>
xor    %edi,%edi
call   1060 <exit@plt>

x86 assembly (clang, -O2):
<precision_test>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x418,%rsp
movl   $0x40,0x4(%rsp)
movl   $0x1,0xc(%rsp)
movl   $0xffffffff,0x8(%rsp)
ff
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
xor    %eax,%eax
call   1030 <dnormal@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x310(%rsp,%rbx,4)
movss  %xmm0,0x110(%rsp,%rbx,4)
xor    %eax,%eax
call   1030 <dnormal@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x210(%rsp,%rbx,4)
movss  %xmm0,0x10(%rsp,%rbx,4)
add    $0x1,%rbx
cmp    $0x40,%rbx
jne    1170 <precision_test+0x30>
lea    0x110(%rsp),%r12
lea    0x10(%rsp),%r14
lea    0x4(%rsp),%r15
lea    0xc(%rsp),%rbx
lea    0x8(%rsp),%r9
mov    %r12,%rdi
mov    %r14,%rsi
mov    %rbx,%rdx
mov    %r15,%rcx
mov    %rbx,%r8
call   1050 <fft_@plt>
mov    %r12,%rdi
mov    %r14,%rsi
mov    %rbx,%rdx
mov    %r15,%rcx
mov    %rbx,%r8
mov    %rbx,%r9
call   1050 <fft_@plt>
cmpl   $0x0,0x4(%rsp)
jle    1257 <precision_test+0x117>
xor    %ebx,%ebx
lea    0xdf8(%rip),%r14
nopl   0x0(%rax,%rax,1)
movss  0x110(%rsp,%rbx,4),%xmm0
movss  0x10(%rsp,%rbx,4),%xmm1
movss  0x310(%rsp,%rbx,4),%xmm2
movss  0x210(%rsp,%rbx,4),%xmm3
movaps %xmm2,%xmm4
subss  %xmm0,%xmm4
movaps %xmm3,%xmm5
subss  %xmm1,%xmm5
mov    %r14,%rdi
mov    %ebx,%esi
call   1040 <printf@plt>
add    $0x1,%rbx
movslq 0x4(%rsp),%rax
cmp    %rax,%rbx
jl     1210 <precision_test+0xd0>
xor    %edi,%edi
call   1060 <exit@plt>

LLVM IR (-O0):
@.str = [42 x i8] c"%10d%10.6f%10.6f%10.6f%10.6f%10.6f%10.6f\0A\00"
define void @precision_test() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca [64 x float]
%8 = alloca [64 x float]
%9 = alloca [64 x float]
%10 = alloca [64 x float]
store i32 64, i32* %1
store i32 1, i32* %2
store i32 -1, i32* %3
store i32 0, i32* %4
br label %11
11: ; preds = %32, %0
%12 = load i32, i32* %4
%13 = load i32, i32* %1
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %35
15: ; preds = %11
%16 = call double (...) @dnormal()
%17 = fptrunc double %16 to float
%18 = load i32, i32* %4
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %19
store float %17, float* %20
%21 = load i32, i32* %4
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds [64 x float], [64 x float]* %7, i64 0, i64 %22
store float %17, float* %23
%24 = call double (...) @dnormal()
%25 = fptrunc double %24 to float
%26 = load i32, i32* %4
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds [64 x float], [64 x float]* %10, i64 0, i64 %27
store float %25, float* %28
%29 = load i32, i32* %4
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %30
store float %25, float* %31
br label %32
32: ; preds = %15
%33 = load i32, i32* %4
%34 = add nsw i32 %33, 1
store i32 %34, i32* %4
br label %11
35: ; preds = %11
%36 = getelementptr inbounds [64 x float], [64 x float]* %7, i64 0, i64 0
%37 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 0
%38 = call i32 @fft_(float* %36, float* %37, i32* %2, i32* %1, i32* %2, i32* %3)
%39 = getelementptr inbounds [64 x float], [64 x float]* %7, i64 0, i64 0
%40 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 0
%41 = call i32 @fft_(float* %39, float* %40, i32* %2, i32* %1, i32* %2, i32* %2)
store i32 0, i32* %4
br label %42
42: ; preds = %83, %35
%43 = load i32, i32* %4
%44 = load i32, i32* %1
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %86
46: ; preds = %42
%47 = load i32, i32* %4
%48 = load i32, i32* %4
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds [64 x float], [64 x float]* %7, i64 0, i64 %49
%51 = load float, float* %50
%52 = load i32, i32* %4
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %53
%55 = load float, float* %54
%56 = load i32, i32* %4
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %57
%59 = load float, float* %58
%60 = load i32, i32* %4
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds [64 x float], [64 x float]* %10, i64 0, i64 %61
%63 = load float, float* %62
%64 = load i32, i32* %4
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds [64 x float], [64 x float]* %9, i64 0, i64 %65
%67 = load float, float* %66
%68 = load i32, i32* %4
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds [64 x float], [64 x float]* %7, i64 0, i64 %69
%71 = load float, float* %70
%72 = fsub float %67, %71
%73 = load i32, i32* %4
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds [64 x float], [64 x float]* %10, i64 0, i64 %74
%76 = load float, float* %75
%77 = load i32, i32* %4
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds [64 x float], [64 x float]* %8, i64 0, i64 %78
%80 = load float, float* %79
%81 = fsub float %76, %80
%82 = call i32 @printf(i8* getelementptr inbounds ([42 x i8], [42 x i8]* @.str, i64 0, i64 0), i32 %47, float %51, float %55, float %59, float %63, float %72, float %81)
br label %83
83: ; preds = %46
%84 = load i32, i32* %4
%85 = add nsw i32 %84, 1
store i32 %85, i32* %4
br label %42
86: ; preds = %42
%87 = call i32 @exit(i32 0)
unreachable
}

LLVM IR (-O2):
@.str = [42 x i8] c"%10d%10.6f%10.6f%10.6f%10.6f%10.6f%10.6f\0A\00"
define void @precision_test() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca i32
%4 = alloca [64 x float]
%5 = alloca [64 x float]
%6 = alloca [64 x float]
%7 = alloca [64 x float]
%8 = bitcast i32* %1 to i8*
store i32 64, i32* %1
%9 = bitcast i32* %2 to i8*
store i32 1, i32* %2
%10 = bitcast i32* %3 to i8*
store i32 -1, i32* %3
%11 = bitcast [64 x float]* %4 to i8*
%12 = bitcast [64 x float]* %5 to i8*
%13 = bitcast [64 x float]* %6 to i8*
%14 = bitcast [64 x float]* %7 to i8*
br label %15
15: ; preds = %0, %15
%16 = phi i64 [ 0, %0 ], [ %25, %15 ]
%17 = tail call double (...) @dnormal()
%18 = fptrunc double %17 to float
%19 = getelementptr inbounds [64 x float], [64 x float]* %6, i64 0, i64 %16
store float %18, float* %19
%20 = getelementptr inbounds [64 x float], [64 x float]* %4, i64 0, i64 %16
store float %18, float* %20
%21 = tail call double (...) @dnormal()
%22 = fptrunc double %21 to float
%23 = getelementptr inbounds [64 x float], [64 x float]* %7, i64 0, i64 %16
store float %22, float* %23
%24 = getelementptr inbounds [64 x float], [64 x float]* %5, i64 0, i64 %16
store float %22, float* %24
%25 = add nuw nsw i64 %16, 1
%26 = icmp eq i64 %25, 64
br i1 %26, label %27, label %15
27: ; preds = %15
%28 = getelementptr inbounds [64 x float], [64 x float]* %4, i64 0, i64 0
%29 = getelementptr inbounds [64 x float], [64 x float]* %5, i64 0, i64 0
%30 = call i32 @fft_(float* %28, float* %29, i32* %2, i32* %1, i32* %2, i32* %3)
%31 = call i32 @fft_(float* %28, float* %29, i32* %2, i32* %1, i32* %2, i32* %2)
%32 = load i32, i32* %1
%33 = icmp sgt i32 %32, 0
br i1 %33, label %34, label %52
34: ; preds = %27, %34
%35 = phi i64 [ %48, %34 ], [ 0, %27 ]
%36 = getelementptr inbounds [64 x float], [64 x float]* %4, i64 0, i64 %35
%37 = load float, float* %36
%38 = getelementptr inbounds [64 x float], [64 x float]* %5, i64 0, i64 %35
%39 = load float, float* %38
%40 = getelementptr inbounds [64 x float], [64 x float]* %6, i64 0, i64 %35
%41 = load float, float* %40
%42 = getelementptr inbounds [64 x float], [64 x float]* %7, i64 0, i64 %35
%43 = load float, float* %42
%44 = fsub float %41, %37
%45 = fsub float %43, %39
%46 = trunc i64 %35 to i32
%47 = call i32 @printf(i8* getelementptr inbounds ([42 x i8], [42 x i8]* @.str, i64 0, i64 0), i32 %46, float %37, float %39, float %41, float %43, float %44, float %45)
%48 = add nuw nsw i64 %35, 1
%49 = load i32, i32* %1
%50 = sext i32 %49 to i64
%51 = icmp slt i64 %48, %50
br i1 %51, label %34, label %52
52: ; preds = %34, %27
%53 = call i32 @exit(i32 0)
unreachable
}