> Thus I suggest using this listing for the MWC function
> in KISS4691:
> ------------------------------------------------

>   unsigned long MWC(void)
>   {unsigned long t,x; static c=0,j=4691;

> and recommend that it be the pattern for implementations in
> PL/I, asm, Fortran or others.  For signed integers, flip the
> sign bits so that t<x becomes (t^m)<(x^m) or equivalents.
> (C versions using signed integers should also avoid the
> problem of right shifts, replacing (x>>19) by
> ((unsigned)x>>19) or ((unsigned)xs>>17) in the XS component.)

> In case you do not have, or don't want to fetch, the original,
> here is the entire listing, with the recommended form for MWC():

> ------------------------------------------------------------
> static unsigned long xs=521288629,xcng=362436069,Q[4691];

> unsigned long MWC(void)
> {unsigned long t,x; static c=0,j=4691;
>    j=(j<4690)? j+1:0;
>    x=Q[j]; t=(x<<13)+c;
>      if(t<c) {c=(x>>19)+1; t=t+x;}
>      else {t=t+x; c=(x>>19)+(t<x);}
>    return (Q[j]=t);
> }

> #define CNG ( xcng=69069*xcng+123 )
> #define XS ( xs^=(xs<<13), xs^=(xs>>17), xs^=(xs<<5) )
> #define KISS ( MWC()+CNG+XS )

> #include<stdio.h>
> int main()
> {unsigned long i,x;
>   for(i=0;i<4691;i++)  Q[i]=CNG+XS;
> printf("Does  MWC result=3740121002 ?\n");
>   for(i=0;i<1000000000;i++) x=MWC();
>     printf("%27u\n",x);
>     printf("Does KISS result=2224631993 ?\n");
>   for(i=0;i<1000000000;i++) x=KISS;
>     printf("%27u\n",x);
> }