void vertical_sum_cpu(float *m, float *s) {
  int i, j;
  for (j = 0; j < 2000; j++) {
    for (i = 0; i < 800; i++) {
      if (j == 0)
        s[i + j * 800] = 0.0;
      else
        s[i + j * 800] = s[i + (j - 1) * 800] + m[i + (j - 1) * 800];
    }
  }
}


x86 assembly (clang, -O0):
<vertical_sum_cpu>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x18(%rbp)
cmpl   $0x7d0,-0x18(%rbp)
jge    11d6 <vertical_sum_cpu+0xd6>
movl   $0x0,-0x14(%rbp)
cmpl   $0x320,-0x14(%rbp)
jge    11c3 <vertical_sum_cpu+0xc3>
cmpl   $0x0,-0x18(%rbp)
jne    115e <vertical_sum_cpu+0x5e>
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   $0x320,-0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11b0 <vertical_sum_cpu+0xb0>
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
sub    $0x1,%edx
imul   $0x320,%edx,%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
mov    -0x18(%rbp),%edx
sub    $0x1,%edx
imul   $0x320,%edx,%edx
add    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   $0x320,-0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    11b5 <vertical_sum_cpu+0xb5>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1127 <vertical_sum_cpu+0x27>
jmp    11c8 <vertical_sum_cpu+0xc8>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1113 <vertical_sum_cpu+0x13>
pop    %rbp
ret

x86 assembly (clang, -O2):
<vertical_sum_cpu>:
add    $0xfffffffffffff384,%rdi
lea    -0xc7c(%rsi),%r10
mov    $0x1,%r8d
xor    %edx,%edx
jmp    1142 <vertical_sum_cpu+0x42>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdx
add    $0x320,%r8
add    $0xc80,%rdi
add    $0xc80,%r10
cmp    $0x7d0,%rdx
je     11a3 <vertical_sum_cpu+0xa3>
imul   $0x320,%rdx,%r9
xor    %ecx,%ecx
jmp    1171 <vertical_sum_cpu+0x71>
nopl   (%rax)
lea    (%r8,%rcx,1),%rax
movss  (%r10,%rcx,4),%xmm0
addss  (%rdi,%rcx,4),%xmm0
movss  %xmm0,(%rsi,%rax,4)
add    $0x2,%rcx
cmp    $0x320,%rcx
je     1120 <vertical_sum_cpu+0x20>
xorps  %xmm0,%xmm0
mov    %rcx,%rax
xorps  %xmm1,%xmm1
test   %rdx,%rdx
je     1190 <vertical_sum_cpu+0x90>
lea    (%rcx,%r9,1),%rax
movss  -0x4(%r10,%rcx,4),%xmm1
addss  -0x4(%rdi,%rcx,4),%xmm1
movss  %xmm1,(%rsi,%rax,4)
test   %rdx,%rdx
jne    1150 <vertical_sum_cpu+0x50>
mov    %rcx,%rax
or     $0x1,%rax
jmp    115f <vertical_sum_cpu+0x5f>
ret

LLVM IR (-O0):
define void @vertical_sum_cpu(float* %0, float* %1) {
%3 = alloca float*
%4 = alloca float*
%5 = alloca i32
%6 = alloca i32
store float* %0, float** %3
store float* %1, float** %4
store i32 0, i32* %6
br label %7
7: ; preds = %57, %2
%8 = load i32, i32* %6
%9 = icmp slt i32 %8, 2000
br i1 %9, label %10, label %60
10: ; preds = %7
store i32 0, i32* %5
br label %11
11: ; preds = %53, %10
%12 = load i32, i32* %5
%13 = icmp slt i32 %12, 800
br i1 %13, label %14, label %56
14: ; preds = %11
%15 = load i32, i32* %6
%16 = icmp eq i32 %15, 0
br i1 %16, label %17, label %25
17: ; preds = %14
%18 = load float*, float** %4
%19 = load i32, i32* %5
%20 = load i32, i32* %6
%21 = mul nsw i32 %20, 800
%22 = add nsw i32 %19, %21
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds float, float* %18, i64 %23
store float 0.000000e+00, float* %24
br label %52
25: ; preds = %14
%26 = load float*, float** %4
%27 = load i32, i32* %5
%28 = load i32, i32* %6
%29 = sub nsw i32 %28, 1
%30 = mul nsw i32 %29, 800
%31 = add nsw i32 %27, %30
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds float, float* %26, i64 %32
%34 = load float, float* %33
%35 = load float*, float** %3
%36 = load i32, i32* %5
%37 = load i32, i32* %6
%38 = sub nsw i32 %37, 1
%39 = mul nsw i32 %38, 800
%40 = add nsw i32 %36, %39
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float, float* %35, i64 %41
%43 = load float, float* %42
%44 = fadd float %34, %43
%45 = load float*, float** %4
%46 = load i32, i32* %5
%47 = load i32, i32* %6
%48 = mul nsw i32 %47, 800
%49 = add nsw i32 %46, %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds float, float* %45, i64 %50
store float %44, float* %51
br label %52
52: ; preds = %25, %17
br label %53
53: ; preds = %52
%54 = load i32, i32* %5
%55 = add nsw i32 %54, 1
store i32 %55, i32* %5
br label %11
56: ; preds = %11
br label %57
57: ; preds = %56
%58 = load i32, i32* %6
%59 = add nsw i32 %58, 1
store i32 %59, i32* %6
br label %7
60: ; preds = %7
ret void
}

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