#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 my_mean1(const double data[], size_t stride, const size_t size) {
  long double mean = 0;
  size_t i;
  for (i = 0; i < size; i++) {
    mean += (data[i * stride] - mean) / (i + 1);
  }
  return mean;
}


x86 assembly (clang, -O0):
<my_mean1>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
fldz
fstpt  -0x30(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x18(%rbp),%rax
jae    117c <my_mean1+0x7c>
mov    -0x8(%rbp),%rcx
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rsi
mov    %rax,%rdx
imul   %rsi,%rdx
fldl   (%rcx,%rdx,8)
fldt   -0x30(%rbp)
fsubr  %st,%st(1)
inc    %rax
mov    %rax,-0x48(%rbp)
sets   %al
movzbl %al,%eax
mov    %eax,%ecx
fildll -0x48(%rbp)
lea    0xea1(%rip),%rax
flds   (%rax,%rcx,4)
faddp  %st,%st(1)
fdivrp %st,%st(2)
faddp  %st,%st(1)
fstpt  -0x30(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    111d <my_mean1+0x1d>
fldt   -0x30(%rbp)
fstpl  -0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<my_mean1>:
test   %rdx,%rdx
je     1113 <my_mean1+0x13>
push   %rbx
fldz
cmp    $0x1,%rdx
jne    1117 <my_mean1+0x17>
xor    %r11d,%r11d
jmp    1187 <my_mean1+0x87>
xorps  %xmm0,%xmm0
ret
fstp   %st(0)
mov    %rdx,%r8
and    $0xfffffffffffffffe,%r8
mov    %rsi,%r10
shl    $0x4,%r10
fldz
xor    %r11d,%r11d
lea    0xecd(%rip),%r9
mov    %rdi,%rcx
cs nopw 0x0(%rax,%rax,1)
fld    %st(0)
fsubrl (%rcx)
xor    %eax,%eax
mov    %r11,%rbx
inc    %rbx
mov    %rbx,-0x8(%rsp)
sets   %al
fildll -0x8(%rsp)
fadds  (%r9,%rax,4)
fdivrp %st,%st(1)
faddp  %st,%st(1)
fld    %st(0)
fsubrl (%rcx,%rsi,8)
xor    %eax,%eax
add    $0x2,%r11
mov    %r11,-0x10(%rsp)
sets   %al
fildll -0x10(%rsp)
fadds  (%r9,%rax,4)
fdivrp %st,%st(1)
faddp  %st,%st(1)
add    %r10,%rcx
cmp    %r11,%r8
jne    1140 <my_mean1+0x40>
test   $0x1,%dl
je     11b4 <my_mean1+0xb4>
imul   %r11,%rsi
fld    %st(0)
fsubrl (%rdi,%rsi,8)
xor    %eax,%eax
inc    %r11
mov    %r11,-0x18(%rsp)
sets   %al
fildll -0x18(%rsp)
lea    0xe53(%rip),%rcx
fadds  (%rcx,%rax,4)
fdivrp %st,%st(1)
faddp  %st,%st(1)
fstpl  -0x20(%rsp)
movsd  -0x20(%rsp),%xmm0
pop    %rbx
ret

LLVM IR (-O0):
define double @my_mean1(double* %0, i64 %1, i64 %2) {
%4 = alloca double*
%5 = alloca i64
%6 = alloca i64
%7 = alloca x86_fp80
%8 = alloca i64
store double* %0, double** %4
store i64 %1, i64* %5
store i64 %2, i64* %6
store x86_fp80 0xK00000000000000000000, x86_fp80* %7
store i64 0, i64* %8
br label %9
9: ; preds = %29, %3
%10 = load i64, i64* %8
%11 = load i64, i64* %6
%12 = icmp ult i64 %10, %11
br i1 %12, label %13, label %32
13: ; preds = %9
%14 = load double*, double** %4
%15 = load i64, i64* %8
%16 = load i64, i64* %5
%17 = mul i64 %15, %16
%18 = getelementptr inbounds double, double* %14, i64 %17
%19 = load double, double* %18
%20 = fpext double %19 to x86_fp80
%21 = load x86_fp80, x86_fp80* %7
%22 = fsub x86_fp80 %20, %21
%23 = load i64, i64* %8
%24 = add i64 %23, 1
%25 = uitofp i64 %24 to x86_fp80
%26 = fdiv x86_fp80 %22, %25
%27 = load x86_fp80, x86_fp80* %7
%28 = fadd x86_fp80 %27, %26
store x86_fp80 %28, x86_fp80* %7
br label %29
29: ; preds = %13
%30 = load i64, i64* %8
%31 = add i64 %30, 1
store i64 %31, i64* %8
br label %9
32: ; preds = %9
%33 = load x86_fp80, x86_fp80* %7
%34 = fptrunc x86_fp80 %33 to double
ret double %34
}

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