#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 pow(int, double);
double quantile(double *, int, double);
double smoothPar(double *x, int n) {
  int i;
  double M = 0;
  double S = 0;
  double Q[2];
  for (i = 0; i < n; i++)
    M += x[i];
  M /= n;
  for (i = 0; i < n; i++)
    S += pow(x[i] - M, 2);
  S /= (n - 1);
  S = pow(S, 0.5);
  Q[0] = quantile(x, n, 0.25);
  Q[1] = quantile(x, n, 0.75);
  Q[0] = (Q[1] - Q[0]) / 1.349;
  if (S < Q[0])
    Q[0] = S;
  Q[0] = (0.51 / pow(n, 0.2)) * Q[0];
  return (Q[0]);
}


x86 assembly (clang, -O0):
<smoothPar>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1177 <smoothPar+0x57>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1146 <smoothPar+0x26>
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11d8 <smoothPar+0xb8>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x18(%rbp),%xmm0
cvttsd2si %xmm0,%edi
movsd  0xe75(%rip),%xmm0
call   1030 <pow@plt>
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1191 <smoothPar+0x71>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  -0x20(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
cvttsd2si -0x20(%rbp),%edi
movsd  0xe1b(%rip),%xmm0
call   1030 <pow@plt>
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
movsd  0xdfa(%rip),%xmm0
call   1040 <quantile@plt>
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
movsd  0xdd9(%rip),%xmm0
call   1040 <quantile@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  0xdb5(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1272 <smoothPar+0x152>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0xc(%rbp),%edi
movsd  0xdab(%rip),%xmm0
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xd93(%rip),%xmm0
divsd  %xmm1,%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<smoothPar>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %esi,%r14d
mov    %rdi,%r13
test   %esi,%esi
jle    1219 <smoothPar+0xf9>
mov    %r14d,%r15d
lea    -0x1(%r15),%rax
mov    %r15d,%ecx
and    $0x7,%ecx
cmp    $0x7,%rax
jae    1156 <smoothPar+0x36>
xorpd  %xmm0,%xmm0
xor    %edx,%edx
jmp    11b1 <smoothPar+0x91>
mov    %r15d,%esi
and    $0xfffffff8,%esi
xorpd  %xmm0,%xmm0
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addsd  0x0(%r13,%rdx,8),%xmm0
addsd  0x8(%r13,%rdx,8),%xmm0
addsd  0x10(%r13,%rdx,8),%xmm0
addsd  0x18(%r13,%rdx,8),%xmm0
addsd  0x20(%r13,%rdx,8),%xmm0
addsd  0x28(%r13,%rdx,8),%xmm0
addsd  0x30(%r13,%rdx,8),%xmm0
addsd  0x38(%r13,%rdx,8),%xmm0
add    $0x8,%rdx
cmp    %rdx,%rsi
jne    1170 <smoothPar+0x50>
test   %rcx,%rcx
je     11de <smoothPar+0xbe>
lea    0x0(,%rdx,8),%rdx
add    %r13,%rdx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rdx,%rsi,8),%xmm0
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    11d0 <smoothPar+0xb0>
test   %r14d,%r14d
jle    1219 <smoothPar+0xf9>
movapd %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
mov    %r15d,%r12d
and    $0x3,%r12d
cmp    $0x3,%rax
movsd  %xmm1,0x8(%rsp)
jae    1222 <smoothPar+0x102>
xorpd  %xmm2,%xmm2
xor    %ebx,%ebx
movapd %xmm0,%xmm1
jmp    12e5 <smoothPar+0x1c5>
xorpd  %xmm2,%xmm2
jmp    1338 <smoothPar+0x218>
and    $0xfffffffc,%r15d
xorpd  %xmm2,%xmm2
xor    %ebx,%ebx
movapd %xmm0,%xmm1
movsd  %xmm2,(%rsp)
movsd  0x0(%r13,%rbx,8),%xmm0
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
movsd  0xdb4(%rip),%xmm0
call   1030 <pow@plt>
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%r13,%rbx,8),%xmm0
subsd  0x8(%rsp),%xmm0
cvttsd2si %xmm0,%edi
movsd  0xd8c(%rip),%xmm0
call   1030 <pow@plt>
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x10(%r13,%rbx,8),%xmm0
subsd  0x8(%rsp),%xmm0
cvttsd2si %xmm0,%edi
movsd  0xd64(%rip),%xmm0
call   1030 <pow@plt>
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x18(%r13,%rbx,8),%xmm0
subsd  0x8(%rsp),%xmm0
cvttsd2si %xmm0,%edi
movsd  0xd3c(%rip),%xmm0
call   1030 <pow@plt>
movsd  0x8(%rsp),%xmm1
movapd %xmm0,%xmm2
addsd  (%rsp),%xmm2
add    $0x4,%rbx
cmp    %rbx,%r15
jne    1230 <smoothPar+0x110>
test   %r12,%r12
je     1338 <smoothPar+0x218>
lea    0x0(,%rbx,8),%r15
add    %r13,%r15
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movsd  %xmm2,(%rsp)
movsd  (%r15,%rbx,8),%xmm0
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
movsd  0xce5(%rip),%xmm0
call   1030 <pow@plt>
movsd  (%rsp),%xmm2
movsd  0x8(%rsp),%xmm1
addsd  %xmm0,%xmm2
add    $0x1,%rbx
cmp    %rbx,%r12
jne    1300 <smoothPar+0x1e0>
lea    -0x1(%r14),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
divsd  %xmm0,%xmm2
cvttsd2si %xmm2,%edi
movsd  0xcb5(%rip),%xmm0
call   1030 <pow@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0xcaa(%rip),%xmm0
mov    %r13,%rdi
mov    %r14d,%esi
call   1040 <quantile@plt>
movsd  %xmm0,(%rsp)
movsd  0xc9a(%rip),%xmm0
mov    %r13,%rdi
mov    %r14d,%esi
call   1040 <quantile@plt>
subsd  (%rsp),%xmm0
divsd  0xc8a(%rip),%xmm0
movsd  0x8(%rsp),%xmm1
minsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  0xc7a(%rip),%xmm0
mov    %r14d,%edi
call   1030 <pow@plt>
movsd  0xc72(%rip),%xmm1
divsd  %xmm0,%xmm1
mulsd  0x8(%rsp),%xmm1
movapd %xmm1,%xmm0
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define double @smoothPar(double* %0, i32 %1) {
%3 = alloca double*
%4 = alloca i32
%5 = alloca i32
%6 = alloca double
%7 = alloca double
%8 = alloca [2 x double]
store double* %0, double** %3
store i32 %1, i32* %4
store double 0.000000e+00, double* %6
store double 0.000000e+00, double* %7
store i32 0, i32* %5
br label %9
9: ; preds = %21, %2
%10 = load i32, i32* %5
%11 = load i32, i32* %4
%12 = icmp slt i32 %10, %11
br i1 %12, label %13, label %24
13: ; preds = %9
%14 = load double*, double** %3
%15 = load i32, i32* %5
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds double, double* %14, i64 %16
%18 = load double, double* %17
%19 = load double, double* %6
%20 = fadd double %19, %18
store double %20, double* %6
br label %21
21: ; preds = %13
%22 = load i32, i32* %5
%23 = add nsw i32 %22, 1
store i32 %23, i32* %5
br label %9
24: ; preds = %9
%25 = load i32, i32* %4
%26 = sitofp i32 %25 to double
%27 = load double, double* %6
%28 = fdiv double %27, %26
store double %28, double* %6
store i32 0, i32* %5
br label %29
29: ; preds = %45, %24
%30 = load i32, i32* %5
%31 = load i32, i32* %4
%32 = icmp slt i32 %30, %31
br i1 %32, label %33, label %48
33: ; preds = %29
%34 = load double*, double** %3
%35 = load i32, i32* %5
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
%39 = load double, double* %6
%40 = fsub double %38, %39
%41 = fptosi double %40 to i32
%42 = call double @pow(i32 %41, double 2.000000e+00)
%43 = load double, double* %7
%44 = fadd double %43, %42
store double %44, double* %7
br label %45
45: ; preds = %33
%46 = load i32, i32* %5
%47 = add nsw i32 %46, 1
store i32 %47, i32* %5
br label %29
48: ; preds = %29
%49 = load i32, i32* %4
%50 = sub nsw i32 %49, 1
%51 = sitofp i32 %50 to double
%52 = load double, double* %7
%53 = fdiv double %52, %51
store double %53, double* %7
%54 = load double, double* %7
%55 = fptosi double %54 to i32
%56 = call double @pow(i32 %55, double 5.000000e-01)
store double %56, double* %7
%57 = load double*, double** %3
%58 = load i32, i32* %4
%59 = call double @quantile(double* %57, i32 %58, double 2.500000e-01)
%60 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
store double %59, double* %60
%61 = load double*, double** %3
%62 = load i32, i32* %4
%63 = call double @quantile(double* %61, i32 %62, double 7.500000e-01)
%64 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 1
store double %63, double* %64
%65 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 1
%66 = load double, double* %65
%67 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
%68 = load double, double* %67
%69 = fsub double %66, %68
%70 = fdiv double %69, 1.349000e+00
%71 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
store double %70, double* %71
%72 = load double, double* %7
%73 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
%74 = load double, double* %73
%75 = fcmp olt double %72, %74
br i1 %75, label %76, label %79
76: ; preds = %48
%77 = load double, double* %7
%78 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
store double %77, double* %78
br label %79
79: ; preds = %76, %48
%80 = load i32, i32* %4
%81 = call double @pow(i32 %80, double 2.000000e-01)
%82 = fdiv double 5.100000e-01, %81
%83 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
%84 = load double, double* %83
%85 = fmul double %82, %84
%86 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
store double %85, double* %86
%87 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
%88 = load double, double* %87
ret double %88
}

LLVM IR (-O2):
define double @smoothPar(double* %0, i32 %1) {
%3 = icmp sgt i32 %1, 0
br i1 %3, label %4, label %125
4: ; preds = %2
%5 = zext i32 %1 to i64
%6 = add nsw i64 %5, -1
%7 = and i64 %5, 7
%8 = icmp ult i64 %6, 7
br i1 %8, label %49, label %9
9: ; preds = %4
%10 = and i64 %5, 4294967288
br label %11
11: ; preds = %11, %9
%12 = phi i64 [ 0, %9 ], [ %46, %11 ]
%13 = phi double [ 0.000000e+00, %9 ], [ %45, %11 ]
%14 = phi i64 [ 0, %9 ], [ %47, %11 ]
%15 = getelementptr inbounds double, double* %0, i64 %12
%16 = load double, double* %15
%17 = fadd double %13, %16
%18 = or i64 %12, 1
%19 = getelementptr inbounds double, double* %0, i64 %18
%20 = load double, double* %19
%21 = fadd double %17, %20
%22 = or i64 %12, 2
%23 = getelementptr inbounds double, double* %0, i64 %22
%24 = load double, double* %23
%25 = fadd double %21, %24
%26 = or i64 %12, 3
%27 = getelementptr inbounds double, double* %0, i64 %26
%28 = load double, double* %27
%29 = fadd double %25, %28
%30 = or i64 %12, 4
%31 = getelementptr inbounds double, double* %0, i64 %30
%32 = load double, double* %31
%33 = fadd double %29, %32
%34 = or i64 %12, 5
%35 = getelementptr inbounds double, double* %0, i64 %34
%36 = load double, double* %35
%37 = fadd double %33, %36
%38 = or i64 %12, 6
%39 = getelementptr inbounds double, double* %0, i64 %38
%40 = load double, double* %39
%41 = fadd double %37, %40
%42 = or i64 %12, 7
%43 = getelementptr inbounds double, double* %0, i64 %42
%44 = load double, double* %43
%45 = fadd double %41, %44
%46 = add nuw nsw i64 %12, 8
%47 = add i64 %14, 8
%48 = icmp eq i64 %47, %10
br i1 %48, label %49, label %11
49: ; preds = %11, %4
%50 = phi double [ undef, %4 ], [ %45, %11 ]
%51 = phi i64 [ 0, %4 ], [ %46, %11 ]
%52 = phi double [ 0.000000e+00, %4 ], [ %45, %11 ]
%53 = icmp eq i64 %7, 0
br i1 %53, label %64, label %54
54: ; preds = %49, %54
%55 = phi i64 [ %61, %54 ], [ %51, %49 ]
%56 = phi double [ %60, %54 ], [ %52, %49 ]
%57 = phi i64 [ %62, %54 ], [ 0, %49 ]
%58 = getelementptr inbounds double, double* %0, i64 %55
%59 = load double, double* %58
%60 = fadd double %56, %59
%61 = add nuw nsw i64 %55, 1
%62 = add i64 %57, 1
%63 = icmp eq i64 %62, %7
br i1 %63, label %64, label %54
64: ; preds = %54, %49
%65 = phi double [ %50, %49 ], [ %60, %54 ]
%66 = sitofp i32 %1 to double
%67 = fdiv double %65, %66
br i1 %3, label %68, label %125
68: ; preds = %64
%69 = and i64 %5, 3
%70 = icmp ult i64 %6, 3
br i1 %70, label %107, label %71
71: ; preds = %68
%72 = and i64 %5, 4294967292
br label %73
73: ; preds = %73, %71
%74 = phi i64 [ 0, %71 ], [ %104, %73 ]
%75 = phi double [ 0.000000e+00, %71 ], [ %103, %73 ]
%76 = phi i64 [ 0, %71 ], [ %105, %73 ]
%77 = getelementptr inbounds double, double* %0, i64 %74
%78 = load double, double* %77
%79 = fsub double %78, %67
%80 = fptosi double %79 to i32
%81 = tail call double @pow(i32 %80, double 2.000000e+00)
%82 = fadd double %75, %81
%83 = or i64 %74, 1
%84 = getelementptr inbounds double, double* %0, i64 %83
%85 = load double, double* %84
%86 = fsub double %85, %67
%87 = fptosi double %86 to i32
%88 = tail call double @pow(i32 %87, double 2.000000e+00)
%89 = fadd double %82, %88
%90 = or i64 %74, 2
%91 = getelementptr inbounds double, double* %0, i64 %90
%92 = load double, double* %91
%93 = fsub double %92, %67
%94 = fptosi double %93 to i32
%95 = tail call double @pow(i32 %94, double 2.000000e+00)
%96 = fadd double %89, %95
%97 = or i64 %74, 3
%98 = getelementptr inbounds double, double* %0, i64 %97
%99 = load double, double* %98
%100 = fsub double %99, %67
%101 = fptosi double %100 to i32
%102 = tail call double @pow(i32 %101, double 2.000000e+00)
%103 = fadd double %96, %102
%104 = add nuw nsw i64 %74, 4
%105 = add i64 %76, 4
%106 = icmp eq i64 %105, %72
br i1 %106, label %107, label %73
107: ; preds = %73, %68
%108 = phi double [ undef, %68 ], [ %103, %73 ]
%109 = phi i64 [ 0, %68 ], [ %104, %73 ]
%110 = phi double [ 0.000000e+00, %68 ], [ %103, %73 ]
%111 = icmp eq i64 %69, 0
br i1 %111, label %125, label %112
112: ; preds = %107, %112
%113 = phi i64 [ %122, %112 ], [ %109, %107 ]
%114 = phi double [ %121, %112 ], [ %110, %107 ]
%115 = phi i64 [ %123, %112 ], [ 0, %107 ]
%116 = getelementptr inbounds double, double* %0, i64 %113
%117 = load double, double* %116
%118 = fsub double %117, %67
%119 = fptosi double %118 to i32
%120 = tail call double @pow(i32 %119, double 2.000000e+00)
%121 = fadd double %114, %120
%122 = add nuw nsw i64 %113, 1
%123 = add i64 %115, 1
%124 = icmp eq i64 %123, %69
br i1 %124, label %125, label %112
125: ; preds = %107, %112, %2, %64
%126 = phi double [ 0.000000e+00, %64 ], [ 0.000000e+00, %2 ], [ %108, %107 ], [ %121, %112 ]
%127 = add nsw i32 %1, -1
%128 = sitofp i32 %127 to double
%129 = fdiv double %126, %128
%130 = fptosi double %129 to i32
%131 = tail call double @pow(i32 %130, double 5.000000e-01)
%132 = tail call double @quantile(double* %0, i32 %1, double 2.500000e-01)
%133 = tail call double @quantile(double* %0, i32 %1, double 7.500000e-01)
%134 = fsub double %133, %132
%135 = fdiv double %134, 1.349000e+00
%136 = fcmp olt double %131, %135
%137 = select i1 %136, double %131, double %135
%138 = tail call double @pow(i32 %1, double 2.000000e-01)
%139 = fdiv double 5.100000e-01, %138
%140 = fmul double %139, %137
ret double %140
}