void DVaxpy22(int n, double y0[], double y1[], double alpha[], double x0[],
              double x1[]) {
  register double a00, a01, a10, a11;
  register double x0i, x1i;
  int ii;
  a00 = alpha[0];
  a01 = alpha[1];
  a10 = alpha[2];
  a11 = alpha[3];
  for (ii = 0; ii < n; ii++) {
    x0i = x0[ii];
    x1i = x1[ii];
    y0[ii] += a00 * x0i + a01 * x1i;
    y1[ii] += a10 * x0i + a11 * x1i;
  }
  return;
}


x86 assembly (clang, -O0):
<DVaxpy22>:
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    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11f3 <DVaxpy22+0xf3>
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),%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%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    1159 <DVaxpy22+0x59>
pop    %rbp
ret

x86 assembly (clang, -O2):
<DVaxpy22>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
jle    1299 <DVaxpy22+0x199>
movsd  (%rcx),%xmm0
movsd  0x8(%rcx),%xmm1
movsd  0x10(%rcx),%xmm2
movsd  0x18(%rcx),%xmm3
mov    %edi,%eax
cmp    $0x1,%edi
jne    1133 <DVaxpy22+0x33>
xor    %ecx,%ecx
jmp    1250 <DVaxpy22+0x150>
lea    (%rdx,%rax,8),%rcx
cmp    %rsi,%rcx
seta   -0x1(%rsp)
lea    (%rsi,%rax,8),%rbx
cmp    %rdx,%rbx
seta   %r12b
lea    (%r8,%rax,8),%rbp
cmp    %rsi,%rbp
seta   %dil
cmp    %r8,%rbx
seta   %r15b
lea    (%r9,%rax,8),%r14
cmp    %rsi,%r14
seta   %r10b
cmp    %r9,%rbx
seta   %r13b
cmp    %rdx,%rbp
seta   %r11b
cmp    %r8,%rcx
seta   %bpl
cmp    %rdx,%r14
seta   %bl
cmp    %r9,%rcx
seta   %r14b
xor    %ecx,%ecx
test   %r12b,-0x1(%rsp)
jne    1250 <DVaxpy22+0x150>
and    %r15b,%dil
jne    1250 <DVaxpy22+0x150>
and    %r13b,%r10b
jne    1250 <DVaxpy22+0x150>
and    %bpl,%r11b
jne    1250 <DVaxpy22+0x150>
and    %r14b,%bl
jne    1250 <DVaxpy22+0x150>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movapd %xmm1,%xmm8
unpcklpd %xmm1,%xmm8
movapd %xmm0,%xmm9
unpcklpd %xmm0,%xmm9
movapd %xmm3,%xmm10
unpcklpd %xmm3,%xmm10
movapd %xmm2,%xmm11
unpcklpd %xmm2,%xmm11
xor    %edi,%edi
nopl   0x0(%rax)
movupd (%r8,%rdi,8),%xmm4
movupd (%r9,%rdi,8),%xmm5
movapd %xmm8,%xmm6
mulpd  %xmm5,%xmm6
movapd %xmm9,%xmm7
mulpd  %xmm4,%xmm7
addpd  %xmm6,%xmm7
movupd (%rsi,%rdi,8),%xmm6
addpd  %xmm7,%xmm6
movupd %xmm6,(%rsi,%rdi,8)
mulpd  %xmm10,%xmm5
mulpd  %xmm11,%xmm4
addpd  %xmm5,%xmm4
movupd (%rdx,%rdi,8),%xmm5
addpd  %xmm4,%xmm5
movupd %xmm5,(%rdx,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rcx
jne    11f0 <DVaxpy22+0xf0>
cmp    %rax,%rcx
je     1299 <DVaxpy22+0x199>
nopw   0x0(%rax,%rax,1)
movsd  (%r8,%rcx,8),%xmm4
movsd  (%r9,%rcx,8),%xmm5
movapd %xmm1,%xmm6
mulsd  %xmm5,%xmm6
movapd %xmm0,%xmm7
mulsd  %xmm4,%xmm7
addsd  %xmm6,%xmm7
addsd  (%rsi,%rcx,8),%xmm7
movsd  %xmm7,(%rsi,%rcx,8)
mulsd  %xmm3,%xmm5
mulsd  %xmm2,%xmm4
addsd  %xmm5,%xmm4
addsd  (%rdx,%rcx,8),%xmm4
movsd  %xmm4,(%rdx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1250 <DVaxpy22+0x150>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @DVaxpy22(i32 %0, double* nocapture %1, double* nocapture %2, double* nocapture %3, double* nocapture %4, double* nocapture %5) {
%7 = load double, double* %3
%8 = getelementptr inbounds double, double* %3, i64 1
%9 = load double, double* %8
%10 = getelementptr inbounds double, double* %3, i64 2
%11 = load double, double* %10
%12 = getelementptr inbounds double, double* %3, i64 3
%13 = load double, double* %12
%14 = icmp sgt i32 %0, 0
br i1 %14, label %15, label %98
15: ; preds = %6
%16 = zext i32 %0 to i64
%17 = icmp eq i32 %0, 1
br i1 %17, label %78, label %18
18: ; preds = %15
%19 = getelementptr double, double* %1, i64 %16
%20 = getelementptr double, double* %2, i64 %16
%21 = getelementptr double, double* %4, i64 %16
%22 = getelementptr double, double* %5, i64 %16
%23 = icmp ugt double* %20, %1
%24 = icmp ugt double* %19, %2
%25 = and i1 %23, %24
%26 = icmp ugt double* %21, %1
%27 = icmp ugt double* %19, %4
%28 = and i1 %26, %27
%29 = or i1 %25, %28
%30 = icmp ugt double* %22, %1
%31 = icmp ugt double* %19, %5
%32 = and i1 %30, %31
%33 = or i1 %29, %32
%34 = icmp ugt double* %21, %2
%35 = icmp ugt double* %20, %4
%36 = and i1 %34, %35
%37 = or i1 %33, %36
%38 = icmp ugt double* %22, %2
%39 = icmp ugt double* %20, %5
%40 = and i1 %38, %39
%41 = or i1 %37, %40
br i1 %41, label %78, label %42
42: ; preds = %18
%43 = and i64 %16, 4294967294
%44 = insertelement <2 x double> poison, double %9, i64 0
%45 = shufflevector <2 x double> %44, <2 x double> poison, <2 x i32> zeroinitializer
%46 = insertelement <2 x double> poison, double %7, i64 0
%47 = shufflevector <2 x double> %46, <2 x double> poison, <2 x i32> zeroinitializer
%48 = insertelement <2 x double> poison, double %13, i64 0
%49 = shufflevector <2 x double> %48, <2 x double> poison, <2 x i32> zeroinitializer
%50 = insertelement <2 x double> poison, double %11, i64 0
%51 = shufflevector <2 x double> %50, <2 x double> poison, <2 x i32> zeroinitializer
br label %52
52: ; preds = %52, %42
%53 = phi i64 [ 0, %42 ], [ %74, %52 ]
%54 = getelementptr inbounds double, double* %4, i64 %53
%55 = bitcast double* %54 to <2 x double>*
%56 = load <2 x double>, <2 x double>* %55
%57 = getelementptr inbounds double, double* %5, i64 %53
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = fmul <2 x double> %45, %59
%61 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %47, <2 x double> %56, <2 x double> %60)
%62 = getelementptr inbounds double, double* %1, i64 %53
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = fadd <2 x double> %64, %61
%66 = bitcast double* %62 to <2 x double>*
store <2 x double> %65, <2 x double>* %66
%67 = fmul <2 x double> %49, %59
%68 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %51, <2 x double> %56, <2 x double> %67)
%69 = getelementptr inbounds double, double* %2, i64 %53
%70 = bitcast double* %69 to <2 x double>*
%71 = load <2 x double>, <2 x double>* %70
%72 = fadd <2 x double> %68, %71
%73 = bitcast double* %69 to <2 x double>*
store <2 x double> %72, <2 x double>* %73
%74 = add nuw i64 %53, 2
%75 = icmp eq i64 %74, %43
br i1 %75, label %76, label %52
76: ; preds = %52
%77 = icmp eq i64 %43, %16
br i1 %77, label %98, label %78
78: ; preds = %18, %15, %76
%79 = phi i64 [ 0, %18 ], [ 0, %15 ], [ %43, %76 ]
br label %80
80: ; preds = %78, %80
%81 = phi i64 [ %96, %80 ], [ %79, %78 ]
%82 = getelementptr inbounds double, double* %4, i64 %81
%83 = load double, double* %82
%84 = getelementptr inbounds double, double* %5, i64 %81
%85 = load double, double* %84
%86 = fmul double %9, %85
%87 = tail call double @llvm.fmuladd.f64(double %7, double %83, double %86)
%88 = getelementptr inbounds double, double* %1, i64 %81
%89 = load double, double* %88
%90 = fadd double %89, %87
store double %90, double* %88
%91 = fmul double %13, %85
%92 = tail call double @llvm.fmuladd.f64(double %11, double %83, double %91)
%93 = getelementptr inbounds double, double* %2, i64 %81
%94 = load double, double* %93
%95 = fadd double %92, %94
store double %95, double* %93
%96 = add nuw nsw i64 %81, 1
%97 = icmp eq i64 %96, %16
br i1 %97, label %98, label %80
98: ; preds = %80, %76, %6
ret void
}