void DVaxpy21(int n, double y0[], double y1[], double alpha[], double x0[]) {
  register double a00, a10;
  register double x0i;
  int ii;
  a00 = alpha[0];
  a10 = alpha[1];
  for (ii = 0; ii < n; ii++) {
    x0i = x0[ii];
    y0[ii] += a00 * x0i;
    y1[ii] += a10 * x0i;
  }
  return;
}


x86 assembly (clang, -O0):
<DVaxpy21>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11ad <DVaxpy21+0xad>
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1139 <DVaxpy21+0x39>
pop    %rbp
ret

x86 assembly (clang, -O2):
<DVaxpy21>:
push   %rbp
push   %r14
push   %rbx
test   %edi,%edi
jle    1304 <DVaxpy21+0x204>
movsd  (%rcx),%xmm0
movsd  0x8(%rcx),%xmm1
mov    %edi,%r10d
cmp    $0x1,%edi
jne    11cf <DVaxpy21+0xcf>
xor    %ecx,%ecx
mov    %rcx,%rax
not    %rax
test   $0x1,%r10b
je     1159 <DVaxpy21+0x59>
movsd  (%r8,%rcx,8),%xmm2
movapd %xmm0,%xmm3
mulsd  %xmm2,%xmm3
addsd  (%rsi,%rcx,8),%xmm3
movsd  %xmm3,(%rsi,%rcx,8)
mulsd  %xmm1,%xmm2
addsd  (%rdx,%rcx,8),%xmm2
movsd  %xmm2,(%rdx,%rcx,8)
or     $0x1,%rcx
add    %r10,%rax
je     1304 <DVaxpy21+0x204>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%r8,%rcx,8),%xmm2
movapd %xmm0,%xmm3
mulsd  %xmm2,%xmm3
addsd  (%rsi,%rcx,8),%xmm3
movsd  %xmm3,(%rsi,%rcx,8)
mulsd  %xmm1,%xmm2
addsd  (%rdx,%rcx,8),%xmm2
movsd  %xmm2,(%rdx,%rcx,8)
movsd  0x8(%r8,%rcx,8),%xmm2
movapd %xmm0,%xmm3
mulsd  %xmm2,%xmm3
addsd  0x8(%rsi,%rcx,8),%xmm3
movsd  %xmm3,0x8(%rsi,%rcx,8)
mulsd  %xmm1,%xmm2
addsd  0x8(%rdx,%rcx,8),%xmm2
movsd  %xmm2,0x8(%rdx,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r10
jne    1170 <DVaxpy21+0x70>
jmp    1304 <DVaxpy21+0x204>
lea    (%rdx,%r10,8),%rcx
cmp    %rsi,%rcx
seta   %r14b
lea    (%rsi,%r10,8),%rdi
cmp    %rdx,%rdi
seta   %bl
lea    (%r8,%r10,8),%rbp
cmp    %rsi,%rbp
seta   %al
cmp    %r8,%rdi
seta   %r9b
cmp    %rdx,%rbp
seta   %dil
cmp    %r8,%rcx
seta   %r11b
xor    %ecx,%ecx
test   %bl,%r14b
jne    1123 <DVaxpy21+0x23>
and    %r9b,%al
jne    1123 <DVaxpy21+0x23>
and    %r11b,%dil
jne    1123 <DVaxpy21+0x23>
mov    %r10d,%ecx
and    $0xfffffffe,%ecx
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
lea    -0x2(%rcx),%rax
mov    %rax,%r9
shr    %r9
add    $0x1,%r9
test   %rax,%rax
je     1309 <DVaxpy21+0x209>
mov    %r9,%rax
and    $0xfffffffffffffffe,%rax
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movupd (%r8,%rdi,8),%xmm4
movupd (%rsi,%rdi,8),%xmm5
movapd %xmm3,%xmm6
mulpd  %xmm4,%xmm6
addpd  %xmm5,%xmm6
movupd %xmm6,(%rsi,%rdi,8)
movupd (%rdx,%rdi,8),%xmm5
mulpd  %xmm2,%xmm4
addpd  %xmm5,%xmm4
movupd %xmm4,(%rdx,%rdi,8)
movupd 0x10(%r8,%rdi,8),%xmm4
movupd 0x10(%rsi,%rdi,8),%xmm5
movapd %xmm3,%xmm6
mulpd  %xmm4,%xmm6
addpd  %xmm5,%xmm6
movupd %xmm6,0x10(%rsi,%rdi,8)
movupd 0x10(%rdx,%rdi,8),%xmm5
mulpd  %xmm2,%xmm4
addpd  %xmm5,%xmm4
movupd %xmm4,0x10(%rdx,%rdi,8)
add    $0x4,%rdi
add    $0xfffffffffffffffe,%rax
jne    1260 <DVaxpy21+0x160>
test   $0x1,%r9b
je     12fb <DVaxpy21+0x1fb>
movupd (%r8,%rdi,8),%xmm4
movupd (%rsi,%rdi,8),%xmm5
mulpd  %xmm4,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm3,(%rsi,%rdi,8)
movupd (%rdx,%rdi,8),%xmm3
mulpd  %xmm4,%xmm2
addpd  %xmm3,%xmm2
movupd %xmm2,(%rdx,%rdi,8)
cmp    %r10,%rcx
jne    1123 <DVaxpy21+0x23>
pop    %rbx
pop    %r14
pop    %rbp
ret
xor    %edi,%edi
test   $0x1,%r9b
jne    12d1 <DVaxpy21+0x1d1>
jmp    12fb <DVaxpy21+0x1fb>

LLVM IR (-O0):
define void @DVaxpy21(i32 %0, double* %1, double* %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca i32
store i32 %0, i32* %6
store double* %1, double** %7
store double* %2, double** %8
store double* %3, double** %9
store double* %4, double** %10
%15 = load double*, double** %9
%16 = getelementptr inbounds double, double* %15, i64 0
%17 = load double, double* %16
store double %17, double* %11
%18 = load double*, double** %9
%19 = getelementptr inbounds double, double* %18, i64 1
%20 = load double, double* %19
store double %20, double* %12
store i32 0, i32* %14
br label %21
21: ; preds = %47, %5
%22 = load i32, i32* %14
%23 = load i32, i32* %6
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %50
25: ; preds = %21
%26 = load double*, double** %10
%27 = load i32, i32* %14
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double, double* %26, i64 %28
%30 = load double, double* %29
store double %30, double* %13
%31 = load double, double* %11
%32 = load double, double* %13
%33 = load double*, double** %7
%34 = load i32, i32* %14
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %33, i64 %35
%37 = load double, double* %36
%38 = call double @llvm.fmuladd.f64(double %31, double %32, double %37)
store double %38, double* %36
%39 = load double, double* %12
%40 = load double, double* %13
%41 = load double*, double** %8
%42 = load i32, i32* %14
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %41, i64 %43
%45 = load double, double* %44
%46 = call double @llvm.fmuladd.f64(double %39, double %40, double %45)
store double %46, double* %44
br label %47
47: ; preds = %25
%48 = load i32, i32* %14
%49 = add nsw i32 %48, 1
store i32 %49, i32* %14
br label %21
50: ; preds = %21
ret void
}

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