void dscal(int n, double da, double *dx, int dx_off, int incx) {
  int i;
  if (n > 0) {
    if (incx != 1) {
      int nincx = n * incx;
      for (i = 0; i < nincx; i += incx)
        dx[i + dx_off] *= da;
    } else {
      for (i = 0; i < n; i++)
        dx[i + dx_off] *= da;
    }
  }
}


x86 assembly (clang, -O0):
<dscal>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
cmpl   $0x0,-0x4(%rbp)
jle    11bd <dscal+0xbd>
cmpl   $0x1,-0x20(%rbp)
je     1176 <dscal+0x76>
mov    -0x4(%rbp),%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    1171 <dscal+0x71>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    113b <dscal+0x3b>
jmp    11b8 <dscal+0xb8>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11b3 <dscal+0xb3>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    117d <dscal+0x7d>
jmp    11b8 <dscal+0xb8>
jmp    11bd <dscal+0xbd>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dscal>:
test   %edi,%edi
jle    1237 <dscal+0x137>
cmp    $0x1,%ecx
jne    111f <dscal+0x1f>
movslq %edx,%r8
mov    %edi,%r10d
cmp    $0x4,%edi
jae    115b <dscal+0x5b>
xor    %ecx,%ecx
jmp    1214 <dscal+0x114>
imul   %ecx,%edi
test   %edi,%edi
jle    1237 <dscal+0x137>
movslq %ecx,%rax
movslq %edx,%rdx
mov    %edi,%ecx
lea    (%rsi,%rdx,8),%rdx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%rsi,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rsi,8)
add    %rax,%rsi
cmp    %rcx,%rsi
jl     1140 <dscal+0x40>
jmp    1237 <dscal+0x137>
mov    %r10d,%ecx
and    $0xfffffffc,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rcx),%rax
mov    %rax,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rax,%rax
je     1238 <dscal+0x138>
mov    %r9,%rdx
and    $0xfffffffffffffffe,%rdx
lea    (%rsi,%r8,8),%rax
add    $0x30,%rax
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd -0x30(%rax,%rdi,8),%xmm2
movupd -0x20(%rax,%rdi,8),%xmm3
movupd -0x10(%rax,%rdi,8),%xmm4
movupd (%rax,%rdi,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,-0x30(%rax,%rdi,8)
movupd %xmm3,-0x20(%rax,%rdi,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,-0x10(%rax,%rdi,8)
movupd %xmm5,(%rax,%rdi,8)
add    $0x8,%rdi
add    $0xfffffffffffffffe,%rdx
jne    11a0 <dscal+0xa0>
test   $0x1,%r9b
je     120f <dscal+0x10f>
add    %r8,%rdi
movupd (%rsi,%rdi,8),%xmm2
movupd 0x10(%rsi,%rdi,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rsi,%rdi,8)
movupd %xmm3,0x10(%rsi,%rdi,8)
cmp    %r10,%rcx
je     1237 <dscal+0x137>
lea    (%rsi,%r8,8),%rax
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r10
jne    1220 <dscal+0x120>
ret
xor    %edi,%edi
test   $0x1,%r9b
jne    11ee <dscal+0xee>
jmp    120f <dscal+0x10f>

LLVM IR (-O0):
define void @dscal(i32 %0, double %1, double* %2, i32 %3, i32 %4) {
%6 = alloca i32
%7 = alloca double
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
store i32 %0, i32* %6
store double %1, double* %7
store double* %2, double** %8
store i32 %3, i32* %9
store i32 %4, i32* %10
%13 = load i32, i32* %6
%14 = icmp sgt i32 %13, 0
br i1 %14, label %15, label %61
15: ; preds = %5
%16 = load i32, i32* %10
%17 = icmp ne i32 %16, 1
br i1 %17, label %18, label %41
18: ; preds = %15
%19 = load i32, i32* %6
%20 = load i32, i32* %10
%21 = mul nsw i32 %19, %20
store i32 %21, i32* %12
store i32 0, i32* %11
br label %22
22: ; preds = %36, %18
%23 = load i32, i32* %11
%24 = load i32, i32* %12
%25 = icmp slt i32 %23, %24
br i1 %25, label %26, label %40
26: ; preds = %22
%27 = load double, double* %7
%28 = load double*, double** %8
%29 = load i32, i32* %11
%30 = load i32, i32* %9
%31 = add nsw i32 %29, %30
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %28, i64 %32
%34 = load double, double* %33
%35 = fmul double %34, %27
store double %35, double* %33
br label %36
36: ; preds = %26
%37 = load i32, i32* %10
%38 = load i32, i32* %11
%39 = add nsw i32 %38, %37
store i32 %39, i32* %11
br label %22
40: ; preds = %22
br label %60
41: ; preds = %15
store i32 0, i32* %11
br label %42
42: ; preds = %56, %41
%43 = load i32, i32* %11
%44 = load i32, i32* %6
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %59
46: ; preds = %42
%47 = load double, double* %7
%48 = load double*, double** %8
%49 = load i32, i32* %11
%50 = load i32, i32* %9
%51 = add nsw i32 %49, %50
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %48, i64 %52
%54 = load double, double* %53
%55 = fmul double %54, %47
store double %55, double* %53
br label %56
56: ; preds = %46
%57 = load i32, i32* %11
%58 = add nsw i32 %57, 1
store i32 %58, i32* %11
br label %42
59: ; preds = %42
br label %60
60: ; preds = %59, %40
br label %61
61: ; preds = %60, %5
ret void
}

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