Re: KISS4691, a potentially top-ranked RNG.

["robin" <robin51@dodo.com.au>]

"Uno" <merrilljensen@q.com> wrote in message news:8bfos5FakcU1@mid.individual.net...
| robin wrote:

| > I have already posted a PL/I version using unsigned arithmetic.
| >
| > Here is another version, this time using signed arithmetic :--
| >
| > (NOSIZE, NOFOFL):
| > RNG: PROCEDURE OPTIONS (MAIN, REORDER);
| >
| >    declare (xs initial (521288629), xcng initial (362436069),
| >             Q(0:4690) ) static fixed binary (31);
| >
| > MWC: procedure () returns (fixed binary (31));
| >    declare (t,x,i) fixed binary (31);
| >    declare (c initial (0), j initial (4691) ) fixed binary (31) static;
| >    declare (t1, t2, t3) fixed binary (31);
| >
| >    if j < hbound(Q,1) then j = j + 1; else j = 0;
| >    x = Q(j);
| >    t = isll(x,13)+c+x;
| >    t1 = iand(x, 3) - iand(t, 3);
| >    t2 = isrl(x, 2) - isrl(t, 2);
| >    if t2 = 0 then t2 = t1;
| >    if t2 > 0 then t3 = 1; else t3 = 0;
| >    c = t3 + isrl(x, 19);
| >    Q(j)=t;
| >    return (t);
| > end MWC;
| >
| > CNG: procedure returns (fixed binary (31));
| >   xcng=bin(69069)*xcng+bin(123);
| >   return (xcng);
| > end CNG;
| >
| > XXS: procedure returns (fixed binary (31));
| >    xs = ieor (xs, isll(xs, 13) );
| >    xs = ieor (xs, isrl(xs, 17) );
| >    xs = ieor (xs, isll(xs,  5) );
| >    return (xs);
| > end XXS;
| >
| > KISS: procedure returns (fixed binary (31));
| >    return ( MWC()+CNG+XXS );
| > end KISS;
| >
| >    declare (i,x) fixed binary (31);
| >    declare y fixed decimal (11);
| >
| >    Q = CNG+XXS; /* Initialize. */
| >    do i = 1 to 1000000000; x=MWC(); end;
| >    put skip edit (" Expected MWC result = 3740121002", 'computed =', x)
| >       (a, skip, x(12), a, f(11));
| >    y = iand(x, 2147483647);
| >    if x < 0 then y = y + 2147483648;
| >    put skip edit (y) (x(11), f(22)); put skip;
| >    do i = 1 to 1000000000; x=KISS; end;
| >    put skip edit ("Expected KISS result = 2224631993", 'computed =', x)
| >       (a, skip, x(12), a, f(11));
| >    y = iand(x, 2147483647);
| >    if x < 0 then y = y + 2147483648;
| >    put skip edit (y) (x(11), f(22));
| >
| > end RNG;

| If you were to comment out the PL/I command line that compiled this,
| what would it be?

???