void get_item_average_matrix(float *rating, double *average_matrix,
                             int start_index) {
  int m, n, i, j;
  int average_index = 0;
  int average_sum = 0;
  int block_m = 0;
  double average_item = 0;
#pragma omp parallel for private(m, n) reduction(+ : average_sum, average_index)
  for (m = start_index;
       m < (1019318 < (start_index + 3500) ? 1019318 : (start_index + 3500));
       m++) {
    average_sum = 0;
    average_index = 0;
    block_m = m - start_index;
    for (n = 0; n < 384546; n++) {
      if (rating[block_m * 384546 + n] != 0) {
        average_sum += rating[block_m * 384546 + n];
        average_index += 1;
      }
    }
    if (average_index != 0) {
      average_matrix[m] = (double)average_sum / (double)average_index;
    } else {
      average_matrix[m] = 0;
    }
  }
#pragma omp parallel for private(m, n)
  for (m = start_index;
       m < (1019318 < (start_index + 3500) ? 1019318 : (start_index + 3500));
       m++) {
    block_m = m - start_index;
    average_item = average_matrix[m];
    for (n = 0; n < 384546; n++) {
      rating[block_m * 384546 + n] -= (float)average_item;
    }
  }
}


x86 assembly (clang, -O0):
<get_item_average_matrix>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x14(%rbp),%ecx
add    $0xdac,%ecx
mov    $0xf8db6,%eax
cmp    %ecx,%eax
jge    115b <get_item_average_matrix+0x5b>
mov    $0xf8db6,%eax
mov    %eax,-0x40(%rbp)
jmp    1166 <get_item_average_matrix+0x66>
mov    -0x14(%rbp),%eax
add    $0xdac,%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x40(%rbp),%ecx
cmp    %ecx,%eax
jge    125b <get_item_average_matrix+0x15b>
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x18(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x30(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x5de22,-0x1c(%rbp)
jge    120e <get_item_average_matrix+0x10e>
mov    -0x8(%rbp),%rax
imul   $0x5de22,-0x30(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    11cc <get_item_average_matrix+0xcc>
jp     11cc <get_item_average_matrix+0xcc>
jmp    11fb <get_item_average_matrix+0xfb>
mov    -0x8(%rbp),%rax
imul   $0x5de22,-0x30(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtsi2ssl -0x2c(%rbp),%xmm0
addss  %xmm1,%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1200 <get_item_average_matrix+0x100>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1192 <get_item_average_matrix+0x92>
cmpl   $0x0,-0x28(%rbp)
je     1238 <get_item_average_matrix+0x138>
cvtsi2sdl -0x2c(%rbp),%xmm0
cvtsi2sdl -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1248 <get_item_average_matrix+0x148>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    124d <get_item_average_matrix+0x14d>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1132 <get_item_average_matrix+0x32>
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x14(%rbp),%ecx
add    $0xdac,%ecx
mov    $0xf8db6,%eax
cmp    %ecx,%eax
jge    128a <get_item_average_matrix+0x18a>
mov    $0xf8db6,%eax
mov    %eax,-0x48(%rbp)
jmp    1295 <get_item_average_matrix+0x195>
mov    -0x14(%rbp),%eax
add    $0xdac,%eax
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x48(%rbp),%ecx
cmp    %ecx,%eax
jge    131b <get_item_average_matrix+0x21b>
mov    -0x18(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x5de22,-0x1c(%rbp)
jge    1308 <get_item_average_matrix+0x208>
movsd  -0x38(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm1
mov    -0x8(%rbp),%rax
imul   $0x5de22,-0x30(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    12c5 <get_item_average_matrix+0x1c5>
jmp    130d <get_item_average_matrix+0x20d>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1261 <get_item_average_matrix+0x161>
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_item_average_matrix>:
push   %rbp
push   %rbx
cmp    $0xf800a,%edx
mov    $0xf800a,%r8d
cmovl  %edx,%r8d
add    $0xdac,%r8d
cmp    %edx,%r8d
jle    12ae <get_item_average_matrix+0x1ae>
movslq %edx,%r10
movslq %r8d,%r9
xorps  %xmm0,%xmm0
mov    %rdi,%r11
jmp    115c <get_item_average_matrix+0x5c>
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,(%rsi,%r10,8)
add    $0x1,%r10
add    $0x177888,%r11
cmp    %r9,%r10
je     11e6 <get_item_average_matrix+0xe6>
xor    %ecx,%ecx
mov    $0x1,%ebx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movss  -0x4(%r11,%rbx,4),%xmm1
movss  (%r11,%rbx,4),%xmm2
xorps  %xmm3,%xmm3
cmpneqss %xmm1,%xmm3
movd   %xmm3,%ebp
sub    %ebp,%eax
ucomiss %xmm1,%xmm0
xorps  %xmm3,%xmm3
cvtsi2ss %ecx,%xmm3
addss  %xmm1,%xmm3
cvttss2si %xmm3,%ebp
cmovne %ebp,%ecx
cmovp  %ebp,%ecx
xorps  %xmm1,%xmm1
cmpneqss %xmm2,%xmm1
movd   %xmm1,%ebp
sub    %ebp,%eax
ucomiss %xmm2,%xmm0
xorps  %xmm1,%xmm1
cvtsi2ss %ecx,%xmm1
addss  %xmm2,%xmm1
cvttss2si %xmm1,%ebp
cmovne %ebp,%ecx
cmovp  %ebp,%ecx
add    $0x2,%rbx
cmp    $0x5de23,%rbx
jne    1170 <get_item_average_matrix+0x70>
test   %eax,%eax
jne    1130 <get_item_average_matrix+0x30>
xorps  %xmm1,%xmm1
jmp    1142 <get_item_average_matrix+0x42>
cmp    %edx,%r8d
jle    12ae <get_item_average_matrix+0x1ae>
movslq %edx,%r8
mov    %rdi,%rcx
mov    %r8,%rdx
nopl   0x0(%rax,%rax,1)
mov    %rdx,%rax
movsd  (%rsi,%rdx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
sub    %r8,%rax
imul   $0x5de22,%rax,%rbx
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
mov    $0xc,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups -0x30(%rcx,%rax,4),%xmm2
movups -0x20(%rcx,%rax,4),%xmm3
movups -0x10(%rcx,%rax,4),%xmm4
movups (%rcx,%rax,4),%xmm5
subps  %xmm1,%xmm2
subps  %xmm1,%xmm3
movups %xmm2,-0x30(%rcx,%rax,4)
movups %xmm3,-0x20(%rcx,%rax,4)
subps  %xmm1,%xmm4
subps  %xmm1,%xmm5
movups %xmm4,-0x10(%rcx,%rax,4)
movups %xmm5,(%rcx,%rax,4)
add    $0x10,%rax
cmp    $0x5de2c,%rax
jne    1230 <get_item_average_matrix+0x130>
movss  0x177880(%rdi,%rbx,4),%xmm1
movss  0x177884(%rdi,%rbx,4),%xmm2
subss  %xmm0,%xmm1
movss  %xmm1,0x177880(%rdi,%rbx,4)
subss  %xmm0,%xmm2
movss  %xmm2,0x177884(%rdi,%rbx,4)
add    $0x1,%rdx
add    $0x177888,%rcx
cmp    %r9,%rdx
jne    1200 <get_item_average_matrix+0x100>
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define void @get_item_average_matrix(float* %0, double* %1, i32 %2) {
%4 = alloca float*
%5 = alloca double*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
store float* %0, float** %4
store double* %1, double** %5
store i32 %2, i32* %6
store i32 0, i32* %11
store i32 0, i32* %12
store i32 0, i32* %13
store double 0.000000e+00, double* %14
%15 = load i32, i32* %6
store i32 %15, i32* %7
br label %16
16: ; preds = %83, %3
%17 = load i32, i32* %7
%18 = load i32, i32* %6
%19 = add nsw i32 %18, 3500
%20 = icmp slt i32 1019318, %19
br i1 %20, label %21, label %22
21: ; preds = %16
br label %25
22: ; preds = %16
%23 = load i32, i32* %6
%24 = add nsw i32 %23, 3500
br label %25
25: ; preds = %22, %21
%26 = phi i32 [ 1019318, %21 ], [ %24, %22 ]
%27 = icmp slt i32 %17, %26
br i1 %27, label %28, label %86
28: ; preds = %25
store i32 0, i32* %12
store i32 0, i32* %11
%29 = load i32, i32* %7
%30 = load i32, i32* %6
%31 = sub nsw i32 %29, %30
store i32 %31, i32* %13
store i32 0, i32* %8
br label %32
32: ; preds = %61, %28
%33 = load i32, i32* %8
%34 = icmp slt i32 %33, 384546
br i1 %34, label %35, label %64
35: ; preds = %32
%36 = load float*, float** %4
%37 = load i32, i32* %13
%38 = mul nsw i32 %37, 384546
%39 = load i32, i32* %8
%40 = add nsw i32 %38, %39
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float, float* %36, i64 %41
%43 = load float, float* %42
%44 = fcmp une float %43, 0.000000e+00
br i1 %44, label %45, label %60
45: ; preds = %35
%46 = load float*, float** %4
%47 = load i32, i32* %13
%48 = mul nsw i32 %47, 384546
%49 = load i32, i32* %8
%50 = add nsw i32 %48, %49
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %46, i64 %51
%53 = load float, float* %52
%54 = load i32, i32* %12
%55 = sitofp i32 %54 to float
%56 = fadd float %55, %53
%57 = fptosi float %56 to i32
store i32 %57, i32* %12
%58 = load i32, i32* %11
%59 = add nsw i32 %58, 1
store i32 %59, i32* %11
br label %60
60: ; preds = %45, %35
br label %61
61: ; preds = %60
%62 = load i32, i32* %8
%63 = add nsw i32 %62, 1
store i32 %63, i32* %8
br label %32
64: ; preds = %32
%65 = load i32, i32* %11
%66 = icmp ne i32 %65, 0
br i1 %66, label %67, label %77
67: ; preds = %64
%68 = load i32, i32* %12
%69 = sitofp i32 %68 to double
%70 = load i32, i32* %11
%71 = sitofp i32 %70 to double
%72 = fdiv double %69, %71
%73 = load double*, double** %5
%74 = load i32, i32* %7
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %73, i64 %75
store double %72, double* %76
br label %82
77: ; preds = %64
%78 = load double*, double** %5
%79 = load i32, i32* %7
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %78, i64 %80
store double 0.000000e+00, double* %81
br label %82
82: ; preds = %77, %67
br label %83
83: ; preds = %82
%84 = load i32, i32* %7
%85 = add nsw i32 %84, 1
store i32 %85, i32* %7
br label %16
86: ; preds = %25
%87 = load i32, i32* %6
store i32 %87, i32* %7
br label %88
88: ; preds = %128, %86
%89 = load i32, i32* %7
%90 = load i32, i32* %6
%91 = add nsw i32 %90, 3500
%92 = icmp slt i32 1019318, %91
br i1 %92, label %93, label %94
93: ; preds = %88
br label %97
94: ; preds = %88
%95 = load i32, i32* %6
%96 = add nsw i32 %95, 3500
br label %97
97: ; preds = %94, %93
%98 = phi i32 [ 1019318, %93 ], [ %96, %94 ]
%99 = icmp slt i32 %89, %98
br i1 %99, label %100, label %131
100: ; preds = %97
%101 = load i32, i32* %7
%102 = load i32, i32* %6
%103 = sub nsw i32 %101, %102
store i32 %103, i32* %13
%104 = load double*, double** %5
%105 = load i32, i32* %7
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds double, double* %104, i64 %106
%108 = load double, double* %107
store double %108, double* %14
store i32 0, i32* %8
br label %109
109: ; preds = %124, %100
%110 = load i32, i32* %8
%111 = icmp slt i32 %110, 384546
br i1 %111, label %112, label %127
112: ; preds = %109
%113 = load double, double* %14
%114 = fptrunc double %113 to float
%115 = load float*, float** %4
%116 = load i32, i32* %13
%117 = mul nsw i32 %116, 384546
%118 = load i32, i32* %8
%119 = add nsw i32 %117, %118
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds float, float* %115, i64 %120
%122 = load float, float* %121
%123 = fsub float %122, %114
store float %123, float* %121
br label %124
124: ; preds = %112
%125 = load i32, i32* %8
%126 = add nsw i32 %125, 1
store i32 %126, i32* %8
br label %109
127: ; preds = %109
br label %128
128: ; preds = %127
%129 = load i32, i32* %7
%130 = add nsw i32 %129, 1
store i32 %130, i32* %7
br label %88
131: ; preds = %97
ret void
}

LLVM IR (-O2):
define void @get_item_average_matrix(float* nocapture %0, double* nocapture %1, i32 %2) {
%4 = icmp slt i32 %2, 1015818
%5 = select i1 %4, i32 %2, i32 1015818
%6 = add nsw i32 %5, 3500
%7 = icmp sgt i32 %6, %2
br i1 %7, label %8, label %106
8: ; preds = %3
%9 = sext i32 %2 to i64
%10 = sext i32 %6 to i64
br label %15
11: ; preds = %52
br i1 %7, label %12, label %106
12: ; preds = %11
%13 = sext i32 %2 to i64
%14 = sext i32 %6 to i64
br label %57
15: ; preds = %8, %52
%16 = phi i64 [ %9, %8 ], [ %55, %52 ]
%17 = sub nsw i64 %16, %9
%18 = mul nuw nsw i64 %17, 384546
br label %19
19: ; preds = %19, %15
%20 = phi i64 [ 0, %15 ], [ %44, %19 ]
%21 = phi i32 [ 0, %15 ], [ %43, %19 ]
%22 = phi i32 [ 0, %15 ], [ %42, %19 ]
%23 = add nuw nsw i64 %20, %18
%24 = getelementptr inbounds float, float* %0, i64 %23
%25 = load float, float* %24
%26 = fcmp une float %25, 0.000000e+00
%27 = sitofp i32 %21 to float
%28 = fadd float %25, %27
%29 = fptosi float %28 to i32
%30 = zext i1 %26 to i32
%31 = add nuw nsw i32 %22, %30
%32 = select i1 %26, i32 %29, i32 %21
%33 = or i64 %20, 1
%34 = add nuw nsw i64 %33, %18
%35 = getelementptr inbounds float, float* %0, i64 %34
%36 = load float, float* %35
%37 = fcmp une float %36, 0.000000e+00
%38 = sitofp i32 %32 to float
%39 = fadd float %36, %38
%40 = fptosi float %39 to i32
%41 = zext i1 %37 to i32
%42 = add nuw nsw i32 %31, %41
%43 = select i1 %37, i32 %40, i32 %32
%44 = add nuw nsw i64 %20, 2
%45 = icmp eq i64 %44, 384546
br i1 %45, label %46, label %19
46: ; preds = %19
%47 = icmp eq i32 %42, 0
br i1 %47, label %52, label %48
48: ; preds = %46
%49 = sitofp i32 %43 to double
%50 = sitofp i32 %42 to double
%51 = fdiv double %49, %50
br label %52
52: ; preds = %46, %48
%53 = phi double [ %51, %48 ], [ 0.000000e+00, %46 ]
%54 = getelementptr inbounds double, double* %1, i64 %16
store double %53, double* %54
%55 = add nsw i64 %16, 1
%56 = icmp eq i64 %55, %10
br i1 %56, label %11, label %15
57: ; preds = %12, %95
%58 = phi i64 [ %13, %12 ], [ %104, %95 ]
%59 = sub nsw i64 %58, %13
%60 = getelementptr inbounds double, double* %1, i64 %58
%61 = load double, double* %60
%62 = fptrunc double %61 to float
%63 = mul nuw nsw i64 %59, 384546
%64 = insertelement <4 x float> poison, float %62, i64 0
%65 = shufflevector <4 x float> %64, <4 x float> poison, <4 x i32> zeroinitializer
%66 = insertelement <4 x float> poison, float %62, i64 0
%67 = shufflevector <4 x float> %66, <4 x float> poison, <4 x i32> zeroinitializer
br label %68
68: ; preds = %68, %57
%69 = phi i64 [ 0, %57 ], [ %93, %68 ]
%70 = add nuw nsw i64 %69, %63
%71 = getelementptr inbounds float, float* %0, i64 %70
%72 = bitcast float* %71 to <4 x float>*
%73 = load <4 x float>, <4 x float>* %72
%74 = getelementptr inbounds float, float* %71, i64 4
%75 = bitcast float* %74 to <4 x float>*
%76 = load <4 x float>, <4 x float>* %75
%77 = fsub <4 x float> %73, %65
%78 = fsub <4 x float> %76, %67
%79 = bitcast float* %71 to <4 x float>*
store <4 x float> %77, <4 x float>* %79
%80 = bitcast float* %74 to <4 x float>*
store <4 x float> %78, <4 x float>* %80
%81 = or i64 %69, 8
%82 = add nuw nsw i64 %81, %63
%83 = getelementptr inbounds float, float* %0, i64 %82
%84 = bitcast float* %83 to <4 x float>*
%85 = load <4 x float>, <4 x float>* %84
%86 = getelementptr inbounds float, float* %83, i64 4
%87 = bitcast float* %86 to <4 x float>*
%88 = load <4 x float>, <4 x float>* %87
%89 = fsub <4 x float> %85, %65
%90 = fsub <4 x float> %88, %67
%91 = bitcast float* %83 to <4 x float>*
store <4 x float> %89, <4 x float>* %91
%92 = bitcast float* %86 to <4 x float>*
store <4 x float> %90, <4 x float>* %92
%93 = add nuw nsw i64 %69, 16
%94 = icmp eq i64 %93, 384544
br i1 %94, label %95, label %68
95: ; preds = %68
%96 = add nuw nsw i64 %63, 384544
%97 = getelementptr inbounds float, float* %0, i64 %96
%98 = load float, float* %97
%99 = fsub float %98, %62
store float %99, float* %97
%100 = add nuw nsw i64 %63, 384545
%101 = getelementptr inbounds float, float* %0, i64 %100
%102 = load float, float* %101
%103 = fsub float %102, %62
store float %103, float* %101
%104 = add nsw i64 %58, 1
%105 = icmp eq i64 %104, %14
br i1 %105, label %106, label %57
106: ; preds = %95, %3, %11
ret void
}