void init_array(int ni, int nj, int nk, int nl, int nm,
                double A[800 + 0][1000 + 0], double B[1000 + 0][900 + 0],
                double C[900 + 0][1200 + 0], double D[1200 + 0][1100 + 0]) {
  int i, j;
  for (i = 0; i < ni; i++)
    for (j = 0; j < nk; j++)
      A[i][j] = (double)(i * j % ni) / (5 * ni);
  for (i = 0; i < nk; i++)
    for (j = 0; j < nj; j++)
      B[i][j] = (double)(i * (j + 1) % nj) / (5 * nj);
  for (i = 0; i < nj; i++)
    for (j = 0; j < nm; j++)
      C[i][j] = (double)(i * (j + 3) % nl) / (5 * nl);
  for (i = 0; i < nm; i++)
    for (j = 0; j < nl; j++)
      D[i][j] = (double)(i * (j + 2) % nk) / (5 * nk);
}


x86 assembly (clang, -O0):
<init_array>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8d,-0x14(%rbp)
mov    %r9,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11a1 <init_array+0xa1>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    118e <init_array+0x8e>
mov    -0x24(%rbp),%eax
imul   -0x28(%rbp),%eax
cltd
idivl  -0x4(%rbp)
cvtsi2sd %edx,%xmm0
imul   $0x5,-0x4(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x1f40,%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    113e <init_array+0x3e>
jmp    1193 <init_array+0x93>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112b <init_array+0x2b>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1223 <init_array+0x123>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1210 <init_array+0x110>
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
cltd
idivl  -0x8(%rbp)
cvtsi2sd %edx,%xmm0
imul   $0x5,-0x8(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x1c20,%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    11bb <init_array+0xbb>
jmp    1215 <init_array+0x115>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11a8 <init_array+0xa8>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12a5 <init_array+0x1a5>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1292 <init_array+0x192>
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
add    $0x3,%ecx
imul   %ecx,%eax
cltd
idivl  -0x10(%rbp)
cvtsi2sd %edx,%xmm0
imul   $0x5,-0x10(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x2580,%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    123d <init_array+0x13d>
jmp    1297 <init_array+0x197>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    122a <init_array+0x12a>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1327 <init_array+0x227>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1314 <init_array+0x214>
mov    -0x24(%rbp),%eax
mov    -0x28(%rbp),%ecx
add    $0x2,%ecx
imul   %ecx,%eax
cltd
idivl  -0xc(%rbp)
cvtsi2sd %edx,%xmm0
imul   $0x5,-0xc(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x2260,%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    12bf <init_array+0x1bf>
jmp    1319 <init_array+0x219>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12ac <init_array+0x1ac>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init_array>:
push   %r15
push   %r14
push   %r12
push   %rbx
mov    %edx,%r10d
mov    %edx,%r14d
test   %edi,%edi
jle    1165 <init_array+0x65>
lea    (%rdi,%rdi,4),%eax
cvtsi2sd %eax,%xmm0
mov    %edi,%r15d
xor    %r12d,%r12d
jmp    1130 <init_array+0x30>
add    $0x1,%r12
add    $0x1f40,%r9
cmp    %r15,%r12
je     1165 <init_array+0x65>
test   %r10d,%r10d
jle    1120 <init_array+0x20>
xor    %r11d,%r11d
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
mov    %r11d,%eax
cltd
idiv   %edi
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r9,%rbx,8)
add    $0x1,%rbx
add    %r12d,%r11d
cmp    %rbx,%r14
jne    1140 <init_array+0x40>
jmp    1120 <init_array+0x20>
mov    %esi,%r9d
test   %r10d,%r10d
jle    11d8 <init_array+0xd8>
mov    0x28(%rsp),%r11
lea    (%rsi,%rsi,4),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xor    %r15d,%r15d
jmp    11a0 <init_array+0xa0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r15
add    $0x1c20,%r11
cmp    %r14,%r15
je     11d8 <init_array+0xd8>
test   %esi,%esi
jle    1190 <init_array+0x90>
mov    %r15d,%edi
xor    %ebx,%ebx
nopl   0x0(%rax)
mov    %edi,%eax
cltd
idiv   %r9d
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r11,%rbx,8)
lea    0x1(%rbx),%rax
add    %r15d,%edi
mov    %rax,%rbx
cmp    %rax,%r9
jne    11b0 <init_array+0xb0>
jmp    1190 <init_array+0x90>
mov    %r8d,%r11d
test   %esi,%esi
jle    1243 <init_array+0x143>
mov    0x30(%rsp),%r15
lea    (%rcx,%rcx,4),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xor    %r14d,%r14d
xor    %edi,%edi
jmp    1214 <init_array+0x114>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rdi
add    $0x3,%r14d
add    $0x2580,%r15
cmp    %r9,%rdi
je     1243 <init_array+0x143>
test   %r8d,%r8d
jle    1200 <init_array+0x100>
mov    %r14d,%esi
xor    %ebx,%ebx
xchg   %ax,%ax
mov    %esi,%eax
cltd
idiv   %ecx
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r15,%rbx,8)
add    $0x1,%rbx
add    %edi,%esi
cmp    %rbx,%r11
jne    1220 <init_array+0x120>
jmp    1200 <init_array+0x100>
test   %r8d,%r8d
jle    12b4 <init_array+0x1b4>
lea    (%r10,%r10,4),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mov    0x38(%rsp),%r14
mov    %ecx,%r9d
xor    %r8d,%r8d
xor    %edi,%edi
jmp    1284 <init_array+0x184>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rdi
add    $0x2,%r8d
add    $0x2260,%r14
cmp    %r11,%rdi
je     12b4 <init_array+0x1b4>
test   %ecx,%ecx
jle    1270 <init_array+0x170>
mov    %r8d,%esi
xor    %ebx,%ebx
nopl   (%rax)
mov    %esi,%eax
cltd
idiv   %r10d
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%rbx,8)
add    $0x1,%rbx
add    %edi,%esi
cmp    %rbx,%r9
jne    1290 <init_array+0x190>
jmp    1270 <init_array+0x170>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @init_array(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, [1000 x double]* %5, [900 x double]* %6, [1200 x double]* %7, [1100 x double]* %8) {
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca [1000 x double]*
%16 = alloca [900 x double]*
%17 = alloca [1200 x double]*
%18 = alloca [1100 x double]*
%19 = alloca i32
%20 = alloca i32
store i32 %0, i32* %10
store i32 %1, i32* %11
store i32 %2, i32* %12
store i32 %3, i32* %13
store i32 %4, i32* %14
store [1000 x double]* %5, [1000 x double]** %15
store [900 x double]* %6, [900 x double]** %16
store [1200 x double]* %7, [1200 x double]** %17
store [1100 x double]* %8, [1100 x double]** %18
store i32 0, i32* %19
br label %21
21: ; preds = %52, %9
%22 = load i32, i32* %19
%23 = load i32, i32* %10
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %55
25: ; preds = %21
store i32 0, i32* %20
br label %26
26: ; preds = %48, %25
%27 = load i32, i32* %20
%28 = load i32, i32* %12
%29 = icmp slt i32 %27, %28
br i1 %29, label %30, label %51
30: ; preds = %26
%31 = load i32, i32* %19
%32 = load i32, i32* %20
%33 = mul nsw i32 %31, %32
%34 = load i32, i32* %10
%35 = srem i32 %33, %34
%36 = sitofp i32 %35 to double
%37 = load i32, i32* %10
%38 = mul nsw i32 5, %37
%39 = sitofp i32 %38 to double
%40 = fdiv double %36, %39
%41 = load [1000 x double]*, [1000 x double]** %15
%42 = load i32, i32* %19
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds [1000 x double], [1000 x double]* %41, i64 %43
%45 = load i32, i32* %20
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds [1000 x double], [1000 x double]* %44, i64 0, i64 %46
store double %40, double* %47
br label %48
48: ; preds = %30
%49 = load i32, i32* %20
%50 = add nsw i32 %49, 1
store i32 %50, i32* %20
br label %26
51: ; preds = %26
br label %52
52: ; preds = %51
%53 = load i32, i32* %19
%54 = add nsw i32 %53, 1
store i32 %54, i32* %19
br label %21
55: ; preds = %21
store i32 0, i32* %19
br label %56
56: ; preds = %88, %55
%57 = load i32, i32* %19
%58 = load i32, i32* %12
%59 = icmp slt i32 %57, %58
br i1 %59, label %60, label %91
60: ; preds = %56
store i32 0, i32* %20
br label %61
61: ; preds = %84, %60
%62 = load i32, i32* %20
%63 = load i32, i32* %11
%64 = icmp slt i32 %62, %63
br i1 %64, label %65, label %87
65: ; preds = %61
%66 = load i32, i32* %19
%67 = load i32, i32* %20
%68 = add nsw i32 %67, 1
%69 = mul nsw i32 %66, %68
%70 = load i32, i32* %11
%71 = srem i32 %69, %70
%72 = sitofp i32 %71 to double
%73 = load i32, i32* %11
%74 = mul nsw i32 5, %73
%75 = sitofp i32 %74 to double
%76 = fdiv double %72, %75
%77 = load [900 x double]*, [900 x double]** %16
%78 = load i32, i32* %19
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds [900 x double], [900 x double]* %77, i64 %79
%81 = load i32, i32* %20
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds [900 x double], [900 x double]* %80, i64 0, i64 %82
store double %76, double* %83
br label %84
84: ; preds = %65
%85 = load i32, i32* %20
%86 = add nsw i32 %85, 1
store i32 %86, i32* %20
br label %61
87: ; preds = %61
br label %88
88: ; preds = %87
%89 = load i32, i32* %19
%90 = add nsw i32 %89, 1
store i32 %90, i32* %19
br label %56
91: ; preds = %56
store i32 0, i32* %19
br label %92
92: ; preds = %124, %91
%93 = load i32, i32* %19
%94 = load i32, i32* %11
%95 = icmp slt i32 %93, %94
br i1 %95, label %96, label %127
96: ; preds = %92
store i32 0, i32* %20
br label %97
97: ; preds = %120, %96
%98 = load i32, i32* %20
%99 = load i32, i32* %14
%100 = icmp slt i32 %98, %99
br i1 %100, label %101, label %123
101: ; preds = %97
%102 = load i32, i32* %19
%103 = load i32, i32* %20
%104 = add nsw i32 %103, 3
%105 = mul nsw i32 %102, %104
%106 = load i32, i32* %13
%107 = srem i32 %105, %106
%108 = sitofp i32 %107 to double
%109 = load i32, i32* %13
%110 = mul nsw i32 5, %109
%111 = sitofp i32 %110 to double
%112 = fdiv double %108, %111
%113 = load [1200 x double]*, [1200 x double]** %17
%114 = load i32, i32* %19
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds [1200 x double], [1200 x double]* %113, i64 %115
%117 = load i32, i32* %20
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds [1200 x double], [1200 x double]* %116, i64 0, i64 %118
store double %112, double* %119
br label %120
120: ; preds = %101
%121 = load i32, i32* %20
%122 = add nsw i32 %121, 1
store i32 %122, i32* %20
br label %97
123: ; preds = %97
br label %124
124: ; preds = %123
%125 = load i32, i32* %19
%126 = add nsw i32 %125, 1
store i32 %126, i32* %19
br label %92
127: ; preds = %92
store i32 0, i32* %19
br label %128
128: ; preds = %160, %127
%129 = load i32, i32* %19
%130 = load i32, i32* %14
%131 = icmp slt i32 %129, %130
br i1 %131, label %132, label %163
132: ; preds = %128
store i32 0, i32* %20
br label %133
133: ; preds = %156, %132
%134 = load i32, i32* %20
%135 = load i32, i32* %13
%136 = icmp slt i32 %134, %135
br i1 %136, label %137, label %159
137: ; preds = %133
%138 = load i32, i32* %19
%139 = load i32, i32* %20
%140 = add nsw i32 %139, 2
%141 = mul nsw i32 %138, %140
%142 = load i32, i32* %12
%143 = srem i32 %141, %142
%144 = sitofp i32 %143 to double
%145 = load i32, i32* %12
%146 = mul nsw i32 5, %145
%147 = sitofp i32 %146 to double
%148 = fdiv double %144, %147
%149 = load [1100 x double]*, [1100 x double]** %18
%150 = load i32, i32* %19
%151 = sext i32 %150 to i64
%152 = getelementptr inbounds [1100 x double], [1100 x double]* %149, i64 %151
%153 = load i32, i32* %20
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds [1100 x double], [1100 x double]* %152, i64 0, i64 %154
store double %148, double* %155
br label %156
156: ; preds = %137
%157 = load i32, i32* %20
%158 = add nsw i32 %157, 1
store i32 %158, i32* %20
br label %133
159: ; preds = %133
br label %160
160: ; preds = %159
%161 = load i32, i32* %19
%162 = add nsw i32 %161, 1
store i32 %162, i32* %19
br label %128
163: ; preds = %128
ret void
}

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