void m_scale(double a, double *A, int nrow, int ncol) {
  int i;
  for (i = nrow * ncol; i > 0; i--)
    *A++ *= a;
}


x86 assembly (clang, -O0):
<m_scale>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jle    1151 <m_scale+0x51>
movsd  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x10(%rbp)
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    111d <m_scale+0x1d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<m_scale>:
imul   %edx,%esi
test   %esi,%esi
jle    11e8 <m_scale+0xe8>
lea    -0x1(%rsi),%r9d
cmp    $0x3,%r9d
jb     11c6 <m_scale+0xc6>
add    $0x1,%r9
mov    %r9,%r8
and    $0xfffffffffffffffc,%r8
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r8),%rcx
mov    %rcx,%rax
shr    $0x2,%rax
add    $0x1,%rax
test   %rcx,%rcx
je     11e9 <m_scale+0xe9>
mov    %rax,%rcx
and    $0xfffffffffffffffe,%rcx
xor    %edx,%edx
nopl   (%rax)
movupd (%rdi,%rdx,8),%xmm2
movupd 0x10(%rdi,%rdx,8),%xmm3
movupd 0x20(%rdi,%rdx,8),%xmm4
movupd 0x30(%rdi,%rdx,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rdi,%rdx,8)
movupd %xmm3,0x10(%rdi,%rdx,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,0x20(%rdi,%rdx,8)
movupd %xmm5,0x30(%rdi,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rcx
jne    1150 <m_scale+0x50>
test   $0x1,%al
je     11ba <m_scale+0xba>
movupd (%rdi,%rdx,8),%xmm2
movupd 0x10(%rdi,%rdx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rdi,%rdx,8)
movupd %xmm3,0x10(%rdi,%rdx,8)
cmp    %r8,%r9
je     11e8 <m_scale+0xe8>
sub    %r8d,%esi
lea    (%rdi,%r8,8),%rdi
add    $0x1,%esi
nopl   0x0(%rax)
movsd  (%rdi),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi)
add    $0x8,%rdi
add    $0xffffffff,%esi
cmp    $0x1,%esi
ja     11d0 <m_scale+0xd0>
ret
xor    %edx,%edx
test   $0x1,%al
jne    119c <m_scale+0x9c>
jmp    11ba <m_scale+0xba>

LLVM IR (-O0):
define void @m_scale(double %0, double* %1, i32 %2, i32 %3) {
%5 = alloca double
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
store double %0, double* %5
store double* %1, double** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%10 = load i32, i32* %7
%11 = load i32, i32* %8
%12 = mul nsw i32 %10, %11
store i32 %12, i32* %9
br label %13
13: ; preds = %22, %4
%14 = load i32, i32* %9
%15 = icmp sgt i32 %14, 0
br i1 %15, label %16, label %25
16: ; preds = %13
%17 = load double, double* %5
%18 = load double*, double** %6
%19 = getelementptr inbounds double, double* %18, i32 1
store double* %19, double** %6
%20 = load double, double* %18
%21 = fmul double %20, %17
store double %21, double* %18
br label %22
22: ; preds = %16
%23 = load i32, i32* %9
%24 = add nsw i32 %23, -1
store i32 %24, i32* %9
br label %13
25: ; preds = %13
ret void
}

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