#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 */
typedef int int64_t;

/* Variables and functions */
double stdlib_strided_dvariancech(const int64_t N, const double correction,
                                  const double *X, const int64_t stride) {
  int64_t ix;
  int64_t i;
  double dN;
  double mu;
  double M2;
  double M;
  double n;
  double d;
  dN = (double)N;
  n = dN - correction;
  if (N <= 0 || n <= 0.0) {
    return 0.0 / 0.0;
  }
  if (N == 1 || stride == 0) {
    return 0.0;
  }
  if (stride < 0) {
    ix = (1 - N) * stride;
  } else {
    ix = 0;
  }
  mu = X[ix];
  ix += stride;
  M2 = 0.0;
  M = 0.0;
  for (i = 1; i < N; i++) {
    d = X[ix] - mu;
    M2 += d * d;
    M += d;
    ix += stride;
  }
  return (M2 / n) - ((M / dN) * (M / n));
}


x86 assembly (clang, -O0):
<stdlib_strided_dvariancech>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %edx,-0x24(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
cmpl   $0x0,-0xc(%rbp)
jle    1144 <stdlib_strided_dvariancech+0x44>
xorps  %xmm0,%xmm0
ucomisd -0x58(%rbp),%xmm0
jb     1156 <stdlib_strided_dvariancech+0x56>
movsd  0xeb4(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1275 <stdlib_strided_dvariancech+0x175>
cmpl   $0x1,-0xc(%rbp)
je     116a <stdlib_strided_dvariancech+0x6a>
cmpl   $0x0,-0x24(%rbp)
jne    1177 <stdlib_strided_dvariancech+0x77>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1275 <stdlib_strided_dvariancech+0x175>
cmpl   $0x0,-0x24(%rbp)
jge    1195 <stdlib_strided_dvariancech+0x95>
mov    $0x1,%eax
sub    -0xc(%rbp),%eax
imul   -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
jmp    119c <stdlib_strided_dvariancech+0x9c>
movl   $0x0,-0x28(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x24(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1233 <stdlib_strided_dvariancech+0x133>
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x24(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ce <stdlib_strided_dvariancech+0xce>
movsd  -0x48(%rbp),%xmm1
divsd  -0x58(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
divsd  -0x58(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<stdlib_strided_dvariancech>:
movapd %xmm0,%xmm1
movsd  0xef4(%rip),%xmm0
test   %edi,%edi
jle    1201 <stdlib_strided_dvariancech+0x101>
cvtsi2sd %edi,%xmm2
movapd %xmm2,%xmm3
subsd  %xmm1,%xmm3
xorpd  %xmm1,%xmm1
ucomisd %xmm3,%xmm1
jae    1201 <stdlib_strided_dvariancech+0x101>
xorpd  %xmm0,%xmm0
cmp    $0x1,%edi
je     1201 <stdlib_strided_dvariancech+0x101>
test   %edx,%edx
je     1201 <stdlib_strided_dvariancech+0x101>
mov    $0x1,%ecx
sub    %edi,%ecx
imul   %edx,%ecx
movslq %edx,%rax
sar    $0x1f,%edx
and    %ecx,%edx
movslq %edx,%rdx
movsd  (%rsi,%rdx,8),%xmm4
lea    -0x1(%rdi),%r8d
cmp    $0x2,%edi
jne    1170 <stdlib_strided_dvariancech+0x70>
xorpd  %xmm1,%xmm1
xorpd  %xmm0,%xmm0
jmp    11cb <stdlib_strided_dvariancech+0xcb>
mov    %r8d,%edi
and    $0xfffffffe,%edi
mov    %rax,%r10
shl    $0x4,%r10
add    %rsi,%r10
lea    (%rsi,%rax,8),%r9
xorpd  %xmm1,%xmm1
mov    %rdx,%rcx
xorpd  %xmm0,%xmm0
nop
lea    (%rcx,%rax,1),%rdx
movsd  (%r9,%rcx,8),%xmm5
subsd  %xmm4,%xmm5
addsd  %xmm5,%xmm1
mulsd  %xmm5,%xmm5
addsd  %xmm0,%xmm5
add    %rax,%rdx
movsd  (%r10,%rcx,8),%xmm0
subsd  %xmm4,%xmm0
addsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm5,%xmm0
mov    %rdx,%rcx
add    $0xfffffffe,%edi
jne    1190 <stdlib_strided_dvariancech+0x90>
test   $0x1,%r8b
je     11e9 <stdlib_strided_dvariancech+0xe9>
add    %rax,%rdx
movsd  (%rsi,%rdx,8),%xmm5
subsd  %xmm4,%xmm5
addsd  %xmm5,%xmm1
mulsd  %xmm5,%xmm5
addsd  %xmm5,%xmm0
divsd  %xmm3,%xmm0
movapd %xmm1,%xmm4
divsd  %xmm3,%xmm4
divsd  %xmm2,%xmm1
mulsd  %xmm4,%xmm1
subsd  %xmm1,%xmm0
ret

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

LLVM IR (-O2):
define double @stdlib_strided_dvariancech(i32 %0, double %1, double* nocapture %2, i32 %3) {
%5 = sitofp i32 %0 to double
%6 = fsub double %5, %1
%7 = icmp slt i32 %0, 1
%8 = fcmp ole double %6, 0.000000e+00
%9 = or i1 %7, %8
br i1 %9, label %69, label %10
10: ; preds = %4
%11 = icmp eq i32 %0, 1
%12 = icmp eq i32 %3, 0
%13 = or i1 %11, %12
br i1 %13, label %69, label %14
14: ; preds = %10
%15 = icmp slt i32 %3, 0
%16 = sub nsw i32 1, %0
%17 = mul nsw i32 %16, %3
%18 = select i1 %15, i32 %17, i32 0
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds double, double* %2, i64 %19
%21 = load double, double* %20
%22 = sext i32 %3 to i64
%23 = add i32 %0, -1
%24 = and i32 %23, 1
%25 = icmp eq i32 %0, 2
br i1 %25, label %47, label %26
26: ; preds = %14
%27 = and i32 %23, -2
br label %28
28: ; preds = %28, %26
%29 = phi i64 [ %19, %26 ], [ %39, %28 ]
%30 = phi double [ 0.000000e+00, %26 ], [ %44, %28 ]
%31 = phi double [ 0.000000e+00, %26 ], [ %43, %28 ]
%32 = phi i32 [ 0, %26 ], [ %45, %28 ]
%33 = add i64 %29, %22
%34 = getelementptr inbounds double, double* %2, i64 %33
%35 = load double, double* %34
%36 = fsub double %35, %21
%37 = tail call double @llvm.fmuladd.f64(double %36, double %36, double %31)
%38 = fadd double %30, %36
%39 = add i64 %33, %22
%40 = getelementptr inbounds double, double* %2, i64 %39
%41 = load double, double* %40
%42 = fsub double %41, %21
%43 = tail call double @llvm.fmuladd.f64(double %42, double %42, double %37)
%44 = fadd double %38, %42
%45 = add i32 %32, 2
%46 = icmp eq i32 %45, %27
br i1 %46, label %47, label %28
47: ; preds = %28, %14
%48 = phi double [ undef, %14 ], [ %43, %28 ]
%49 = phi double [ undef, %14 ], [ %44, %28 ]
%50 = phi i64 [ %19, %14 ], [ %39, %28 ]
%51 = phi double [ 0.000000e+00, %14 ], [ %44, %28 ]
%52 = phi double [ 0.000000e+00, %14 ], [ %43, %28 ]
%53 = icmp eq i32 %24, 0
br i1 %53, label %61, label %54
54: ; preds = %47
%55 = add i64 %50, %22
%56 = getelementptr inbounds double, double* %2, i64 %55
%57 = load double, double* %56
%58 = fsub double %57, %21
%59 = tail call double @llvm.fmuladd.f64(double %58, double %58, double %52)
%60 = fadd double %51, %58
br label %61
61: ; preds = %47, %54
%62 = phi double [ %48, %47 ], [ %59, %54 ]
%63 = phi double [ %49, %47 ], [ %60, %54 ]
%64 = fdiv double %62, %6
%65 = fdiv double %63, %6
%66 = fneg double %63
%67 = fdiv double %66, %5
%68 = tail call double @llvm.fmuladd.f64(double %67, double %65, double %64)
br label %69
69: ; preds = %10, %4, %61
%70 = phi double [ %68, %61 ], [ 0x7FF8000000000000, %4 ], [ 0.000000e+00, %10 ]
ret double %70
}