void calcmean(float *matrix, float *mean) {
  int i, j;
  float sum;
#pragma omp parallel for
  for (i = 0; i < 20530; i++) {
    sum = 0.0;
    for (j = 0; j < 4000; j++) {
      sum += matrix[i * 4000 + j];
    }
    mean[i] = sum / (float)4000;
  }
}


x86 assembly (clang, -O0):
<calcmean>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
cmpl   $0x5032,-0x14(%rbp)
jge    1196 <calcmean+0x96>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x0,-0x18(%rbp)
cmpl   $0xfa0,-0x18(%rbp)
jge    116a <calcmean+0x6a>
mov    -0x8(%rbp),%rax
imul   $0xfa0,-0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    112f <calcmean+0x2f>
movss  -0x1c(%rbp),%xmm0
movss  0xe89(%rip),%xmm1
divss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1113 <calcmean+0x13>
pop    %rbp
ret

x86 assembly (clang, -O2):
<calcmean>:
xor    %eax,%eax
movss  0xef6(%rip),%xmm0
nopw   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
mov    $0x4,%ecx
nopl   0x0(%rax,%rax,1)
addss  -0x10(%rdi,%rcx,4),%xmm1
addss  -0xc(%rdi,%rcx,4),%xmm1
addss  -0x8(%rdi,%rcx,4),%xmm1
addss  -0x4(%rdi,%rcx,4),%xmm1
addss  (%rdi,%rcx,4),%xmm1
add    $0x5,%rcx
cmp    $0xfa4,%rcx
jne    1120 <calcmean+0x20>
divss  %xmm0,%xmm1
movss  %xmm1,(%rsi,%rax,4)
add    $0x1,%rax
add    $0x3e80,%rdi
cmp    $0x5032,%rax
jne    1110 <calcmean+0x10>
ret

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

LLVM IR (-O2):
define void @calcmean(float* nocapture %0, float* nocapture writeonly %1) {
br label %3
3: ; preds = %2, %35
%4 = phi i64 [ 0, %2 ], [ %38, %35 ]
%5 = mul nuw nsw i64 %4, 4000
br label %6
6: ; preds = %6, %3
%7 = phi i64 [ 0, %3 ], [ %33, %6 ]
%8 = phi float [ 0.000000e+00, %3 ], [ %32, %6 ]
%9 = add nuw nsw i64 %7, %5
%10 = getelementptr inbounds float, float* %0, i64 %9
%11 = load float, float* %10
%12 = fadd float %8, %11
%13 = add nuw nsw i64 %7, 1
%14 = add nuw nsw i64 %13, %5
%15 = getelementptr inbounds float, float* %0, i64 %14
%16 = load float, float* %15
%17 = fadd float %12, %16
%18 = add nuw nsw i64 %7, 2
%19 = add nuw nsw i64 %18, %5
%20 = getelementptr inbounds float, float* %0, i64 %19
%21 = load float, float* %20
%22 = fadd float %17, %21
%23 = add nuw nsw i64 %7, 3
%24 = add nuw nsw i64 %23, %5
%25 = getelementptr inbounds float, float* %0, i64 %24
%26 = load float, float* %25
%27 = fadd float %22, %26
%28 = add nuw nsw i64 %7, 4
%29 = add nuw nsw i64 %28, %5
%30 = getelementptr inbounds float, float* %0, i64 %29
%31 = load float, float* %30
%32 = fadd float %27, %31
%33 = add nuw nsw i64 %7, 5
%34 = icmp eq i64 %33, 4000
br i1 %34, label %35, label %6
35: ; preds = %6
%36 = fdiv float %32, 4.000000e+03
%37 = getelementptr inbounds float, float* %1, i64 %4
store float %36, float* %37
%38 = add nuw nsw i64 %4, 1
%39 = icmp eq i64 %38, 20530
br i1 %39, label %40, label %3
40: ; preds = %35
ret void
}