void kernel_dgead_1_0_lib4(int kmax, double alpha, double *A, double *B) {
  if (kmax <= 0)
    return;
  const int bs = 4;
  int k;
  for (k = 0; k < kmax - 3; k += 4) {
    B[0 + bs * 0] += alpha * A[0 + bs * 0];
    B[0 + bs * 1] += alpha * A[0 + bs * 1];
    B[0 + bs * 2] += alpha * A[0 + bs * 2];
    B[0 + bs * 3] += alpha * A[0 + bs * 3];
    A += 16;
    B += 16;
  }
  for (; k < kmax; k++) {
    B[0 + bs * 0] += alpha * A[0 + bs * 0];
    A += 4;
    B += 4;
  }
}


x86 assembly (clang, -O0):
<kernel_dgead_1_0_lib4>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
cmpl   $0x0,-0x4(%rbp)
jg     1123 <kernel_dgead_1_0_lib4+0x23>
jmp    1251 <kernel_dgead_1_0_lib4+0x151>
movl   $0x4,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x3,%ecx
cmp    %ecx,%eax
jge    11f9 <kernel_dgead_1_0_lib4+0xf9>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm2
mov    -0x20(%rbp),%rax
movsd  0x20(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rax)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x40(%rax),%xmm2
mov    -0x20(%rbp),%rax
movsd  0x40(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x40(%rax)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x60(%rax),%xmm2
mov    -0x20(%rbp),%rax
movsd  0x60(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x60(%rax)
mov    -0x18(%rbp),%rax
add    $0x80,%rax
mov    %rax,-0x18(%rbp)
mov    -0x20(%rbp),%rax
add    $0x80,%rax
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x28(%rbp)
jmp    1131 <kernel_dgead_1_0_lib4+0x31>
jmp    11fe <kernel_dgead_1_0_lib4+0xfe>
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1251 <kernel_dgead_1_0_lib4+0x151>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x18(%rbp)
mov    -0x20(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11fe <kernel_dgead_1_0_lib4+0xfe>
pop    %rbp
ret

x86 assembly (clang, -O2):
<kernel_dgead_1_0_lib4>:
test   %edi,%edi
jle    122e <kernel_dgead_1_0_lib4+0x12e>
xor    %r9d,%r9d
cmp    $0x4,%edi
jb     117f <kernel_dgead_1_0_lib4+0x7f>
lea    -0x3(%rdi),%eax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rsi),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdx),%xmm1
movsd  %xmm1,(%rdx)
movsd  0x20(%rsi),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x20(%rdx),%xmm1
movsd  %xmm1,0x20(%rdx)
movsd  0x40(%rsi),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x40(%rdx),%xmm1
movsd  %xmm1,0x40(%rdx)
movsd  0x60(%rsi),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x60(%rdx),%xmm1
movsd  %xmm1,0x60(%rdx)
sub    $0xffffffffffffff80,%rsi
sub    $0xffffffffffffff80,%rdx
add    $0x4,%ecx
cmp    %eax,%ecx
jl     1120 <kernel_dgead_1_0_lib4+0x20>
mov    %edi,%r9d
and    $0xfffffffc,%r9d
cmp    %edi,%r9d
jge    122e <kernel_dgead_1_0_lib4+0x12e>
mov    %r9d,%r8d
not    %r8d
add    %edi,%r8d
mov    %edi,%ecx
and    $0x3,%ecx
je     11c2 <kernel_dgead_1_0_lib4+0xc2>
neg    %ecx
xor    %eax,%eax
nopl   0x0(%rax)
movsd  (%rsi),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdx),%xmm1
movsd  %xmm1,(%rdx)
add    $0x20,%rsi
add    $0x20,%rdx
add    $0xffffffff,%eax
cmp    %eax,%ecx
jne    11a0 <kernel_dgead_1_0_lib4+0xa0>
sub    %eax,%r9d
cmp    $0x3,%r8d
jb     122e <kernel_dgead_1_0_lib4+0x12e>
sub    %r9d,%edi
xor    %eax,%eax
nopl   (%rax)
movsd  (%rsi,%rax,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdx,%rax,1),%xmm1
movsd  %xmm1,(%rdx,%rax,1)
movsd  0x20(%rsi,%rax,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x20(%rdx,%rax,1),%xmm1
movsd  %xmm1,0x20(%rdx,%rax,1)
movsd  0x40(%rsi,%rax,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x40(%rdx,%rax,1),%xmm1
movsd  %xmm1,0x40(%rdx,%rax,1)
movsd  0x60(%rsi,%rax,1),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x60(%rdx,%rax,1),%xmm1
movsd  %xmm1,0x60(%rdx,%rax,1)
sub    $0xffffffffffffff80,%rax
add    $0xfffffffc,%edi
jne    11d0 <kernel_dgead_1_0_lib4+0xd0>
ret

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

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