void DVaxpy13(int n, double y0[], double alpha[], double x0[], double x1[],
              double x2[]) {
  register double a00, a01, a02;
  register double x0i, x1i, x2i;
  int ii;
  a00 = alpha[0];
  a01 = alpha[1];
  a02 = alpha[2];
  for (ii = 0; ii < n; ii++) {
    x0i = x0[ii];
    x1i = x1[ii];
    x2i = x2[ii];
    y0[ii] += a00 * x0i + a01 * x1i + a02 * x2i;
  }
  return;
}


x86 assembly (clang, -O0):
<DVaxpy13>:
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    %r9,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11db <DVaxpy13+0xdb>
mov    -0x20(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    114b <DVaxpy13+0x4b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<DVaxpy13>:
push   %rbp
push   %r14
push   %rbx
test   %edi,%edi
jle    1302 <DVaxpy13+0x202>
movsd  (%rdx),%xmm8
movsd  0x8(%rdx),%xmm9
movsd  0x10(%rdx),%xmm10
mov    %edi,%r10d
cmp    $0x4,%edi
jae    112c <DVaxpy13+0x2c>
xor    %edx,%edx
jmp    123f <DVaxpy13+0x13f>
lea    (%rcx,%r10,8),%rax
cmp    %rsi,%rax
seta   %r14b
lea    (%rsi,%r10,8),%rdx
cmp    %rcx,%rdx
seta   %bl
lea    (%r8,%r10,8),%rax
cmp    %rsi,%rax
seta   %al
lea    (%r9,%r10,8),%rdi
cmp    %r8,%rdx
seta   %r11b
cmp    %rsi,%rdi
seta   %dil
cmp    %r9,%rdx
seta   %bpl
xor    %edx,%edx
test   %bl,%r14b
jne    123f <DVaxpy13+0x13f>
and    %r11b,%al
jne    123f <DVaxpy13+0x13f>
and    %bpl,%dil
jne    123f <DVaxpy13+0x13f>
mov    %r10d,%edx
and    $0xfffffffc,%edx
movapd %xmm9,%xmm11
unpcklpd %xmm9,%xmm11
movapd %xmm8,%xmm4
unpcklpd %xmm8,%xmm4
movapd %xmm10,%xmm5
unpcklpd %xmm10,%xmm5
lea    0x0(,%r10,8),%rdi
and    $0xffffffffffffffe0,%rdi
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd (%rcx,%rax,1),%xmm6
movupd 0x10(%rcx,%rax,1),%xmm7
movupd (%r8,%rax,1),%xmm0
movupd 0x10(%r8,%rax,1),%xmm1
movupd (%r9,%rax,1),%xmm2
movupd 0x10(%r9,%rax,1),%xmm3
mulpd  %xmm11,%xmm0
mulpd  %xmm11,%xmm1
mulpd  %xmm4,%xmm6
addpd  %xmm0,%xmm6
mulpd  %xmm4,%xmm7
addpd  %xmm1,%xmm7
mulpd  %xmm5,%xmm2
addpd  %xmm6,%xmm2
mulpd  %xmm5,%xmm3
addpd  %xmm7,%xmm3
movupd (%rsi,%rax,1),%xmm0
addpd  %xmm2,%xmm0
movupd 0x10(%rsi,%rax,1),%xmm1
addpd  %xmm3,%xmm1
movupd %xmm0,(%rsi,%rax,1)
movupd %xmm1,0x10(%rsi,%rax,1)
add    $0x20,%rax
cmp    %rax,%rdi
jne    11c0 <DVaxpy13+0xc0>
cmp    %r10,%rdx
je     1302 <DVaxpy13+0x202>
mov    %rdx,%rax
not    %rax
test   $0x1,%r10b
je     1281 <DVaxpy13+0x181>
movsd  (%r8,%rdx,8),%xmm0
mulsd  %xmm9,%xmm0
movsd  (%rcx,%rdx,8),%xmm1
mulsd  %xmm8,%xmm1
addsd  %xmm0,%xmm1
movsd  (%r9,%rdx,8),%xmm0
mulsd  %xmm10,%xmm0
addsd  %xmm1,%xmm0
addsd  (%rsi,%rdx,8),%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
or     $0x1,%rdx
add    %r10,%rax
je     1302 <DVaxpy13+0x202>
nopw   0x0(%rax,%rax,1)
movsd  (%r8,%rdx,8),%xmm0
mulsd  %xmm9,%xmm0
movsd  (%rcx,%rdx,8),%xmm1
mulsd  %xmm8,%xmm1
addsd  %xmm0,%xmm1
movsd  (%r9,%rdx,8),%xmm0
mulsd  %xmm10,%xmm0
addsd  %xmm1,%xmm0
addsd  (%rsi,%rdx,8),%xmm0
movsd  %xmm0,(%rsi,%rdx,8)
movsd  0x8(%r8,%rdx,8),%xmm0
mulsd  %xmm9,%xmm0
movsd  0x8(%rcx,%rdx,8),%xmm1
mulsd  %xmm8,%xmm1
addsd  %xmm0,%xmm1
movsd  0x8(%r9,%rdx,8),%xmm0
mulsd  %xmm10,%xmm0
addsd  %xmm1,%xmm0
addsd  0x8(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x8(%rsi,%rdx,8)
add    $0x2,%rdx
cmp    %rdx,%r10
jne    1290 <DVaxpy13+0x190>
pop    %rbx
pop    %r14
pop    %rbp
ret

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

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