randlc, an Octave code which implements a version of the pseudorandom number generator used by the NAS Parallel Benchmarks.
The generator has the form
X(K+1) = A * X(K) mod 2^46
where the suggested valule of the multiplier A is 5^13 = 1220703125.
This scheme generates 2^44 numbers before repeating.
The web site for the NAS Parallel Benchmarks is https://www.nas.nasa.gov/Resources/Software/npb.html.
The information on this web page is distributed under the MIT license.
randlc is available in a C version and a C++ version and a Fortran77 version and a Fortran90 version and a MATLAB version and an Octave version and a Python version.