void divide_array_number(double *out, double *a, double f, int ncol, int nlns) {
  int i, j;
  double *ptr_a, *ptr_out;
  for (i = 0; i < ncol; i++) {
    ptr_out = (&(out[(i) * (nlns) + (0)]));
    ptr_a = (&(a[(i) * (nlns) + (0)]));
    for (j = 0; j < nlns; j++)
      ptr_out[j] = ptr_a[j] / f;
  }
}


x86 assembly (clang, -O0):
<divide_array_number>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11b5 <divide_array_number+0xb5>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x20(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11a2 <divide_array_number+0xa2>
mov    -0x30(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x18(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1169 <divide_array_number+0x69>
jmp    11a7 <divide_array_number+0xa7>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111e <divide_array_number+0x1e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<divide_array_number>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    12ba <divide_array_number+0x1ba>
movslq %ecx,%r15
mov    %edx,%eax
mov    %rax,-0x8(%rsp)
mov    %r15d,%r8d
mov    %r8d,%eax
and    $0xfffffffe,%eax
mov    %rax,-0x20(%rsp)
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x28(%rsp)
mov    %rax,%rdx
shr    %rdx
add    $0x1,%rdx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    %rdx,%rax
mov    %rdx,-0x30(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x38(%rsp)
mov    %r8,%rax
neg    %rax
mov    %rax,-0x10(%rsp)
lea    0x10(%rsi),%r10
lea    0x0(,%r15,8),%r13
lea    0x10(%rdi),%rdx
lea    0x8(%rdi),%rbp
xor    %r14d,%r14d
mov    %rsi,%r12
mov    %r15,-0x18(%rsp)
jmp    11ab <divide_array_number+0xab>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r14
add    %r13,%r10
add    %r13,%rdx
add    %r13,%rbp
add    %r13,%r12
cmp    -0x8(%rsp),%r14
je     12ba <divide_array_number+0x1ba>
test   %ecx,%ecx
jle    1190 <divide_array_number+0x90>
mov    %r14,%rax
imul   %r15,%rax
lea    (%rdi,%rax,8),%r11
lea    (%rsi,%rax,8),%r9
cmp    $0x1,%ecx
je     11d8 <divide_array_number+0xd8>
add    %r8,%rax
lea    (%rsi,%rax,8),%rbx
cmp    %rbx,%r11
jae    123f <divide_array_number+0x13f>
lea    (%rdi,%rax,8),%rax
cmp    %rax,%r9
jae    123f <divide_array_number+0x13f>
xor    %ebx,%ebx
mov    %rbx,%rax
test   $0x1,%r8b
je     11fa <divide_array_number+0xfa>
movsd  (%r9,%rbx,8),%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,(%r11,%rbx,8)
mov    %rbx,%rax
or     $0x1,%rax
not    %rbx
cmp    -0x10(%rsp),%rbx
je     1190 <divide_array_number+0x90>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r12,%rax,8),%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,-0x8(%rbp,%rax,8)
movsd  0x8(%r12,%rax,8),%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,0x0(%rbp,%rax,8)
add    $0x2,%rax
cmp    %rax,%r8
jne    1210 <divide_array_number+0x110>
jmp    1190 <divide_array_number+0x90>
cmpq   $0x0,-0x28(%rsp)
je     12ae <divide_array_number+0x1ae>
mov    -0x38(%rsp),%rax
xor    %r15d,%r15d
nop
movupd -0x10(%r10,%r15,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rdx,%r15,8)
movupd (%r10,%r15,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdx,%r15,8)
add    $0x4,%r15
add    $0xfffffffffffffffe,%rax
jne    1250 <divide_array_number+0x150>
testb  $0x1,-0x30(%rsp)
je     1293 <divide_array_number+0x193>
movupd (%r9,%r15,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%r11,%r15,8)
mov    -0x20(%rsp),%rax
mov    %rax,%rbx
cmp    %r8,%rax
mov    -0x18(%rsp),%r15
je     1190 <divide_array_number+0x90>
jmp    11da <divide_array_number+0xda>
xor    %r15d,%r15d
testb  $0x1,-0x30(%rsp)
jne    1283 <divide_array_number+0x183>
jmp    1293 <divide_array_number+0x193>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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