double fvec_product(const float *v, long n) {
  double p = 1.0;
  while (n--)
    p *= v[n];
  return p;
}


x86 assembly (clang, -O0):
<fvec_product>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  0xeec(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0xffffffffffffffff,%rcx
mov    %rcx,-0x10(%rbp)
cmp    $0x0,%rax
je     1152 <fvec_product+0x52>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    1119 <fvec_product+0x19>
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<fvec_product>:
test   %rsi,%rsi
je     1182 <fvec_product+0x82>
lea    -0x1(%rsi),%rax
mov    %rsi,%rcx
movsd  0xeec(%rip),%xmm0
and    $0x3,%rcx
je     1138 <fvec_product+0x38>
nopw   0x0(%rax,%rax,1)
movss  -0x4(%rdi,%rsi,4),%xmm1
add    $0xffffffffffffffff,%rsi
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
add    $0xffffffffffffffff,%rcx
jne    1120 <fvec_product+0x20>
cmp    $0x3,%rax
jb     1181 <fvec_product+0x81>
xchg   %ax,%ax
movss  -0x4(%rdi,%rsi,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm0,%xmm1
movss  -0x8(%rdi,%rsi,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
movss  -0x10(%rdi,%rsi,4),%xmm1
movss  -0xc(%rdi,%rsi,4),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
mulsd  %xmm2,%xmm0
add    $0xfffffffffffffffc,%rsi
jne    1140 <fvec_product+0x40>
ret
movsd  0xe76(%rip),%xmm0
ret

LLVM IR (-O0):
define double @fvec_product(float* %0, i64 %1) {
%3 = alloca float*
%4 = alloca i64
%5 = alloca double
store float* %0, float** %3
store i64 %1, i64* %4
store double 1.000000e+00, double* %5
br label %6
6: ; preds = %10, %2
%7 = load i64, i64* %4
%8 = add nsw i64 %7, -1
store i64 %8, i64* %4
%9 = icmp ne i64 %7, 0
br i1 %9, label %10, label %18
10: ; preds = %6
%11 = load float*, float** %3
%12 = load i64, i64* %4
%13 = getelementptr inbounds float, float* %11, i64 %12
%14 = load float, float* %13
%15 = fpext float %14 to double
%16 = load double, double* %5
%17 = fmul double %16, %15
store double %17, double* %5
br label %6
18: ; preds = %6
%19 = load double, double* %5
ret double %19
}

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