// park-miller.cc   Park-Miller random number generator

// W. Langdon cs.ucl.ac.uk 5 May 1994

#include <math.h>
#include "park-miller.h"

int intrnd (int& seed) // 1<=seed<=m
{
double const a    = 16807;      //ie 7**5
double const m    = 2147483647; //ie 2**31-1

	double temp = seed * a;
	seed = (int) (temp - m * floor ( temp / m ));
	return seed;
}