#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 */
short median_short(int, short *);
scalar_t__ sqrt(double);
short meddev_short(int n, short *data) {
  short medval1, var2, var3, dn;
  int i;
  if (n <= 1)
    return 0.0;
  medval1 = median_short(n, data);
  var2 = 0.;
  dn = (short)n - 1;
  for (i = 0; i < n; i++) {
    var3 = data[i] - medval1;
    var2 += var3 * var3 / (short)dn;
  }
  return ((short)(sqrt((double)var2)));
}


x86 assembly (clang, -O0):
<meddev_short>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
cmpl   $0x1,-0x8(%rbp)
jg     1144 <meddev_short+0x24>
movw   $0x0,-0x2(%rbp)
jmp    11cb <meddev_short+0xab>
mov    -0x8(%rbp),%edi
mov    -0x10(%rbp),%rsi
call   1040 <median_short@plt>
mov    %ax,-0x12(%rbp)
movw   $0x0,-0x14(%rbp)
mov    -0x8(%rbp),%eax
cwtl
sub    $0x1,%eax
mov    %ax,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11ba <meddev_short+0x9a>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
movswl -0x12(%rbp),%ecx
sub    %ecx,%eax
mov    %ax,-0x16(%rbp)
movswl -0x16(%rbp),%eax
movswl -0x16(%rbp),%ecx
imul   %ecx,%eax
movswl -0x18(%rbp),%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
movswl -0x14(%rbp),%eax
add    %ecx,%eax
mov    %ax,-0x14(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    116c <meddev_short+0x4c>
movswl -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm0
call   1030 <sqrt@plt>
mov    %ax,-0x2(%rbp)
movswl -0x2(%rbp),%eax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<meddev_short>:
cmp    $0x2,%edi
jge    1128 <meddev_short+0x8>
xor    %eax,%eax
ret
push   %r15
push   %r14
push   %rbx
mov    %rsi,%r15
mov    %edi,%r14d
call   1040 <median_short@plt>
mov    %eax,%ecx
lea    -0x1(%r14),%eax
movswl %ax,%esi
mov    %r14d,%r8d
lea    (%r8,%r8,1),%r9
and    $0xfffffffffffffffc,%r9
xor    %ebx,%ebx
xor    %edi,%edi
movzwl (%r15,%rbx,1),%eax
sub    %ecx,%eax
cwtl
imul   %eax,%eax
cltd
idiv   %esi
add    %eax,%edi
movzwl 0x2(%r15,%rbx,1),%eax
sub    %ecx,%eax
cwtl
imul   %eax,%eax
cltd
idiv   %esi
add    %edi,%eax
mov    %eax,%edi
add    $0x4,%rbx
cmp    %rbx,%r9
jne    1150 <meddev_short+0x30>
test   $0x1,%r8b
je     1192 <meddev_short+0x72>
movzwl (%r15,%rbx,1),%eax
sub    %ecx,%eax
cwtl
imul   %eax,%eax
cltd
idiv   %esi
add    %eax,%edi
movswl %di,%eax
cvtsi2sd %eax,%xmm0
call   1030 <sqrt@plt>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define signext i16 @meddev_short(i32 %0, i16* %1) {
%3 = alloca i16
%4 = alloca i32
%5 = alloca i16*
%6 = alloca i16
%7 = alloca i16
%8 = alloca i16
%9 = alloca i16
%10 = alloca i32
store i32 %0, i32* %4
store i16* %1, i16** %5
%11 = load i32, i32* %4
%12 = icmp sle i32 %11, 1
br i1 %12, label %13, label %14
13: ; preds = %2
store i16 0, i16* %3
br label %58
14: ; preds = %2
%15 = load i32, i32* %4
%16 = load i16*, i16** %5
%17 = call signext i16 @median_short(i32 %15, i16* %16)
store i16 %17, i16* %6
store i16 0, i16* %7
%18 = load i32, i32* %4
%19 = trunc i32 %18 to i16
%20 = sext i16 %19 to i32
%21 = sub nsw i32 %20, 1
%22 = trunc i32 %21 to i16
store i16 %22, i16* %9
store i32 0, i32* %10
br label %23
23: ; preds = %50, %14
%24 = load i32, i32* %10
%25 = load i32, i32* %4
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %53
27: ; preds = %23
%28 = load i16*, i16** %5
%29 = load i32, i32* %10
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds i16, i16* %28, i64 %30
%32 = load i16, i16* %31
%33 = sext i16 %32 to i32
%34 = load i16, i16* %6
%35 = sext i16 %34 to i32
%36 = sub nsw i32 %33, %35
%37 = trunc i32 %36 to i16
store i16 %37, i16* %8
%38 = load i16, i16* %8
%39 = sext i16 %38 to i32
%40 = load i16, i16* %8
%41 = sext i16 %40 to i32
%42 = mul nsw i32 %39, %41
%43 = load i16, i16* %9
%44 = sext i16 %43 to i32
%45 = sdiv i32 %42, %44
%46 = load i16, i16* %7
%47 = sext i16 %46 to i32
%48 = add nsw i32 %47, %45
%49 = trunc i32 %48 to i16
store i16 %49, i16* %7
br label %50
50: ; preds = %27
%51 = load i32, i32* %10
%52 = add nsw i32 %51, 1
store i32 %52, i32* %10
br label %23
53: ; preds = %23
%54 = load i16, i16* %7
%55 = sitofp i16 %54 to double
%56 = call i64 @sqrt(double %55)
%57 = trunc i64 %56 to i16
store i16 %57, i16* %3
br label %58
58: ; preds = %53, %13
%59 = load i16, i16* %3
ret i16 %59
}

LLVM IR (-O2):
define signext i16 @meddev_short(i32 %0, i16* %1) {
%3 = icmp slt i32 %0, 2
br i1 %3, label %52, label %4
4: ; preds = %2
%5 = tail call signext i16 @median_short(i32 %0, i16* %1)
%6 = shl i32 %0, 16
%7 = add i32 %6, -65536
%8 = ashr exact i32 %7, 16
%9 = zext i32 %0 to i64
%10 = and i64 %9, 1
%11 = and i64 %9, 4294967294
br label %12
12: ; preds = %12, %4
%13 = phi i64 [ 0, %4 ], [ %33, %12 ]
%14 = phi i16 [ 0, %4 ], [ %32, %12 ]
%15 = phi i64 [ 0, %4 ], [ %34, %12 ]
%16 = getelementptr inbounds i16, i16* %1, i64 %13
%17 = load i16, i16* %16
%18 = sub i16 %17, %5
%19 = sext i16 %18 to i32
%20 = mul nsw i32 %19, %19
%21 = sdiv i32 %20, %8
%22 = trunc i32 %21 to i16
%23 = add i16 %14, %22
%24 = or i64 %13, 1
%25 = getelementptr inbounds i16, i16* %1, i64 %24
%26 = load i16, i16* %25
%27 = sub i16 %26, %5
%28 = sext i16 %27 to i32
%29 = mul nsw i32 %28, %28
%30 = sdiv i32 %29, %8
%31 = trunc i32 %30 to i16
%32 = add i16 %23, %31
%33 = add nuw nsw i64 %13, 2
%34 = add i64 %15, 2
%35 = icmp eq i64 %34, %11
br i1 %35, label %36, label %12
36: ; preds = %12
%37 = icmp eq i64 %10, 0
br i1 %37, label %47, label %38
38: ; preds = %36
%39 = getelementptr inbounds i16, i16* %1, i64 %33
%40 = load i16, i16* %39
%41 = sub i16 %40, %5
%42 = sext i16 %41 to i32
%43 = mul nsw i32 %42, %42
%44 = sdiv i32 %43, %8
%45 = trunc i32 %44 to i16
%46 = add i16 %32, %45
br label %47
47: ; preds = %36, %38
%48 = phi i16 [ %32, %36 ], [ %46, %38 ]
%49 = sitofp i16 %48 to double
%50 = tail call i64 @sqrt(double %49)
%51 = trunc i64 %50 to i16
br label %52
52: ; preds = %2, %47
%53 = phi i16 [ %51, %47 ], [ 0, %2 ]
ret i16 %53
}