void inidat(nx, ny, u1) int nx, ny;
float *u1;
{
  int ix, iy;
  for (ix = 0; ix <= nx - 1; ix++)
    for (iy = 0; iy <= ny - 1; iy++)
      *(u1 + ix * ny + iy) = (float)(ix * (nx - ix - 1) * iy * (ny - iy - 1));
}


x86 assembly (clang, -O0):
<inidat>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     11a0 <inidat+0xa0>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     118d <inidat+0x8d>
mov    -0x14(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
imul   -0x18(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
cvtsi2ss %eax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x8(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    112d <inidat+0x2d>
jmp    1192 <inidat+0x92>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1115 <inidat+0x15>
pop    %rbp
ret

x86 assembly (clang, -O2):
<inidat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdx,-0x38(%rsp)
mov    %rdi,-0x30(%rsp)
test   %edi,%edi
jle    1379 <inidat+0x279>
movslq %esi,%rcx
mov    -0x30(%rsp),%rax
mov    %eax,%r15d
mov    %ecx,%r9d
mov    %r9d,%edx
and    $0xfffffffc,%edx
mov    %rdx,-0x8(%rsp)
add    $0xfffffffffffffffc,%rdx
mov    %rdx,-0x10(%rsp)
mov    %rdx,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
movd   %ecx,%xmm0
pshufd $0x0,%xmm0,%xmm0
mov    %rdi,%rdx
mov    %rdi,-0x18(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,-0x20(%rsp)
mov    -0x38(%rsp),%rdx
lea    0x10(%rdx),%r12
lea    0x0(,%rcx,4),%r13
lea    -0x2(%rax),%r11d
mov    %rcx,-0x28(%rsp)
lea    -0x1(%rcx),%eax
mov    %eax,-0x3c(%rsp)
xor    %ecx,%ecx
movdqa 0xe71(%rip),%xmm8
pcmpeqd %xmm2,%xmm2
movdqa 0xe74(%rip),%xmm9
movdqa 0xe7b(%rip),%xmm10
movdqa 0xe82(%rip),%xmm11
xor    %r8d,%r8d
jmp    11da <inidat+0xda>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r8
add    %r13,%r12
add    %r11d,%ecx
add    $0xfffffffe,%r11d
add    %r13,%rdx
cmp    %r15,%r8
je     1379 <inidat+0x279>
test   %esi,%esi
jle    11c0 <inidat+0xc0>
cmp    $0x4,%esi
jae    11f0 <inidat+0xf0>
xor    %ebp,%ebp
jmp    132c <inidat+0x22c>
nopw   0x0(%rax,%rax,1)
mov    %r8d,%eax
not    %eax
add    -0x30(%rsp),%eax
imul   %r8d,%eax
movd   %eax,%xmm1
pshufd $0x0,%xmm1,%xmm4
pshufd $0xf5,%xmm4,%xmm6
cmpq   $0x0,-0x10(%rsp)
je     1364 <inidat+0x264>
mov    -0x20(%rsp),%r10
xor    %r14d,%r14d
movdqa %xmm8,%xmm5
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
pshufd $0xf5,%xmm5,%xmm3
pmuludq %xmm6,%xmm3
movdqa %xmm4,%xmm7
pmuludq %xmm5,%xmm7
movdqa %xmm5,%xmm1
pxor   %xmm2,%xmm1
paddd  %xmm0,%xmm1
pmuludq %xmm1,%xmm7
pshufd $0xe8,%xmm7,%xmm7
pshufd $0xf5,%xmm1,%xmm1
pmuludq %xmm3,%xmm1
pshufd $0xe8,%xmm1,%xmm1
punpckldq %xmm1,%xmm7
cvtdq2ps %xmm7,%xmm1
movups %xmm1,-0x10(%r12,%r14,4)
movdqa %xmm5,%xmm1
paddd  %xmm9,%xmm1
pshufd $0xf5,%xmm1,%xmm3
pmuludq %xmm6,%xmm3
pmuludq %xmm4,%xmm1
movdqa %xmm0,%xmm7
psubd  %xmm5,%xmm7
paddd  %xmm10,%xmm7
pmuludq %xmm7,%xmm1
pshufd $0xe8,%xmm1,%xmm1
pshufd $0xf5,%xmm7,%xmm7
pmuludq %xmm3,%xmm7
pshufd $0xe8,%xmm7,%xmm3
punpckldq %xmm3,%xmm1
cvtdq2ps %xmm1,%xmm1
movups %xmm1,(%r12,%r14,4)
add    $0x8,%r14
paddd  %xmm11,%xmm5
add    $0xfffffffffffffffe,%r10
jne    1230 <inidat+0x130>
testb  $0x1,-0x18(%rsp)
je     131b <inidat+0x21b>
mov    %r8,%rax
imul   -0x28(%rsp),%rax
mov    -0x38(%rsp),%rdi
lea    (%rdi,%rax,4),%rax
pshufd $0xf5,%xmm5,%xmm1
pmuludq %xmm1,%xmm6
pmuludq %xmm5,%xmm4
pxor   %xmm2,%xmm5
paddd  %xmm0,%xmm5
pmuludq %xmm5,%xmm4
pshufd $0xe8,%xmm4,%xmm1
pshufd $0xf5,%xmm5,%xmm3
pmuludq %xmm6,%xmm3
pshufd $0xe8,%xmm3,%xmm3
punpckldq %xmm3,%xmm1
cvtdq2ps %xmm1,%xmm1
movups %xmm1,(%rax,%r14,4)
mov    -0x8(%rsp),%rax
mov    %rax,%rbp
cmp    %r9,%rax
je     11c0 <inidat+0xc0>
mov    %ecx,%ebx
imul   %ebp,%ebx
mov    -0x3c(%rsp),%eax
sub    %ebp,%eax
nopw   0x0(%rax,%rax,1)
mov    %ebx,%edi
imul   %eax,%edi
xorps  %xmm1,%xmm1
cvtsi2ss %edi,%xmm1
movss  %xmm1,(%rdx,%rbp,4)
add    $0x1,%rbp
add    %ecx,%ebx
add    $0xffffffff,%eax
cmp    %rbp,%r9
jne    1340 <inidat+0x240>
jmp    11c0 <inidat+0xc0>
xor    %r14d,%r14d
movdqa %xmm8,%xmm5
testb  $0x1,-0x18(%rsp)
jne    12d1 <inidat+0x1d1>
jmp    131b <inidat+0x21b>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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