void c_dscal(const int N, const double alpha, double *X, const int stride) {
  int m;
  int i;
  if (N <= 0 || stride <= 0 || alpha == 1.0) {
    return;
  }
  if (stride == 1) {
    m = N % 5;
    if (m > 0) {
      for (i = 0; i < m; i++) {
        X[i] *= alpha;
      }
    }
    if (N < 5) {
      return;
    }
    for (i = m; i < N; i += 5) {
      X[i] *= alpha;
      X[i + 1] *= alpha;
      X[i + 2] *= alpha;
      X[i + 3] *= alpha;
      X[i + 4] *= alpha;
    }
    return;
  }
  for (i = 0; i < N * stride; i += stride) {
    X[i] *= alpha;
  }
  return;
}


x86 assembly (clang, -O0):
<c_dscal>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
cmpl   $0x0,-0x4(%rbp)
jle    1144 <c_dscal+0x44>
cmpl   $0x0,-0x1c(%rbp)
jle    1144 <c_dscal+0x44>
movsd  -0x10(%rbp),%xmm0
movsd  0xecc(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1149 <c_dscal+0x49>
jp     1149 <c_dscal+0x49>
jmp    12a6 <c_dscal+0x1a6>
cmpl   $0x1,-0x1c(%rbp)
jne    1263 <c_dscal+0x163>
mov    -0x4(%rbp),%eax
mov    $0x5,%ecx
cltd
idiv   %ecx
mov    %edx,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jle    11a8 <c_dscal+0xa8>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11a3 <c_dscal+0xa3>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1172 <c_dscal+0x72>
jmp    11a8 <c_dscal+0xa8>
cmpl   $0x5,-0x4(%rbp)
jge    11b7 <c_dscal+0xb7>
jmp    12a6 <c_dscal+0x1a6>
mov    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    125e <c_dscal+0x15e>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x4,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x5,%eax
mov    %eax,-0x24(%rbp)
jmp    11bd <c_dscal+0xbd>
jmp    12a6 <c_dscal+0x1a6>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    12a1 <c_dscal+0x1a1>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    126a <c_dscal+0x16a>
jmp    12a6 <c_dscal+0x1a6>
pop    %rbp
ret

x86 assembly (clang, -O2):
<c_dscal>:
ucomisd 0xef8(%rip),%xmm0
jne    1110 <c_dscal+0x10>
jnp    12c3 <c_dscal+0x1c3>
test   %edi,%edi
jle    12c3 <c_dscal+0x1c3>
test   %edx,%edx
jle    12c3 <c_dscal+0x1c3>
cmp    $0x1,%edx
jne    1154 <c_dscal+0x54>
mov    %edi,%r8d
mov    $0xcccccccd,%eax
imul   %r8,%rax
shr    $0x22,%rax
lea    (%rax,%rax,4),%eax
mov    %edi,%r9d
sub    %eax,%r9d
je     1257 <c_dscal+0x157>
mov    %r9d,%r11d
cmp    $0x4,%r9d
jae    118b <c_dscal+0x8b>
xor    %edx,%edx
jmp    1240 <c_dscal+0x140>
imul   %edx,%edi
test   %edi,%edi
jle    12c3 <c_dscal+0x1c3>
mov    %edx,%eax
mov    %edi,%ecx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rsi,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rdx,8)
add    %rax,%rdx
cmp    %rcx,%rdx
jl     1170 <c_dscal+0x70>
jmp    12c3 <c_dscal+0x1c3>
mov    %r11d,%edx
and    $0x4,%edx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rdx),%rax
mov    %rax,%r10
shr    $0x2,%r10
add    $0x1,%r10
test   %rax,%rax
je     12c4 <c_dscal+0x1c4>
mov    %r10,%rax
and    $0xfffffffffffffffe,%rax
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movupd (%rsi,%rcx,8),%xmm2
movupd 0x10(%rsi,%rcx,8),%xmm3
movupd 0x20(%rsi,%rcx,8),%xmm4
movupd 0x30(%rsi,%rcx,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rsi,%rcx,8)
movupd %xmm3,0x10(%rsi,%rcx,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,0x20(%rsi,%rcx,8)
movupd %xmm5,0x30(%rsi,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rax
jne    11c0 <c_dscal+0xc0>
test   $0x1,%r10b
je     122c <c_dscal+0x12c>
movupd (%rsi,%rcx,8),%xmm2
movupd 0x10(%rsi,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rsi,%rcx,8)
movupd %xmm3,0x10(%rsi,%rcx,8)
cmp    %r11,%rdx
je     1257 <c_dscal+0x157>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%r11
jne    1240 <c_dscal+0x140>
cmp    $0x5,%edi
jl     12c3 <c_dscal+0x1c3>
cmp    %edi,%r9d
jge    12c3 <c_dscal+0x1c3>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    $0xcccccccd,%ecx
imul   %r8,%rcx
shr    $0x22,%rcx
lea    (%rcx,%rcx,4),%rax
sub    %rax,%r8
lea    (%rsi,%r8,8),%rax
add    $0x20,%rax
lea    (%rcx,%rcx,4),%edx
mov    %edi,%ecx
sub    %edx,%ecx
nopl   0x0(%rax)
movupd -0x20(%rax),%xmm2
movupd -0x10(%rax),%xmm3
mulpd  %xmm1,%xmm2
movupd %xmm2,-0x20(%rax)
mulpd  %xmm1,%xmm3
movupd %xmm3,-0x10(%rax)
movsd  (%rax),%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,(%rax)
add    $0x28,%rax
add    $0x5,%ecx
cmp    %edi,%ecx
jl     1290 <c_dscal+0x190>
ret
xor    %ecx,%ecx
test   $0x1,%r10b
jne    120e <c_dscal+0x10e>
jmp    122c <c_dscal+0x12c>

LLVM IR (-O0):
define void @c_dscal(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 i32, i32* %5
%12 = icmp sle i32 %11, 0
br i1 %12, label %19, label %13
13: ; preds = %4
%14 = load i32, i32* %8
%15 = icmp sle i32 %14, 0
br i1 %15, label %19, label %16
16: ; preds = %13
%17 = load double, double* %6
%18 = fcmp oeq double %17, 1.000000e+00
br i1 %18, label %19, label %20
19: ; preds = %16, %13, %4
br label %119
20: ; preds = %16
%21 = load i32, i32* %8
%22 = icmp eq i32 %21, 1
br i1 %22, label %23, label %99
23: ; preds = %20
%24 = load i32, i32* %5
%25 = srem i32 %24, 5
store i32 %25, i32* %9
%26 = load i32, i32* %9
%27 = icmp sgt i32 %26, 0
br i1 %27, label %28, label %45
28: ; preds = %23
store i32 0, i32* %10
br label %29
29: ; preds = %41, %28
%30 = load i32, i32* %10
%31 = load i32, i32* %9
%32 = icmp slt i32 %30, %31
br i1 %32, label %33, label %44
33: ; preds = %29
%34 = load double, double* %6
%35 = load double*, double** %7
%36 = load i32, i32* %10
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %35, i64 %37
%39 = load double, double* %38
%40 = fmul double %39, %34
store double %40, double* %38
br label %41
41: ; preds = %33
%42 = load i32, i32* %10
%43 = add nsw i32 %42, 1
store i32 %43, i32* %10
br label %29
44: ; preds = %29
br label %45
45: ; preds = %44, %23
%46 = load i32, i32* %5
%47 = icmp slt i32 %46, 5
br i1 %47, label %48, label %49
48: ; preds = %45
br label %119
49: ; preds = %45
%50 = load i32, i32* %9
store i32 %50, i32* %10
br label %51
51: ; preds = %95, %49
%52 = load i32, i32* %10
%53 = load i32, i32* %5
%54 = icmp slt i32 %52, %53
br i1 %54, label %55, label %98
55: ; preds = %51
%56 = load double, double* %6
%57 = load double*, double** %7
%58 = load i32, i32* %10
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %57, i64 %59
%61 = load double, double* %60
%62 = fmul double %61, %56
store double %62, double* %60
%63 = load double, double* %6
%64 = load double*, double** %7
%65 = load i32, i32* %10
%66 = add nsw i32 %65, 1
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %64, i64 %67
%69 = load double, double* %68
%70 = fmul double %69, %63
store double %70, double* %68
%71 = load double, double* %6
%72 = load double*, double** %7
%73 = load i32, i32* %10
%74 = add nsw i32 %73, 2
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %72, i64 %75
%77 = load double, double* %76
%78 = fmul double %77, %71
store double %78, double* %76
%79 = load double, double* %6
%80 = load double*, double** %7
%81 = load i32, i32* %10
%82 = add nsw i32 %81, 3
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds double, double* %80, i64 %83
%85 = load double, double* %84
%86 = fmul double %85, %79
store double %86, double* %84
%87 = load double, double* %6
%88 = load double*, double** %7
%89 = load i32, i32* %10
%90 = add nsw i32 %89, 4
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %88, i64 %91
%93 = load double, double* %92
%94 = fmul double %93, %87
store double %94, double* %92
br label %95
95: ; preds = %55
%96 = load i32, i32* %10
%97 = add nsw i32 %96, 5
store i32 %97, i32* %10
br label %51
98: ; preds = %51
br label %119
99: ; preds = %20
store i32 0, i32* %10
br label %100
100: ; preds = %114, %99
%101 = load i32, i32* %10
%102 = load i32, i32* %5
%103 = load i32, i32* %8
%104 = mul nsw i32 %102, %103
%105 = icmp slt i32 %101, %104
br i1 %105, label %106, label %118
106: ; preds = %100
%107 = load double, double* %6
%108 = load double*, double** %7
%109 = load i32, i32* %10
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %108, i64 %110
%112 = load double, double* %111
%113 = fmul double %112, %107
store double %113, double* %111
br label %114
114: ; preds = %106
%115 = load i32, i32* %8
%116 = load i32, i32* %10
%117 = add nsw i32 %116, %115
store i32 %117, i32* %10
br label %100
118: ; preds = %100
br label %119
119: ; preds = %118, %98, %48, %19
ret void
}

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