#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_snanmeanwd(const int64_t N, const float *X,
                                const int64_t stride) {
  int64_t ix;
  int64_t i;
  double n;
  float mu;
  float v;
  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++) {
    v = X[ix];
    if (v == v) {
      n += 1.0;
      mu += (float)((double)(v - mu) / n);
    }
    ix += stride;
  }
  if (n == 0.0) {
    return 0.0f / 0.0f;
  }
  return mu;
}


x86 assembly (clang, -O0):
<stdlib_strided_snanmeanwd>:
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_snanmeanwd+0x2d>
movss  0xedd(%rip),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1250 <stdlib_strided_snanmeanwd+0x150>
cmpq   $0x1,-0x10(%rbp)
je     1143 <stdlib_strided_snanmeanwd+0x43>
cmpq   $0x0,-0x20(%rbp)
jne    1155 <stdlib_strided_snanmeanwd+0x55>
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1250 <stdlib_strided_snanmeanwd+0x150>
cmpq   $0x0,-0x20(%rbp)
jge    1177 <stdlib_strided_snanmeanwd+0x77>
mov    $0x1,%eax
sub    -0x10(%rbp),%rax
imul   -0x20(%rbp),%rax
mov    %rax,-0x28(%rbp)
jmp    117f <stdlib_strided_snanmeanwd+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    121c <stdlib_strided_snanmeanwd+0x11c>
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
ucomiss -0x40(%rbp),%xmm0
jne    11ff <stdlib_strided_snanmeanwd+0xff>
jp     11ff <stdlib_strided_snanmeanwd+0xff>
movsd  0xe34(%rip),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movss  -0x40(%rbp),%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_snanmeanwd+0x97>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1246 <stdlib_strided_snanmeanwd+0x146>
jp     1246 <stdlib_strided_snanmeanwd+0x146>
movss  0xdc4(%rip),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1250 <stdlib_strided_snanmeanwd+0x150>
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<stdlib_strided_snanmeanwd>:
test   %rdi,%rdi
jle    1176 <stdlib_strided_snanmeanwd+0x76>
cmp    $0x1,%rdi
je     117f <stdlib_strided_snanmeanwd+0x7f>
test   %rdx,%rdx
je     117f <stdlib_strided_snanmeanwd+0x7f>
mov    $0x1,%eax
sub    %rdi,%rax
imul   %rdx,%rax
mov    %rdx,%rcx
sar    $0x3f,%rcx
and    %rax,%rcx
lea    (%rsi,%rcx,4),%rax
shl    $0x2,%rdx
xorpd  %xmm1,%xmm1
xorps  %xmm0,%xmm0
movsd  0xecb(%rip),%xmm2
nopl   (%rax)
movss  (%rax),%xmm3
ucomiss %xmm3,%xmm3
jp     1161 <stdlib_strided_snanmeanwd+0x61>
addsd  %xmm2,%xmm1
subss  %xmm0,%xmm3
cvtss2sd %xmm3,%xmm3
divsd  %xmm1,%xmm3
cvtsd2ss %xmm3,%xmm3
addss  %xmm3,%xmm0
add    %rdx,%rax
add    $0xffffffffffffffff,%rdi
jne    1140 <stdlib_strided_snanmeanwd+0x40>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    117e <stdlib_strided_snanmeanwd+0x7e>
jp     117e <stdlib_strided_snanmeanwd+0x7e>
movss  0xe82(%rip),%xmm0
ret
movss  (%rsi),%xmm0
ret

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

LLVM IR (-O2):
define float @stdlib_strided_snanmeanwd(i64 %0, float* nocapture %1, i64 %2) {
%4 = icmp slt i64 %0, 1
br i1 %4, label %40, 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 %40
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
br label %16
16: ; preds = %11, %31
%17 = phi float [ %33, %31 ], [ 0.000000e+00, %11 ]
%18 = phi double [ %32, %31 ], [ 0.000000e+00, %11 ]
%19 = phi i64 [ %35, %31 ], [ 0, %11 ]
%20 = phi i64 [ %34, %31 ], [ %15, %11 ]
%21 = getelementptr inbounds float, float* %1, i64 %20
%22 = load float, float* %21
%23 = fcmp ord float %22, 0.000000e+00
br i1 %23, label %24, label %31
24: ; preds = %16
%25 = fadd double %18, 1.000000e+00
%26 = fsub float %22, %17
%27 = fpext float %26 to double
%28 = fdiv double %27, %25
%29 = fptrunc double %28 to float
%30 = fadd float %17, %29
br label %31
31: ; preds = %24, %16
%32 = phi double [ %25, %24 ], [ %18, %16 ]
%33 = phi float [ %30, %24 ], [ %17, %16 ]
%34 = add nsw i64 %20, %2
%35 = add nuw nsw i64 %19, 1
%36 = icmp eq i64 %35, %0
br i1 %36, label %37, label %16
37: ; preds = %31
%38 = fcmp oeq double %32, 0.000000e+00
%39 = select i1 %38, float 0x7FF8000000000000, float %33
br label %40
40: ; preds = %37, %3, %9
%41 = phi float [ %10, %9 ], [ 0x7FF8000000000000, %3 ], [ %39, %37 ]
ret float %41
}