/// <summary> /// Run example /// </summary> public void Run() { // 1. Generate 20 samples of the function f(x) = x on interval [-5, 5] Console.WriteLine(@"1. Generate 20 samples of the function f(x) = x on interval [-5, 5]"); double[] points; var values = SignalGenerator.EquidistantInterval(TargetFunction, -5, 5, 20, out points); Console.WriteLine(); // 2. Create a burlish stoer rational interpolation based on arbitrary points var method = Interpolate.RationalWithPoles(points, values); Console.WriteLine(@"2. Create a burlish stoer rational interpolation based on arbitrary points"); Console.WriteLine(); // 3. Check if interpolation support integration Console.WriteLine(@"3. Support integration = {0}", method.SupportsIntegration); Console.WriteLine(); // 4. Check if interpolation support differentiation Console.WriteLine(@"4. Support differentiation = {0}", method.SupportsDifferentiation); Console.WriteLine(); // 5. Interpolate ten random points and compare to function results Console.WriteLine(@"5. Interpolate ten random points and compare to function results"); var rng = new MersenneTwister(1); for (var i = 0; i < 10; i++) { // Generate random value from [0, 5] var point = rng.Next(0, 5); Console.WriteLine(@"Interpolate at {0} = {1}. Function({0}) = {2}", point.ToString("N05"), method.Interpolate(point).ToString("N05"), TargetFunction(point).ToString("N05")); } Console.WriteLine(); }
/// <summary> /// Run example /// </summary> /// <seealso cref="http://en.wikipedia.org/wiki/Random_number_generation">Random number generation</seealso> /// <seealso cref="http://en.wikipedia.org/wiki/Linear_congruential_generator">Linear congruential generator</seealso> /// <seealso cref="http://en.wikipedia.org/wiki/Mersenne_twister">Mersenne twister</seealso> /// <seealso cref="http://en.wikipedia.org/wiki/Lagged_Fibonacci_generator">Lagged Fibonacci generator</seealso> /// <seealso cref="http://en.wikipedia.org/wiki/Xorshift">Xorshift</seealso> public void Run() { // All RNG classes in MathNet have next counstructors: // - RNG(int seed, bool threadSafe): initializes a new instance with specific seed value and thread safe property // - RNG(int seed): iуууnitializes a new instance with specific seed value. Thread safe property is set to Control.ThreadSafeRandomNumberGenerators // - RNG(bool threadSafe) : initializes a new instance with the seed value set to DateTime.Now.Ticks and specific thread safe property // - RNG(bool threadSafe) : initializes a new instance with the seed value set to DateTime.Now.Ticks and thread safe property set to Control.ThreadSafeRandomNumberGenerators // All RNG classes in MathNet have next methods to produce random values: // - double[] NextDouble(int n): returns an "n"-size array of uniformly distributed random doubles in the interval [0.0,1.0]; // - int Next(): returns a nonnegative random number; // - int Next(int maxValue): returns a random number less then a specified maximum; // - int Next(int minValue, int maxValue): returns a random number within a specified range; // - void NextBytes(byte[] buffer): fills the elements of a specified array of bytes with random numbers; // All RNG classes in MathNet have next extension methods to produce random values: // - long NextInt64(): returns a nonnegative random number less than "Int64.MaxValue"; // - int NextFullRangeInt32(): returns a random number of the full Int32 range; // - long NextFullRangeInt64(): returns a random number of the full Int64 range; // - decimal NextDecimal(): returns a nonnegative decimal floating point random number less than 1.0; // 1. Multiplicative congruential generator using a modulus of 2^31-1 and a multiplier of 1132489760 var mcg31M1 = new Mcg31m1(1); Console.WriteLine(@"1. Generate 10 random double values using Multiplicative congruential generator with a modulus of 2^31-1 and a multiplier of 1132489760"); var randomValues = mcg31M1.NextDouble(10); for (var i = 0; i < randomValues.Length; i++) { Console.Write(randomValues[i].ToString("N") + @" "); } Console.WriteLine(); Console.WriteLine(); // 2. Multiplicative congruential generator using a modulus of 2^59 and a multiplier of 13^13 var mcg59 = new Mcg59(1); Console.WriteLine(@"2. Generate 10 random integer values using Multiplicative congruential generator with a modulus of 2^59 and a multiplier of 13^13"); for (var i = 0; i < 10; i++) { Console.Write(mcg59.Next() + @" "); } Console.WriteLine(); Console.WriteLine(); // 3. Random number generator using Mersenne Twister 19937 algorithm var mersenneTwister = new MersenneTwister(1); Console.WriteLine(@"3. Generate 10 random integer values less then 100 using Mersenne Twister 19937 algorithm"); for (var i = 0; i < 10; i++) { Console.Write(mersenneTwister.Next(100) + @" "); } Console.WriteLine(); Console.WriteLine(); // 4. Multiple recursive generator with 2 components of order 3 var mrg32K3A = new Mrg32k3a(1); Console.WriteLine(@"4. Generate 10 random integer values in range [50;100] using multiple recursive generator with 2 components of order 3"); for (var i = 0; i < 10; i++) { Console.Write(mrg32K3A.Next(50, 100) + @" "); } Console.WriteLine(); Console.WriteLine(); // 5. Parallel Additive Lagged Fibonacci pseudo-random number generator var palf = new Palf(1); Console.WriteLine(@"5. Generate 10 random bytes using Parallel Additive Lagged Fibonacci pseudo-random number generator"); var bytes = new byte[10]; palf.NextBytes(bytes); for (var i = 0; i < bytes.Length; i++) { Console.Write(bytes[i] + @" "); } Console.WriteLine(); Console.WriteLine(); // 6. A random number generator based on the "System.Security.Cryptography.RandomNumberGenerator" class in the .NET library var systemCryptoRandomNumberGenerator = new SystemCryptoRandomNumberGenerator(); Console.WriteLine(@"6. Generate 10 random decimal values using RNG based on the 'System.Security.Cryptography.RandomNumberGenerator'"); for (var i = 0; i < 10; i++) { Console.Write(systemCryptoRandomNumberGenerator.NextDecimal().ToString("N") + @" "); } Console.WriteLine(); Console.WriteLine(); // 7. Wichmann-Hill’s 1982 combined multiplicative congruential generator var rngWh1982 = new WH1982(); Console.WriteLine(@"7. Generate 10 random full Int32 range values using Wichmann-Hill’s 1982 combined multiplicative congruential generator"); for (var i = 0; i < 10; i++) { Console.Write(rngWh1982.NextFullRangeInt32() + @" "); } Console.WriteLine(); Console.WriteLine(); // 8. Wichmann-Hill’s 2006 combined multiplicative congruential generator. var rngWh2006 = new WH2006(); Console.WriteLine(@"8. Generate 10 random full Int64 range values using Wichmann-Hill’s 2006 combined multiplicative congruential generator"); for (var i = 0; i < 10; i++) { Console.Write(rngWh2006.NextFullRangeInt32() + @" "); } Console.WriteLine(); Console.WriteLine(); // 9. Multiply-with-carry Xorshift pseudo random number generator var xorshift = new Xorshift(); Console.WriteLine(@"9. Generate 10 random nonnegative values less than Int64.MaxValue using Multiply-with-carry Xorshift pseudo random number generator"); for (var i = 0; i < 10; i++) { Console.Write(xorshift.NextInt64() + @" "); } Console.WriteLine(); }
private DreamRadarFrame GeneratePokemon(BWRng pidrng2, MersenneTwister ivrng2) { var pidrng = new BWRng(pidrng2.Seed); var ivrng = new MersenneTwister(ivrng2); var frame = new DreamRadarFrame(); frame.Pid = GeneratePID(pidrng); // two unknown advances pidrng.GetNext64BitNumber(); pidrng.GetNext64BitNumber(); frame.Nature = pidrng.GetNext32BitNumber(25); // IVs frame.Hp = ivrng.Next() >> 27; frame.Atk = ivrng.Next() >> 27; frame.Def = ivrng.Next() >> 27; frame.Spa = ivrng.Next() >> 27; frame.Spd = ivrng.Next() >> 27; frame.Spe = ivrng.Next() >> 27; return frame; }
private void Advance(BWRng pidrng, MersenneTwister ivrng, List<DreamRadarFrame.Spin> spins) { // first PIDRNG advance = spin spins.Add((DreamRadarFrame.Spin) pidrng.GetNext32BitNumber(8)); pidrng.GetNext64BitNumber(); ivrng.Next(); ivrng.Next(); }
// todo: remove this private static void SlotAdvances(BWRng pidrng, MersenneTwister ivrng) { pidrng.GetNext64BitNumber(); pidrng.GetNext64BitNumber(); pidrng.GetNext64BitNumber(); pidrng.GetNext64BitNumber(); pidrng.GetNext64BitNumber(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); ivrng.Next(); }
public List<DreamRadarFrame> Generate(ulong seed, Profile profile) { var frames = new List<DreamRadarFrame>(); // Build the PIDRNG uint initialFrame = Functions.initialPIDRNG(seed, profile); var pidrng = new BWRng(seed); pidrng.Advance(initialFrame); // Build the MTRNG // todo: use fast MTRNG when available var ivrng = new MersenneTwister((uint) (seed >> 32)); // advance 8 frames for BW2 for (uint i = 0; i < 10; ++i) ivrng.Next(); // one single advancement for entering the menu pidrng.GetNext64BitNumber(); var spins = new List<DreamRadarFrame.Spin>(); // initial advances for (uint i = 0; i < initialFrame; ++i) Advance(pidrng, ivrng, spins); // slot advances // we're always doing the slot 1 advance here pidrng.GetNext64BitNumber(); for (uint i = 1; i < TargetSlot; ++i) SlotAdvances(pidrng, ivrng); for (uint i = initialFrame; i <= MaxFrame; ++i) { DreamRadarFrame frame = GeneratePokemon(pidrng, ivrng); var arrSpins = new DreamRadarFrame.Spin[spins.Count]; spins.CopyTo(arrSpins); frame.Spins = arrSpins; // add checks/comparisons on the frame here // nature/IVs frames.Add(frame); Advance(pidrng, ivrng, spins); } return frames; }