#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*/ memcpy(float *, float *, int);
float spectral_median(float *x, int n) {
  float temp;
  int i, j;
  float tmp[n + 1];
  memcpy(tmp, x, sizeof(float) * (n + 1));
  for (i = 0; i < n; i++) {
    for (j = i + 1; j <= n; j++) {
      if (tmp[j] < tmp[i]) {
        temp = tmp[i];
        tmp[i] = tmp[j];
        tmp[j] = temp;
      }
    }
  }
  if (n % 2 == 0) {
    return ((tmp[n / 2] + tmp[n / 2 - 1]) / 2.f);
  } else {
    return tmp[n / 2];
  }
}


x86 assembly (clang, -O0):
<spectral_median>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %rsp,%rcx
mov    %rcx,-0x28(%rbp)
lea    0xf(,%rax,4),%rcx
and    $0xfffffffffffffff0,%rcx
mov    %rsp,%rdi
sub    %rcx,%rdi
mov    %rdi,-0x40(%rbp)
mov    %rdi,%rsp
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rsi
mov    -0x14(%rbp),%eax
add    $0x1,%eax
cltq
shl    $0x2,%rax
mov    %eax,%edx
call   1030 <memcpy@plt>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11ff <spectral_median+0xef>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     11ec <spectral_median+0xdc>
mov    -0x40(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss %xmm1,%xmm0
jbe    11d9 <spectral_median+0xc9>
mov    -0x40(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x18(%rbp)
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x18(%rbp),%xmm0
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    11de <spectral_median+0xce>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    117c <spectral_median+0x6c>
jmp    11f1 <spectral_median+0xe1>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1167 <spectral_median+0x57>
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    1265 <spectral_median+0x155>
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  0xdb0(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
movl   $0x1,-0x34(%rbp)
jmp    128a <spectral_median+0x17a>
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x40(%rbp),%rax
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x4(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,%rsp
movss  -0x4(%rbp),%xmm0
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<spectral_median>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %esi,%r14d
mov    %rdi,%rsi
lea    0x1(%r14),%r15d
mov    %rsp,%rbx
lea    0xf(,%r15,4),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%rbx
mov    %rbx,%rsp
lea    0x4(,%r14,4),%edx
mov    %rbx,%rdi
call   1030 <memcpy@plt>
test   %r14d,%r14d
jle    11a4 <spectral_median+0x94>
mov    %r14d,%eax
xor    %edx,%edx
mov    %rbx,%rcx
jmp    116d <spectral_median+0x5d>
nopl   0x0(%rax,%rax,1)
add    $0x4,%rcx
add    $0xffffffffffffffff,%r15
cmp    %rax,%rdx
je     11a4 <spectral_median+0x94>
mov    %rdx,%rsi
add    $0x1,%rdx
mov    $0x1,%edi
jmp    1189 <spectral_median+0x79>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    %rdi,%r15
je     1160 <spectral_median+0x50>
movss  (%rcx,%rdi,4),%xmm0
movss  (%rbx,%rsi,4),%xmm1
ucomiss %xmm0,%xmm1
jbe    1180 <spectral_median+0x70>
movss  %xmm0,(%rbx,%rsi,4)
movss  %xmm1,(%rcx,%rdi,4)
jmp    1180 <spectral_median+0x70>
mov    %r14d,%eax
shr    $0x1f,%eax
add    %r14d,%eax
sar    %eax
cltq
movss  (%rbx,%rax,4),%xmm0
test   $0x1,%r14b
jne    11ca <spectral_median+0xba>
addss  -0x4(%rbx,%rax,4),%xmm0
mulss  0xe36(%rip),%xmm0
lea    -0x18(%rbp),%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define float @spectral_median(float* %0, i32 %1) {
%3 = add i32 %1, 1
%4 = zext i32 %3 to i64
%5 = alloca float, i64 %4
%6 = shl i32 %3, 2
%7 = call i32 @memcpy(float* %5, float* %0, i32 %6)
%8 = icmp sgt i32 %1, 0
br i1 %8, label %9, label %29
9: ; preds = %2
%10 = zext i32 %1 to i64
br label %14
11: ; preds = %26
%12 = add nuw nsw i64 %16, 1
%13 = icmp eq i64 %17, %10
br i1 %13, label %29, label %14
14: ; preds = %11, %9
%15 = phi i64 [ 0, %9 ], [ %17, %11 ]
%16 = phi i64 [ 1, %9 ], [ %12, %11 ]
%17 = add nuw nsw i64 %15, 1
%18 = getelementptr inbounds float, float* %5, i64 %15
br label %19
19: ; preds = %14, %26
%20 = phi i64 [ %16, %14 ], [ %27, %26 ]
%21 = getelementptr inbounds float, float* %5, i64 %20
%22 = load float, float* %21
%23 = load float, float* %18
%24 = fcmp olt float %22, %23
br i1 %24, label %25, label %26
25: ; preds = %19
store float %22, float* %18
store float %23, float* %21
br label %26
26: ; preds = %19, %25
%27 = add nuw nsw i64 %20, 1
%28 = icmp eq i64 %27, %4
br i1 %28, label %11, label %19
29: ; preds = %11, %2
%30 = and i32 %1, 1
%31 = icmp eq i32 %30, 0
%32 = sdiv i32 %1, 2
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %5, i64 %33
%35 = load float, float* %34
br i1 %31, label %36, label %43
36: ; preds = %29
%37 = add nsw i32 %32, -1
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds float, float* %5, i64 %38
%40 = load float, float* %39
%41 = fadd float %35, %40
%42 = fmul float %41, 5.000000e-01
br label %43
43: ; preds = %29, %36
%44 = phi float [ %42, %36 ], [ %35, %29 ]
ret float %44
}