void blas_scal(const int N, const double ALPHA, double *X, const int INCX) {
  if (ALPHA != 0) {
    for (int i = 0; i < N; ++i) {
      X[i * INCX] *= ALPHA;
    }
  } else {
    for (int i = 0; i < N; ++i) {
      X[i * INCX] = 0;
    }
  }
}


x86 assembly (clang, -O0):
<blas_scal>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1130 <blas_scal+0x30>
jp     1130 <blas_scal+0x30>
jmp    1173 <blas_scal+0x73>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    116e <blas_scal+0x6e>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1137 <blas_scal+0x37>
jmp    11af <blas_scal+0xaf>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11aa <blas_scal+0xaa>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    117a <blas_scal+0x7a>
jmp    11af <blas_scal+0xaf>
pop    %rbp
ret

x86 assembly (clang, -O2):
<blas_scal>:
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    110c <blas_scal+0xc>
jnp    1131 <blas_scal+0x31>
test   %edi,%edi
jle    1270 <blas_scal+0x170>
movslq %edx,%r10
mov    %edi,%r8d
lea    -0x1(%r8),%rax
mov    %r8d,%ecx
and    $0x3,%ecx
cmp    $0x3,%rax
jae    115a <blas_scal+0x5a>
xor    %edi,%edi
jmp    11c6 <blas_scal+0xc6>
test   %edi,%edi
jle    1270 <blas_scal+0x170>
movslq %edx,%rax
mov    %edi,%r8d
lea    -0x1(%r8),%rdx
mov    %r8d,%ecx
and    $0x3,%ecx
cmp    $0x3,%rdx
jae    11f7 <blas_scal+0xf7>
xor    %edi,%edi
jmp    124b <blas_scal+0x14b>
and    $0xfffffffc,%r8d
mov    %r10,%r11
shl    $0x4,%r11
mov    %r10,%r9
shl    $0x5,%r9
lea    0x0(,%r10,8),%rax
lea    (%rax,%rax,2),%rdx
xor    %edi,%edi
mov    %rsi,%rax
nopl   (%rax)
movsd  (%rax),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax)
movsd  (%rax,%r10,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%r10,8)
movsd  (%rax,%r11,1),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%r11,1)
movsd  (%rax,%rdx,1),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rdx,1)
add    $0x4,%rdi
add    %r9,%rax
cmp    %rdi,%r8
jne    1180 <blas_scal+0x80>
test   %rcx,%rcx
je     1270 <blas_scal+0x170>
imul   %r10,%rdi
lea    (%rsi,%rdi,8),%rax
shl    $0x3,%r10
nopl   0x0(%rax,%rax,1)
movsd  (%rax),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax)
add    %r10,%rax
add    $0xffffffffffffffff,%rcx
jne    11e0 <blas_scal+0xe0>
jmp    1270 <blas_scal+0x170>
and    $0xfffffffc,%r8d
mov    %rax,%r9
shl    $0x4,%r9
mov    %rax,%r11
shl    $0x5,%r11
lea    0x0(,%rax,8),%rdx
lea    (%rdx,%rdx,2),%r10
xor    %edi,%edi
mov    %rsi,%rdx
nopw   0x0(%rax,%rax,1)
movq   $0x0,(%rdx)
movq   $0x0,(%rdx,%rax,8)
movq   $0x0,(%rdx,%r9,1)
movq   $0x0,(%rdx,%r10,1)
add    $0x4,%rdi
add    %r11,%rdx
cmp    %rdi,%r8
jne    1220 <blas_scal+0x120>
test   %rcx,%rcx
je     1270 <blas_scal+0x170>
imul   %rax,%rdi
lea    (%rsi,%rdi,8),%rdx
shl    $0x3,%rax
nopl   0x0(%rax)
movq   $0x0,(%rdx)
add    %rax,%rdx
add    $0xffffffffffffffff,%rcx
jne    1260 <blas_scal+0x160>
ret

LLVM IR (-O0):
define void @blas_scal(i32 %0, double %1, double* %2, i32 %3) {
%5 = alloca i32
%6 = alloca double
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i32 %0, i32* %5
store double %1, double* %6
store double* %2, double** %7
store i32 %3, i32* %8
%11 = load double, double* %6
%12 = fcmp une double %11, 0.000000e+00
br i1 %12, label %13, label %32
13: ; preds = %4
store i32 0, i32* %9
br label %14
14: ; preds = %28, %13
%15 = load i32, i32* %9
%16 = load i32, i32* %5
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %31
18: ; preds = %14
%19 = load double, double* %6
%20 = load double*, double** %7
%21 = load i32, i32* %9
%22 = load i32, i32* %8
%23 = mul nsw i32 %21, %22
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %20, i64 %24
%26 = load double, double* %25
%27 = fmul double %26, %19
store double %27, double* %25
br label %28
28: ; preds = %18
%29 = load i32, i32* %9
%30 = add nsw i32 %29, 1
store i32 %30, i32* %9
br label %14
31: ; preds = %14
br label %48
32: ; preds = %4
store i32 0, i32* %10
br label %33
33: ; preds = %44, %32
%34 = load i32, i32* %10
%35 = load i32, i32* %5
%36 = icmp slt i32 %34, %35
br i1 %36, label %37, label %47
37: ; preds = %33
%38 = load double*, double** %7
%39 = load i32, i32* %10
%40 = load i32, i32* %8
%41 = mul nsw i32 %39, %40
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %38, i64 %42
store double 0.000000e+00, double* %43
br label %44
44: ; preds = %37
%45 = load i32, i32* %10
%46 = add nsw i32 %45, 1
store i32 %46, i32* %10
br label %33
47: ; preds = %33
br label %48
48: ; preds = %47, %31
ret void
}

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