void c_sscal(int n, double sa, double *sx) {
  register int i, m;
  if (n <= 0) {
    return;
  }
  m = n % 4;
  if (m != 0) {
    for (i = 1; i <= m; i++) {
      sx[i - 1] *= sa;
    }
  }
  for (i = m + 1; i <= n; i += 4) {
    sx[i - 1] *= sa;
    sx[i + 1 - 1] *= sa;
    sx[i + 2 - 1] *= sa;
    sx[i + 3 - 1] *= sa;
  }
  return;
}


x86 assembly (clang, -O0):
<c_sscal>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
cmpl   $0x0,-0x4(%rbp)
jg     111f <c_sscal+0x1f>
jmp    121a <c_sscal+0x11a>
mov    -0x4(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %edx,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
je     1179 <c_sscal+0x79>
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     1174 <c_sscal+0x74>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    113e <c_sscal+0x3e>
jmp    1179 <c_sscal+0x79>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1215 <c_sscal+0x115>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x2,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x3,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x1c(%rbp)
jmp    1182 <c_sscal+0x82>
jmp    121a <c_sscal+0x11a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<c_sscal>:
test   %edi,%edi
jle    1259 <c_sscal+0x159>
mov    %edi,%ecx
and    $0x3,%ecx
je     1144 <c_sscal+0x44>
lea    0x1(%rcx),%eax
movsd  (%rsi),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi)
cmp    $0x2,%eax
je     1144 <c_sscal+0x44>
movsd  0x8(%rsi),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsi)
cmp    $0x3,%eax
je     1144 <c_sscal+0x44>
movsd  0x10(%rsi),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsi)
cmp    $0x4,%edi
jl     1259 <c_sscal+0x159>
lea    0x1(%rcx),%r10
mov    %edi,%r11d
add    $0x5,%rcx
lea    0x1(%r11),%rax
cmp    %rax,%rcx
cmovbe %rax,%rcx
mov    %r10,%r9
not    %r9
add    %rcx,%r9
cmp    $0x4,%r9
jae    117a <c_sscal+0x7a>
mov    %r10,%rcx
jmp    1221 <c_sscal+0x121>
shr    $0x2,%r9
add    $0x1,%r9
mov    %r9,%r8
and    $0xfffffffffffffffe,%r8
lea    (%r10,%r8,4),%rcx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
xor    %eax,%eax
mov    %r8,%rdx
nopw   0x0(%rax,%rax,1)
mov    %rax,%rdi
or     %r10,%rdi
movupd -0x8(%rsi,%rdi,8),%xmm2
movupd 0x8(%rsi,%rdi,8),%xmm3
movupd 0x18(%rsi,%rdi,8),%xmm4
movupd 0x28(%rsi,%rdi,8),%xmm5
movapd %xmm2,%xmm6
unpcklpd %xmm4,%xmm6
unpckhpd %xmm4,%xmm2
movapd %xmm3,%xmm4
unpcklpd %xmm5,%xmm4
unpckhpd %xmm5,%xmm3
mulpd  %xmm1,%xmm6
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm3
movapd %xmm4,%xmm5
unpcklpd %xmm3,%xmm5
movapd %xmm6,%xmm7
unpcklpd %xmm2,%xmm7
movhlps %xmm4,%xmm3
movhlps %xmm6,%xmm2
movups %xmm2,0x18(%rsi,%rdi,8)
movups %xmm3,0x28(%rsi,%rdi,8)
movupd %xmm7,-0x8(%rsi,%rdi,8)
movupd %xmm5,0x8(%rsi,%rdi,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rdx
jne    11a0 <c_sscal+0xa0>
cmp    %r8,%r9
je     1259 <c_sscal+0x159>
unpcklpd %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movupd -0x8(%rsi,%rcx,8),%xmm1
movupd 0x8(%rsi,%rcx,8),%xmm2
mulpd  %xmm0,%xmm1
movupd %xmm1,-0x8(%rsi,%rcx,8)
mulpd  %xmm0,%xmm2
movupd %xmm2,0x8(%rsi,%rcx,8)
add    $0x4,%rcx
cmp    %r11,%rcx
jbe    1230 <c_sscal+0x130>
ret

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

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