void kernel_dsyttmm_lu_nt_2x2_lib4(int kmax, double *A, double *C) {
  const int lda = 4;
  const int ldc = 4;
  int k;
  double a_0, a_1, c_00 = 0, c_10 = 0, c_11 = 0;
  for (k = 0; k < kmax - 4; k += 4) {
    a_0 = A[0 + lda * 0];
    a_1 = A[1 + lda * 0];
    c_00 += a_0 * a_0;
    c_10 += a_1 * a_0;
    c_11 += a_1 * a_1;
    a_0 = A[0 + lda * 1];
    a_1 = A[1 + lda * 1];
    c_00 += a_0 * a_0;
    c_10 += a_1 * a_0;
    c_11 += a_1 * a_1;
    a_0 = A[0 + lda * 2];
    a_1 = A[1 + lda * 2];
    c_00 += a_0 * a_0;
    c_10 += a_1 * a_0;
    c_11 += a_1 * a_1;
    a_0 = A[0 + lda * 3];
    a_1 = A[1 + lda * 3];
    c_00 += a_0 * a_0;
    c_10 += a_1 * a_0;
    c_11 += a_1 * a_1;
    A += 16;
  }
  for (; k < kmax - 2; k += 2) {
    a_0 = A[0 + lda * 0];
    a_1 = A[1 + lda * 0];
    c_00 += a_0 * a_0;
    c_10 += a_1 * a_0;
    c_11 += a_1 * a_1;
    a_0 = A[0 + lda * 1];
    a_1 = A[1 + lda * 1];
    c_00 += a_0 * a_0;
    c_10 += a_1 * a_0;
    c_11 += a_1 * a_1;
    A += 8;
  }
  a_0 = A[0 + lda * 0];
  a_1 = A[1 + lda * 0];
  c_00 += a_0 * a_0;
  c_10 += a_1 * a_0;
  c_11 += a_1 * a_1;
  a_1 = A[1 + lda * 1];
  c_11 += a_1 * a_1;
  C[0 + ldc * 0] = c_00;
  C[1 + ldc * 0] = c_10;
  C[1 + ldc * 1] = c_11;
}


x86 assembly (clang, -O0):
<kernel_dsyttmm_lu_nt_2x2_lib4>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x4,-0x1c(%rbp)
movl   $0x4,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x4,%ecx
cmp    %ecx,%eax
jge    1328 <kernel_dsyttmm_lu_nt_2x2_lib4+0x228>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x20(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x40(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x48(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x60(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x68(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
add    $0x80,%rax
mov    %rax,-0x10(%rbp)
mov    -0x24(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x24(%rbp)
jmp    113c <kernel_dsyttmm_lu_nt_2x2_lib4+0x3c>
jmp    132d <kernel_dsyttmm_lu_nt_2x2_lib4+0x22d>
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    1437 <kernel_dsyttmm_lu_nt_2x2_lib4+0x337>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x20(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
add    $0x40,%rax
mov    %rax,-0x10(%rbp)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    132d <kernel_dsyttmm_lu_nt_2x2_lib4+0x22d>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,0x28(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<kernel_dsyttmm_lu_nt_2x2_lib4>:
cmp    $0x5,%edi
jl     11c2 <kernel_dsyttmm_lu_nt_2x2_lib4+0xc2>
lea    -0x4(%rdi),%eax
xorpd  %xmm3,%xmm3
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
movupd (%rsi),%xmm4
movupd 0x20(%rsi),%xmm1
movupd 0x40(%rsi),%xmm2
movapd %xmm4,%xmm5
unpcklpd %xmm4,%xmm5
mulpd  %xmm4,%xmm5
unpckhpd %xmm4,%xmm4
mulsd  %xmm4,%xmm4
addsd  %xmm0,%xmm4
movupd 0x60(%rsi),%xmm0
addpd  %xmm3,%xmm5
movapd %xmm1,%xmm3
unpcklpd %xmm1,%xmm3
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
unpckhpd %xmm1,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm4,%xmm1
movapd %xmm2,%xmm4
unpcklpd %xmm2,%xmm4
mulpd  %xmm2,%xmm4
addpd  %xmm3,%xmm4
unpckhpd %xmm2,%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm1,%xmm3
addpd  %xmm4,%xmm3
unpckhpd %xmm0,%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
sub    $0xffffffffffffff80,%rsi
add    $0x4,%ecx
cmp    %eax,%ecx
jl     1120 <kernel_dsyttmm_lu_nt_2x2_lib4+0x20>
addpd  %xmm4,%xmm1
lea    -0x1(%rdi),%eax
and    $0xfffffffc,%eax
jmp    11cc <kernel_dsyttmm_lu_nt_2x2_lib4+0xcc>
xorpd  %xmm1,%xmm1
xorpd  %xmm0,%xmm0
xor    %eax,%eax
add    $0xfffffffe,%edi
cmp    %edi,%eax
jge    1238 <kernel_dsyttmm_lu_nt_2x2_lib4+0x138>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd (%rsi),%xmm4
movapd %xmm4,%xmm2
unpcklpd %xmm4,%xmm2
mulpd  %xmm4,%xmm2
unpckhpd %xmm4,%xmm4
mulsd  %xmm4,%xmm4
addsd  %xmm0,%xmm4
movupd 0x20(%rsi),%xmm0
addpd  %xmm1,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
mulpd  %xmm0,%xmm3
movapd %xmm3,%xmm1
addpd  %xmm2,%xmm1
unpckhpd %xmm0,%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm4,%xmm0
add    $0x40,%rsi
add    $0x2,%eax
cmp    %edi,%eax
jl     11e0 <kernel_dsyttmm_lu_nt_2x2_lib4+0xe0>
addpd  %xmm2,%xmm3
movapd %xmm3,%xmm1
movupd (%rsi),%xmm2
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
mulpd  %xmm2,%xmm3
addpd  %xmm1,%xmm3
unpckhpd %xmm2,%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
movsd  0x28(%rsi),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
movupd %xmm3,(%rdx)
movsd  %xmm0,0x28(%rdx)
ret

LLVM IR (-O0):
define void @kernel_dsyttmm_lu_nt_2x2_lib4(i32 %0, double* %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
store i32 %0, i32* %4
store double* %1, double** %5
store double* %2, double** %6
store i32 4, i32* %7
store i32 4, i32* %8
store double 0.000000e+00, double* %12
store double 0.000000e+00, double* %13
store double 0.000000e+00, double* %14
store i32 0, i32* %9
br label %15
15: ; preds = %95, %3
%16 = load i32, i32* %9
%17 = load i32, i32* %4
%18 = sub nsw i32 %17, 4
%19 = icmp slt i32 %16, %18
br i1 %19, label %20, label %98
20: ; preds = %15
%21 = load double*, double** %5
%22 = getelementptr inbounds double, double* %21, i64 0
%23 = load double, double* %22
store double %23, double* %10
%24 = load double*, double** %5
%25 = getelementptr inbounds double, double* %24, i64 1
%26 = load double, double* %25
store double %26, double* %11
%27 = load double, double* %10
%28 = load double, double* %10
%29 = load double, double* %12
%30 = call double @llvm.fmuladd.f64(double %27, double %28, double %29)
store double %30, double* %12
%31 = load double, double* %11
%32 = load double, double* %10
%33 = load double, double* %13
%34 = call double @llvm.fmuladd.f64(double %31, double %32, double %33)
store double %34, double* %13
%35 = load double, double* %11
%36 = load double, double* %11
%37 = load double, double* %14
%38 = call double @llvm.fmuladd.f64(double %35, double %36, double %37)
store double %38, double* %14
%39 = load double*, double** %5
%40 = getelementptr inbounds double, double* %39, i64 4
%41 = load double, double* %40
store double %41, double* %10
%42 = load double*, double** %5
%43 = getelementptr inbounds double, double* %42, i64 5
%44 = load double, double* %43
store double %44, double* %11
%45 = load double, double* %10
%46 = load double, double* %10
%47 = load double, double* %12
%48 = call double @llvm.fmuladd.f64(double %45, double %46, double %47)
store double %48, double* %12
%49 = load double, double* %11
%50 = load double, double* %10
%51 = load double, double* %13
%52 = call double @llvm.fmuladd.f64(double %49, double %50, double %51)
store double %52, double* %13
%53 = load double, double* %11
%54 = load double, double* %11
%55 = load double, double* %14
%56 = call double @llvm.fmuladd.f64(double %53, double %54, double %55)
store double %56, double* %14
%57 = load double*, double** %5
%58 = getelementptr inbounds double, double* %57, i64 8
%59 = load double, double* %58
store double %59, double* %10
%60 = load double*, double** %5
%61 = getelementptr inbounds double, double* %60, i64 9
%62 = load double, double* %61
store double %62, double* %11
%63 = load double, double* %10
%64 = load double, double* %10
%65 = load double, double* %12
%66 = call double @llvm.fmuladd.f64(double %63, double %64, double %65)
store double %66, double* %12
%67 = load double, double* %11
%68 = load double, double* %10
%69 = load double, double* %13
%70 = call double @llvm.fmuladd.f64(double %67, double %68, double %69)
store double %70, double* %13
%71 = load double, double* %11
%72 = load double, double* %11
%73 = load double, double* %14
%74 = call double @llvm.fmuladd.f64(double %71, double %72, double %73)
store double %74, double* %14
%75 = load double*, double** %5
%76 = getelementptr inbounds double, double* %75, i64 12
%77 = load double, double* %76
store double %77, double* %10
%78 = load double*, double** %5
%79 = getelementptr inbounds double, double* %78, i64 13
%80 = load double, double* %79
store double %80, double* %11
%81 = load double, double* %10
%82 = load double, double* %10
%83 = load double, double* %12
%84 = call double @llvm.fmuladd.f64(double %81, double %82, double %83)
store double %84, double* %12
%85 = load double, double* %11
%86 = load double, double* %10
%87 = load double, double* %13
%88 = call double @llvm.fmuladd.f64(double %85, double %86, double %87)
store double %88, double* %13
%89 = load double, double* %11
%90 = load double, double* %11
%91 = load double, double* %14
%92 = call double @llvm.fmuladd.f64(double %89, double %90, double %91)
store double %92, double* %14
%93 = load double*, double** %5
%94 = getelementptr inbounds double, double* %93, i64 16
store double* %94, double** %5
br label %95
95: ; preds = %20
%96 = load i32, i32* %9
%97 = add nsw i32 %96, 4
store i32 %97, i32* %9
br label %15
98: ; preds = %15
br label %99
99: ; preds = %143, %98
%100 = load i32, i32* %9
%101 = load i32, i32* %4
%102 = sub nsw i32 %101, 2
%103 = icmp slt i32 %100, %102
br i1 %103, label %104, label %146
104: ; preds = %99
%105 = load double*, double** %5
%106 = getelementptr inbounds double, double* %105, i64 0
%107 = load double, double* %106
store double %107, double* %10
%108 = load double*, double** %5
%109 = getelementptr inbounds double, double* %108, i64 1
%110 = load double, double* %109
store double %110, double* %11
%111 = load double, double* %10
%112 = load double, double* %10
%113 = load double, double* %12
%114 = call double @llvm.fmuladd.f64(double %111, double %112, double %113)
store double %114, double* %12
%115 = load double, double* %11
%116 = load double, double* %10
%117 = load double, double* %13
%118 = call double @llvm.fmuladd.f64(double %115, double %116, double %117)
store double %118, double* %13
%119 = load double, double* %11
%120 = load double, double* %11
%121 = load double, double* %14
%122 = call double @llvm.fmuladd.f64(double %119, double %120, double %121)
store double %122, double* %14
%123 = load double*, double** %5
%124 = getelementptr inbounds double, double* %123, i64 4
%125 = load double, double* %124
store double %125, double* %10
%126 = load double*, double** %5
%127 = getelementptr inbounds double, double* %126, i64 5
%128 = load double, double* %127
store double %128, double* %11
%129 = load double, double* %10
%130 = load double, double* %10
%131 = load double, double* %12
%132 = call double @llvm.fmuladd.f64(double %129, double %130, double %131)
store double %132, double* %12
%133 = load double, double* %11
%134 = load double, double* %10
%135 = load double, double* %13
%136 = call double @llvm.fmuladd.f64(double %133, double %134, double %135)
store double %136, double* %13
%137 = load double, double* %11
%138 = load double, double* %11
%139 = load double, double* %14
%140 = call double @llvm.fmuladd.f64(double %137, double %138, double %139)
store double %140, double* %14
%141 = load double*, double** %5
%142 = getelementptr inbounds double, double* %141, i64 8
store double* %142, double** %5
br label %143
143: ; preds = %104
%144 = load i32, i32* %9
%145 = add nsw i32 %144, 2
store i32 %145, i32* %9
br label %99
146: ; preds = %99
%147 = load double*, double** %5
%148 = getelementptr inbounds double, double* %147, i64 0
%149 = load double, double* %148
store double %149, double* %10
%150 = load double*, double** %5
%151 = getelementptr inbounds double, double* %150, i64 1
%152 = load double, double* %151
store double %152, double* %11
%153 = load double, double* %10
%154 = load double, double* %10
%155 = load double, double* %12
%156 = call double @llvm.fmuladd.f64(double %153, double %154, double %155)
store double %156, double* %12
%157 = load double, double* %11
%158 = load double, double* %10
%159 = load double, double* %13
%160 = call double @llvm.fmuladd.f64(double %157, double %158, double %159)
store double %160, double* %13
%161 = load double, double* %11
%162 = load double, double* %11
%163 = load double, double* %14
%164 = call double @llvm.fmuladd.f64(double %161, double %162, double %163)
store double %164, double* %14
%165 = load double*, double** %5
%166 = getelementptr inbounds double, double* %165, i64 5
%167 = load double, double* %166
store double %167, double* %11
%168 = load double, double* %11
%169 = load double, double* %11
%170 = load double, double* %14
%171 = call double @llvm.fmuladd.f64(double %168, double %169, double %170)
store double %171, double* %14
%172 = load double, double* %12
%173 = load double*, double** %6
%174 = getelementptr inbounds double, double* %173, i64 0
store double %172, double* %174
%175 = load double, double* %13
%176 = load double*, double** %6
%177 = getelementptr inbounds double, double* %176, i64 1
store double %175, double* %177
%178 = load double, double* %14
%179 = load double*, double** %6
%180 = getelementptr inbounds double, double* %179, i64 5
store double %178, double* %180
ret void
}

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