double svd_ddot(long n, double *dx, long incx, double *dy, long incy) {
  long i;
  double dot_product;
  if (n <= 0 || incx == 0 || incy == 0)
    return (0.0);
  dot_product = 0.0;
  if (incx == 1 && incy == 1)
    for (i = 0; i < n; i++)
      dot_product += (*dx++) * (*dy++);
  else {
    if (incx < 0)
      dx += (-n + 1) * incx;
    if (incy < 0)
      dy += (-n + 1) * incy;
    for (i = 0; i < n; i++) {
      dot_product += (*dx) * (*dy);
      dx += incx;
      dy += incy;
    }
  }
  return (dot_product);
}


x86 assembly (clang, -O0):
<svd_ddot>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
cmpq   $0x0,-0x10(%rbp)
jle    1139 <svd_ddot+0x39>
cmpq   $0x0,-0x20(%rbp)
je     1139 <svd_ddot+0x39>
cmpq   $0x0,-0x30(%rbp)
jne    1146 <svd_ddot+0x46>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    128c <svd_ddot+0x18c>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
cmpq   $0x1,-0x20(%rbp)
jne    11c8 <svd_ddot+0xc8>
cmpq   $0x1,-0x30(%rbp)
jne    11c8 <svd_ddot+0xc8>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11c3 <svd_ddot+0xc3>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x18(%rbp)
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x28(%rbp)
movsd  (%rax),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    116c <svd_ddot+0x6c>
jmp    1282 <svd_ddot+0x182>
cmpq   $0x0,-0x20(%rbp)
jge    11ee <svd_ddot+0xee>
xor    %eax,%eax
sub    -0x10(%rbp),%rax
add    $0x1,%rax
imul   -0x20(%rbp),%rax
shl    $0x3,%rax
add    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
cmpq   $0x0,-0x30(%rbp)
jge    1214 <svd_ddot+0x114>
xor    %eax,%eax
sub    -0x10(%rbp),%rax
add    $0x1,%rax
imul   -0x30(%rbp),%rax
shl    $0x3,%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x28(%rbp)
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    127d <svd_ddot+0x17d>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
shl    $0x3,%rax
add    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
mov    -0x30(%rbp),%rax
shl    $0x3,%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    121c <svd_ddot+0x11c>
jmp    1282 <svd_ddot+0x182>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<svd_ddot>:
push   %r15
push   %r14
push   %rbx
xorpd  %xmm0,%xmm0
test   %rdi,%rdi
jle    1308 <svd_ddot+0x208>
test   %rdx,%rdx
je     1308 <svd_ddot+0x208>
test   %r8,%r8
je     1308 <svd_ddot+0x208>
cmp    $0x1,%rdx
jne    1151 <svd_ddot+0x51>
cmp    $0x1,%r8
jne    1151 <svd_ddot+0x51>
lea    -0x1(%rdi),%rdx
mov    %edi,%eax
and    $0x3,%eax
cmp    $0x3,%rdx
jae    11a7 <svd_ddot+0xa7>
xorpd  %xmm0,%xmm0
test   %rax,%rax
jne    121a <svd_ddot+0x11a>
jmp    1308 <svd_ddot+0x208>
mov    $0x1,%eax
sub    %rdi,%rax
mov    %rax,%r9
imul   %r8,%r9
mov    %r8,%rbx
sar    $0x3f,%rbx
and    %r9,%rbx
lea    (%rcx,%rbx,8),%rcx
imul   %rdx,%rax
mov    %rdx,%rbx
sar    $0x3f,%rbx
and    %rax,%rbx
lea    (%rsi,%rbx,8),%rsi
lea    -0x1(%rdi),%rax
mov    %edi,%r11d
and    $0x3,%r11d
cmp    $0x3,%rax
jae    123c <svd_ddot+0x13c>
xorpd  %xmm0,%xmm0
test   %r11,%r11
jne    12dd <svd_ddot+0x1dd>
jmp    1308 <svd_ddot+0x208>
and    $0xfffffffffffffffc,%rdi
xorpd  %xmm0,%xmm0
xor    %edx,%edx
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rsi,%rbx,8),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
movsd  0x8(%rsi,%rbx,8),%xmm2
addsd  %xmm0,%xmm1
mulsd  0x8(%rcx,%rbx,8),%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rsi,%rbx,8),%xmm1
mulsd  0x10(%rcx,%rbx,8),%xmm1
addsd  %xmm2,%xmm1
movsd  0x18(%rsi,%rbx,8),%xmm0
mulsd  0x18(%rcx,%rbx,8),%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rbx
add    $0xffffffffffffffe0,%rdx
cmp    %rbx,%rdi
jne    11c0 <svd_ddot+0xc0>
sub    %rdx,%rcx
sub    %rdx,%rsi
test   %rax,%rax
je     1308 <svd_ddot+0x208>
xor    %edx,%edx
nopl   0x0(%rax)
movsd  (%rsi,%rdx,8),%xmm1
mulsd  (%rcx,%rdx,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1220 <svd_ddot+0x120>
jmp    1308 <svd_ddot+0x208>
and    $0xfffffffffffffffc,%rdi
mov    %r8,%r9
shl    $0x5,%r9
mov    %rdx,%r10
shl    $0x5,%r10
lea    0x0(,%rdx,8),%rax
lea    0x0(,%r8,8),%rbx
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %rcx,%r14
mov    %rsi,%r15
movsd  (%rsi),%xmm1
mulsd  (%rcx),%xmm1
addsd  %xmm0,%xmm1
add    %rax,%rsi
movsd  (%r15,%rax,1),%xmm0
add    %rbx,%rcx
mulsd  (%r14,%rbx,1),%xmm0
addsd  %xmm1,%xmm0
movsd  (%rax,%rsi,1),%xmm1
mulsd  (%rbx,%rcx,1),%xmm1
add    %rax,%rsi
add    %rbx,%rcx
addsd  %xmm0,%xmm1
movsd  (%rax,%rsi,1),%xmm0
add    %rax,%rsi
mulsd  (%rbx,%rcx,1),%xmm0
add    %rbx,%rcx
addsd  %xmm1,%xmm0
add    %rax,%rsi
add    %rbx,%rcx
add    $0xfffffffffffffffc,%rdi
jne    1270 <svd_ddot+0x170>
add    %r9,%r14
add    %r10,%r15
mov    %r14,%rcx
mov    %r15,%rsi
test   %r11,%r11
je     1308 <svd_ddot+0x208>
shl    $0x3,%r8
shl    $0x3,%rdx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rsi),%xmm1
mulsd  (%rcx),%xmm1
addsd  %xmm1,%xmm0
add    %r8,%rcx
add    %rdx,%rsi
add    $0xffffffffffffffff,%r11
jne    12f0 <svd_ddot+0x1f0>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define double @svd_ddot(i64 %0, double* %1, i64 %2, double* %3, i64 %4) {
%6 = alloca double
%7 = alloca i64
%8 = alloca double*
%9 = alloca i64
%10 = alloca double*
%11 = alloca i64
%12 = alloca i64
%13 = alloca double
store i64 %0, i64* %7
store double* %1, double** %8
store i64 %2, i64* %9
store double* %3, double** %10
store i64 %4, i64* %11
%14 = load i64, i64* %7
%15 = icmp sle i64 %14, 0
br i1 %15, label %22, label %16
16: ; preds = %5
%17 = load i64, i64* %9
%18 = icmp eq i64 %17, 0
br i1 %18, label %22, label %19
19: ; preds = %16
%20 = load i64, i64* %11
%21 = icmp eq i64 %20, 0
br i1 %21, label %22, label %23
22: ; preds = %19, %16, %5
store double 0.000000e+00, double* %6
br label %93
23: ; preds = %19
store double 0.000000e+00, double* %13
%24 = load i64, i64* %9
%25 = icmp eq i64 %24, 1
br i1 %25, label %26, label %47
26: ; preds = %23
%27 = load i64, i64* %11
%28 = icmp eq i64 %27, 1
br i1 %28, label %29, label %47
29: ; preds = %26
store i64 0, i64* %12
br label %30
30: ; preds = %43, %29
%31 = load i64, i64* %12
%32 = load i64, i64* %7
%33 = icmp slt i64 %31, %32
br i1 %33, label %34, label %46
34: ; preds = %30
%35 = load double*, double** %8
%36 = getelementptr inbounds double, double* %35, i32 1
store double* %36, double** %8
%37 = load double, double* %35
%38 = load double*, double** %10
%39 = getelementptr inbounds double, double* %38, i32 1
store double* %39, double** %10
%40 = load double, double* %38
%41 = load double, double* %13
%42 = call double @llvm.fmuladd.f64(double %37, double %40, double %41)
store double %42, double* %13
br label %43
43: ; preds = %34
%44 = load i64, i64* %12
%45 = add nsw i64 %44, 1
store i64 %45, i64* %12
br label %30
46: ; preds = %30
br label %91
47: ; preds = %26, %23
%48 = load i64, i64* %9
%49 = icmp slt i64 %48, 0
br i1 %49, label %50, label %58
50: ; preds = %47
%51 = load i64, i64* %7
%52 = sub nsw i64 0, %51
%53 = add nsw i64 %52, 1
%54 = load i64, i64* %9
%55 = mul nsw i64 %53, %54
%56 = load double*, double** %8
%57 = getelementptr inbounds double, double* %56, i64 %55
store double* %57, double** %8
br label %58
58: ; preds = %50, %47
%59 = load i64, i64* %11
%60 = icmp slt i64 %59, 0
br i1 %60, label %61, label %69
61: ; preds = %58
%62 = load i64, i64* %7
%63 = sub nsw i64 0, %62
%64 = add nsw i64 %63, 1
%65 = load i64, i64* %11
%66 = mul nsw i64 %64, %65
%67 = load double*, double** %10
%68 = getelementptr inbounds double, double* %67, i64 %66
store double* %68, double** %10
br label %69
69: ; preds = %61, %58
store i64 0, i64* %12
br label %70
70: ; preds = %87, %69
%71 = load i64, i64* %12
%72 = load i64, i64* %7
%73 = icmp slt i64 %71, %72
br i1 %73, label %74, label %90
74: ; preds = %70
%75 = load double*, double** %8
%76 = load double, double* %75
%77 = load double*, double** %10
%78 = load double, double* %77
%79 = load double, double* %13
%80 = call double @llvm.fmuladd.f64(double %76, double %78, double %79)
store double %80, double* %13
%81 = load i64, i64* %9
%82 = load double*, double** %8
%83 = getelementptr inbounds double, double* %82, i64 %81
store double* %83, double** %8
%84 = load i64, i64* %11
%85 = load double*, double** %10
%86 = getelementptr inbounds double, double* %85, i64 %84
store double* %86, double** %10
br label %87
87: ; preds = %74
%88 = load i64, i64* %12
%89 = add nsw i64 %88, 1
store i64 %89, i64* %12
br label %70
90: ; preds = %70
br label %91
91: ; preds = %90, %46
%92 = load double, double* %13
store double %92, double* %6
br label %93
93: ; preds = %91, %22
%94 = load double, double* %6
ret double %94
}

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