<table border="1" cellspacing="0" cellpadding="8">
    <tr>
        <th>Issue</th>
        <td>
            <a href=https://github.com/llvm/llvm-project/issues/106221>106221</a>
        </td>
    </tr>

    <tr>
        <th>Summary</th>
        <td>
            [flang][runtime] random_init collision rate is high (>1%) on some platforms and time dependent
        </td>
    </tr>

    <tr>
      <th>Labels</th>
      <td>
            enhancement,
            flang:runtime
      </td>
    </tr>

    <tr>
      <th>Assignees</th>
      <td>
      </td>
    </tr>

    <tr>
      <th>Reporter</th>
      <td>
          jeanPerier
      </td>
    </tr>
</table>

<pre>
    Two calls to random_init close to one another may init random to the same seed with a high probability (>20%) on certain platform at a given time (the collision probability seems to cycle with the time in periods of about 1 minute). This likely shows a lack or randomness in f18 random_init implementation.

This was initially caught by gfortran random_init_2.f90 test in LLVM test suite that sporadically failed: https://github.com/llvm/llvm-test-suite/blob/f6d8e5692b65387703a687a556c2f77111cbfaed/Fortran/gfortran/regression/random_init_2.f90

I went on and made a sampling program 
[sampling.f90](https://github.com/user-attachments/files/16762635/sampling.f90.txt)
 to check random_init collisions between two random_init  calls next to one another.

Attached is a graph that shows the collision rate (in %) in function of the runtime (in seconds).


f18 random_init clearly has a cyclic behavior with spikes where the collision rate goes above 20% for 20 seconds every 2 minutes or so, while it can be zero (no collisions out of 100000 runs) for periods of 20 seconds too. Attached is gnuplot graph and the program to sample it.

The same graph produced by a gfortran compiled exec on the same is mostly flat and never goes above 0.005% collision rate (most samples are 0., with sporadic 4.99999989E-03% samples). Interestingly enough, gfortran is 40x slower in executing random intrinsic calls.

</pre>
<img width="1px" height="1px" alt="" src="http://email.email.llvm.org/o/eJysVkuP2zYQ_jX0ZWBDot4HH3azMRAgBXoIeg1G1EhiliIFklqv--uLkeyNdxsURVEf9CLnwW---cYYgh4s0VEUj6J42uESR-ePPwjt7-Q1-V3rusvx29mBQmMCRAcebeem79rqCMq4QPzRWQK0Lo7kYcILrKvbTl6OI0HAiSAQdXDWcQSEUQ8jzN612Gqj4wWErEX2WSZCFkI24Cwo8hG1hdlg7J2fACMgDPqFLEQ9EZuwb-WM0UE7-85fIJrWlNVFGdrC8u7Vkr2S164L4HrA1i0RUpi0XSIJ2Rzg26gDGP1M5gJhdOcACAbVMzh_PZmlENhPn9bvUNHTbGgiGzFqZw8ieRLJw3ZdnZ4xrABpNOYCCpdhjNBeYOidjx7tvbPv8tA3CUQKkUN9_frHb9tLWHQkiCNGCLPz2Gm1uutRG-pE9gBjjHMQ2YOQJyFPg47j0h6Um4Q8GfNyu-3Z2371JuSpNa4V8tSXXU1F2ci2LLK6qpIMy7rCoiiV7KsqTVPV9kidkKfTljNH6N8ePQ2eAheEXz6e5h6QL3AmG7nWaDuYsCNApspstB24moPHCa4WxeNtZXVTPAlZ_8Mpl0B-jzGiGrkagQ-mDfE9LatSllkh5One5SG-RiGbLdzKnJHU83vK36gWoKV4Jmbi-X1XXHvF0mv80BvvuPCwZkYdaKbW4HEer_Vc2fae1x7jynZt4doeTLzFKqYYM5i3-8XeukJbCKSc7QKT-T6skKkoHvu0_n6X83_BEkOg9SHP-qbu6nrfSlns8xzbfV3JbF9UeZ8kXZdWdf6G6nb92DLKEHpzgREZC-5XraClEV-081vnhlk_U4DzSJ5-Bc7gKHAfvxCsEgK98yCTGwxAL-QvIK8tHriNgxPyE5xHbQg4CbTQEvxJ3jGE1t0Xm_XB9ZAm_GOkGdg1xp2M3IWLzh3gvsSDXWbj4rXQzHY-xI3h0W2s50Q-KMZVOje72btuUdSxXuBPxVBumrnvgV5JcTe9Ka4OMLkQWRkMq6ftwDIU93glhyQpGLK_841tr5kFQM97V8y2imy6A_mhWX9183mfZOzoarEK6RcbyVOI2g7mAmTdMozs4y15HSBPXiEYdybPvOZDLLz_NkK0jV7boNXWWr_g883Z_0bqLEOsZZrsa9X0-xxz2jd1S_s27ZpE9pgXKhGy2XXHrGuyBnd0TCuZV2WRJtVuPJZ5UzZ1L1vV59hh2XeqUyUWZVaVbV03O32UicyTWlZpJoskP6Qqaztqi6bCloo0F3lCE2pzYJk-OD_sdAgLHdOklDLdGWzJhHVyS0l2RKvWoSOkFPKTkLI3aAeRPVxVgb8XTzt_XFW_XYYg8sToEMPPAFFHs_4Z2GyLJ1E83syLp1_r4EYVHbZ5vs3w9OcID26it_kdNt6zRnU0k-3Ixt3izfFfzqrZux-kopCnFYpVyTc0Xo7yrwAAAP__uiLulQ">