KISS07

Javascript

  function KISS07() {
    return (function(args) {
      // George Marsaglia, 2007-06-23
      //http://groups.google.com/group/comp.lang.fortran/msg/6edb8ad6ec5421a5
      var x = 123456789;
      var y = 362436069;
      var z =  21288629;
      var w =  14921776;
      var c = 0;
 
      if (args.length == 0) {
        args = [+new Date];
      }
      var mash = Mash();
      for (var i = 0; i < args.length; i++) {
        x ^= mash(args[i]) * 0x100000000; // 2^32
        y ^= mash(args[i]) * 0x100000000;
        z ^= mash(args[i]) * 0x100000000;
        w ^= mash(args[i]) * 0x100000000;
      }
      if (y === 0) {
        y = 1;
      }
      c ^= z >>> 31;
      z &= 0x7fffffff;
      if ((z % 7559) === 0) {
        z++;
      }
      w &= 0x7fffffff;
      if ((w % 7559) === 0) {
        w++;
      }
      mash = null;
 
      var uint32 = function() {
        var t;
 
        x += 545925293;
        x >>>= 0;
 
        y ^= y << 13;
        y ^= y >>> 17;
        y ^= y << 5;
 
        t = z + w + c;
        z = w;
        c = t >>> 31;
        w = t & 0x7fffffff;
 
        return x + y + w >>> 0;
      };
 
      var random = function() {
        return uint32() * 2.3283064365386963e-10; // 2^-32
      };
      random.uint32 = uint32;
      random.fract53 = function() {
        return random() +
          (uint32() & 0x1fffff) * 1.1102230246251565e-16; // 2^-53
      };
      random.args = args;
      random.version = 'KISS07 0.9';
 
      return random;
    } (Array.prototype.slice.call(arguments)));
  };

C

#include <stdio.h>
#include <stdint.h>
#include "mash.h"
 
// George Marsaglia, 2007-06-23
//http://groups.google.com/group/comp.lang.fortran/msg/6edb8ad6ec5421a5
 
const double norm32 = 2.3283064365386963e-10;   // 2^-32
 
static uint32_t x = 123456789;
static uint32_t y = 362436069;
static uint32_t z = 21288629;
static uint32_t w = 14921776;
static uint32_t c = 0;
 
static inline double kiss07(void) {
  uint32_t t;
 
  x += 545925293;
 
  y ^= y << 13;
  y ^= y >> 17;
  y ^= y << 5;
 
  t = z + w + c;
  z = w;
  c = t >> 31;
  w = t & 0x7fffffff;
 
  return (x + y + w) * norm32;  // 2^-32;
}
 
int main(int argc, unsigned char *argv[]) {
  double buffer[256];
  int i;
 
  for (i = 1; i < argc; i++) {
    x ^= (uint32_t) (mash(argv[i]) * 0x100000000ULL);   // 2^32
    y ^= (uint32_t) (mash(argv[i]) * 0x100000000ULL);   // 2^32
    z ^= (uint32_t) (mash(argv[i]) * 0x100000000ULL);   // 2^32
    w ^= (uint32_t) (mash(argv[i]) * 0x100000000ULL);   // 2^32
  }
  if (y == 0) {
    y = 1;
  }
  c ^= z >> 31;
  z &= 0x7fffffff;
  if ((z % 7559) == 0) {
    z++;
  }
  w &= 0x7fffffff;
  if ((w % 7559) == 0) {
    w++;
  }
 
  while (1) {
    for (i = 0; i < 256; i++) {
      buffer[i] = kiss07();
    }
    fwrite(buffer, sizeof buffer, 1, stdout);
  }
}