#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 fabs(float);
int logf(int);
int /*<<< orphan*/ median(float *, int, float *);
float sqrt(int);
void denoise_donohoHard_DWTAnlz(float *ip, float *op, int n) {
  float med, trshld;
  for (int j = 1; j < n; j *= 2) {
    median(ip + j, j, &med);
    trshld = med / 0.6745 * sqrt(2 * logf(j));
    for (int i = j; i < 2 * j; i++) {
      if (fabs(ip[i]) > trshld)
        op[i] = ip[i];
      else
        op[i] = 0;
    }
  }
}


x86 assembly (clang, -O0):
<denoise_donohoHard_DWTAnlz>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    124e <denoise_donohoHard_DWTAnlz+0x10e>
mov    -0x8(%rbp),%rdi
movslq -0x20(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x20(%rbp),%esi
lea    -0x18(%rbp),%rdx
call   1030 <median@plt>
movss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe6e(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%edi
call   1040 <logf@plt>
mov    %eax,%edi
shl    $0x1,%edi
call   1060 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x30(%rbp),%xmm0
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x20(%rbp),%ecx
shl    $0x1,%ecx
cmp    %ecx,%eax
jge    123b <denoise_donohoHard_DWTAnlz+0xfb>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
call   1050 <fabs@plt>
ucomiss -0x1c(%rbp),%xmm0
jbe    1218 <denoise_donohoHard_DWTAnlz+0xd8>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1228 <denoise_donohoHard_DWTAnlz+0xe8>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    122d <denoise_donohoHard_DWTAnlz+0xed>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11cc <denoise_donohoHard_DWTAnlz+0x8c>
jmp    1240 <denoise_donohoHard_DWTAnlz+0x100>
mov    -0x20(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    115a <denoise_donohoHard_DWTAnlz+0x1a>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<denoise_donohoHard_DWTAnlz>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
cmp    $0x2,%edx
jl     1329 <denoise_donohoHard_DWTAnlz+0x1f9>
mov    %edx,%r14d
mov    %rsi,%r12
mov    %rdi,%rbp
lea    0x10(%rdi),%rax
mov    %rax,0x8(%rsp)
lea    0x10(%rsi),%r13
mov    $0x1,%r8d
lea    0x4(%rsp),%r15
jmp    1179 <denoise_donohoHard_DWTAnlz+0x49>
nopw   0x0(%rax,%rax,1)
cmp    %r14d,%r8d
jge    1329 <denoise_donohoHard_DWTAnlz+0x1f9>
mov    %r8d,%ebx
mov    %r8d,%eax
lea    0x0(,%rax,4),%rdi
add    %rbp,%rdi
mov    %r8d,%esi
mov    %r15,%rdx
call   1030 <median@plt>
movss  0x4(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
divsd  0xe59(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
mov    %ebx,%edi
call   1040 <logf@plt>
lea    (%rax,%rax,1),%edi
call   1050 <sqrt@plt>
movaps 0xe4d(%rip),%xmm7
lea    (%rbx,%rbx,1),%r8d
cmp    %r8d,%ebx
jge    1170 <denoise_donohoHard_DWTAnlz+0x40>
cvtss2sd %xmm0,%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movslq %ebx,%rax
movslq %r8d,%rcx
mov    %rcx,%r9
sub    %rax,%r9
cmp    $0x8,%r9
jb     1223 <denoise_donohoHard_DWTAnlz+0xf3>
lea    (%r12,%rax,4),%rdx
lea    0x0(,%rcx,4),%rsi
add    %rbp,%rsi
cmp    %rsi,%rdx
jae    12a9 <denoise_donohoHard_DWTAnlz+0x179>
lea    (%r12,%rcx,4),%rdx
lea    0x0(,%rax,4),%rsi
add    %rbp,%rsi
cmp    %rdx,%rsi
jae    12a9 <denoise_donohoHard_DWTAnlz+0x179>
mov    %rax,%rdx
mov    %rdx,%rsi
test   $0x1,%dl
je     124f <denoise_donohoHard_DWTAnlz+0x11f>
movss  0x0(%rbp,%rdx,4),%xmm1
movaps %xmm1,%xmm2
andps  %xmm7,%xmm2
movaps %xmm0,%xmm3
cmpltss %xmm2,%xmm3
andps  %xmm1,%xmm3
movss  %xmm3,(%r12,%rdx,4)
lea    0x1(%rdx),%rsi
not    %rdx
add    %rcx,%rdx
je     1170 <denoise_donohoHard_DWTAnlz+0x40>
nopl   0x0(%rax,%rax,1)
movss  0x0(%rbp,%rsi,4),%xmm1
movaps %xmm1,%xmm2
andps  %xmm7,%xmm2
movaps %xmm0,%xmm3
cmpltss %xmm2,%xmm3
andps  %xmm1,%xmm3
movss  %xmm3,(%r12,%rsi,4)
movss  0x4(%rbp,%rsi,4),%xmm1
movaps %xmm1,%xmm2
andps  %xmm7,%xmm2
movaps %xmm0,%xmm3
cmpltss %xmm2,%xmm3
andps  %xmm1,%xmm3
movss  %xmm3,0x4(%r12,%rsi,4)
add    $0x2,%rsi
cmp    %rsi,%rcx
jne    1260 <denoise_donohoHard_DWTAnlz+0x130>
jmp    1170 <denoise_donohoHard_DWTAnlz+0x40>
mov    %r9,%rdi
and    $0xfffffffffffffff8,%rdi
lea    (%rdi,%rax,1),%rdx
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
mov    0x8(%rsp),%rsi
lea    (%rsi,%rax,4),%rbx
lea    0x0(,%rax,4),%rax
add    %r13,%rax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movups -0x10(%rbx,%rsi,4),%xmm2
movups (%rbx,%rsi,4),%xmm3
movaps %xmm2,%xmm4
andps  %xmm7,%xmm4
movaps %xmm3,%xmm5
andps  %xmm7,%xmm5
movaps %xmm1,%xmm6
cmpltps %xmm4,%xmm6
movaps %xmm1,%xmm4
cmpltps %xmm5,%xmm4
andps  %xmm2,%xmm6
andps  %xmm3,%xmm4
movups %xmm6,-0x10(%rax,%rsi,4)
movups %xmm4,(%rax,%rsi,4)
add    $0x8,%rsi
cmp    %rsi,%rdi
jne    12e0 <denoise_donohoHard_DWTAnlz+0x1b0>
cmp    %rdi,%r9
je     1170 <denoise_donohoHard_DWTAnlz+0x40>
jmp    1226 <denoise_donohoHard_DWTAnlz+0xf6>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @denoise_donohoHard_DWTAnlz(float* %0, float* nocapture writeonly %1, i32 %2) {
%4 = alloca float
%5 = bitcast float* %4 to i8*
%6 = icmp sgt i32 %2, 1
br i1 %6, label %10, label %9
7: ; preds = %82, %86, %67, %10
%8 = icmp slt i32 %24, %2
br i1 %8, label %10, label %9
9: ; preds = %7, %3
ret void
10: ; preds = %3, %7
%11 = phi i32 [ %24, %7 ], [ 1, %3 ]
%12 = zext i32 %11 to i64
%13 = getelementptr inbounds float, float* %0, i64 %12
%14 = call i32 @median(float* %13, i32 %11, float* %4)
%15 = load float, float* %4
%16 = fpext float %15 to double
%17 = fdiv double %16, 6.745000e-01
%18 = call i32 @logf(i32 %11)
%19 = shl nsw i32 %18, 1
%20 = call float @sqrt(i32 %19)
%21 = fpext float %20 to double
%22 = fmul double %17, %21
%23 = fptrunc double %22 to float
%24 = shl i32 %11, 1
%25 = icmp slt i32 %11, %24
br i1 %25, label %26, label %7
26: ; preds = %10
%27 = sext i32 %11 to i64
%28 = sext i32 %24 to i64
%29 = sub nsw i64 %28, %27
%30 = icmp ult i64 %29, 8
br i1 %30, label %69, label %31
31: ; preds = %26
%32 = getelementptr float, float* %1, i64 %27
%33 = getelementptr float, float* %1, i64 %28
%34 = getelementptr float, float* %0, i64 %27
%35 = getelementptr float, float* %0, i64 %28
%36 = icmp ult float* %32, %35
%37 = icmp ult float* %34, %33
%38 = and i1 %36, %37
br i1 %38, label %69, label %39
39: ; preds = %31
%40 = and i64 %29, -8
%41 = add nsw i64 %40, %27
%42 = insertelement <4 x float> poison, float %23, i64 0
%43 = shufflevector <4 x float> %42, <4 x float> poison, <4 x i32> zeroinitializer
%44 = insertelement <4 x float> poison, float %23, i64 0
%45 = shufflevector <4 x float> %44, <4 x float> poison, <4 x i32> zeroinitializer
br label %46
46: ; preds = %46, %39
%47 = phi i64 [ 0, %39 ], [ %65, %46 ]
%48 = add i64 %47, %27
%49 = getelementptr inbounds float, float* %0, i64 %48
%50 = bitcast float* %49 to <4 x float>*
%51 = load <4 x float>, <4 x float>* %50
%52 = getelementptr inbounds float, float* %49, i64 4
%53 = bitcast float* %52 to <4 x float>*
%54 = load <4 x float>, <4 x float>* %53
%55 = call <4 x float> @llvm.fabs.v4f32(<4 x float> %51)
%56 = call <4 x float> @llvm.fabs.v4f32(<4 x float> %54)
%57 = fcmp ogt <4 x float> %55, %43
%58 = fcmp ogt <4 x float> %56, %45
%59 = select <4 x i1> %57, <4 x float> %51, <4 x float> zeroinitializer
%60 = select <4 x i1> %58, <4 x float> %54, <4 x float> zeroinitializer
%61 = getelementptr inbounds float, float* %1, i64 %48
%62 = bitcast float* %61 to <4 x float>*
store <4 x float> %59, <4 x float>* %62
%63 = getelementptr inbounds float, float* %61, i64 4
%64 = bitcast float* %63 to <4 x float>*
store <4 x float> %60, <4 x float>* %64
%65 = add nuw i64 %47, 8
%66 = icmp eq i64 %65, %40
br i1 %66, label %67, label %46
67: ; preds = %46
%68 = icmp eq i64 %29, %40
br i1 %68, label %7, label %69
69: ; preds = %31, %26, %67
%70 = phi i64 [ %27, %31 ], [ %27, %26 ], [ %41, %67 ]
%71 = xor i64 %70, -1
%72 = and i64 %70, 1
%73 = icmp eq i64 %72, 0
br i1 %73, label %82, label %74
74: ; preds = %69
%75 = getelementptr inbounds float, float* %0, i64 %70
%76 = load float, float* %75
%77 = call float @llvm.fabs.f32(float %76)
%78 = fcmp ogt float %77, %23
%79 = select i1 %78, float %76, float 0.000000e+00
%80 = getelementptr inbounds float, float* %1, i64 %70
store float %79, float* %80
%81 = add nsw i64 %70, 1
br label %82
82: ; preds = %74, %69
%83 = phi i64 [ %70, %69 ], [ %81, %74 ]
%84 = sub nsw i64 0, %28
%85 = icmp eq i64 %71, %84
br i1 %85, label %7, label %86
86: ; preds = %82, %86
%87 = phi i64 [ %101, %86 ], [ %83, %82 ]
%88 = getelementptr inbounds float, float* %0, i64 %87
%89 = load float, float* %88
%90 = call float @llvm.fabs.f32(float %89)
%91 = fcmp ogt float %90, %23
%92 = select i1 %91, float %89, float 0.000000e+00
%93 = getelementptr inbounds float, float* %1, i64 %87
store float %92, float* %93
%94 = add nsw i64 %87, 1
%95 = getelementptr inbounds float, float* %0, i64 %94
%96 = load float, float* %95
%97 = call float @llvm.fabs.f32(float %96)
%98 = fcmp ogt float %97, %23
%99 = select i1 %98, float %96, float 0.000000e+00
%100 = getelementptr inbounds float, float* %1, i64 %94
store float %99, float* %100
%101 = add nsw i64 %87, 2
%102 = icmp eq i64 %101, %28
br i1 %102, label %7, label %86
}