void dscal(int n, double sa, double sx[], int incx) {
  int i, j;
  if (incx == 1) {
#pragma omp for
    for (i = 0; i < n; i++)
      sx[i] *= sa;
  } else {
    j = 0;
    for (i = 0; i < n; i++) {
      sx[j] *= sa;
      j += incx;
    }
  }
}


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)
cmpl   $0x1,-0x1c(%rbp)
jne    115a <dscal+0x5a>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1155 <dscal+0x55>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1124 <dscal+0x24>
jmp    11a7 <dscal+0xa7>
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11a2 <dscal+0xa2>
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)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1168 <dscal+0x68>
jmp    11a7 <dscal+0xa7>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dscal>:
push   %r14
push   %rbx
cmp    $0x1,%edx
jne    1122 <dscal+0x22>
test   %edi,%edi
jle    13a1 <dscal+0x2a1>
mov    %edi,%eax
cmp    $0x4,%edi
jae    1222 <dscal+0x122>
xor    %ecx,%ecx
jmp    12d0 <dscal+0x1d0>
test   %edi,%edi
jle    13a1 <dscal+0x2a1>
movslq %edx,%r10
lea    -0x1(%rdi),%r8d
xor    %ecx,%ecx
cmp    $0x3,%r8d
jb     12ec <dscal+0x1ec>
mov    $0x0,%r9d
cmp    $0x1,%edx
jne    12ef <dscal+0x1ef>
add    $0x1,%r8
mov    %r8,%r9
and    $0xfffffffffffffffc,%r9
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r9),%rax
mov    %rax,%r11
shr    $0x2,%r11
add    $0x1,%r11
test   %rax,%rax
je     13b5 <dscal+0x2b5>
mov    %r11,%rax
and    $0xfffffffffffffffe,%rax
mov    %r10,%rcx
shl    $0x5,%rcx
lea    0x10(%rsi),%rbx
mov    %r10,%r14
shl    $0x6,%r14
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd -0x10(%rbx),%xmm2
movupd (%rbx),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,-0x10(%rbx)
movupd %xmm3,(%rbx)
movupd -0x10(%rbx,%rcx,1),%xmm2
movupd (%rbx,%rcx,1),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,-0x10(%rbx,%rcx,1)
movupd %xmm3,(%rbx,%rcx,1)
add    $0x8,%rdx
add    %r14,%rbx
add    $0xfffffffffffffffe,%rax
jne    11a0 <dscal+0xa0>
test   $0x1,%r11b
je     120d <dscal+0x10d>
imul   %r10,%rdx
movupd (%rsi,%rdx,8),%xmm2
movupd 0x10(%rsi,%rdx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rsi,%rdx,8)
movupd %xmm3,0x10(%rsi,%rdx,8)
cmp    %r9,%r8
je     13a1 <dscal+0x2a1>
mov    %r9,%rcx
imul   %r10,%rcx
jmp    12ef <dscal+0x1ef>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%rcx),%rdx
mov    %rdx,%rbx
shr    $0x2,%rbx
add    $0x1,%rbx
test   %rdx,%rdx
je     13a5 <dscal+0x2a5>
mov    %rbx,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %edx,%edx
movupd (%rsi,%rdx,8),%xmm2
movupd 0x10(%rsi,%rdx,8),%xmm3
movupd 0x20(%rsi,%rdx,8),%xmm4
movupd 0x30(%rsi,%rdx,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rsi,%rdx,8)
movupd %xmm3,0x10(%rsi,%rdx,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,0x20(%rsi,%rdx,8)
movupd %xmm5,0x30(%rsi,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1250 <dscal+0x150>
test   $0x1,%bl
je     12bb <dscal+0x1bb>
movupd (%rsi,%rdx,8),%xmm2
movupd 0x10(%rsi,%rdx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rsi,%rdx,8)
movupd %xmm3,0x10(%rsi,%rdx,8)
cmp    %rax,%rcx
je     13a1 <dscal+0x2a1>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    12d0 <dscal+0x1d0>
jmp    13a1 <dscal+0x2a1>
xor    %r9d,%r9d
mov    %edi,%edx
sub    %r9d,%edx
mov    %r9d,%eax
not    %eax
add    %edi,%eax
and    $0x3,%edx
je     132b <dscal+0x22b>
neg    %edx
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
add    %r10,%rcx
add    $0xffffffff,%ebx
cmp    %ebx,%edx
jne    1310 <dscal+0x210>
sub    %ebx,%r9d
cmp    $0x3,%eax
jb     13a1 <dscal+0x2a1>
lea    (%rsi,%rcx,8),%rdx
mov    %r10,%r8
shl    $0x5,%r8
lea    0x0(,%r10,8),%rax
lea    (%rax,%rax,2),%rax
add    %r10,%rcx
shl    $0x4,%r10
lea    (%rsi,%rcx,8),%rcx
sub    %r9d,%edi
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx)
movsd  (%rcx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx)
movsd  (%rdx,%r10,1),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%r10,1)
movsd  (%rdx,%rax,1),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rax,1)
add    %r8,%rdx
add    %r8,%rcx
add    $0xfffffffc,%edi
jne    1360 <dscal+0x260>
pop    %rbx
pop    %r14
ret
xor    %edx,%edx
test   $0x1,%bl
jne    129d <dscal+0x19d>
jmp    12bb <dscal+0x1bb>
xor    %edx,%edx
test   $0x1,%r11b
jne    11eb <dscal+0xeb>
jmp    120d <dscal+0x10d>

LLVM IR (-O0):
define void @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* %8
%12 = icmp eq i32 %11, 1
br i1 %12, label %13, label %30
13: ; preds = %4
store i32 0, i32* %9
br label %14
14: ; preds = %26, %13
%15 = load i32, i32* %9
%16 = load i32, i32* %5
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %29
18: ; preds = %14
%19 = load double, double* %6
%20 = load double*, double** %7
%21 = load i32, i32* %9
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %20, i64 %22
%24 = load double, double* %23
%25 = fmul double %24, %19
store double %25, double* %23
br label %26
26: ; preds = %18
%27 = load i32, i32* %9
%28 = add nsw i32 %27, 1
store i32 %28, i32* %9
br label %14
29: ; preds = %14
br label %50
30: ; preds = %4
store i32 0, i32* %10
store i32 0, i32* %9
br label %31
31: ; preds = %46, %30
%32 = load i32, i32* %9
%33 = load i32, i32* %5
%34 = icmp slt i32 %32, %33
br i1 %34, label %35, label %49
35: ; preds = %31
%36 = load double, double* %6
%37 = load double*, double** %7
%38 = load i32, i32* %10
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %37, i64 %39
%41 = load double, double* %40
%42 = fmul double %41, %36
store double %42, double* %40
%43 = load i32, i32* %8
%44 = load i32, i32* %10
%45 = add nsw i32 %44, %43
store i32 %45, i32* %10
br label %46
46: ; preds = %35
%47 = load i32, i32* %9
%48 = add nsw i32 %47, 1
store i32 %48, i32* %9
br label %31
49: ; preds = %31
br label %50
50: ; preds = %49, %29
ret void
}

LLVM IR (-O2):
define void @dscal(i32 %0, double %1, double* nocapture %2, i32 %3) {
%5 = icmp eq i32 %3, 1
%6 = icmp sgt i32 %0, 0
br i1 %5, label %100, label %7
7: ; preds = %4
br i1 %6, label %8, label %190
8: ; preds = %7
%9 = sext i32 %3 to i64
%10 = add i32 %0, -1
%11 = zext i32 %10 to i64
%12 = add nuw nsw i64 %11, 1
%13 = icmp ugt i32 %10, 2
%14 = icmp eq i32 %3, 1
%15 = and i1 %13, %14
br i1 %15, label %16, label %77
16: ; preds = %8
%17 = and i64 %12, 8589934588
%18 = mul i64 %17, %9
%19 = trunc i64 %17 to i32
%20 = insertelement <2 x double> poison, double %1, i64 0
%21 = shufflevector <2 x double> %20, <2 x double> poison, <2 x i32> zeroinitializer
%22 = insertelement <2 x double> poison, double %1, i64 0
%23 = shufflevector <2 x double> %22, <2 x double> poison, <2 x i32> zeroinitializer
%24 = add nsw i64 %17, -4
%25 = lshr exact i64 %24, 2
%26 = add nuw nsw i64 %25, 1
%27 = and i64 %26, 1
%28 = icmp eq i64 %24, 0
br i1 %28, label %60, label %29
29: ; preds = %16
%30 = and i64 %26, 9223372036854775806
br label %31
31: ; preds = %31, %29
%32 = phi i64 [ 0, %29 ], [ %57, %31 ]
%33 = phi i64 [ 0, %29 ], [ %58, %31 ]
%34 = mul i64 %32, %9
%35 = getelementptr inbounds double, double* %2, i64 %34
%36 = bitcast double* %35 to <2 x double>*
%37 = load <2 x double>, <2 x double>* %36
%38 = getelementptr inbounds double, double* %35, i64 2
%39 = bitcast double* %38 to <2 x double>*
%40 = load <2 x double>, <2 x double>* %39
%41 = fmul <2 x double> %37, %21
%42 = fmul <2 x double> %40, %23
%43 = bitcast double* %35 to <2 x double>*
store <2 x double> %41, <2 x double>* %43
%44 = bitcast double* %38 to <2 x double>*
store <2 x double> %42, <2 x double>* %44
%45 = or i64 %32, 4
%46 = mul i64 %45, %9
%47 = getelementptr inbounds double, double* %2, i64 %46
%48 = bitcast double* %47 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = getelementptr inbounds double, double* %47, i64 2
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = fmul <2 x double> %49, %21
%54 = fmul <2 x double> %52, %23
%55 = bitcast double* %47 to <2 x double>*
store <2 x double> %53, <2 x double>* %55
%56 = bitcast double* %50 to <2 x double>*
store <2 x double> %54, <2 x double>* %56
%57 = add nuw i64 %32, 8
%58 = add i64 %33, 2
%59 = icmp eq i64 %58, %30
br i1 %59, label %60, label %31
60: ; preds = %31, %16
%61 = phi i64 [ 0, %16 ], [ %57, %31 ]
%62 = icmp eq i64 %27, 0
br i1 %62, label %75, label %63
63: ; preds = %60
%64 = mul i64 %61, %9
%65 = getelementptr inbounds double, double* %2, i64 %64
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = getelementptr inbounds double, double* %65, i64 2
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = fmul <2 x double> %67, %21
%72 = fmul <2 x double> %70, %23
%73 = bitcast double* %65 to <2 x double>*
store <2 x double> %71, <2 x double>* %73
%74 = bitcast double* %68 to <2 x double>*
store <2 x double> %72, <2 x double>* %74
br label %75
75: ; preds = %60, %63
%76 = icmp eq i64 %12, %17
br i1 %76, label %190, label %77
77: ; preds = %8, %75
%78 = phi i64 [ 0, %8 ], [ %18, %75 ]
%79 = phi i32 [ 0, %8 ], [ %19, %75 ]
%80 = sub i32 %0, %79
%81 = xor i32 %79, -1
%82 = add i32 %81, %0
%83 = and i32 %80, 3
%84 = icmp eq i32 %83, 0
br i1 %84, label %96, label %85
85: ; preds = %77, %85
%86 = phi i64 [ %92, %85 ], [ %78, %77 ]
%87 = phi i32 [ %93, %85 ], [ %79, %77 ]
%88 = phi i32 [ %94, %85 ], [ 0, %77 ]
%89 = getelementptr inbounds double, double* %2, i64 %86
%90 = load double, double* %89
%91 = fmul double %90, %1
store double %91, double* %89
%92 = add nsw i64 %86, %9
%93 = add nuw nsw i32 %87, 1
%94 = add i32 %88, 1
%95 = icmp eq i32 %94, %83
br i1 %95, label %96, label %85
96: ; preds = %85, %77
%97 = phi i64 [ %78, %77 ], [ %92, %85 ]
%98 = phi i32 [ %79, %77 ], [ %93, %85 ]
%99 = icmp ult i32 %82, 3
br i1 %99, label %190, label %169
100: ; preds = %4
br i1 %6, label %101, label %190
101: ; preds = %100
%102 = zext i32 %0 to i64
%103 = icmp ult i32 %0, 4
br i1 %103, label %160, label %104
104: ; preds = %101
%105 = and i64 %102, 4294967292
%106 = insertelement <2 x double> poison, double %1, i64 0
%107 = shufflevector <2 x double> %106, <2 x double> poison, <2 x i32> zeroinitializer
%108 = insertelement <2 x double> poison, double %1, i64 0
%109 = shufflevector <2 x double> %108, <2 x double> poison, <2 x i32> zeroinitializer
%110 = add nsw i64 %105, -4
%111 = lshr exact i64 %110, 2
%112 = add nuw nsw i64 %111, 1
%113 = and i64 %112, 1
%114 = icmp eq i64 %110, 0
br i1 %114, label %144, label %115
115: ; preds = %104
%116 = and i64 %112, 9223372036854775806
br label %117
117: ; preds = %117, %115
%118 = phi i64 [ 0, %115 ], [ %141, %117 ]
%119 = phi i64 [ 0, %115 ], [ %142, %117 ]
%120 = getelementptr inbounds double, double* %2, i64 %118
%121 = bitcast double* %120 to <2 x double>*
%122 = load <2 x double>, <2 x double>* %121
%123 = getelementptr inbounds double, double* %120, i64 2
%124 = bitcast double* %123 to <2 x double>*
%125 = load <2 x double>, <2 x double>* %124
%126 = fmul <2 x double> %122, %107
%127 = fmul <2 x double> %125, %109
%128 = bitcast double* %120 to <2 x double>*
store <2 x double> %126, <2 x double>* %128
%129 = bitcast double* %123 to <2 x double>*
store <2 x double> %127, <2 x double>* %129
%130 = or i64 %118, 4
%131 = getelementptr inbounds double, double* %2, i64 %130
%132 = bitcast double* %131 to <2 x double>*
%133 = load <2 x double>, <2 x double>* %132
%134 = getelementptr inbounds double, double* %131, i64 2
%135 = bitcast double* %134 to <2 x double>*
%136 = load <2 x double>, <2 x double>* %135
%137 = fmul <2 x double> %133, %107
%138 = fmul <2 x double> %136, %109
%139 = bitcast double* %131 to <2 x double>*
store <2 x double> %137, <2 x double>* %139
%140 = bitcast double* %134 to <2 x double>*
store <2 x double> %138, <2 x double>* %140
%141 = add nuw i64 %118, 8
%142 = add i64 %119, 2
%143 = icmp eq i64 %142, %116
br i1 %143, label %144, label %117
144: ; preds = %117, %104
%145 = phi i64 [ 0, %104 ], [ %141, %117 ]
%146 = icmp eq i64 %113, 0
br i1 %146, label %158, label %147
147: ; preds = %144
%148 = getelementptr inbounds double, double* %2, i64 %145
%149 = bitcast double* %148 to <2 x double>*
%150 = load <2 x double>, <2 x double>* %149
%151 = getelementptr inbounds double, double* %148, i64 2
%152 = bitcast double* %151 to <2 x double>*
%153 = load <2 x double>, <2 x double>* %152
%154 = fmul <2 x double> %150, %107
%155 = fmul <2 x double> %153, %109
%156 = bitcast double* %148 to <2 x double>*
store <2 x double> %154, <2 x double>* %156
%157 = bitcast double* %151 to <2 x double>*
store <2 x double> %155, <2 x double>* %157
br label %158
158: ; preds = %144, %147
%159 = icmp eq i64 %105, %102
br i1 %159, label %190, label %160
160: ; preds = %101, %158
%161 = phi i64 [ 0, %101 ], [ %105, %158 ]
br label %162
162: ; preds = %160, %162
%163 = phi i64 [ %167, %162 ], [ %161, %160 ]
%164 = getelementptr inbounds double, double* %2, i64 %163
%165 = load double, double* %164
%166 = fmul double %165, %1
store double %166, double* %164
%167 = add nuw nsw i64 %163, 1
%168 = icmp eq i64 %167, %102
br i1 %168, label %190, label %162
169: ; preds = %96, %169
%170 = phi i64 [ %187, %169 ], [ %97, %96 ]
%171 = phi i32 [ %188, %169 ], [ %98, %96 ]
%172 = getelementptr inbounds double, double* %2, i64 %170
%173 = load double, double* %172
%174 = fmul double %173, %1
store double %174, double* %172
%175 = add nsw i64 %170, %9
%176 = getelementptr inbounds double, double* %2, i64 %175
%177 = load double, double* %176
%178 = fmul double %177, %1
store double %178, double* %176
%179 = add nsw i64 %175, %9
%180 = getelementptr inbounds double, double* %2, i64 %179
%181 = load double, double* %180
%182 = fmul double %181, %1
store double %182, double* %180
%183 = add nsw i64 %179, %9
%184 = getelementptr inbounds double, double* %2, i64 %183
%185 = load double, double* %184
%186 = fmul double %185, %1
store double %186, double* %184
%187 = add nsw i64 %183, %9
%188 = add nuw nsw i32 %171, 4
%189 = icmp eq i32 %188, %0
br i1 %189, label %190, label %169
190: ; preds = %96, %169, %162, %75, %158, %7, %100
ret void
}