#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 */
double floatMAD(int, float *);
float floatmedian(int, float *);
float Tukey(int ndat, float *dat) {
  int niter = 10;
  int n, i;
  float s;
  float t, u, c;
  double sdenom, snumer;
  s = floatMAD(ndat, dat) / 0.6745;
  t = floatmedian(ndat, dat);
  c = 6.0;
  for (n = 0; n < niter; n++) {
    snumer = 0;
    sdenom = 0;
    for (i = 0; i < ndat; i++) {
      u = (dat[i] - t) / (c * s);
      if (u < 1 && u > -1) {
        snumer += dat[i] * (1. - u * u) * (1. - u * u);
        sdenom += (1. - u * u) * (1. - u * u);
      }
    }
    if (sdenom != 0) {
      t = snumer / sdenom;
    } else {
      c *= 1.1;
    }
  }
  return t;
}


x86 assembly (clang, -O0):
<Tukey>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0xa,-0x14(%rbp)
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
call   1040 <floatMAD@plt>
movsd  0xec6(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
call   1030 <floatmedian@plt>
movss  %xmm0,-0x24(%rbp)
movss  0xe90(%rip),%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1312 <Tukey+0x1f2>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12ac <Tukey+0x18c>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x24(%rbp),%xmm0
movss  -0x2c(%rbp),%xmm1
mulss  -0x20(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  0xe2c(%rip),%xmm0
ucomiss -0x28(%rbp),%xmm0
jbe    1299 <Tukey+0x179>
movss  -0x28(%rbp),%xmm0
movss  0xe19(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    1299 <Tukey+0x179>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x28(%rbp),%xmm1
mulss  -0x28(%rbp),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xe01(%rip),%xmm1
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movss  -0x28(%rbp),%xmm1
mulss  -0x28(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xde3(%rip),%xmm2
subsd  %xmm1,%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movss  -0x28(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdb7(%rip),%xmm0
subsd  %xmm1,%xmm0
movss  -0x28(%rbp),%xmm1
mulss  -0x28(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xd9d(%rip),%xmm2
subsd  %xmm1,%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    129e <Tukey+0x17e>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    119f <Tukey+0x7f>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12c9 <Tukey+0x1a9>
jp     12c9 <Tukey+0x1a9>
jmp    12e1 <Tukey+0x1c1>
movsd  -0x40(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
jmp    12ff <Tukey+0x1df>
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xd26(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
jmp    1304 <Tukey+0x1e4>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    117c <Tukey+0x5c>
movss  -0x24(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Tukey>:
push   %rbp
push   %rbx
push   %rax
mov    %rsi,%rbx
mov    %edi,%ebp
call   1040 <floatMAD@plt>
divsd  0xecb(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x4(%rsp)
mov    %ebp,%edi
mov    %rbx,%rsi
call   1030 <floatmedian@plt>
movss  0x4(%rsp),%xmm12
mov    %ebp,%eax
xor    %ecx,%ecx
movss  0xebb(%rip),%xmm10
movss  0xeb7(%rip),%xmm1
movss  0xeb3(%rip),%xmm2
movsd  0xe92(%rip),%xmm11
xorpd  %xmm9,%xmm9
movsd  0xe8c(%rip),%xmm8
jmp    11b2 <Tukey+0x92>
cs nopw 0x0(%rax,%rax,1)
xorps  %xmm3,%xmm3
cvtss2sd %xmm10,%xmm3
mulsd  %xmm8,%xmm3
xorps  %xmm10,%xmm10
cvtsd2ss %xmm3,%xmm10
add    $0x1,%ecx
cmp    $0xa,%ecx
je     1252 <Tukey+0x132>
test   %ebp,%ebp
jle    1190 <Tukey+0x70>
movaps %xmm10,%xmm4
mulss  %xmm12,%xmm4
xorpd  %xmm7,%xmm7
xor    %edx,%edx
jmp    11d9 <Tukey+0xb9>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rdx
cmp    %rdx,%rax
je     1230 <Tukey+0x110>
movss  (%rbx,%rdx,4),%xmm6
movaps %xmm6,%xmm5
subss  %xmm0,%xmm5
divss  %xmm4,%xmm5
ucomiss %xmm5,%xmm1
jbe    11d0 <Tukey+0xb0>
ucomiss %xmm2,%xmm5
jbe    11d0 <Tukey+0xb0>
cvtss2sd %xmm6,%xmm6
mulss  %xmm5,%xmm5
cvtss2sd %xmm5,%xmm5
movapd %xmm11,%xmm3
subsd  %xmm5,%xmm3
mulsd  %xmm3,%xmm6
movapd %xmm3,%xmm5
unpcklpd %xmm6,%xmm5
unpcklpd %xmm3,%xmm3
mulpd  %xmm5,%xmm3
addpd  %xmm3,%xmm7
jmp    11d0 <Tukey+0xb0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
ucomisd %xmm9,%xmm7
jne    123d <Tukey+0x11d>
jnp    1190 <Tukey+0x70>
movapd %xmm7,%xmm0
unpckhpd %xmm7,%xmm0
divsd  %xmm7,%xmm0
cvtsd2ss %xmm0,%xmm0
jmp    11a6 <Tukey+0x86>
add    $0x8,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define float @Tukey(i32 %0, float* %1) {
%3 = alloca i32
%4 = alloca float*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca float
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca double
%13 = alloca double
store i32 %0, i32* %3
store float* %1, float** %4
store i32 10, i32* %5
%14 = load i32, i32* %3
%15 = load float*, float** %4
%16 = call double @floatMAD(i32 %14, float* %15)
%17 = fdiv double %16, 6.745000e-01
%18 = fptrunc double %17 to float
store float %18, float* %8
%19 = load i32, i32* %3
%20 = load float*, float** %4
%21 = call float @floatmedian(i32 %19, float* %20)
store float %21, float* %9
store float 6.000000e+00, float* %11
store i32 0, i32* %6
br label %22
22: ; preds = %98, %2
%23 = load i32, i32* %6
%24 = load i32, i32* %5
%25 = icmp slt i32 %23, %24
br i1 %25, label %26, label %101
26: ; preds = %22
store double 0.000000e+00, double* %13
store double 0.000000e+00, double* %12
store i32 0, i32* %7
br label %27
27: ; preds = %81, %26
%28 = load i32, i32* %7
%29 = load i32, i32* %3
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %84
31: ; preds = %27
%32 = load float*, float** %4
%33 = load i32, i32* %7
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds float, float* %32, i64 %34
%36 = load float, float* %35
%37 = load float, float* %9
%38 = fsub float %36, %37
%39 = load float, float* %11
%40 = load float, float* %8
%41 = fmul float %39, %40
%42 = fdiv float %38, %41
store float %42, float* %10
%43 = load float, float* %10
%44 = fcmp olt float %43, 1.000000e+00
br i1 %44, label %45, label %80
45: ; preds = %31
%46 = load float, float* %10
%47 = fcmp ogt float %46, -1.000000e+00
br i1 %47, label %48, label %80
48: ; preds = %45
%49 = load float*, float** %4
%50 = load i32, i32* %7
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %49, i64 %51
%53 = load float, float* %52
%54 = fpext float %53 to double
%55 = load float, float* %10
%56 = load float, float* %10
%57 = fmul float %55, %56
%58 = fpext float %57 to double
%59 = fsub double 1.000000e+00, %58
%60 = fmul double %54, %59
%61 = load float, float* %10
%62 = load float, float* %10
%63 = fmul float %61, %62
%64 = fpext float %63 to double
%65 = fsub double 1.000000e+00, %64
%66 = load double, double* %13
%67 = call double @llvm.fmuladd.f64(double %60, double %65, double %66)
store double %67, double* %13
%68 = load float, float* %10
%69 = load float, float* %10
%70 = fmul float %68, %69
%71 = fpext float %70 to double
%72 = fsub double 1.000000e+00, %71
%73 = load float, float* %10
%74 = load float, float* %10
%75 = fmul float %73, %74
%76 = fpext float %75 to double
%77 = fsub double 1.000000e+00, %76
%78 = load double, double* %12
%79 = call double @llvm.fmuladd.f64(double %72, double %77, double %78)
store double %79, double* %12
br label %80
80: ; preds = %48, %45, %31
br label %81
81: ; preds = %80
%82 = load i32, i32* %7
%83 = add nsw i32 %82, 1
store i32 %83, i32* %7
br label %27
84: ; preds = %27
%85 = load double, double* %12
%86 = fcmp une double %85, 0.000000e+00
br i1 %86, label %87, label %92
87: ; preds = %84
%88 = load double, double* %13
%89 = load double, double* %12
%90 = fdiv double %88, %89
%91 = fptrunc double %90 to float
store float %91, float* %9
br label %97
92: ; preds = %84
%93 = load float, float* %11
%94 = fpext float %93 to double
%95 = fmul double %94, 1.100000e+00
%96 = fptrunc double %95 to float
store float %96, float* %11
br label %97
97: ; preds = %92, %87
br label %98
98: ; preds = %97
%99 = load i32, i32* %6
%100 = add nsw i32 %99, 1
store i32 %100, i32* %6
br label %22
101: ; preds = %22
%102 = load float, float* %9
ret float %102
}

LLVM IR (-O2):
define float @Tukey(i32 %0, float* %1) {
%3 = tail call double @floatMAD(i32 %0, float* %1)
%4 = fdiv double %3, 6.745000e-01
%5 = fptrunc double %4 to float
%6 = tail call float @floatmedian(i32 %0, float* %1)
%7 = icmp sgt i32 %0, 0
%8 = zext i32 %0 to i64
br label %9
9: ; preds = %2, %50
%10 = phi float [ 6.000000e+00, %2 ], [ %52, %50 ]
%11 = phi i32 [ 0, %2 ], [ %53, %50 ]
%12 = phi float [ %6, %2 ], [ %51, %50 ]
br i1 %7, label %13, label %46
13: ; preds = %9
%14 = fmul float %10, %5
br label %15
15: ; preds = %13, %35
%16 = phi i64 [ 0, %13 ], [ %37, %35 ]
%17 = phi <2 x double> [ zeroinitializer, %13 ], [ %36, %35 ]
%18 = getelementptr inbounds float, float* %1, i64 %16
%19 = load float, float* %18
%20 = fsub float %19, %12
%21 = fdiv float %20, %14
%22 = fcmp olt float %21, 1.000000e+00
%23 = fcmp ogt float %21, -1.000000e+00
%24 = and i1 %22, %23
br i1 %24, label %25, label %35
25: ; preds = %15
%26 = fpext float %19 to double
%27 = fmul float %21, %21
%28 = fpext float %27 to double
%29 = fsub double 1.000000e+00, %28
%30 = fmul double %29, %26
%31 = insertelement <2 x double> poison, double %29, i64 0
%32 = insertelement <2 x double> %31, double %30, i64 1
%33 = shufflevector <2 x double> %31, <2 x double> poison, <2 x i32> zeroinitializer
%34 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %32, <2 x double> %33, <2 x double> %17)
br label %35
35: ; preds = %15, %25
%36 = phi <2 x double> [ %34, %25 ], [ %17, %15 ]
%37 = add nuw nsw i64 %16, 1
%38 = icmp eq i64 %37, %8
br i1 %38, label %39, label %15
39: ; preds = %35
%40 = extractelement <2 x double> %36, i64 0
%41 = fcmp une double %40, 0.000000e+00
br i1 %41, label %42, label %46
42: ; preds = %39
%43 = extractelement <2 x double> %36, i64 1
%44 = fdiv double %43, %40
%45 = fptrunc double %44 to float
br label %50
46: ; preds = %9, %39
%47 = fpext float %10 to double
%48 = fmul double %47, 1.100000e+00
%49 = fptrunc double %48 to float
br label %50
50: ; preds = %42, %46
%51 = phi float [ %45, %42 ], [ %12, %46 ]
%52 = phi float [ %10, %42 ], [ %49, %46 ]
%53 = add nuw nsw i32 %11, 1
%54 = icmp eq i32 %53, 10
br i1 %54, label %55, label %9
55: ; preds = %50
ret float %51
}