double fvec_distance_L2sqr(const float *v1, const float *v2, long n) {
  long i;
  double dis = 0, a;
  for (i = 0; i < n; i++) {
    a = (double)v1[i] - v2[i];
    dis += a * a;
  }
  return dis;
}


x86 assembly (clang, -O0):
<fvec_distance_L2sqr>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movq   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    1186 <fvec_distance_L2sqr+0x86>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    1120 <fvec_distance_L2sqr+0x20>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<fvec_distance_L2sqr>:
test   %rdx,%rdx
jle    1113 <fvec_distance_L2sqr+0x13>
cmp    $0x1,%rdx
jne    1117 <fvec_distance_L2sqr+0x17>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    117d <fvec_distance_L2sqr+0x7d>
xorps  %xmm0,%xmm0
ret
mov    %rdx,%rcx
and    $0xfffffffffffffffe,%rcx
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rdi,%rax,4),%xmm1
movss  0x4(%rdi,%rax,4),%xmm2
cvtss2sd %xmm1,%xmm1
movss  (%rsi,%rax,4),%xmm3
movss  0x4(%rsi,%rax,4),%xmm4
cvtss2sd %xmm3,%xmm3
subsd  %xmm3,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
xorps  %xmm2,%xmm2
cvtss2sd %xmm4,%xmm2
subsd  %xmm2,%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x2,%rax
cmp    %rax,%rcx
jne    1130 <fvec_distance_L2sqr+0x30>
test   $0x1,%dl
je     11a4 <fvec_distance_L2sqr+0xa4>
movss  (%rdi,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
movss  (%rsi,%rax,4),%xmm2
cvtss2sd %xmm2,%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @fvec_distance_L2sqr(float* %0, float* %1, i64 %2) {
%4 = alloca float*
%5 = alloca float*
%6 = alloca i64
%7 = alloca i64
%8 = alloca double
%9 = alloca double
store float* %0, float** %4
store float* %1, float** %5
store i64 %2, i64* %6
store double 0.000000e+00, double* %8
store i64 0, i64* %7
br label %10
10: ; preds = %30, %3
%11 = load i64, i64* %7
%12 = load i64, i64* %6
%13 = icmp slt i64 %11, %12
br i1 %13, label %14, label %33
14: ; preds = %10
%15 = load float*, float** %4
%16 = load i64, i64* %7
%17 = getelementptr inbounds float, float* %15, i64 %16
%18 = load float, float* %17
%19 = fpext float %18 to double
%20 = load float*, float** %5
%21 = load i64, i64* %7
%22 = getelementptr inbounds float, float* %20, i64 %21
%23 = load float, float* %22
%24 = fpext float %23 to double
%25 = fsub double %19, %24
store double %25, double* %9
%26 = load double, double* %9
%27 = load double, double* %9
%28 = load double, double* %8
%29 = call double @llvm.fmuladd.f64(double %26, double %27, double %28)
store double %29, double* %8
br label %30
30: ; preds = %14
%31 = load i64, i64* %7
%32 = add nsw i64 %31, 1
store i64 %32, i64* %7
br label %10
33: ; preds = %10
%34 = load double, double* %8
ret double %34
}

LLVM IR (-O2):
define double @fvec_distance_L2sqr(float* nocapture %0, float* nocapture %1, i64 %2) {
%4 = icmp sgt i64 %2, 0
br i1 %4, label %5, label %48
5: ; preds = %3
%6 = and i64 %2, 1
%7 = icmp eq i64 %2, 1
br i1 %7, label %34, label %8
8: ; preds = %5
%9 = and i64 %2, -2
br label %10
10: ; preds = %10, %8
%11 = phi i64 [ 0, %8 ], [ %31, %10 ]
%12 = phi double [ 0.000000e+00, %8 ], [ %30, %10 ]
%13 = phi i64 [ 0, %8 ], [ %32, %10 ]
%14 = getelementptr inbounds float, float* %0, i64 %11
%15 = load float, float* %14
%16 = fpext float %15 to double
%17 = getelementptr inbounds float, float* %1, i64 %11
%18 = load float, float* %17
%19 = fpext float %18 to double
%20 = fsub double %16, %19
%21 = tail call double @llvm.fmuladd.f64(double %20, double %20, double %12)
%22 = or i64 %11, 1
%23 = getelementptr inbounds float, float* %0, i64 %22
%24 = load float, float* %23
%25 = fpext float %24 to double
%26 = getelementptr inbounds float, float* %1, i64 %22
%27 = load float, float* %26
%28 = fpext float %27 to double
%29 = fsub double %25, %28
%30 = tail call double @llvm.fmuladd.f64(double %29, double %29, double %21)
%31 = add nuw nsw i64 %11, 2
%32 = add i64 %13, 2
%33 = icmp eq i64 %32, %9
br i1 %33, label %34, label %10
34: ; preds = %10, %5
%35 = phi double [ undef, %5 ], [ %30, %10 ]
%36 = phi i64 [ 0, %5 ], [ %31, %10 ]
%37 = phi double [ 0.000000e+00, %5 ], [ %30, %10 ]
%38 = icmp eq i64 %6, 0
br i1 %38, label %48, label %39
39: ; preds = %34
%40 = getelementptr inbounds float, float* %0, i64 %36
%41 = load float, float* %40
%42 = fpext float %41 to double
%43 = getelementptr inbounds float, float* %1, i64 %36
%44 = load float, float* %43
%45 = fpext float %44 to double
%46 = fsub double %42, %45
%47 = tail call double @llvm.fmuladd.f64(double %46, double %46, double %37)
br label %48
48: ; preds = %39, %34, %3
%49 = phi double [ 0.000000e+00, %3 ], [ %35, %34 ], [ %47, %39 ]
ret double %49
}