double my_ddot(int *_n, double *dx, int *_incx, double *dy, int *_incy) {
  register double dtemp;
  int i, ix, iy;
  int n = *_n, incx = *_incx, incy = *_incy;
  dtemp = 0.0;
  if (n <= 0)
    return ((double)dtemp);
  dx--;
  dy--;
  ix = 1;
  iy = 1;
  if (incx < 0)
    ix = (-n + 1) * incx + 1;
  if (incy < 0)
    iy = (-n + 1) * incy + 1;
  {
    double *xptr, *yptr;
    for (i = 1, xptr = dx + ix, yptr = dy + iy; i <= n;
         i++, xptr += incx, yptr += incy)
      dtemp += (*yptr) * (*xptr);
    return (dtemp);
  }
}


x86 assembly (clang, -O0):
<my_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)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
cmpl   $0x0,-0x48(%rbp)
jg     1154 <my_ddot+0x54>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    124f <my_ddot+0x14f>
mov    -0x18(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x18(%rbp)
mov    -0x28(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x28(%rbp)
movl   $0x1,-0x40(%rbp)
movl   $0x1,-0x44(%rbp)
cmpl   $0x0,-0x4c(%rbp)
jge    1196 <my_ddot+0x96>
xor    %eax,%eax
sub    -0x48(%rbp),%eax
add    $0x1,%eax
imul   -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x0,-0x50(%rbp)
jge    11b2 <my_ddot+0xb2>
xor    %eax,%eax
sub    -0x48(%rbp),%eax
add    $0x1,%eax
imul   -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
movl   $0x1,-0x3c(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x48(%rbp),%eax
jg     1245 <my_ddot+0x145>
mov    -0x60(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x4c(%rbp),%ecx
mov    -0x58(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x50(%rbp),%ecx
mov    -0x60(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
jmp    11df <my_ddot+0xdf>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<my_ddot>:
push   %rbp
push   %r14
push   %rbx
mov    (%rdi),%eax
test   %eax,%eax
jle    116c <my_ddot+0x6c>
movslq (%rdx),%r9
movslq (%r8),%r11
mov    $0x1,%edx
mov    $0x1,%edi
sub    %eax,%edi
mov    %r11d,%ebp
imul   %edi,%ebp
imul   %r9d,%edi
add    $0x1,%edi
test   %r9d,%r9d
cmovns %edx,%edi
movslq %edi,%rdi
add    $0x1,%ebp
test   %r11d,%r11d
cmovns %edx,%ebp
movslq %ebp,%rdx
lea    -0x8(%rsi,%rdi,8),%rdi
lea    (%rcx,%rdx,8),%rsi
add    $0xfffffffffffffff8,%rsi
lea    -0x1(%rax),%ecx
mov    %eax,%r14d
and    $0x3,%r14d
cmp    $0x3,%ecx
jae    1175 <my_ddot+0x75>
xorpd  %xmm0,%xmm0
test   %r14d,%r14d
jne    120b <my_ddot+0x10b>
jmp    1238 <my_ddot+0x138>
xorpd  %xmm0,%xmm0
jmp    1238 <my_ddot+0x138>
and    $0xfffffffc,%eax
mov    %r11,%r8
shl    $0x5,%r8
mov    %r9,%r10
shl    $0x5,%r10
lea    0x0(,%r11,8),%rdx
lea    0x0(,%r9,8),%rbx
xorpd  %xmm0,%xmm0
nopw   0x0(%rax,%rax,1)
mov    %rsi,%rcx
mov    %rdi,%rbp
movsd  (%rsi),%xmm1
mulsd  (%rdi),%xmm1
addsd  %xmm0,%xmm1
add    %rdx,%rsi
movsd  (%rcx,%rdx,1),%xmm0
add    %rbx,%rdi
mulsd  0x0(%rbp,%rbx,1),%xmm0
addsd  %xmm1,%xmm0
movsd  (%rdx,%rsi,1),%xmm1
mulsd  (%rbx,%rdi,1),%xmm1
add    %rdx,%rsi
add    %rbx,%rdi
addsd  %xmm0,%xmm1
movsd  (%rdx,%rsi,1),%xmm0
add    %rdx,%rsi
mulsd  (%rbx,%rdi,1),%xmm0
add    %rbx,%rdi
addsd  %xmm1,%xmm0
add    %rdx,%rsi
add    %rbx,%rdi
add    $0xfffffffc,%eax
jne    11a0 <my_ddot+0xa0>
add    %r8,%rcx
add    %r10,%rbp
mov    %rcx,%rsi
mov    %rbp,%rdi
test   %r14d,%r14d
je     1238 <my_ddot+0x138>
shl    $0x3,%r11
shl    $0x3,%r9
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rsi),%xmm1
mulsd  (%rdi),%xmm1
addsd  %xmm1,%xmm0
add    %r11,%rsi
add    %r9,%rdi
add    $0xffffffff,%r14d
jne    1220 <my_ddot+0x120>
pop    %rbx
pop    %r14
pop    %rbp
ret

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

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