void svd_dscal(long n, double da, double *dx, long incx) {
  long i;
  if (n <= 0 || incx == 0)
    return;
  if (incx < 0)
    dx += (-n + 1) * incx;
  for (i = 0; i < n; i++) {
    *dx *= da;
    dx += incx;
  }
  return;
}


x86 assembly (clang, -O0):
<svd_dscal>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
cmpq   $0x0,-0x8(%rbp)
jle    112b <svd_dscal+0x2b>
cmpq   $0x0,-0x20(%rbp)
jne    1130 <svd_dscal+0x30>
jmp    11a3 <svd_dscal+0xa3>
cmpq   $0x0,-0x20(%rbp)
jge    1156 <svd_dscal+0x56>
xor    %eax,%eax
sub    -0x8(%rbp),%rax
add    $0x1,%rax
imul   -0x20(%rbp),%rax
shl    $0x3,%rax
add    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x8(%rbp),%rax
jge    119e <svd_dscal+0x9e>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
shl    $0x3,%rax
add    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    115e <svd_dscal+0x5e>
jmp    11a3 <svd_dscal+0xa3>
pop    %rbp
ret

x86 assembly (clang, -O2):
<svd_dscal>:
test   %rdi,%rdi
jle    11c5 <svd_dscal+0xc5>
test   %rdx,%rdx
je     11c5 <svd_dscal+0xc5>
mov    $0x1,%eax
sub    %rdi,%rax
imul   %rdx,%rax
mov    %rdx,%rcx
sar    $0x3f,%rcx
and    %rax,%rcx
lea    (%rsi,%rcx,8),%rcx
lea    -0x1(%rdi),%rax
mov    %edi,%r9d
and    $0x3,%r9d
cmp    $0x3,%rax
jb     11a5 <svd_dscal+0xa5>
and    $0xfffffffffffffffc,%rdi
mov    %rdx,%r8
shl    $0x5,%r8
lea    0x0(,%rdx,8),%rsi
mov    %rcx,%rax
movsd  (%rcx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx)
movsd  (%rcx,%rsi,1),%xmm1
add    %rsi,%rcx
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rsi,1)
movsd  (%rsi,%rcx,1),%xmm1
lea    (%rcx,%rsi,1),%r10
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%rcx,1)
movsd  (%rsi,%r10,1),%xmm1
lea    (%r10,%rsi,1),%rcx
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%r10,1)
add    %rsi,%rcx
add    $0xfffffffffffffffc,%rdi
jne    1150 <svd_dscal+0x50>
add    %r8,%rax
mov    %rax,%rcx
test   %r9,%r9
je     11c5 <svd_dscal+0xc5>
shl    $0x3,%rdx
xchg   %ax,%ax
movsd  (%rcx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx)
add    %rdx,%rcx
add    $0xffffffffffffffff,%r9
jne    11b0 <svd_dscal+0xb0>
ret

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

LLVM IR (-O2):
define void @svd_dscal(i64 %0, double %1, double* nocapture %2, i64 %3) {
%5 = icmp slt i64 %0, 1
%6 = icmp eq i64 %3, 0
%7 = or i1 %5, %6
br i1 %7, label %47, label %8
8: ; preds = %4
%9 = icmp slt i64 %3, 0
%10 = sub nsw i64 1, %0
%11 = mul nsw i64 %10, %3
%12 = select i1 %9, i64 %11, i64 0
%13 = getelementptr double, double* %2, i64 %12
%14 = add i64 %0, -1
%15 = and i64 %0, 3
%16 = icmp ult i64 %14, 3
br i1 %16, label %36, label %17
17: ; preds = %8
%18 = and i64 %0, -4
br label %19
19: ; preds = %19, %17
%20 = phi double* [ %13, %17 ], [ %33, %19 ]
%21 = phi i64 [ 0, %17 ], [ %34, %19 ]
%22 = load double, double* %20
%23 = fmul double %22, %1
store double %23, double* %20
%24 = getelementptr inbounds double, double* %20, i64 %3
%25 = load double, double* %24
%26 = fmul double %25, %1
store double %26, double* %24
%27 = getelementptr inbounds double, double* %24, i64 %3
%28 = load double, double* %27
%29 = fmul double %28, %1
store double %29, double* %27
%30 = getelementptr inbounds double, double* %27, i64 %3
%31 = load double, double* %30
%32 = fmul double %31, %1
store double %32, double* %30
%33 = getelementptr inbounds double, double* %30, i64 %3
%34 = add i64 %21, 4
%35 = icmp eq i64 %34, %18
br i1 %35, label %36, label %19
36: ; preds = %19, %8
%37 = phi double* [ %13, %8 ], [ %33, %19 ]
%38 = icmp eq i64 %15, 0
br i1 %38, label %47, label %39
39: ; preds = %36, %39
%40 = phi double* [ %44, %39 ], [ %37, %36 ]
%41 = phi i64 [ %45, %39 ], [ 0, %36 ]
%42 = load double, double* %40
%43 = fmul double %42, %1
store double %43, double* %40
%44 = getelementptr inbounds double, double* %40, i64 %3
%45 = add i64 %41, 1
%46 = icmp eq i64 %45, %15
br i1 %46, label %47, label %39
47: ; preds = %36, %39, %4
ret void
}