double fvec_norm2sqr(const float *v, long n) {
  double s = 0;
  long i;
  for (i = 0; i < n; i++) {
    s += v[i] * v[i];
  }
  return s;
}


x86 assembly (clang, -O0):
<fvec_norm2sqr>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movq   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    1163 <fvec_norm2sqr+0x63>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
jmp    111c <fvec_norm2sqr+0x1c>
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<fvec_norm2sqr>:
test   %rsi,%rsi
jle    111c <fvec_norm2sqr+0x1c>
lea    -0x1(%rsi),%rcx
mov    %esi,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    1120 <fvec_norm2sqr+0x20>
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
jmp    1183 <fvec_norm2sqr+0x83>
xorps  %xmm0,%xmm0
ret
and    $0xfffffffffffffffc,%rsi
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movss  (%rdi,%rcx,4),%xmm1
movss  0x4(%rdi,%rcx,4),%xmm2
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm0,%xmm1
mulss  %xmm2,%xmm2
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
addsd  %xmm1,%xmm0
movss  0x8(%rdi,%rcx,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm0,%xmm1
movss  0xc(%rdi,%rcx,4),%xmm0
mulss  %xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rcx
cmp    %rcx,%rsi
jne    1130 <fvec_norm2sqr+0x30>
test   %rax,%rax
je     11aa <fvec_norm2sqr+0xaa>
lea    (%rdi,%rcx,4),%rcx
xor    %edx,%edx
xchg   %ax,%ax
movss  (%rcx,%rdx,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1190 <fvec_norm2sqr+0x90>
ret

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

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