#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(double, double);
double sqrt(double);
double *omnibus_moments(double *x, int n) {
  double diff, mean = 0., fssm, tssm, sum_cube = 0., sum_four = 0., sum_sq = 0.;
  static double y[2];
  int i;
  for (i = 0; i < n; ++i)
    mean += x[i];
  mean /= n;
  for (i = 0; i < n; ++i) {
    diff = x[i] - mean;
    sum_sq += diff * diff;
    sum_cube += diff * diff * diff;
    sum_four += diff * diff * diff * diff;
  }
  tssm = sqrt((double)n) * sum_cube / pow(sum_sq, 1.5);
  fssm = n * sum_four / (sum_sq * sum_sq);
  y[0] = tssm;
  y[1] = fssm;
  return y;
}


x86 assembly (clang, -O0):
<omnibus_moments>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1187 <omnibus_moments+0x67>
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1156 <omnibus_moments+0x36>
cvtsi2sdl -0xc(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1235 <omnibus_moments+0x115>
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x18(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11a1 <omnibus_moments+0x81>
cvtsi2sdl -0xc(%rbp),%xmm0
call   1040 <sqrt@plt>
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  0xdaa(%rip),%xmm1
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0x58(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,0x2daa(%rip)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,0x2da5(%rip)
lea    0x2d96(%rip),%rax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<omnibus_moments>:
sub    $0x48,%rsp
test   %esi,%esi
jle    1143 <omnibus_moments+0x23>
mov    %esi,%r8d
lea    -0x1(%r8),%r9
mov    %r8d,%edx
and    $0x7,%edx
cmp    $0x7,%r9
jae    114c <omnibus_moments+0x2c>
xorpd  %xmm5,%xmm5
xor    %eax,%eax
jmp    1198 <omnibus_moments+0x78>
cvtsi2sd %esi,%xmm0
jmp    1264 <omnibus_moments+0x144>
mov    %r8d,%ecx
and    $0xfffffff8,%ecx
xorpd  %xmm5,%xmm5
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
addsd  (%rdi,%rax,8),%xmm5
addsd  0x8(%rdi,%rax,8),%xmm5
addsd  0x10(%rdi,%rax,8),%xmm5
addsd  0x18(%rdi,%rax,8),%xmm5
addsd  0x20(%rdi,%rax,8),%xmm5
addsd  0x28(%rdi,%rax,8),%xmm5
addsd  0x30(%rdi,%rax,8),%xmm5
addsd  0x38(%rdi,%rax,8),%xmm5
add    $0x8,%rax
cmp    %rax,%rcx
jne    1160 <omnibus_moments+0x40>
test   %rdx,%rdx
je     11be <omnibus_moments+0x9e>
lea    (%rdi,%rax,8),%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rax,%rcx,8),%xmm5
add    $0x1,%rcx
cmp    %rcx,%rdx
jne    11b0 <omnibus_moments+0x90>
cvtsi2sd %esi,%xmm0
test   %esi,%esi
jle    1264 <omnibus_moments+0x144>
divsd  %xmm0,%xmm5
test   %r9,%r9
je     126e <omnibus_moments+0x14e>
mov    %r8d,%edx
and    $0xfffffffe,%edx
xorpd  %xmm1,%xmm1
xorpd  %xmm6,%xmm6
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rcx,8),%xmm2
movsd  0x8(%rdi,%rcx,8),%xmm7
subsd  %xmm5,%xmm2
movapd %xmm2,%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm3,%xmm6
movapd %xmm2,%xmm4
mulsd  %xmm3,%xmm4
unpcklpd %xmm4,%xmm3
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
subsd  %xmm5,%xmm7
movapd %xmm7,%xmm1
mulsd  %xmm7,%xmm1
addsd  %xmm1,%xmm6
movapd %xmm7,%xmm3
mulsd  %xmm1,%xmm3
unpcklpd %xmm3,%xmm1
unpcklpd %xmm7,%xmm7
mulpd  %xmm1,%xmm7
movapd %xmm7,%xmm1
addpd  %xmm2,%xmm1
add    $0x2,%rcx
cmp    %rcx,%rdx
jne    11f0 <omnibus_moments+0xd0>
addpd  %xmm2,%xmm7
test   $0x1,%r8b
jne    127e <omnibus_moments+0x15e>
jmp    12ab <omnibus_moments+0x18b>
xorpd  %xmm7,%xmm7
xorpd  %xmm6,%xmm6
jmp    12ab <omnibus_moments+0x18b>
xorpd  %xmm1,%xmm1
xorpd  %xmm6,%xmm6
xor    %ecx,%ecx
test   $0x1,%r8b
je     12ab <omnibus_moments+0x18b>
movsd  (%rdi,%rcx,8),%xmm7
subsd  %xmm5,%xmm7
movapd %xmm7,%xmm3
mulsd  %xmm7,%xmm3
addsd  %xmm3,%xmm6
movapd %xmm7,%xmm2
mulsd  %xmm3,%xmm2
unpcklpd %xmm2,%xmm3
unpcklpd %xmm7,%xmm7
mulpd  %xmm3,%xmm7
addpd  %xmm1,%xmm7
movapd %xmm7,0x30(%rsp)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
movapd %xmm6,0x10(%rsp)
movapd %xmm0,0x20(%rsp)
jb     12d5 <omnibus_moments+0x1b5>
xorps  %xmm1,%xmm1
sqrtsd %xmm0,%xmm1
movapd %xmm1,(%rsp)
jmp    12e5 <omnibus_moments+0x1c5>
call   1040 <sqrt@plt>
movapd 0x10(%rsp),%xmm6
movapd %xmm0,(%rsp)
movsd  0xd13(%rip),%xmm1
movapd %xmm6,%xmm0
call   1030 <pow@plt>
movapd (%rsp),%xmm2
unpcklpd 0x20(%rsp),%xmm2
movapd 0x30(%rsp),%xmm1
mulpd  %xmm2,%xmm1
movapd 0x10(%rsp),%xmm2
mulsd  %xmm2,%xmm2
unpcklpd %xmm2,%xmm0
divpd  %xmm0,%xmm1
movapd %xmm1,0x2d1b(%rip)
lea    0x2d14(%rip),%rax
add    $0x48,%rsp
ret

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

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