Imported from 5.0/lrandom.tar.gz.

Makefile

@@ -0,0 +1,59 @@
+# makefile for random library for Lua
+
+# change these to reflect your Lua installation
+LUA= /tmp/lhf/lua-5.0
+LUAINC= $(LUA)/include
+LUALIB= $(LUA)/lib
+LUABIN= $(LUA)/bin
+
+# no need to change anything below here
+CFLAGS= $(INCS) $(WARN) -O2 $G
+WARN= -ansi -pedantic -Wall
+INCS= -I$(LUAINC)
+
+MYNAME= random
+MYLIB= l$(MYNAME)
+T= $(MYLIB).so
+OBJS= $(MYLIB).o
+TEST= test.lua
+
+all:	test
+
+test:	$T
+	$(LUABIN)/lua -l$(MYNAME) $(TEST)
+
+o:	$(MYLIB).o
+
+so:	$T
+
+$T:	$(OBJS)
+	$(CC) -o $@ -shared $(OBJS)
+
+$(OBJS): random.c
+
+clean:
+	rm -f $(OBJS) $T core core.* a.out
+
+doc:
+	@echo "$(MYNAME) library:"
+	@fgrep '/**' $(MYLIB).c | cut -f2 -d/ | tr -d '*' | sort | column
+
+# distribution
+
+FTP= $(HOME)/public/ftp/lua/5.0
+D= $(MYNAME)
+A= $(MYLIB).tar.gz
+TOTAR= Makefile,README,$(MYLIB).c,$(MYNAME).lua,test.lua,random.c
+
+tar:	clean
+	tar zcvf $A -C .. $D/{$(TOTAR)}
+
+distr:	tar
+	touch -r $A .stamp
+	mv $A $(FTP)
+
+diff:	clean
+	tar zxf $(FTP)/$A
+	diff $D .
+
+# eof

README

@@ -0,0 +1,25 @@
+This is a random-number library for Lua 5.0. It is based on the Mersenne
+Twister random number generator available at
+	http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/MT2002/emt19937ar.html
+
+To try the library, just edit Makefile to reflect your installation of Lua and
+then run make. This will build the library and run a simple test. For detailed
+installation instructions, see
+	http://www.tecgraf.puc-rio.br/~lhf/ftp/lua/install.html
+
+There is no manual but the library is simple and intuitive; see the summary
+below. Read also test.lua, which shows the library in action.
+
+This code is hereby placed in the public domain, except random.c, which is
+subject to the BSD license.
+
+Please send comments, suggestions, and bug reports to [email protected] .
+
+-------------------------------------------------------------------------------
+
+random library:
+ clone(c)                tostring(c)             valuex(c) 
+ new([seed])             value(c)                version 
+ seed(c,[seed])          valuei(c,a,[b]) 
+
+-------------------------------------------------------------------------------

lrandom.c

@@ -0,0 +1,125 @@
+/*
+* lrandom.c
+* random interface for Lua 5.0
+* Luiz Henrique de Figueiredo <[email protected]>
+* 09 Jun 2006 21:09:17
+* This code is hereby placed in the public domain.
+*/
+
+#include <math.h>
+#include <stdio.h>
+
+#include "lua.h"
+#include "lauxlib.h"
+
+#define MYNAME		"random"
+#define MYVERSION	MYNAME " library for " LUA_VERSION " / Jun 2006"
+#define MYTYPE		MYNAME " handle"
+
+#define SEED		2004UL
+#include "random.c"
+
+static MT *Pget(lua_State *L, int i)
+{
+ if (luaL_checkudata(L,i,MYTYPE)==NULL) luaL_typerror(L,i,MYTYPE);
+ return lua_touserdata(L,i);
+}
+
+static MT *Pnew(lua_State *L)
+{
+ MT *c=lua_newuserdata(L,sizeof(MT));
+ luaL_getmetatable(L,MYTYPE);
+ lua_setmetatable(L,-2);
+ return c;
+}
+
+static int Lnew(lua_State *L)			/** new([seed]) */
+{
+ MT *c=Pnew(L);
+ init_genrand(c,luaL_optlong(L,1,SEED));
+ return 1;
+}
+
+static int Lclone(lua_State *L)			/** clone(c) */
+{
+ MT *c=Pget(L,1);
+ MT *d=Pnew(L);
+ *d=*c;
+ return 1;
+}
+
+static int Lseed(lua_State *L)			/** seed(c,[seed]) */
+{
+ MT *c=Pget(L,1);
+ init_genrand(c,luaL_optlong(L,2,SEED));
+ return 0;
+}
+
+static int Lvalue(lua_State *L)			/** value(c) */
+{
+ MT *c=Pget(L,1);
+ lua_pushnumber(L,genrand_real1(c));
+ return 1;
+}
+
+static int Lvaluei(lua_State *L)		/** valuei(c,a,[b]) */
+{
+ MT *c=Pget(L,1);
+ int a,b;
+ if (lua_gettop(L)==2)
+ {
+  a=1;
+  b=luaL_checkint(L,2);
+ }
+ else
+ {
+  a=luaL_checkint(L,2);
+  b=luaL_checkint(L,3);
+ }
+ lua_pushnumber(L,floor(a+genrand_real2(c)*(b-a+1)));
+ return 1;
+}
+
+static int Lvaluex(lua_State *L)		/** valuex(c) */
+{
+ MT *c=Pget(L,1);
+ lua_pushnumber(L,genrand_res53(c));
+ return 1;
+}
+
+static int Ltostring(lua_State *L)		/** tostring(c) */
+{
+ MT *c=Pget(L,1);
+ char s[64];
+ sprintf(s,"%s %p",MYTYPE,(void*)c);
+ lua_pushstring(L,s);
+ return 1;
+}
+
+static const luaL_reg R[] =
+{
+	{ "__tostring",	Ltostring	},
+	{ "clone",	Lclone		},
+	{ "new",	Lnew		},
+	{ "seed",	Lseed		},
+	{ "tostring",	Ltostring	},
+	{ "value",	Lvalue		},
+	{ "valuei",	Lvaluei		},
+	{ "valuex",	Lvaluex		},
+	{ NULL,		NULL		}
+};
+
+LUALIB_API int luaopen_random(lua_State *L)
+{
+ lua_pushliteral(L,MYNAME);
+ luaL_newmetatable(L,MYTYPE);
+ luaL_openlib(L,NULL,R,0);
+ lua_pushliteral(L,"version");			/** version */
+ lua_pushliteral(L,MYVERSION);
+ lua_settable(L,-3);
+ lua_pushliteral(L,"__index");
+ lua_pushvalue(L,-2);
+ lua_settable(L,-3);
+ lua_rawset(L,LUA_GLOBALSINDEX);
+ return 1;
+}

random.c

@@ -0,0 +1,142 @@
+/*
+* random.c
+* Mersenne Twister random number generator
+* Luiz Henrique de Figueiredo <[email protected]>
+* 26 Jun 2004 10:14:59
+* slightly modified from mt19937ar.c available at
+* http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/MT2002/emt19937ar.html
+*/
+
+/* 
+   A C-program for MT19937, with initialization improved 2002/1/26.
+   Coded by Takuji Nishimura and Makoto Matsumoto.
+
+   Before using, initialize the state by using init_genrand(seed)  
+   or init_by_array(init_key, key_length).
+
+   Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura,
+   All rights reserved.                          
+
+   Redistribution and use in source and binary forms, with or without
+   modification, are permitted provided that the following conditions
+   are met:
+
+     1. Redistributions of source code must retain the above copyright
+        notice, this list of conditions and the following disclaimer.
+
+     2. Redistributions in binary form must reproduce the above copyright
+        notice, this list of conditions and the following disclaimer in the
+        documentation and/or other materials provided with the distribution.
+
+     3. The names of its contributors may not be used to endorse or promote 
+        products derived from this software without specific prior written 
+        permission.
+
+   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
+   "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
+   LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
+   A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR
+   CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
+   EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
+   PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
+   PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
+   LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
+   NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
+   SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+
+
+   Any feedback is very welcome.
+   http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html
+   email: m-mat @ math.sci.hiroshima-u.ac.jp (remove space)
+*/
+
+/* Period parameters */  
+#define N 624
+#define M 397
+#define MATRIX_A 0x9908b0dfUL   /* constant vector a */
+#define UPPER_MASK 0x80000000UL /* most significant w-r bits */
+#define LOWER_MASK 0x7fffffffUL /* least significant r bits */
+
+typedef struct {
+    unsigned long v[N]; /* the array for the state vector  */
+    int i;
+} MT;
+
+#define mt	(o->v)
+#define mti	(o->i)
+
+/* initializes mt[N] with a seed */
+static void init_genrand(MT *o, unsigned long s)
+{
+    mt[0]= s & 0xffffffffUL;
+    for (mti=1; mti<N; mti++) {
+        mt[mti] = 
+	    (1812433253UL * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti); 
+        /* See Knuth TAOCP Vol2. 3rd Ed. P.106 for multiplier. */
+        /* In the previous versions, MSBs of the seed affect   */
+        /* only MSBs of the array mt[].                        */
+        /* 2002/01/09 modified by Makoto Matsumoto             */
+        mt[mti] &= 0xffffffffUL;
+        /* for >32 bit machines */
+    }
+}
+
+/* generates a random number on [0,0xffffffff]-interval */
+static unsigned long genrand_int32(MT *o)
+{
+    unsigned long y;
+    static unsigned long mag01[2]={0x0UL, MATRIX_A};
+    /* mag01[x] = x * MATRIX_A  for x=0,1 */
+
+    if (mti >= N) { /* generate N words at one time */
+        int kk;
+
+        if (mti == N+1)   /* if init_genrand() has not been called, */
+            init_genrand(o,5489UL); /* a default initial seed is used */
+
+        for (kk=0;kk<N-M;kk++) {
+            y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
+            mt[kk] = mt[kk+M] ^ (y >> 1) ^ mag01[y & 0x1UL];
+        }
+        for (;kk<N-1;kk++) {
+            y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK);
+            mt[kk] = mt[kk+(M-N)] ^ (y >> 1) ^ mag01[y & 0x1UL];
+        }
+        y = (mt[N-1]&UPPER_MASK)|(mt[0]&LOWER_MASK);
+        mt[N-1] = mt[M-1] ^ (y >> 1) ^ mag01[y & 0x1UL];
+
+        mti = 0;
+    }
+
+    y = mt[mti++];
+
+    /* Tempering */
+    y ^= (y >> 11);
+    y ^= (y << 7) & 0x9d2c5680UL;
+    y ^= (y << 15) & 0xefc60000UL;
+    y ^= (y >> 18);
+
+    return y;
+}
+
+/* generates a random number on [0,1]-real-interval */
+static double genrand_real1(MT *o)
+{
+    return genrand_int32(o)*(1.0/4294967295.0); 
+    /* divided by 2^32-1 */ 
+}
+
+/* generates a random number on [0,1)-real-interval */
+static double genrand_real2(MT *o)
+{
+    return genrand_int32(o)*(1.0/4294967296.0); 
+    /* divided by 2^32 */
+}
+
+/* generates a random number on [0,1) with 53-bit resolution*/
+static double genrand_res53(MT *o) 
+{ 
+    unsigned long a=genrand_int32(o)>>5, b=genrand_int32(o)>>6; 
+    return(a*67108864.0+b)*(1.0/9007199254740992.0); 
+} 
+/* These real versions are due to Isaku Wada, 2002/01/09 added */

random.lua

@@ -0,0 +1,10 @@
+-- random.lua
+-- support code for random library
+-- usage lua -lrandom ...
+
+local function so(x)
+	local SOPATH= os.getenv"LUA_SOPATH" or "./"
+	assert(loadlib(SOPATH.."l"..x..".so","luaopen_"..x))()
+end
+
+so"random"