#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 */
float FLT_MIN;
void wiener_subtraction(int fft_size_2, float *spectrum,
                        float *noise_thresholds, float *Gk) {
  int k;
  for (k = 0; k <= fft_size_2; k++) {
    if (noise_thresholds[k] > FLT_MIN) {
      if (spectrum[k] > noise_thresholds[k]) {
        Gk[k] = (spectrum[k] - noise_thresholds[k]) / spectrum[k];
      } else {
        Gk[k] = 0.f;
      }
    } else {
      Gk[k] = 1.f;
    }
  }
  for (k = 1; k < fft_size_2; k++) {
    Gk[(2 * fft_size_2) - k] = Gk[k];
  }
}


x86 assembly (clang, -O0):
<wiener_subtraction>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     11d8 <wiener_subtraction+0xd8>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x2eae(%rip),%rax
ucomiss (%rax),%xmm0
jbe    11b0 <wiener_subtraction+0xb0>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
jbe    119b <wiener_subtraction+0x9b>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    11ab <wiener_subtraction+0xab>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11c5 <wiener_subtraction+0xc5>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  0xe40(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11ca <wiener_subtraction+0xca>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111a <wiener_subtraction+0x1a>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    121b <wiener_subtraction+0x11b>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x4(%rbp),%ecx
shl    $0x1,%ecx
sub    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11df <wiener_subtraction+0xdf>
pop    %rbp
ret

x86 assembly (clang, -O2):
<wiener_subtraction>:
test   %edi,%edi
js     1216 <wiener_subtraction+0x116>
lea    0x1(%rdi),%r8d
xor    %eax,%eax
mov    0x2ed3(%rip),%r9
movss  0xee3(%rip),%xmm0
jmp    112e <wiener_subtraction+0x2e>
nop
movss  %xmm1,(%rcx,%rax,4)
add    $0x1,%rax
cmp    %rax,%r8
je     1156 <wiener_subtraction+0x56>
movss  (%rdx,%rax,4),%xmm2
ucomiss (%r9),%xmm2
movaps %xmm0,%xmm1
jbe    1120 <wiener_subtraction+0x20>
movss  (%rsi,%rax,4),%xmm3
xorps  %xmm1,%xmm1
ucomiss %xmm2,%xmm3
jbe    1120 <wiener_subtraction+0x20>
movaps %xmm3,%xmm1
subss  %xmm2,%xmm1
divss  %xmm3,%xmm1
jmp    1120 <wiener_subtraction+0x20>
cmp    $0x2,%edi
jl     1216 <wiener_subtraction+0x116>
lea    (%rdi,%rdi,1),%eax
movslq %eax,%r8
mov    %edi,%eax
lea    -0x1(%rax),%r10
add    $0xfffffffffffffffe,%rax
mov    %r10d,%r9d
and    $0x3,%r9d
mov    $0x1,%esi
cmp    $0x3,%rax
jb     11e3 <wiener_subtraction+0xe3>
and    $0xfffffffffffffffc,%r10
lea    (%rcx,%r8,4),%rdx
add    $0xfffffffffffffffc,%rdx
lea    0x10(%rcx),%rax
mov    $0xffffffffffffffff,%rsi
nopl   0x0(%rax,%rax,1)
movss  -0xc(%rax),%xmm0
movss  %xmm0,0x4(%rdx,%rsi,4)
movss  -0x8(%rax),%xmm0
movss  %xmm0,(%rdx,%rsi,4)
movss  -0x4(%rax),%xmm0
movss  %xmm0,-0x4(%rdx,%rsi,4)
movss  (%rax),%xmm0
movss  %xmm0,-0x8(%rdx,%rsi,4)
add    $0x10,%rax
lea    (%r10,%rsi,1),%rdi
add    $0xfffffffffffffffc,%rdi
add    $0xfffffffffffffffc,%rsi
cmp    $0xffffffffffffffff,%rdi
jne    11a0 <wiener_subtraction+0xa0>
neg    %rsi
test   %r9,%r9
je     1216 <wiener_subtraction+0x116>
sub    %rsi,%r8
lea    (%rcx,%r8,4),%rax
lea    (%rcx,%rsi,4),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rcx,%rdx,4),%xmm0
movss  %xmm0,(%rax)
add    $0x1,%rdx
add    $0xfffffffffffffffc,%rax
cmp    %rdx,%r9
jne    1200 <wiener_subtraction+0x100>
ret

LLVM IR (-O0):
@FLT_MIN = global float 0.000000e+00
define void @wiener_subtraction(i32 %0, float* %1, float* %2, float* %3) {
%5 = alloca i32
%6 = alloca float*
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32
store i32 %0, i32* %5
store float* %1, float** %6
store float* %2, float** %7
store float* %3, float** %8
store i32 0, i32* %9
br label %10
10: ; preds = %68, %4
%11 = load i32, i32* %9
%12 = load i32, i32* %5
%13 = icmp sle i32 %11, %12
br i1 %13, label %14, label %71
14: ; preds = %10
%15 = load float*, float** %7
%16 = load i32, i32* %9
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds float, float* %15, i64 %17
%19 = load float, float* %18
%20 = load float, float* @FLT_MIN
%21 = fcmp ogt float %19, %20
br i1 %21, label %22, label %62
22: ; preds = %14
%23 = load float*, float** %6
%24 = load i32, i32* %9
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds float, float* %23, i64 %25
%27 = load float, float* %26
%28 = load float*, float** %7
%29 = load i32, i32* %9
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %28, i64 %30
%32 = load float, float* %31
%33 = fcmp ogt float %27, %32
br i1 %33, label %34, label %56
34: ; preds = %22
%35 = load float*, float** %6
%36 = load i32, i32* %9
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds float, float* %35, i64 %37
%39 = load float, float* %38
%40 = load float*, float** %7
%41 = load i32, i32* %9
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds float, float* %40, i64 %42
%44 = load float, float* %43
%45 = fsub float %39, %44
%46 = load float*, float** %6
%47 = load i32, i32* %9
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %46, i64 %48
%50 = load float, float* %49
%51 = fdiv float %45, %50
%52 = load float*, float** %8
%53 = load i32, i32* %9
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float, float* %52, i64 %54
store float %51, float* %55
br label %61
56: ; preds = %22
%57 = load float*, float** %8
%58 = load i32, i32* %9
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %57, i64 %59
store float 0.000000e+00, float* %60
br label %61
61: ; preds = %56, %34
br label %67
62: ; preds = %14
%63 = load float*, float** %8
%64 = load i32, i32* %9
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds float, float* %63, i64 %65
store float 1.000000e+00, float* %66
br label %67
67: ; preds = %62, %61
br label %68
68: ; preds = %67
%69 = load i32, i32* %9
%70 = add nsw i32 %69, 1
store i32 %70, i32* %9
br label %10
71: ; preds = %10
store i32 1, i32* %9
br label %72
72: ; preds = %89, %71
%73 = load i32, i32* %9
%74 = load i32, i32* %5
%75 = icmp slt i32 %73, %74
br i1 %75, label %76, label %92
76: ; preds = %72
%77 = load float*, float** %8
%78 = load i32, i32* %9
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds float, float* %77, i64 %79
%81 = load float, float* %80
%82 = load float*, float** %8
%83 = load i32, i32* %5
%84 = mul nsw i32 2, %83
%85 = load i32, i32* %9
%86 = sub nsw i32 %84, %85
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds float, float* %82, i64 %87
store float %81, float* %88
br label %89
89: ; preds = %76
%90 = load i32, i32* %9
%91 = add nsw i32 %90, 1
store i32 %91, i32* %9
br label %72
92: ; preds = %72
ret void
}

LLVM IR (-O2):
@FLT_MIN = global float 0.000000e+00
define void @wiener_subtraction(i32 %0, float* nocapture %1, float* nocapture %2, float* nocapture %3) {
%5 = icmp slt i32 %0, 0
br i1 %5, label %77, label %6
6: ; preds = %4
%7 = add nuw i32 %0, 1
%8 = zext i32 %7 to i64
br label %21
9: ; preds = %34
%10 = icmp sgt i32 %0, 1
br i1 %10, label %11, label %77
11: ; preds = %9
%12 = shl nuw nsw i32 %0, 1
%13 = sext i32 %12 to i64
%14 = zext i32 %0 to i64
%15 = add nsw i64 %14, -1
%16 = add nsw i64 %14, -2
%17 = and i64 %15, 3
%18 = icmp ult i64 %16, 3
br i1 %18, label %64, label %19
19: ; preds = %11
%20 = and i64 %15, -4
br label %39
21: ; preds = %6, %34
%22 = phi i64 [ 0, %6 ], [ %37, %34 ]
%23 = getelementptr inbounds float, float* %2, i64 %22
%24 = load float, float* %23
%25 = load float, float* @FLT_MIN
%26 = fcmp ogt float %24, %25
br i1 %26, label %27, label %34
27: ; preds = %21
%28 = getelementptr inbounds float, float* %1, i64 %22
%29 = load float, float* %28
%30 = fcmp ogt float %29, %24
br i1 %30, label %31, label %34
31: ; preds = %27
%32 = fsub float %29, %24
%33 = fdiv float %32, %29
br label %34
34: ; preds = %21, %27, %31
%35 = phi float [ %33, %31 ], [ 0.000000e+00, %27 ], [ 1.000000e+00, %21 ]
%36 = getelementptr inbounds float, float* %3, i64 %22
store float %35, float* %36
%37 = add nuw nsw i64 %22, 1
%38 = icmp eq i64 %37, %8
br i1 %38, label %9, label %21
39: ; preds = %39, %19
%40 = phi i64 [ 1, %19 ], [ %61, %39 ]
%41 = phi i64 [ 0, %19 ], [ %62, %39 ]
%42 = getelementptr inbounds float, float* %3, i64 %40
%43 = load float, float* %42
%44 = sub nsw i64 %13, %40
%45 = getelementptr inbounds float, float* %3, i64 %44
store float %43, float* %45
%46 = add nuw nsw i64 %40, 1
%47 = getelementptr inbounds float, float* %3, i64 %46
%48 = load float, float* %47
%49 = sub nsw i64 %13, %46
%50 = getelementptr inbounds float, float* %3, i64 %49
store float %48, float* %50
%51 = add nuw nsw i64 %40, 2
%52 = getelementptr inbounds float, float* %3, i64 %51
%53 = load float, float* %52
%54 = sub nsw i64 %13, %51
%55 = getelementptr inbounds float, float* %3, i64 %54
store float %53, float* %55
%56 = add nuw nsw i64 %40, 3
%57 = getelementptr inbounds float, float* %3, i64 %56
%58 = load float, float* %57
%59 = sub nsw i64 %13, %56
%60 = getelementptr inbounds float, float* %3, i64 %59
store float %58, float* %60
%61 = add nuw nsw i64 %40, 4
%62 = add nuw nsw i64 %41, 4
%63 = icmp eq i64 %62, %20
br i1 %63, label %64, label %39
64: ; preds = %39, %11
%65 = phi i64 [ 1, %11 ], [ %61, %39 ]
%66 = icmp eq i64 %17, 0
br i1 %66, label %77, label %67
67: ; preds = %64, %67
%68 = phi i64 [ %74, %67 ], [ %65, %64 ]
%69 = phi i64 [ %75, %67 ], [ 0, %64 ]
%70 = getelementptr inbounds float, float* %3, i64 %68
%71 = load float, float* %70
%72 = sub nsw i64 %13, %68
%73 = getelementptr inbounds float, float* %3, i64 %72
store float %71, float* %73
%74 = add nuw nsw i64 %68, 1
%75 = add i64 %69, 1
%76 = icmp eq i64 %75, %17
br i1 %76, label %77, label %67
77: ; preds = %64, %67, %4, %9
ret void
}