void d_set_mat(int row, int col, double alpha, double *A, int lda) {
  int i, j;
  for (j = 0; j < col; j++) {
    i = 0;
    for (; i < row - 3; i += 4) {
      A[i + 0 + j * lda] = alpha;
      A[i + 1 + j * lda] = alpha;
      A[i + 2 + j * lda] = alpha;
      A[i + 3 + j * lda] = alpha;
    }
    for (; i < row; i++) {
      A[i + j * lda] = alpha;
    }
  }
}


x86 assembly (clang, -O0):
<d_set_mat>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    121e <d_set_mat+0x11e>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x3,%ecx
cmp    %ecx,%eax
jge    11cf <d_set_mat+0xcf>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x0,%ecx
mov    -0x24(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
mov    -0x24(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x2,%ecx
mov    -0x24(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x3,%ecx
mov    -0x24(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x20(%rbp)
jmp    1130 <d_set_mat+0x30>
jmp    11d4 <d_set_mat+0xd4>
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    120b <d_set_mat+0x10b>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    -0x24(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11d4 <d_set_mat+0xd4>
jmp    1210 <d_set_mat+0x110>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111d <d_set_mat+0x1d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<d_set_mat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %esi,%esi
jle    1365 <d_set_mat+0x265>
movabs $0x400000000,%r12
lea    -0x3(%rdi),%eax
movslq %eax,%r8
movslq %ecx,%rcx
mov    %esi,%eax
mov    %rax,-0x8(%rsp)
mov    %edi,%r15d
movaps %xmm0,%xmm1
movlhps %xmm0,%xmm1
mov    %rcx,%rax
shl    $0x20,%rax
mov    %rax,-0x10(%rsp)
shl    $0x3,%rcx
mov    %rcx,-0x18(%rsp)
lea    0x70(%rdx),%r10
lea    0x10(%rdx),%rsi
movabs $0x300000000,%rcx
movabs $0x200000000,%rbx
movabs $0x100000000,%rax
mov    %rax,-0x30(%rsp)
xor    %r14d,%r14d
mov    %rdx,%r13
mov    %rdi,-0x28(%rsp)
jmp    11c0 <d_set_mat+0xc0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r14
mov    -0x10(%rsp),%rax
add    %rax,%rcx
add    %rax,%rbx
add    %rax,-0x30(%rsp)
mov    -0x18(%rsp),%rax
add    %rax,%r13
add    %rax,%r10
add    %rax,%rsi
cmp    -0x8(%rsp),%r14
mov    %r11,%rdi
je     1365 <d_set_mat+0x265>
mov    $0x0,%ebp
cmp    $0x4,%edi
jl     1240 <d_set_mat+0x140>
mov    %rdi,%r11
mov    %rsi,-0x40(%rsp)
mov    %r10,-0x38(%rsp)
mov    -0x30(%rsp),%rax
mov    %rbx,%r9
mov    %rcx,%rdi
mov    %rcx,%rsi
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
movsd  %xmm0,0x0(%r13,%rbp,8)
mov    %rax,%rcx
sar    $0x1d,%rcx
movsd  %xmm0,(%rdx,%rcx,1)
mov    %rbx,%rcx
sar    $0x1d,%rcx
movsd  %xmm0,(%rdx,%rcx,1)
mov    %rsi,%rcx
sar    $0x1d,%rcx
movsd  %xmm0,(%rdx,%rcx,1)
add    $0x4,%rbp
add    %r12,%rsi
add    %r12,%rbx
add    %r12,%rax
cmp    %r8,%rbp
jl     11f0 <d_set_mat+0xf0>
mov    -0x38(%rsp),%r10
mov    -0x40(%rsp),%rsi
mov    %rdi,%rcx
mov    %r9,%rbx
mov    %r11,%rdi
mov    %rdi,%r11
cmp    %edi,%ebp
jge    1190 <d_set_mat+0x90>
mov    %rcx,%rax
mov    %ebp,%r9d
mov    %r15,%rcx
sub    %r9,%rcx
cmp    $0x4,%rcx
jae    1270 <d_set_mat+0x170>
mov    %rax,%rcx
jmp    1350 <d_set_mat+0x250>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %rbx,-0x40(%rsp)
mov    %rsi,%rbp
mov    %rcx,%rsi
and    $0xfffffffffffffffc,%rsi
mov    %rsi,-0x20(%rsp)
add    $0xfffffffffffffffc,%rsi
mov    %rsi,%r11
shr    $0x2,%r11
add    $0x1,%r11
mov    %r11d,%edi
and    $0x3,%edi
cmp    $0xc,%rsi
mov    %r10,-0x38(%rsp)
jae    12a8 <d_set_mat+0x1a8>
xor    %esi,%esi
jmp    12f1 <d_set_mat+0x1f1>
and    $0xfffffffffffffffc,%r11
lea    (%r10,%r9,8),%rbx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups %xmm1,-0x70(%rbx,%rsi,8)
movups %xmm1,-0x60(%rbx,%rsi,8)
movups %xmm1,-0x50(%rbx,%rsi,8)
movups %xmm1,-0x40(%rbx,%rsi,8)
movups %xmm1,-0x30(%rbx,%rsi,8)
movups %xmm1,-0x20(%rbx,%rsi,8)
movups %xmm1,-0x10(%rbx,%rsi,8)
movups %xmm1,(%rbx,%rsi,8)
add    $0x10,%rsi
add    $0xfffffffffffffffc,%r11
jne    12c0 <d_set_mat+0x1c0>
test   %rdi,%rdi
je     1322 <d_set_mat+0x222>
add    %r9,%rsi
lea    0x0(,%rsi,8),%rsi
add    %rbp,%rsi
shl    $0x5,%rdi
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movups %xmm1,-0x10(%rsi,%rbx,1)
movups %xmm1,(%rsi,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%rdi
jne    1310 <d_set_mat+0x210>
mov    -0x20(%rsp),%rdi
cmp    %rdi,%rcx
mov    -0x28(%rsp),%r11
mov    -0x38(%rsp),%r10
mov    %rbp,%rsi
mov    %rax,%rcx
mov    -0x40(%rsp),%rbx
je     1190 <d_set_mat+0x90>
add    %rdi,%r9
nopl   0x0(%rax,%rax,1)
movsd  %xmm0,0x0(%r13,%r9,8)
add    $0x1,%r9
cmp    %r9,%r15
jne    1350 <d_set_mat+0x250>
jmp    1190 <d_set_mat+0x90>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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