void linearMean(int n, int p, int j, double *data, double *mean) {
  double temp;
  int i;
  temp = 0.0;
  for (i = 0; i < n; i++) {
    temp += data[i * p + j];
  }
  *mean = temp / n;
}


x86 assembly (clang, -O0):
<linearMean>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    115e <linearMean+0x5e>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1124 <linearMean+0x24>
movsd  -0x28(%rbp),%xmm0
cvtsi2sdl -0x4(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<linearMean>:
push   %r15
push   %r14
push   %r12
push   %rbx
test   %edi,%edi
jle    112c <linearMean+0x2c>
movslq %esi,%r10
movslq %edx,%r9
mov    %edi,%r11d
lea    -0x1(%r11),%rax
mov    %r11d,%edx
and    $0x3,%edx
cmp    $0x3,%rax
jae    1135 <linearMean+0x35>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    1193 <linearMean+0x93>
xorpd  %xmm0,%xmm0
jmp    11bd <linearMean+0xbd>
and    $0xfffffffc,%r11d
mov    %r10,%r14
shl    $0x4,%r14
add    %rcx,%r14
lea    0x0(,%r9,8),%rbx
mov    %r10,%rsi
shl    $0x5,%rsi
lea    (%rcx,%r10,8),%r15
lea    (%r10,%r10,2),%rax
lea    (%rcx,%rax,8),%r12
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addsd  (%rcx,%rbx,1),%xmm0
addsd  (%r15,%rbx,1),%xmm0
addsd  (%r14,%rbx,1),%xmm0
addsd  (%r12,%rbx,1),%xmm0
add    $0x4,%rax
add    %rsi,%rbx
cmp    %rax,%r11
jne    1170 <linearMean+0x70>
test   %rdx,%rdx
je     11bd <linearMean+0xbd>
imul   %r10,%rax
add    %r9,%rax
lea    (%rcx,%rax,8),%rax
shl    $0x3,%r10
nopw   0x0(%rax,%rax,1)
addsd  (%rax),%xmm0
add    %r10,%rax
add    $0xffffffffffffffff,%rdx
jne    11b0 <linearMean+0xb0>
cvtsi2sd %edi,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%r8)
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @linearMean(i32 %0, i32 %1, i32 %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca double*
%11 = alloca double
%12 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store double* %3, double** %9
store double* %4, double** %10
store double 0.000000e+00, double* %11
store i32 0, i32* %12
br label %13
13: ; preds = %29, %5
%14 = load i32, i32* %12
%15 = load i32, i32* %6
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %32
17: ; preds = %13
%18 = load double*, double** %9
%19 = load i32, i32* %12
%20 = load i32, i32* %7
%21 = mul nsw i32 %19, %20
%22 = load i32, i32* %8
%23 = add nsw i32 %21, %22
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %18, i64 %24
%26 = load double, double* %25
%27 = load double, double* %11
%28 = fadd double %27, %26
store double %28, double* %11
br label %29
29: ; preds = %17
%30 = load i32, i32* %12
%31 = add nsw i32 %30, 1
store i32 %31, i32* %12
br label %13
32: ; preds = %13
%33 = load double, double* %11
%34 = load i32, i32* %6
%35 = sitofp i32 %34 to double
%36 = fdiv double %33, %35
%37 = load double*, double** %10
store double %36, double* %37
ret void
}

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