/* -*-c-*-
*
- * $Id: lcrand.c,v 1.1 1999/12/10 23:15:27 mdw Exp $
+ * $Id: lcrand.c,v 1.4 2000/12/06 20:31:06 mdw Exp $
*
* Simple linear congruential generator
*
/*----- Revision history --------------------------------------------------*
*
* $Log: lcrand.c,v $
+ * Revision 1.4 2000/12/06 20:31:06 mdw
+ * Simplify uniform range transformation.
+ *
+ * Revision 1.3 2000/06/17 11:29:03 mdw
+ * Add the flags word to the generic generator.
+ *
+ * Revision 1.2 1999/12/13 15:34:01 mdw
+ * Add support for seeding from a generic pseudorandom source.
+ *
* Revision 1.1 1999/12/10 23:15:27 mdw
* Noncryptographic random number generator.
*
uint32 r = P - P % m;
do xx = lcrand(xx); while (xx >= r);
*x = xx;
- return (xx / (r / m));
+ return (xx % m);
}
/*----- Generic interface -------------------------------------------------*/
case GRAND_CHECK:
case GRAND_SEEDINT:
case GRAND_SEEDUINT32:
+ case GRAND_SEEDRAND:
rc = 1;
break;
default:
case GRAND_SEEDUINT32:
g->x = va_arg(ap, uint32);
break;
+ case GRAND_SEEDRAND: {
+ grand *rr = va_arg(ap, grand *);
+ uint32 x;
+ do x = rr->ops->word(rr); while (x >= P || x == LCRAND_FIXEDPT);
+ g->x = x;
+ } break;
default:
GRAND_BADOP;
break;
static const grand_ops gops = {
"lcrand",
- LCRAND_P,
+ LCRAND_P, 0,
gmisc, gdestroy,
graw, gbyte, grand_word, grange, grand_fill
};