#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 scalar_t__ int64_t;

/* Variables and functions */
float stdlib_strided_smeanwd(const int64_t N, const float *X,
                             const int64_t stride) {
  int64_t ix;
  int64_t i;
  double n;
  float mu;
  if (N <= 0) {
    return 0.0f / 0.0f;
  }
  if (N == 1 || stride == 0) {
    return X[0];
  }
  if (stride < 0) {
    ix = (1 - N) * stride;
  } else {
    ix = 0;
  }
  mu = 0.0f;
  n = 0.0;
  for (i = 0; i < N; i++) {
    n += 1.0;
    mu += (float)((double)(X[ix] - mu) / n);
    ix += stride;
  }
  return mu;
}


x86 assembly (clang, -O0):
<stdlib_strided_smeanwd>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
cmpq   $0x0,-0x10(%rbp)
jg     112d <stdlib_strided_smeanwd+0x2d>
movss  0xee5(%rip),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1207 <stdlib_strided_smeanwd+0x107>
cmpq   $0x1,-0x10(%rbp)
je     1143 <stdlib_strided_smeanwd+0x43>
cmpq   $0x0,-0x20(%rbp)
jne    1155 <stdlib_strided_smeanwd+0x55>
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1207 <stdlib_strided_smeanwd+0x107>
cmpq   $0x0,-0x20(%rbp)
jge    1177 <stdlib_strided_smeanwd+0x77>
mov    $0x1,%eax
sub    -0x10(%rbp),%rax
imul   -0x20(%rbp),%rax
mov    %rax,-0x28(%rbp)
jmp    117f <stdlib_strided_smeanwd+0x7f>
movq   $0x0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11fd <stdlib_strided_smeanwd+0xfd>
movsd  0xe53(%rip),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x3c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
divsd  -0x38(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
addss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x20(%rbp),%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1197 <stdlib_strided_smeanwd+0x97>
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<stdlib_strided_smeanwd>:
test   %rdi,%rdi
jle    11d0 <stdlib_strided_smeanwd+0xd0>
cmp    $0x1,%rdi
je     11d9 <stdlib_strided_smeanwd+0xd9>
test   %rdx,%rdx
je     11d9 <stdlib_strided_smeanwd+0xd9>
mov    $0x1,%eax
sub    %rdi,%rax
imul   %rdx,%rax
mov    %rdx,%r8
sar    $0x3f,%r8
and    %rax,%r8
mov    %rdi,%rcx
and    $0xfffffffffffffffe,%rcx
lea    0x0(,%r8,4),%r9
lea    0x0(,%rdx,8),%rax
shl    $0x2,%rdx
lea    (%rdx,%r8,4),%rdx
xorpd  %xmm1,%xmm1
xorps  %xmm0,%xmm0
movsd  0xea8(%rip),%xmm2
addsd  %xmm2,%xmm1
movss  (%rsi,%r9,1),%xmm3
subss  %xmm0,%xmm3
cvtss2sd %xmm3,%xmm3
divsd  %xmm1,%xmm3
cvtsd2ss %xmm3,%xmm3
addss  %xmm0,%xmm3
addsd  %xmm2,%xmm1
movss  (%rsi,%rdx,1),%xmm0
subss  %xmm3,%xmm0
cvtss2sd %xmm0,%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
addss  %xmm3,%xmm0
add    %rax,%rsi
add    $0xfffffffffffffffe,%rcx
jne    1160 <stdlib_strided_smeanwd+0x60>
test   $0x1,%dil
je     11cf <stdlib_strided_smeanwd+0xcf>
addsd  0xe56(%rip),%xmm1
movss  (%rsi,%r8,4),%xmm2
subss  %xmm0,%xmm2
cvtss2sd %xmm2,%xmm2
divsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsd2ss %xmm2,%xmm1
addss  %xmm1,%xmm0
ret
movss  0xe28(%rip),%xmm0
ret
movss  (%rsi),%xmm0
ret

LLVM IR (-O0):
define float @stdlib_strided_smeanwd(i64 %0, float* %1, i64 %2) {
%4 = alloca float
%5 = alloca i64
%6 = alloca float*
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca double
%11 = alloca float
store i64 %0, i64* %5
store float* %1, float** %6
store i64 %2, i64* %7
%12 = load i64, i64* %5
%13 = icmp sle i64 %12, 0
br i1 %13, label %14, label %15
14: ; preds = %3
store float 0x7FF8000000000000, float* %4
br label %62
15: ; preds = %3
%16 = load i64, i64* %5
%17 = icmp eq i64 %16, 1
br i1 %17, label %21, label %18
18: ; preds = %15
%19 = load i64, i64* %7
%20 = icmp eq i64 %19, 0
br i1 %20, label %21, label %25
21: ; preds = %18, %15
%22 = load float*, float** %6
%23 = getelementptr inbounds float, float* %22, i64 0
%24 = load float, float* %23
store float %24, float* %4
br label %62
25: ; preds = %18
%26 = load i64, i64* %7
%27 = icmp slt i64 %26, 0
br i1 %27, label %28, label %33
28: ; preds = %25
%29 = load i64, i64* %5
%30 = sub nsw i64 1, %29
%31 = load i64, i64* %7
%32 = mul nsw i64 %30, %31
store i64 %32, i64* %8
br label %34
33: ; preds = %25
store i64 0, i64* %8
br label %34
34: ; preds = %33, %28
store float 0.000000e+00, float* %11
store double 0.000000e+00, double* %10
store i64 0, i64* %9
br label %35
35: ; preds = %57, %34
%36 = load i64, i64* %9
%37 = load i64, i64* %5
%38 = icmp slt i64 %36, %37
br i1 %38, label %39, label %60
39: ; preds = %35
%40 = load double, double* %10
%41 = fadd double %40, 1.000000e+00
store double %41, double* %10
%42 = load float*, float** %6
%43 = load i64, i64* %8
%44 = getelementptr inbounds float, float* %42, i64 %43
%45 = load float, float* %44
%46 = load float, float* %11
%47 = fsub float %45, %46
%48 = fpext float %47 to double
%49 = load double, double* %10
%50 = fdiv double %48, %49
%51 = fptrunc double %50 to float
%52 = load float, float* %11
%53 = fadd float %52, %51
store float %53, float* %11
%54 = load i64, i64* %7
%55 = load i64, i64* %8
%56 = add nsw i64 %55, %54
store i64 %56, i64* %8
br label %57
57: ; preds = %39
%58 = load i64, i64* %9
%59 = add nsw i64 %58, 1
store i64 %59, i64* %9
br label %35
60: ; preds = %35
%61 = load float, float* %11
store float %61, float* %4
br label %62
62: ; preds = %60, %21, %14
%63 = load float, float* %4
ret float %63
}

LLVM IR (-O2):
define float @stdlib_strided_smeanwd(i64 %0, float* nocapture %1, i64 %2) {
%4 = icmp slt i64 %0, 1
br i1 %4, label %54, label %5
5: ; preds = %3
%6 = icmp eq i64 %0, 1
%7 = icmp eq i64 %2, 0
%8 = or i1 %6, %7
br i1 %8, label %9, label %11
9: ; preds = %5
%10 = load float, float* %1
br label %54
11: ; preds = %5
%12 = icmp slt i64 %2, 0
%13 = sub nsw i64 1, %0
%14 = mul nsw i64 %13, %2
%15 = select i1 %12, i64 %14, i64 0
%16 = and i64 %0, 1
%17 = and i64 %0, -2
br label %18
18: ; preds = %18, %11
%19 = phi float [ 0.000000e+00, %11 ], [ %39, %18 ]
%20 = phi double [ 0.000000e+00, %11 ], [ %32, %18 ]
%21 = phi i64 [ %15, %11 ], [ %40, %18 ]
%22 = phi i64 [ 0, %11 ], [ %41, %18 ]
%23 = fadd double %20, 1.000000e+00
%24 = getelementptr inbounds float, float* %1, i64 %21
%25 = load float, float* %24
%26 = fsub float %25, %19
%27 = fpext float %26 to double
%28 = fdiv double %27, %23
%29 = fptrunc double %28 to float
%30 = fadd float %19, %29
%31 = add nsw i64 %21, %2
%32 = fadd double %23, 1.000000e+00
%33 = getelementptr inbounds float, float* %1, i64 %31
%34 = load float, float* %33
%35 = fsub float %34, %30
%36 = fpext float %35 to double
%37 = fdiv double %36, %32
%38 = fptrunc double %37 to float
%39 = fadd float %30, %38
%40 = add nsw i64 %31, %2
%41 = add i64 %22, 2
%42 = icmp eq i64 %41, %17
br i1 %42, label %43, label %18
43: ; preds = %18
%44 = fadd double %32, 1.000000e+00
%45 = icmp eq i64 %16, 0
br i1 %45, label %54, label %46
46: ; preds = %43
%47 = getelementptr inbounds float, float* %1, i64 %40
%48 = load float, float* %47
%49 = fsub float %48, %39
%50 = fpext float %49 to double
%51 = fdiv double %50, %44
%52 = fptrunc double %51 to float
%53 = fadd float %39, %52
br label %54
54: ; preds = %46, %43, %3, %9
%55 = phi float [ %10, %9 ], [ 0x7FF8000000000000, %3 ], [ %39, %43 ], [ %53, %46 ]
ret float %55
}