float difference(v1, v2, pic_x, pic_y)
float v1[675][675][2], v2[675][675][2];
int pic_x, pic_y;
{
  int i, j, n;
  float sum, t1, t2;
  sum = 0.0;
  n = 0;
  for (i = 0; i < pic_x; i++)
    for (j = 0; j < pic_y; j++) {
      t1 = v1[i][j][0] - v2[i][j][0];
      t2 = v1[i][j][1] - v2[i][j][1];
      sum += sqrt(t1 * t1 + t2 * t2);
      n++;
    }
  sum = sum / n;
  return (sum);
}


x86 assembly (clang, -O0):
<difference>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    123b <difference+0x12b>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1228 <difference+0x118>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x1518,%rcx,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x1518,%rcx,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
subss  (%rax,%rcx,8),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x1518,%rcx,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  0x4(%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
imul   $0x1518,%rcx,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
subss  0x4(%rax,%rcx,8),%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x2c(%rbp),%xmm0
movss  -0x2c(%rbp),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  -0x30(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    114f <difference+0x3f>
jmp    122d <difference+0x11d>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    113c <difference+0x2c>
movss  -0x28(%rbp),%xmm0
cvtsi2ssl -0x24(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<difference>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
test   %edx,%edx
jle    120d <difference+0xfd>
mov    %ecx,%r15d
mov    %rsi,%rbx
mov    %rdi,%rbp
xor    %r13d,%r13d
test   %ecx,%ecx
mov    $0x0,%eax
cmovg  %ecx,%eax
mov    %eax,0x4(%rsp)
mov    %edx,0x8(%rsp)
mov    %edx,%eax
mov    %rax,0x10(%rsp)
mov    %ecx,%r12d
add    $0x4,%rbp
add    $0x4,%rbx
xorps  %xmm2,%xmm2
xorpd  %xmm3,%xmm3
jmp    1179 <difference+0x69>
nop
add    $0x1,%r13
add    $0x1518,%rbp
add    $0x1518,%rbx
cmp    0x10(%rsp),%r13
je     11f7 <difference+0xe7>
test   %r15d,%r15d
jle    1160 <difference+0x50>
xor    %r14d,%r14d
jmp    11c0 <difference+0xb0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  %xmm2,0xc(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm3,%xmm3
movss  0xc(%rsp),%xmm2
xorps  %xmm1,%xmm1
cvtss2sd %xmm2,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm2,%xmm2
cvtsd2ss %xmm1,%xmm2
add    $0x1,%r14
cmp    %r14,%r12
je     1160 <difference+0x50>
movss  -0x4(%rbp,%r14,8),%xmm0
movss  0x0(%rbp,%r14,8),%xmm1
subss  -0x4(%rbx,%r14,8),%xmm0
subss  (%rbx,%r14,8),%xmm1
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
ucomisd %xmm3,%xmm0
jb     1190 <difference+0x80>
sqrtsd %xmm0,%xmm0
jmp    11a5 <difference+0x95>
mov    0x4(%rsp),%eax
imul   0x8(%rsp),%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
divss  %xmm0,%xmm2
jmp    1215 <difference+0x105>
movss  0xdeb(%rip),%xmm2
movaps %xmm2,%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define float @difference([675 x [2 x float]]* %0, [675 x [2 x float]]* %1, i32 %2, i32 %3) {
%5 = alloca [675 x [2 x float]]*
%6 = alloca [675 x [2 x float]]*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca float
%13 = alloca float
%14 = alloca float
store [675 x [2 x float]]* %0, [675 x [2 x float]]** %5
store [675 x [2 x float]]* %1, [675 x [2 x float]]** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
store float 0.000000e+00, float* %12
store i32 0, i32* %11
store i32 0, i32* %9
br label %15
15: ; preds = %81, %4
%16 = load i32, i32* %9
%17 = load i32, i32* %7
%18 = icmp slt i32 %16, %17
br i1 %18, label %19, label %84
19: ; preds = %15
store i32 0, i32* %10
br label %20
20: ; preds = %77, %19
%21 = load i32, i32* %10
%22 = load i32, i32* %8
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %80
24: ; preds = %20
%25 = load [675 x [2 x float]]*, [675 x [2 x float]]** %5
%26 = load i32, i32* %9
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %25, i64 %27
%29 = load i32, i32* %10
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %28, i64 0, i64 %30
%32 = getelementptr inbounds [2 x float], [2 x float]* %31, i64 0, i64 0
%33 = load float, float* %32
%34 = load [675 x [2 x float]]*, [675 x [2 x float]]** %6
%35 = load i32, i32* %9
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %34, i64 %36
%38 = load i32, i32* %10
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %37, i64 0, i64 %39
%41 = getelementptr inbounds [2 x float], [2 x float]* %40, i64 0, i64 0
%42 = load float, float* %41
%43 = fsub float %33, %42
store float %43, float* %13
%44 = load [675 x [2 x float]]*, [675 x [2 x float]]** %5
%45 = load i32, i32* %9
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %44, i64 %46
%48 = load i32, i32* %10
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %47, i64 0, i64 %49
%51 = getelementptr inbounds [2 x float], [2 x float]* %50, i64 0, i64 1
%52 = load float, float* %51
%53 = load [675 x [2 x float]]*, [675 x [2 x float]]** %6
%54 = load i32, i32* %9
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %53, i64 %55
%57 = load i32, i32* %10
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %56, i64 0, i64 %58
%60 = getelementptr inbounds [2 x float], [2 x float]* %59, i64 0, i64 1
%61 = load float, float* %60
%62 = fsub float %52, %61
store float %62, float* %14
%63 = load float, float* %13
%64 = load float, float* %13
%65 = load float, float* %14
%66 = load float, float* %14
%67 = fmul float %65, %66
%68 = call float @llvm.fmuladd.f32(float %63, float %64, float %67)
%69 = fpext float %68 to double
%70 = call double @sqrt(double %69)
%71 = load float, float* %12
%72 = fpext float %71 to double
%73 = fadd double %72, %70
%74 = fptrunc double %73 to float
store float %74, float* %12
%75 = load i32, i32* %11
%76 = add nsw i32 %75, 1
store i32 %76, i32* %11
br label %77
77: ; preds = %24
%78 = load i32, i32* %10
%79 = add nsw i32 %78, 1
store i32 %79, i32* %10
br label %20
80: ; preds = %20
br label %81
81: ; preds = %80
%82 = load i32, i32* %9
%83 = add nsw i32 %82, 1
store i32 %83, i32* %9
br label %15
84: ; preds = %15
%85 = load float, float* %12
%86 = load i32, i32* %11
%87 = sitofp i32 %86 to float
%88 = fdiv float %85, %87
store float %88, float* %12
%89 = load float, float* %12
ret float %89
}

LLVM IR (-O2):
define float @difference([675 x [2 x float]]* nocapture %0, [675 x [2 x float]]* nocapture %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %44
6: ; preds = %4
%7 = icmp sgt i32 %3, 0
%8 = call i32 @llvm.smax.i32(i32 %3, i32 0)
%9 = zext i32 %2 to i64
%10 = zext i32 %3 to i64
br label %11
11: ; preds = %6, %36
%12 = phi i64 [ 0, %6 ], [ %38, %36 ]
%13 = phi float [ 0.000000e+00, %6 ], [ %37, %36 ]
br i1 %7, label %14, label %36
14: ; preds = %11, %14
%15 = phi i64 [ %34, %14 ], [ 0, %11 ]
%16 = phi float [ %33, %14 ], [ %13, %11 ]
%17 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %0, i64 %12, i64 %15, i64 0
%18 = load float, float* %17
%19 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %1, i64 %12, i64 %15, i64 0
%20 = load float, float* %19
%21 = fsub float %18, %20
%22 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %0, i64 %12, i64 %15, i64 1
%23 = load float, float* %22
%24 = getelementptr inbounds [675 x [2 x float]], [675 x [2 x float]]* %1, i64 %12, i64 %15, i64 1
%25 = load float, float* %24
%26 = fsub float %23, %25
%27 = fmul float %26, %26
%28 = tail call float @llvm.fmuladd.f32(float %21, float %21, float %27)
%29 = fpext float %28 to double
%30 = tail call double @sqrt(double %29)
%31 = fpext float %16 to double
%32 = fadd double %30, %31
%33 = fptrunc double %32 to float
%34 = add nuw nsw i64 %15, 1
%35 = icmp eq i64 %34, %10
br i1 %35, label %36, label %14
36: ; preds = %14, %11
%37 = phi float [ %13, %11 ], [ %33, %14 ]
%38 = add nuw nsw i64 %12, 1
%39 = icmp eq i64 %38, %9
br i1 %39, label %40, label %11
40: ; preds = %36
%41 = mul i32 %8, %2
%42 = sitofp i32 %41 to float
%43 = fdiv float %37, %42
br label %44
44: ; preds = %40, %4
%45 = phi float [ %43, %40 ], [ 0x7FF8000000000000, %4 ]
ret float %45
}