Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.lang.forth > #132332

Re: KISS 64-bit pseudo-random number generator

Show all headers | View raw

On 9/19/24 06:45, Krishna Myneni wrote:
> On 9/19/24 03:57, albert@spenarnc.xs4all.nl wrote:
>> In article <vcgok8$gol7$1@dont-email.me>,
>> Krishna Myneni  <krishna.myneni@ccreweb.org> wrote:
>> <SNIP>
>>> Moments of speed
>>>   N       <v> (m/s)    <v^2> (m/s)^2    <v^3> (m/s)^3
>>> 10^2     1181.0956     1656472.7       2604709063.
>>> 10^3     1293.3130     1952149.7       3300955817.
>>> 10^4     1259.3279     1862988.3       3108515117.
>>> 10^5     1260.5577     1872157.8       3147664636.
>>> 10^6     1259.4425     1868918.9       3139487337.
>>> 10^7     1259.6136     1869145.0       3139092438.
>>
>> I think for a Monte Carlo simulation at least three tests
>> must be done with different seeds.
> 
> Good point. For a meaningful comparison of errors between PRNGs at a 
> specific N, the statistical variation of the <v^n> need to be measured 
> for different seed values.
> 
> I can add some code to measure this sigma at each N, with 32 seeds 
> uniformly spaced between 0 and UMAX.
> 

I've calculated the statistical variation in the moments for each set of 
N, using 16 different seeds (spaced apart over the interval for UMAX). 
The standard dev. for the 16 <v^i>, computed for N trials is comparable 
to the relative error between the moment and its theoretical value. 
Thus, the relative errors are indeed a meaningful comparison between the 
two prngs tested here, and I think this implies that for N > 10^5 the 
LCG PRNG  (RANDOM) gives more accurate answers than the KISS 64 bit PRNG 
(RAN-KISS), for this problem. The LCG PRNG is faster than the KISS 
64-bit PRNG.

minforth stated earlier that he would prefer to use diehard tests to 
decide between which of these two PRNGs to use for computing these 
results from random trials. It will be interesting to see if diehard 
tests are consistent with what I find from actually using the PRNGs and 
comparing the results to the expected results (for large N and ideal PRNG).

--
Krishna

Back to comp.lang.forth | Previous | Next — Previous in thread | Next in thread | Find similar

Thread

KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-08 22:09 -0500
  Re: KISS 64-bit pseudo-random number generator Lars Brinkhoff <lars.spam@nocrew.org> - 2024-09-09 06:55 +0000
    Re: KISS 64-bit pseudo-random number generator mhx@iae.nl (mhx) - 2024-09-09 07:41 +0000
      Re: KISS 64-bit pseudo-random number generator anton@mips.complang.tuwien.ac.at (Anton Ertl) - 2024-09-09 08:55 +0000
        Re: KISS 64-bit pseudo-random number generator mhx@iae.nl (mhx) - 2024-09-09 10:04 +0000
        Re: KISS 64-bit pseudo-random number generator albert@spenarnc.xs4all.nl - 2024-09-09 12:28 +0200
          Re: KISS 64-bit pseudo-random number generator anton@mips.complang.tuwien.ac.at (Anton Ertl) - 2024-09-09 15:26 +0000
            Re: KISS 64-bit pseudo-random number generator albert@spenarnc.xs4all.nl - 2024-09-10 12:00 +0200
        Re: KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-10 19:30 -0500
          Re: KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-12 21:10 -0500
            Re: KISS 64-bit pseudo-random number generator Paul Rubin <no.email@nospam.invalid> - 2024-09-12 19:15 -0700
              Re: KISS 64-bit pseudo-random number generator mhx@iae.nl (mhx) - 2024-09-13 06:07 +0000
                Re: KISS 64-bit pseudo-random number generator Paul Rubin <no.email@nospam.invalid> - 2024-09-13 03:46 -0700
              Re: KISS 64-bit pseudo-random number generator minforth@gmx.net (minforth) - 2024-09-13 06:56 +0000
            Re: KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-18 19:28 -0500
              Re: KISS 64-bit pseudo-random number generator mhx@iae.nl (mhx) - 2024-09-19 06:33 +0000
                Re: KISS 64-bit pseudo-random number generator minforth@gmx.net (minforth) - 2024-09-19 08:18 +0000
                Re: KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-19 03:50 -0500
                Re: KISS 64-bit pseudo-random number generator albert@spenarnc.xs4all.nl - 2024-09-19 10:57 +0200
                Re: KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-19 06:45 -0500
                Re: KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-25 18:15 -0500
                Re: KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-25 20:45 -0500
    Re: KISS 64-bit pseudo-random number generator Krishna Myneni <krishna.myneni@ccreweb.org> - 2024-09-09 18:36 -0500

csiph-web