double c_sdot(int n, double *sx, double *sy) {
  register int i, m;
  double ans;
  ans = 0.;
  if (n <= 0) {
    return ans;
  }
  m = n % 4;
  if (m != 0) {
    for (i = 1; i <= m; i++) {
      ans += sx[i - 1] * sy[i - 1];
    }
  }
  for (i = m + 1; i <= n; i += 4) {
    ans += sx[i - 1] * sy[i - 1] + sx[i + 1 - 1] * sy[i + 1 - 1] +
           sx[i + 2 - 1] * sy[i + 2 - 1] + sx[i + 3 - 1] * sy[i + 3 - 1];
  }
  return ans;
}


x86 assembly (clang, -O0):
<c_sdot>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
cmpl   $0x0,-0xc(%rbp)
jg     1130 <c_sdot+0x30>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1295 <c_sdot+0x195>
mov    -0xc(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %edx,-0x28(%rbp)
cmpl   $0x0,-0x28(%rbp)
je     11a4 <c_sdot+0xa4>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
jg     119f <c_sdot+0x9f>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    114f <c_sdot+0x4f>
jmp    11a4 <c_sdot+0xa4>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     128b <c_sdot+0x18b>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x2,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x2,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x3,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x3,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x24(%rbp)
jmp    11ad <c_sdot+0xad>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<c_sdot>:
test   %edi,%edi
jle    1146 <c_sdot+0x46>
mov    %edi,%eax
and    $0x3,%eax
je     114a <c_sdot+0x4a>
movsd  (%rsi),%xmm1
mulsd  (%rdx),%xmm1
lea    0x1(%rax),%ecx
xorpd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
cmp    $0x2,%ecx
je     114e <c_sdot+0x4e>
movsd  0x8(%rsi),%xmm1
mulsd  0x8(%rdx),%xmm1
addsd  %xmm1,%xmm0
cmp    $0x3,%ecx
je     114e <c_sdot+0x4e>
movsd  0x10(%rsi),%xmm1
mulsd  0x10(%rdx),%xmm1
addsd  %xmm1,%xmm0
jmp    114e <c_sdot+0x4e>
xorps  %xmm0,%xmm0
ret
xorpd  %xmm0,%xmm0
cmp    $0x4,%edi
jl     11a7 <c_sdot+0xa7>
add    $0x1,%eax
mov    %edi,%ecx
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rsi,%rax,8),%xmm1
movsd  (%rsi,%rax,8),%xmm2
mulsd  (%rdx,%rax,8),%xmm2
mulsd  -0x8(%rdx,%rax,8),%xmm1
addsd  %xmm2,%xmm1
movsd  0x8(%rsi,%rax,8),%xmm2
mulsd  0x8(%rdx,%rax,8),%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rsi,%rax,8),%xmm1
mulsd  0x10(%rdx,%rax,8),%xmm1
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
add    $0x4,%rax
cmp    %rcx,%rax
jbe    1160 <c_sdot+0x60>
ret

LLVM IR (-O0):
define double @c_sdot(i32 %0, double* %1, double* %2) {
%4 = alloca double
%5 = alloca i32
%6 = alloca double*
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
store i32 %0, i32* %5
store double* %1, double** %6
store double* %2, double** %7
store double 0.000000e+00, double* %10
%11 = load i32, i32* %5
%12 = icmp sle i32 %11, 0
br i1 %12, label %13, label %15
13: ; preds = %3
%14 = load double, double* %10
store double %14, double* %4
br label %117
15: ; preds = %3
%16 = load i32, i32* %5
%17 = srem i32 %16, 4
store i32 %17, i32* %9
%18 = load i32, i32* %9
%19 = icmp ne i32 %18, 0
br i1 %19, label %20, label %44
20: ; preds = %15
store i32 1, i32* %8
br label %21
21: ; preds = %40, %20
%22 = load i32, i32* %8
%23 = load i32, i32* %9
%24 = icmp sle i32 %22, %23
br i1 %24, label %25, label %43
25: ; preds = %21
%26 = load double*, double** %6
%27 = load i32, i32* %8
%28 = sub nsw i32 %27, 1
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %26, i64 %29
%31 = load double, double* %30
%32 = load double*, double** %7
%33 = load i32, i32* %8
%34 = sub nsw i32 %33, 1
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %32, i64 %35
%37 = load double, double* %36
%38 = load double, double* %10
%39 = call double @llvm.fmuladd.f64(double %31, double %37, double %38)
store double %39, double* %10
br label %40
40: ; preds = %25
%41 = load i32, i32* %8
%42 = add nsw i32 %41, 1
store i32 %42, i32* %8
br label %21
43: ; preds = %21
br label %44
44: ; preds = %43, %15
%45 = load i32, i32* %9
%46 = add nsw i32 %45, 1
store i32 %46, i32* %8
br label %47
47: ; preds = %112, %44
%48 = load i32, i32* %8
%49 = load i32, i32* %5
%50 = icmp sle i32 %48, %49
br i1 %50, label %51, label %115
51: ; preds = %47
%52 = load double*, double** %6
%53 = load i32, i32* %8
%54 = sub nsw i32 %53, 1
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %52, i64 %55
%57 = load double, double* %56
%58 = load double*, double** %7
%59 = load i32, i32* %8
%60 = sub nsw i32 %59, 1
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %58, i64 %61
%63 = load double, double* %62
%64 = load double*, double** %6
%65 = load i32, i32* %8
%66 = add nsw i32 %65, 1
%67 = sub nsw i32 %66, 1
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %64, i64 %68
%70 = load double, double* %69
%71 = load double*, double** %7
%72 = load i32, i32* %8
%73 = add nsw i32 %72, 1
%74 = sub nsw i32 %73, 1
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %71, i64 %75
%77 = load double, double* %76
%78 = fmul double %70, %77
%79 = call double @llvm.fmuladd.f64(double %57, double %63, double %78)
%80 = load double*, double** %6
%81 = load i32, i32* %8
%82 = add nsw i32 %81, 2
%83 = sub nsw i32 %82, 1
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %80, i64 %84
%86 = load double, double* %85
%87 = load double*, double** %7
%88 = load i32, i32* %8
%89 = add nsw i32 %88, 2
%90 = sub nsw i32 %89, 1
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %87, i64 %91
%93 = load double, double* %92
%94 = call double @llvm.fmuladd.f64(double %86, double %93, double %79)
%95 = load double*, double** %6
%96 = load i32, i32* %8
%97 = add nsw i32 %96, 3
%98 = sub nsw i32 %97, 1
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds double, double* %95, i64 %99
%101 = load double, double* %100
%102 = load double*, double** %7
%103 = load i32, i32* %8
%104 = add nsw i32 %103, 3
%105 = sub nsw i32 %104, 1
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds double, double* %102, i64 %106
%108 = load double, double* %107
%109 = call double @llvm.fmuladd.f64(double %101, double %108, double %94)
%110 = load double, double* %10
%111 = fadd double %110, %109
store double %111, double* %10
br label %112
112: ; preds = %51
%113 = load i32, i32* %8
%114 = add nsw i32 %113, 4
store i32 %114, i32* %8
br label %47
115: ; preds = %47
%116 = load double, double* %10
store double %116, double* %4
br label %117
117: ; preds = %115, %13
%118 = load double, double* %4
ret double %118
}

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