void init_array(int m, int n, double A[2100 + 0][1900 + 0], double r[2100 + 0],
                double p[1900 + 0]) {
  int i, j;
  for (i = 0; i < m; i++)
    p[i] = (double)(i % m) / m;
  for (i = 0; i < n; i++) {
    r[i] = (double)(i % n) / n;
    for (j = 0; j < m; j++)
      A[i][j] = (double)(i * (j + 1) % n) / n;
  }
}


x86 assembly (clang, -O0):
<init_array>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1158 <init_array+0x58>
mov    -0x24(%rbp),%eax
cltd
idivl  -0x4(%rbp)
cvtsi2sd %edx,%xmm0
cvtsi2sdl -0x4(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111d <init_array+0x1d>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11f8 <init_array+0xf8>
mov    -0x24(%rbp),%eax
cltd
idivl  -0x8(%rbp)
cvtsi2sd %edx,%xmm0
cvtsi2sdl -0x8(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11e5 <init_array+0xe5>
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
cltd
idivl  -0x8(%rbp)
cvtsi2sd %edx,%xmm0
cvtsi2sdl -0x8(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x3b60,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1193 <init_array+0x93>
jmp    11ea <init_array+0xea>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    115f <init_array+0x5f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init_array>:
push   %r14
push   %rbx
mov    %rdx,%r9
mov    %edi,%r10d
test   %edi,%edi
jle    11da <init_array+0xda>
cvtsi2sd %edi,%xmm0
cmp    $0x1,%edi
jne    1121 <init_array+0x21>
xor    %eax,%eax
jmp    11c0 <init_array+0xc0>
mov    %r10d,%eax
and    $0xfffffffe,%eax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rax),%rdx
mov    %rdx,%r11
shr    %r11
add    $0x1,%r11
test   %rdx,%rdx
je     124c <init_array+0x14c>
mov    %r11,%rbx
and    $0xfffffffffffffffe,%rbx
movdqa 0xebb(%rip),%xmm2
xor    %edx,%edx
movdqa 0xec1(%rip),%xmm3
movdqa 0xec9(%rip),%xmm4
nopw   0x0(%rax,%rax,1)
cvtdq2pd %xmm2,%xmm5
divpd  %xmm1,%xmm5
movupd %xmm5,(%r8,%rdx,8)
movdqa %xmm2,%xmm5
paddd  %xmm3,%xmm5
cvtdq2pd %xmm5,%xmm5
divpd  %xmm1,%xmm5
movupd %xmm5,0x10(%r8,%rdx,8)
add    $0x4,%rdx
paddd  %xmm4,%xmm2
add    $0xfffffffffffffffe,%rbx
jne    1170 <init_array+0x70>
cvtdq2pd %xmm2,%xmm2
test   $0x1,%r11b
je     11b7 <init_array+0xb7>
divpd  %xmm1,%xmm2
movupd %xmm2,(%r8,%rdx,8)
cmp    %r10,%rax
je     11da <init_array+0xda>
nopl   0x0(%rax)
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r8,%rax,8)
add    $0x1,%rax
cmp    %rax,%r10
jne    11c0 <init_array+0xc0>
test   %esi,%esi
jle    1248 <init_array+0x148>
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
mov    %esi,%r11d
xor    %r14d,%r14d
jmp    1200 <init_array+0x100>
nopl   (%rax)
add    $0x1,%r14
add    $0x3b60,%r9
cmp    %r11,%r14
je     1248 <init_array+0x148>
xorps  %xmm1,%xmm1
cvtsi2sd %r14d,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx,%r14,8)
test   %edi,%edi
jle    11f0 <init_array+0xf0>
mov    %r14d,%r8d
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
mov    %r8d,%eax
cltd
idiv   %esi
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r9,%rbx,8)
lea    0x1(%rbx),%rax
add    %r14d,%r8d
mov    %rax,%rbx
cmp    %rax,%r10
jne    1220 <init_array+0x120>
jmp    11f0 <init_array+0xf0>
pop    %rbx
pop    %r14
ret
xorpd  %xmm2,%xmm2
movhpd 0xda8(%rip),%xmm2
xor    %edx,%edx
test   $0x1,%r11b
jne    11ad <init_array+0xad>
jmp    11b7 <init_array+0xb7>

LLVM IR (-O0):
define void @init_array(i32 %0, i32 %1, [1900 x double]* %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca [1900 x double]*
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32
%12 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store [1900 x double]* %2, [1900 x double]** %8
store double* %3, double** %9
store double* %4, double** %10
store i32 0, i32* %11
br label %13
13: ; preds = %29, %5
%14 = load i32, i32* %11
%15 = load i32, i32* %6
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %32
17: ; preds = %13
%18 = load i32, i32* %11
%19 = load i32, i32* %6
%20 = srem i32 %18, %19
%21 = sitofp i32 %20 to double
%22 = load i32, i32* %6
%23 = sitofp i32 %22 to double
%24 = fdiv double %21, %23
%25 = load double*, double** %10
%26 = load i32, i32* %11
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
store double %24, double* %28
br label %29
29: ; preds = %17
%30 = load i32, i32* %11
%31 = add nsw i32 %30, 1
store i32 %31, i32* %11
br label %13
32: ; preds = %13
store i32 0, i32* %11
br label %33
33: ; preds = %75, %32
%34 = load i32, i32* %11
%35 = load i32, i32* %7
%36 = icmp slt i32 %34, %35
br i1 %36, label %37, label %78
37: ; preds = %33
%38 = load i32, i32* %11
%39 = load i32, i32* %7
%40 = srem i32 %38, %39
%41 = sitofp i32 %40 to double
%42 = load i32, i32* %7
%43 = sitofp i32 %42 to double
%44 = fdiv double %41, %43
%45 = load double*, double** %9
%46 = load i32, i32* %11
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %45, i64 %47
store double %44, double* %48
store i32 0, i32* %12
br label %49
49: ; preds = %71, %37
%50 = load i32, i32* %12
%51 = load i32, i32* %6
%52 = icmp slt i32 %50, %51
br i1 %52, label %53, label %74
53: ; preds = %49
%54 = load i32, i32* %11
%55 = load i32, i32* %12
%56 = add nsw i32 %55, 1
%57 = mul nsw i32 %54, %56
%58 = load i32, i32* %7
%59 = srem i32 %57, %58
%60 = sitofp i32 %59 to double
%61 = load i32, i32* %7
%62 = sitofp i32 %61 to double
%63 = fdiv double %60, %62
%64 = load [1900 x double]*, [1900 x double]** %8
%65 = load i32, i32* %11
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds [1900 x double], [1900 x double]* %64, i64 %66
%68 = load i32, i32* %12
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds [1900 x double], [1900 x double]* %67, i64 0, i64 %69
store double %63, double* %70
br label %71
71: ; preds = %53
%72 = load i32, i32* %12
%73 = add nsw i32 %72, 1
store i32 %73, i32* %12
br label %49
74: ; preds = %49
br label %75
75: ; preds = %74
%76 = load i32, i32* %11
%77 = add nsw i32 %76, 1
store i32 %77, i32* %11
br label %33
78: ; preds = %33
ret void
}

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