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


x86 assembly (clang, -O0):
<my_daxpy>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x54(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x58(%rbp)
mov    -0x30(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x5c(%rbp)
cmpl   $0x0,-0x54(%rbp)
jg     1153 <my_daxpy+0x53>
jmp    126f <my_daxpy+0x16f>
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1170 <my_daxpy+0x70>
jp     1170 <my_daxpy+0x70>
jmp    126f <my_daxpy+0x16f>
mov    -0x18(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x18(%rbp)
mov    -0x28(%rbp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x28(%rbp)
movl   $0x1,-0x38(%rbp)
movl   $0x1,-0x3c(%rbp)
cmpl   $0x0,-0x58(%rbp)
jge    11b2 <my_daxpy+0xb2>
xor    %eax,%eax
sub    -0x54(%rbp),%eax
add    $0x1,%eax
imul   -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x5c(%rbp)
jge    11ce <my_daxpy+0xce>
xor    %eax,%eax
sub    -0x54(%rbp),%eax
add    $0x1,%eax
imul   -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x68(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x70(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x54(%rbp),%eax
jg     126a <my_daxpy+0x16a>
movsd  -0x48(%rbp),%xmm0
mov    -0x68(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x70(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x58(%rbp),%ecx
mov    -0x68(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x68(%rbp)
mov    -0x5c(%rbp),%ecx
mov    -0x70(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x70(%rbp)
jmp    1205 <my_daxpy+0x105>
jmp    126f <my_daxpy+0x16f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<my_daxpy>:
push   %rbp
push   %r15
push   %r14
push   %rbx
mov    (%rdi),%r14d
test   %r14d,%r14d
jle    124c <my_daxpy+0x14c>
movsd  (%rsi),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1126 <my_daxpy+0x26>
jnp    124c <my_daxpy+0x14c>
movslq (%rcx),%r10
movslq (%r9),%r11
mov    $0x1,%eax
mov    $0x1,%ecx
sub    %r14d,%ecx
mov    %r11d,%esi
imul   %ecx,%esi
imul   %r10d,%ecx
add    $0x1,%ecx
test   %r10d,%r10d
cmovns %eax,%ecx
movslq %ecx,%rcx
add    $0x1,%esi
test   %r11d,%r11d
cmovns %eax,%esi
movslq %esi,%rax
lea    -0x8(%rdx,%rcx,8),%rsi
lea    (%r8,%rax,8),%rax
add    $0xfffffffffffffff8,%rax
lea    -0x1(%r14),%ecx
mov    %r14d,%r15d
and    $0x3,%r15d
cmp    $0x3,%ecx
jb     121b <my_daxpy+0x11b>
and    $0xfffffffc,%r14d
mov    %r11,%r8
shl    $0x5,%r8
mov    %r10,%r9
shl    $0x5,%r9
lea    0x0(,%r10,8),%rcx
lea    0x0(,%r11,8),%rdx
xchg   %ax,%ax
mov    %rax,%rbx
mov    %rsi,%rdi
movsd  (%rsi),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rax),%xmm1
movsd  %xmm1,(%rax)
add    %rcx,%rsi
add    %rdx,%rax
movsd  (%rdi,%rcx,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rbx,%rdx,1),%xmm1
movsd  %xmm1,(%rbx,%rdx,1)
lea    (%rax,%rdx,1),%rbp
movsd  (%rcx,%rsi,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdx,%rax,1),%xmm1
add    %rcx,%rsi
movsd  %xmm1,(%rdx,%rax,1)
lea    (%rdx,%rbp,1),%rax
movsd  (%rcx,%rsi,1),%xmm1
mulsd  %xmm0,%xmm1
add    %rcx,%rsi
addsd  (%rdx,%rbp,1),%xmm1
movsd  %xmm1,(%rdx,%rbp,1)
add    %rcx,%rsi
add    %rdx,%rax
add    $0xfffffffc,%r14d
jne    11a0 <my_daxpy+0xa0>
add    %r8,%rbx
add    %r9,%rdi
mov    %rbx,%rax
mov    %rdi,%rsi
test   %r15d,%r15d
je     124c <my_daxpy+0x14c>
shl    $0x3,%r11
shl    $0x3,%r10
nopl   0x0(%rax,%rax,1)
movsd  (%rsi),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rax),%xmm1
movsd  %xmm1,(%rax)
add    %r11,%rax
add    %r10,%rsi
add    $0xffffffff,%r15d
jne    1230 <my_daxpy+0x130>
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

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