void putdata2d(int n1, int n2, double **a) {
  int j1, j2, seed = 0;
  for (j1 = 0; j1 <= n1 - 1; j1++) {
    for (j2 = 0; j2 <= n2 - 1; j2++) {
      a[j1][j2] =
          ((*(&seed) = (*(&seed) * 7141 + 54773) % 259200) * (1.0 / 259200));
    }
  }
}


x86 assembly (clang, -O0):
<putdata2d>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     11a2 <putdata2d+0xa2>
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     118f <putdata2d+0x8f>
imul   $0x1be5,-0x1c(%rbp),%eax
add    $0xd5f5,%eax
mov    $0x3f480,%ecx
cltd
idiv   %ecx
mov    %edx,-0x1c(%rbp)
cvtsi2sd %edx,%xmm0
movsd  0xe98(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1134 <putdata2d+0x34>
jmp    1194 <putdata2d+0x94>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    111c <putdata2d+0x1c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<putdata2d>:
push   %r14
push   %rbx
test   %edi,%edi
jle    123f <putdata2d+0x13f>
mov    %edi,%r9d
mov    %esi,%r8d
mov    %r8d,%r11d
and    $0xfffffffe,%r11d
xor    %r10d,%r10d
movsd  0xedd(%rip),%xmm0
xor    %edi,%edi
jmp    113d <putdata2d+0x3d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r10
cmp    %r9,%r10
je     123f <putdata2d+0x13f>
test   %esi,%esi
jle    1130 <putdata2d+0x30>
mov    (%rdx,%r10,8),%r14
xor    %ecx,%ecx
cmp    $0x1,%esi
je     11ea <putdata2d+0xea>
imul   $0x1be5,%edi,%edi
mov    %edi,%ebx
add    $0xd5f5,%ebx
movslq %ebx,%rbx
imul   $0xffffffff81742e05,%rbx,%rax
shr    $0x20,%rax
add    %edi,%eax
add    $0xd5f5,%eax
mov    %eax,%edi
shr    $0x1f,%edi
sar    $0x11,%eax
add    %edi,%eax
imul   $0x3f480,%eax,%eax
sub    %eax,%ebx
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%rcx,8)
imul   $0x1be5,%ebx,%eax
mov    %eax,%edi
add    $0xd5f5,%edi
movslq %edi,%rdi
imul   $0xffffffff81742e05,%rdi,%rbx
shr    $0x20,%rbx
add    %ebx,%eax
add    $0xd5f5,%eax
mov    %eax,%ebx
shr    $0x1f,%ebx
sar    $0x11,%eax
add    %ebx,%eax
imul   $0x3f480,%eax,%eax
sub    %eax,%edi
xorps  %xmm1,%xmm1
cvtsi2sd %edi,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%r14,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r11
jne    1150 <putdata2d+0x50>
test   $0x1,%r8b
je     1130 <putdata2d+0x30>
imul   $0x1be5,%edi,%eax
mov    %eax,%edi
add    $0xd5f5,%edi
movslq %edi,%rdi
imul   $0xffffffff81742e05,%rdi,%rbx
shr    $0x20,%rbx
add    %ebx,%eax
add    $0xd5f5,%eax
mov    %eax,%ebx
shr    $0x1f,%ebx
sar    $0x11,%eax
add    %ebx,%eax
imul   $0x3f480,%eax,%eax
sub    %eax,%edi
xorps  %xmm1,%xmm1
cvtsi2sd %edi,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%rcx,8)
jmp    1130 <putdata2d+0x30>
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define void @putdata2d(i32 %0, i32 %1, double** nocapture %2) {
%4 = icmp sgt i32 %0, 0
br i1 %4, label %5, label %54
5: ; preds = %3
%6 = icmp sgt i32 %1, 0
%7 = zext i32 %0 to i64
%8 = zext i32 %1 to i64
%9 = and i64 %8, 1
%10 = icmp eq i32 %1, 1
%11 = and i64 %8, 4294967294
%12 = icmp eq i64 %9, 0
br label %13
13: ; preds = %5, %50
%14 = phi i64 [ 0, %5 ], [ %52, %50 ]
%15 = phi i32 [ 0, %5 ], [ %51, %50 ]
br i1 %6, label %16, label %50
16: ; preds = %13
%17 = getelementptr inbounds double*, double** %2, i64 %14
%18 = load double*, double** %17
br i1 %10, label %39, label %19
19: ; preds = %16, %19
%20 = phi i64 [ %36, %19 ], [ 0, %16 ]
%21 = phi i32 [ %32, %19 ], [ %15, %16 ]
%22 = phi i64 [ %37, %19 ], [ 0, %16 ]
%23 = mul nsw i32 %21, 7141
%24 = add nsw i32 %23, 54773
%25 = srem i32 %24, 259200
%26 = sitofp i32 %25 to double
%27 = fmul double %26, 0x3ED02E85C0898B71
%28 = getelementptr inbounds double, double* %18, i64 %20
store double %27, double* %28
%29 = or i64 %20, 1
%30 = mul nsw i32 %25, 7141
%31 = add nsw i32 %30, 54773
%32 = srem i32 %31, 259200
%33 = sitofp i32 %32 to double
%34 = fmul double %33, 0x3ED02E85C0898B71
%35 = getelementptr inbounds double, double* %18, i64 %29
store double %34, double* %35
%36 = add nuw nsw i64 %20, 2
%37 = add i64 %22, 2
%38 = icmp eq i64 %37, %11
br i1 %38, label %39, label %19
39: ; preds = %19, %16
%40 = phi i32 [ undef, %16 ], [ %32, %19 ]
%41 = phi i64 [ 0, %16 ], [ %36, %19 ]
%42 = phi i32 [ %15, %16 ], [ %32, %19 ]
br i1 %12, label %50, label %43
43: ; preds = %39
%44 = mul nsw i32 %42, 7141
%45 = add nsw i32 %44, 54773
%46 = srem i32 %45, 259200
%47 = sitofp i32 %46 to double
%48 = fmul double %47, 0x3ED02E85C0898B71
%49 = getelementptr inbounds double, double* %18, i64 %41
store double %48, double* %49
br label %50
50: ; preds = %43, %39, %13
%51 = phi i32 [ %15, %13 ], [ %40, %39 ], [ %46, %43 ]
%52 = add nuw nsw i64 %14, 1
%53 = icmp eq i64 %52, %7
br i1 %53, label %54, label %13
54: ; preds = %50, %3
ret void
}