void are_blocks_equal(float ref_block[32][32], float to_comp[32][32],
                      float *mse) {
  int i, j;
  float diff;
  *mse = 0.0;
  for (i = 0; i < 32; i++)
    for (j = 0; j < 32; j++) {
      diff = ref_block[i][j] - to_comp[i][j];
      *mse += diff * diff;
    }
}


x86 assembly (clang, -O0):
<are_blocks_equal>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x20,-0x1c(%rbp)
jge    11b1 <are_blocks_equal+0xb1>
movl   $0x0,-0x20(%rbp)
cmpl   $0x20,-0x20(%rbp)
jge    119e <are_blocks_equal+0x9e>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
movss  -0x24(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1133 <are_blocks_equal+0x33>
jmp    11a3 <are_blocks_equal+0xa3>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1122 <are_blocks_equal+0x22>
pop    %rbp
ret

x86 assembly (clang, -O2):
<are_blocks_equal>:
movl   $0x0,(%rdx)
add    $0x4,%rdi
add    $0x4,%rsi
xorps  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  -0x4(%rdi,%rcx,4),%xmm1
subss  -0x4(%rsi,%rcx,4),%xmm1
mulss  %xmm1,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rdx)
movss  (%rdi,%rcx,4),%xmm0
subss  (%rsi,%rcx,4),%xmm0
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rdx)
add    $0x2,%rcx
cmp    $0x20,%rcx
jne    1130 <are_blocks_equal+0x30>
add    $0x1,%rax
sub    $0xffffffffffffff80,%rdi
sub    $0xffffffffffffff80,%rsi
cmp    $0x20,%rax
jne    1120 <are_blocks_equal+0x20>
ret

LLVM IR (-O0):
define void @are_blocks_equal([32 x float]* %0, [32 x float]* %1, float* %2) {
%4 = alloca [32 x float]*
%5 = alloca [32 x float]*
%6 = alloca float*
%7 = alloca i32
%8 = alloca i32
%9 = alloca float
store [32 x float]* %0, [32 x float]** %4
store [32 x float]* %1, [32 x float]** %5
store float* %2, float** %6
%10 = load float*, float** %6
store float 0.000000e+00, float* %10
store i32 0, i32* %7
br label %11
11: ; preds = %45, %3
%12 = load i32, i32* %7
%13 = icmp slt i32 %12, 32
br i1 %13, label %14, label %48
14: ; preds = %11
store i32 0, i32* %8
br label %15
15: ; preds = %41, %14
%16 = load i32, i32* %8
%17 = icmp slt i32 %16, 32
br i1 %17, label %18, label %44
18: ; preds = %15
%19 = load [32 x float]*, [32 x float]** %4
%20 = load i32, i32* %7
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds [32 x float], [32 x float]* %19, i64 %21
%23 = load i32, i32* %8
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds [32 x float], [32 x float]* %22, i64 0, i64 %24
%26 = load float, float* %25
%27 = load [32 x float]*, [32 x float]** %5
%28 = load i32, i32* %7
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds [32 x float], [32 x float]* %27, i64 %29
%31 = load i32, i32* %8
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds [32 x float], [32 x float]* %30, i64 0, i64 %32
%34 = load float, float* %33
%35 = fsub float %26, %34
store float %35, float* %9
%36 = load float, float* %9
%37 = load float, float* %9
%38 = load float*, float** %6
%39 = load float, float* %38
%40 = call float @llvm.fmuladd.f32(float %36, float %37, float %39)
store float %40, float* %38
br label %41
41: ; preds = %18
%42 = load i32, i32* %8
%43 = add nsw i32 %42, 1
store i32 %43, i32* %8
br label %15
44: ; preds = %15
br label %45
45: ; preds = %44
%46 = load i32, i32* %7
%47 = add nsw i32 %46, 1
store i32 %47, i32* %7
br label %11
48: ; preds = %11
ret void
}

LLVM IR (-O2):
define void @are_blocks_equal([32 x float]* nocapture %0, [32 x float]* nocapture %1, float* nocapture %2) {
store float 0.000000e+00, float* %2
br label %4
4: ; preds = %3, %25
%5 = phi float [ 0.000000e+00, %3 ], [ %22, %25 ]
%6 = phi i64 [ 0, %3 ], [ %26, %25 ]
br label %7
7: ; preds = %7, %4
%8 = phi float [ %5, %4 ], [ %22, %7 ]
%9 = phi i64 [ 0, %4 ], [ %23, %7 ]
%10 = getelementptr inbounds [32 x float], [32 x float]* %0, i64 %6, i64 %9
%11 = load float, float* %10
%12 = getelementptr inbounds [32 x float], [32 x float]* %1, i64 %6, i64 %9
%13 = load float, float* %12
%14 = fsub float %11, %13
%15 = tail call float @llvm.fmuladd.f32(float %14, float %14, float %8)
store float %15, float* %2
%16 = or i64 %9, 1
%17 = getelementptr inbounds [32 x float], [32 x float]* %0, i64 %6, i64 %16
%18 = load float, float* %17
%19 = getelementptr inbounds [32 x float], [32 x float]* %1, i64 %6, i64 %16
%20 = load float, float* %19
%21 = fsub float %18, %20
%22 = tail call float @llvm.fmuladd.f32(float %21, float %21, float %15)
store float %22, float* %2
%23 = add nuw nsw i64 %9, 2
%24 = icmp eq i64 %23, 32
br i1 %24, label %25, label %7
25: ; preds = %7
%26 = add nuw nsw i64 %6, 1
%27 = icmp eq i64 %26, 32
br i1 %27, label %28, label %4
28: ; preds = %25
ret void
}