added c program to generate lists of random chunks from GSL random number generators, these can serve as comparison to test correctness of implementation

Author: dcatteeu

tests/gsl-rng-dump.c

@@ -0,0 +1,68 @@
+#include <stdio.h>
+#include <string.h>
+#include <gsl/gsl_rng.h>
+
+/* Return the gsl_rng_type for the GSL random number generator with
+   the provided name. */
+const gsl_rng_type* type_for_name (const char* name) 
+{
+    const gsl_rng_type **t, **t0;
+    t0 = gsl_rng_types_setup ();
+    for (t = t0; *t != 0; t++){
+	if (strcmp((*t)->name, name) == 0) {
+	    return *t;
+	}
+    }
+    return 0;
+}
+
+/* Given the name of a GSL random number generator and (optionally) a
+   seed, print the name, the seed, and the first n random chunks to
+   stdout. */
+int main (int argc, char** argv)
+{
+    const char * name;
+    unsigned long seed;
+    const gsl_rng_type * T;
+    gsl_rng * r;
+    
+    int i, n = 10;
+    
+    if (argc < 2) {
+	fprintf (stderr, 
+		 "Expects one or two parameters. The first, a string, is the name of a random\n"
+		 "number generator and is mandatory. The second, an unsigned long, is a seed\n"
+		 "and is optional, if it is not supplied, the default seed is used.\n");
+	exit (1);
+    } else if (argc < 3) {
+	seed = 0; // In GSL, 0 is replaced by the rng's default seed.
+    } else {
+	seed = atol(argv[2]);
+    }
+    name = argv[1];
+    
+    gsl_rng_env_setup();
+    
+    T = type_for_name (name);
+
+    if (T == 0) {
+	fprintf (stderr, "\"%s\" is not the name of a gsl_rng_type.\n", name);
+	exit (2);
+    }
+
+    r = gsl_rng_alloc (T);
+    gsl_rng_set(r, seed);
+    
+    /* Print name of rng, the seed used to initialized it, and its
+       first n random chunks. Items are separated by a <space> and
+       terminated with a <newline>. */
+    printf("\"%s\" %lu", name, seed);
+    for (i = 0; i < n; i++) {
+	unsigned long u = gsl_rng_get (r);
+	printf (" %lu", u);
+    }
+    printf("\n");
+    
+    gsl_rng_free (r);
+    return 0;
+}

tests/gsl-rng-dump.csv

@@ -0,0 +1,11 @@
+"randu", 0, 65539, 393225, 1769499, 7077969, 26542323, 95552217, 334432395, 1146624417, 1722371299, 14608041
+"minstd", 0, 16807, 282475249, 1622650073, 984943658, 1144108930, 470211272, 101027544, 1457850878, 1458777923, 2007237709
+"mt19937", 0, 4293858116, 699692587, 1213834231, 4068197670, 994957275, 2082945813, 4112332215, 3196767107, 2319469851, 3178073856
+"mrg", 0, 572361259, 521023500, 563045572, 393759085, 1080953451, 130004609, 893178225, 1206078822, 397595998, 1497657786
+"cmrg", 0, 240037626, 2059795007, 1807165044, 1987289342, 591431996, 1523218864, 1189992209, 1891279779, 1282347090, 1614637363
+"transputer", 0, 1664525, 389569705, 2940799637, 158984081, 2862450781, 3211393721, 1851289957, 3934847009, 2184914861, 246739401
+"borosh13", 0, 1812433253, 88293849, 1790253981, 42330609, 3130934549, 318537289, 3711034317, 2922550497, 3385337285, 3700446137
+"fishman18", 0, 62089911, 847344462, 1061653656, 1954074819, 226824280, 953102500, 1452288378, 50913524, 2133871779, 1843965925
+"fishman20", 0, 48271, 182605794, 1291394886, 1914720637, 2078669041, 407355683, 1105902161, 854716505, 564586691, 1596680831
+"lecuyer21", 0, 40692, 1655838864, 2103410263, 1872071452, 652912057, 1780294415, 535353314, 525453832, 1422611300, 1336516156
+"waterman14", 0, 1566083941, 2203506137, 1324822941, 1986974193, 2643373845, 1922267721, 3584877005, 1128752353, 3108566981, 1963571129

tests/make-gsl-rng-dump.sh

@@ -0,0 +1,2 @@
+gcc -Wall -I/usr/local/include -c gsl-rng-dump.c
+gcc -L/usr/local/lib gsl-rng-dump.o -lgsl -lm -o gsl-rng-dump