long double cvector_dot_l(const double *a, const double *b, int n) {
  long double sum = 0.0;
  for (int i = 0; i < n; i++) {
    sum += a[i] * b[i];
  }
  return sum;
}


x86 assembly (clang, -O0):
<cvector_dot_l>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
fldz
fstpt  -0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    115f <cvector_dot_l+0x5f>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
fldl   -0x40(%rbp)
fldt   -0x30(%rbp)
faddp  %st,%st(1)
fstpt  -0x30(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    111b <cvector_dot_l+0x1b>
fldt   -0x30(%rbp)
pop    %rbp
ret

x86 assembly (clang, -O2):
<cvector_dot_l>:
fldz
test   %edx,%edx
jle    11bd <cvector_dot_l+0xbd>
fstp   %st(0)
mov    %edx,%ecx
lea    -0x1(%rcx),%rdx
mov    %ecx,%eax
and    $0x3,%eax
fldz
cmp    $0x3,%rdx
jae    1123 <cvector_dot_l+0x23>
xor    %edx,%edx
jmp    118f <cvector_dot_l+0x8f>
fstp   %st(0)
and    $0xfffffffc,%ecx
fldz
xor    %edx,%edx
nopl   0x0(%rax)
movsd  (%rdi,%rdx,8),%xmm0
movsd  0x8(%rdi,%rdx,8),%xmm1
mulsd  (%rsi,%rdx,8),%xmm0
movsd  %xmm0,-0x8(%rsp)
faddl  -0x8(%rsp)
mulsd  0x8(%rsi,%rdx,8),%xmm1
movsd  %xmm1,-0x20(%rsp)
faddl  -0x20(%rsp)
movsd  0x10(%rdi,%rdx,8),%xmm0
mulsd  0x10(%rsi,%rdx,8),%xmm0
movsd  %xmm0,-0x18(%rsp)
faddl  -0x18(%rsp)
movsd  0x18(%rdi,%rdx,8),%xmm0
mulsd  0x18(%rsi,%rdx,8),%xmm0
movsd  %xmm0,-0x10(%rsp)
faddl  -0x10(%rsp)
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    1130 <cvector_dot_l+0x30>
test   %rax,%rax
je     11bd <cvector_dot_l+0xbd>
lea    (%rsi,%rdx,8),%rcx
lea    (%rdi,%rdx,8),%rdx
xor    %esi,%esi
xchg   %ax,%ax
movsd  (%rdx,%rsi,8),%xmm0
mulsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,-0x28(%rsp)
faddl  -0x28(%rsp)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    11a0 <cvector_dot_l+0xa0>
ret

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

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