test.it.unimi.dsi.util.XoShiRo256PlusRandomGeneratorTest Maven / Gradle / Ivy

Go to download
Show more of this group Show more artifacts with this name
Show all versions of dsiutils Show documentation
The DSI utilities are a mishmash of classes accumulated during the last twenty years in projects developed at the DSI (Dipartimento di Scienze dell'Informazione, i.e., Information Sciences Department), now DI (Dipartimento di Informatica, i.e., Informatics Department), of the Universita` degli Studi di Milano.
There is a newer version: 2.7.3
Show newest version
/*
 * DSI utilities
 *
 * Copyright (C) 2014-2020 Sebastiano Vigna
 *
 *  This library is free software; you can redistribute it and/or modify it
 *  under the terms of the GNU Lesser General Public License as published by the Free
 *  Software Foundation; either version 3 of the License, or (at your option)
 *  any later version.
 *
 *  This library is distributed in the hope that it will be useful, but
 *  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
 *  or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public License
 *  for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with this program; if not, see .
 *
 */

package it.unimi.dsi.util;

import static org.junit.Assert.assertEquals;
import static org.junit.Assert.assertTrue;

import org.junit.Test;


public class XoShiRo256PlusRandomGeneratorTest {
	private final static long seeds[] = { 0, 1, 1024, 0x5555555555555555L };

	@Test
	public void testNextFloat() {
		for (final long seed : seeds) {
			final XoShiRo256PlusRandomGenerator xoroshiro = new XoShiRo256PlusRandomGenerator(seed);
			double avg = 0;
			for (int i = 1000000; i-- != 0;) {
				final float d = xoroshiro.nextFloat();
				assertTrue(Float.toString(d), d < 1);
				assertTrue(Float.toString(d), d >= 0);
				avg += d;
			}

			assertEquals(500000, avg, 1000);
		}
	}

	@Test
	public void testNextDouble() {
		for (final long seed : seeds) {
			final XoShiRo256PlusRandomGenerator xoroshiro = new XoShiRo256PlusRandomGenerator(seed);
			double avg = 0;
			for (int i = 1000000; i-- != 0;) {
				final double d = xoroshiro.nextDouble();
				assertTrue(d < 1);
				assertTrue(d >= 0);
				avg += d;
			}

			assertEquals(500000, avg, 1000);
		}
	}

	@Test
	public void testNextDoubleFast() {
		for (final long seed : seeds) {
			final XoShiRo256PlusRandomGenerator xoroshiro = new XoShiRo256PlusRandomGenerator(seed);
			double avg = 0;
			for (int i = 1000000; i-- != 0;) {
				final double d = xoroshiro.nextDoubleFast();
				assertTrue(d < 1);
				assertTrue(d >= 0);
				avg += d;
			}

			assertEquals(500000, avg, 1000);
		}
	}

	@Test
	public void testNextInt() {
		for (final long seed : seeds) {
			final XoShiRo256PlusRandomGenerator xoroshiro = new XoShiRo256PlusRandomGenerator(seed);
			double avg = 0;
			for (int i = 100000000; i-- != 0;) {
				final int d = xoroshiro.nextInt(101);
				assertTrue(d <= 100);
				assertTrue(d >= 0);
				avg += d;
			}

			assertEquals(5000000000L, avg, 1000000);
		}
	}

	@Test
	public void testNextInt2() {
		for (final long seed : seeds) {
			final XoShiRo256PlusRandomGenerator xoroshiro = new XoShiRo256PlusRandomGenerator(seed);
			final int[] count = new int[32];
			long change = 0;
			int prev = 0;
			for (int i = 1000000; i-- != 0;) {
				final int d = xoroshiro.nextInt();
				change += Long.bitCount(d ^ prev);
				for (int b = 32; b-- != 0;)
					if ((d & (1 << b)) != 0) count[b]++;
				prev = d;
			}

			assertEquals(32 * 1000000L, change, 33000);
			for (int b = 32; b-- != 0;) assertEquals(500000, count[b], 2000);
		}
	}

	@Test
	public void testNextInt3() {
		final XoShiRo256PlusRandomGenerator xoroshiro = new XoShiRo256PlusRandomGenerator(0);
		for(int i = 0; i < 100; i++) assertTrue(xoroshiro.nextInt(16) < 16);
	}

	@Test
	public void testNextInt4() {
		final XoShiRo256PlusRandomGenerator xoroshiro = new XoShiRo256PlusRandomGenerator(0);
		assertEquals(0, xoroshiro.nextInt(1));
	}

	@Test
	public void testNextLong() {
		for (final long seed : seeds) {
			final XoShiRo256PlusRandomGenerator xoroshiro = new XoShiRo256PlusRandomGenerator(seed);
			final int[] count = new int[64];
			long change = 0;
			long prev = 0;
			for (int i = 1000000; i-- != 0;) {
				final long d = xoroshiro.nextLong();
				change += Long.bitCount(d ^ prev);
				for (int b = 64; b-- != 0;)
					if ((d & (1L << b)) != 0) count[b]++;
				prev = d;
			}

			assertEquals(32 * 1000000L, change, 7000);
			for (int b = 64; b-- != 0;) assertEquals(500000, count[b], 2000);
		}
	}

	@Test
	public void testSameAsRandom() {
		final XoShiRo256PlusRandom xoshiroPlusRandom = new XoShiRo256PlusRandom(0);
		final XoShiRo256PlusRandomGenerator xoshiroPlus = new XoShiRo256PlusRandomGenerator(0);
		for(int i = 1000; i-- != 0;) {
			assertEquals(xoshiroPlus.nextLong(), xoshiroPlusRandom.nextLong());
			assertEquals(0, xoshiroPlus.nextDouble(), xoshiroPlusRandom.nextDouble());
			assertEquals(xoshiroPlus.nextInt(), xoshiroPlusRandom.nextInt());
			assertEquals(xoshiroPlus.nextInt(99), xoshiroPlusRandom.nextInt(99));
			assertEquals(xoshiroPlus.nextInt(128), xoshiroPlusRandom.nextInt(128));
		}
		xoshiroPlus.jump();
		xoshiroPlusRandom.jump();
		for(int i = 1000; i-- != 0;) {
			assertEquals(xoshiroPlus.nextLong(), xoshiroPlusRandom.nextLong());
			assertEquals(0, xoshiroPlus.nextDouble(), xoshiroPlusRandom.nextDouble());
			assertEquals(xoshiroPlus.nextInt(), xoshiroPlusRandom.nextInt());
			assertEquals(xoshiroPlus.nextInt(99), xoshiroPlusRandom.nextInt(99));
			assertEquals(xoshiroPlus.nextInt(128), xoshiroPlusRandom.nextInt(128));
		}
		xoshiroPlus.longJump();
		xoshiroPlusRandom.longJump();
		for(int i = 1000; i-- != 0;) {
			assertEquals(xoshiroPlus.nextLong(), xoshiroPlusRandom.nextLong());
			assertEquals(0, xoshiroPlus.nextDouble(), xoshiroPlusRandom.nextDouble());
			assertEquals(xoshiroPlus.nextInt(), xoshiroPlusRandom.nextInt());
			assertEquals(xoshiroPlus.nextInt(99), xoshiroPlusRandom.nextInt(99));
			assertEquals(xoshiroPlus.nextInt(128), xoshiroPlusRandom.nextInt(128));
		}
	}
}