void init_array(int w, int h, float *alpha, float imgIn[4096 + 0][2160 + 0],
                float imgOut[4096 + 0][2160 + 0]) {
  int i, j;
  *alpha = 0.25;
  for (i = 0; i < w; i++)
    for (j = 0; j < h; j++)
      imgIn[i][j] = (float)((313 * i + 991 * j) % 65536) / 65535.0f;
}


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)
mov    -0x10(%rbp),%rax
movss  0xede(%rip),%xmm0
movss  %xmm0,(%rax)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11b0 <init_array+0xb0>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    119d <init_array+0x9d>
imul   $0x139,-0x24(%rbp),%eax
imul   $0x3df,-0x28(%rbp),%ecx
add    %ecx,%eax
mov    $0x10000,%ecx
cltd
idiv   %ecx
cvtsi2ss %edx,%xmm0
movss  0xe94(%rip),%xmm1
divss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0x21c0,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1140 <init_array+0x40>
jmp    11a2 <init_array+0xa2>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112d <init_array+0x2d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init_array>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
movl   $0x3e800000,(%rdx)
test   %edi,%edi
jle    12ea <init_array+0x1ea>
mov    %edi,%r14d
mov    %esi,%r13d
mov    %r13d,%r11d
and    $0xfffffffc,%r11d
lea    -0x4(%r11),%r8
mov    %r8,%r9
shr    $0x2,%r9
add    $0x1,%r9
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
lea    0x10(%rcx),%rdi
xor    %r15d,%r15d
movss  0xf16(%rip),%xmm0
movdqa 0xead(%rip),%xmm8
movdqa 0xeb5(%rip),%xmm2
movdqa 0xebd(%rip),%xmm3
movaps 0xec6(%rip),%xmm4
movdqa 0xecd(%rip),%xmm9
movdqa 0xed4(%rip),%xmm10
mov    %rcx,%rax
xor    %r12d,%r12d
jmp    11b1 <init_array+0xb1>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r12
add    $0x21c0,%rdi
add    $0x139,%r15d
add    $0x21c0,%rax
cmp    %r14,%r12
je     12ea <init_array+0x1ea>
test   %esi,%esi
jle    1190 <init_array+0x90>
cmp    $0x4,%esi
jae    11d0 <init_array+0xd0>
xor    %edx,%edx
jmp    12a1 <init_array+0x1a1>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
imul   $0x139,%r12d,%edx
movd   %edx,%xmm1
pshufd $0x0,%xmm1,%xmm7
test   %r8,%r8
je     12d7 <init_array+0x1d7>
mov    %r10,%rdx
xor    %ebx,%ebx
movdqa %xmm8,%xmm5
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movdqa %xmm5,%xmm1
pmuludq %xmm2,%xmm1
pshufd $0xe8,%xmm1,%xmm1
pshufd $0xf5,%xmm5,%xmm6
pmuludq %xmm2,%xmm6
pshufd $0xe8,%xmm6,%xmm6
punpckldq %xmm6,%xmm1
paddd  %xmm7,%xmm1
movdqa %xmm1,%xmm6
pand   %xmm3,%xmm6
cvtdq2ps %xmm6,%xmm6
divps  %xmm4,%xmm6
movups %xmm6,-0x10(%rdi,%rbx,4)
paddd  %xmm9,%xmm1
pand   %xmm3,%xmm1
cvtdq2ps %xmm1,%xmm1
divps  %xmm4,%xmm1
movups %xmm1,(%rdi,%rbx,4)
add    $0x8,%rbx
paddd  %xmm10,%xmm5
add    $0xfffffffffffffffe,%rdx
jne    1200 <init_array+0x100>
test   $0x1,%r9b
je     1295 <init_array+0x195>
pshufd $0xf5,%xmm5,%xmm1
pmuludq %xmm2,%xmm5
pshufd $0xe8,%xmm5,%xmm5
pmuludq %xmm2,%xmm1
pshufd $0xe8,%xmm1,%xmm1
punpckldq %xmm1,%xmm5
paddd  %xmm7,%xmm5
pand   %xmm3,%xmm5
cvtdq2ps %xmm5,%xmm1
divps  %xmm4,%xmm1
imul   $0x21c0,%r12,%rdx
add    %rcx,%rdx
movups %xmm1,(%rdx,%rbx,4)
mov    %r11,%rdx
cmp    %r13,%r11
je     1190 <init_array+0x90>
imul   $0x3df,%edx,%ebx
add    %r15d,%ebx
nopw   0x0(%rax,%rax,1)
movzwl %bx,%ebp
xorps  %xmm1,%xmm1
cvtsi2ss %ebp,%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,(%rax,%rdx,4)
add    $0x1,%rdx
add    $0x3df,%ebx
cmp    %rdx,%r13
jne    12b0 <init_array+0x1b0>
jmp    1190 <init_array+0x90>
xor    %ebx,%ebx
movdqa %xmm8,%xmm5
test   $0x1,%r9b
jne    125e <init_array+0x15e>
jmp    1295 <init_array+0x195>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @init_array(i32 %0, i32 %1, float* %2, [2160 x float]* %3, [2160 x float]* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca float*
%9 = alloca [2160 x float]*
%10 = alloca [2160 x float]*
%11 = alloca i32
%12 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store float* %2, float** %8
store [2160 x float]* %3, [2160 x float]** %9
store [2160 x float]* %4, [2160 x float]** %10
%13 = load float*, float** %8
store float 2.500000e-01, float* %13
store i32 0, i32* %11
br label %14
14: ; preds = %43, %5
%15 = load i32, i32* %11
%16 = load i32, i32* %6
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %46
18: ; preds = %14
store i32 0, i32* %12
br label %19
19: ; preds = %39, %18
%20 = load i32, i32* %12
%21 = load i32, i32* %7
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %42
23: ; preds = %19
%24 = load i32, i32* %11
%25 = mul nsw i32 313, %24
%26 = load i32, i32* %12
%27 = mul nsw i32 991, %26
%28 = add nsw i32 %25, %27
%29 = srem i32 %28, 65536
%30 = sitofp i32 %29 to float
%31 = fdiv float %30, 6.553500e+04
%32 = load [2160 x float]*, [2160 x float]** %9
%33 = load i32, i32* %11
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds [2160 x float], [2160 x float]* %32, i64 %34
%36 = load i32, i32* %12
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds [2160 x float], [2160 x float]* %35, i64 0, i64 %37
store float %31, float* %38
br label %39
39: ; preds = %23
%40 = load i32, i32* %12
%41 = add nsw i32 %40, 1
store i32 %41, i32* %12
br label %19
42: ; preds = %19
br label %43
43: ; preds = %42
%44 = load i32, i32* %11
%45 = add nsw i32 %44, 1
store i32 %45, i32* %11
br label %14
46: ; preds = %14
ret void
}

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