diff --git a/src/UnitTests/DistributionTests/CommonDistributionTests.cs b/src/UnitTests/DistributionTests/CommonDistributionTests.cs
index ca28005e..7bc90dab 100644
--- a/src/UnitTests/DistributionTests/CommonDistributionTests.cs
+++ b/src/UnitTests/DistributionTests/CommonDistributionTests.cs
@@ -151,12 +151,12 @@ namespace MathNet.Numerics.UnitTests.DistributionTests
         {
             foreach (var dd in _discreteDistributions)
             {
-                dd.RandomSource = new Random();
+                dd.RandomSource = MersenneTwister.Default;
             }
 
             foreach (var cd in _continuousDistributions)
             {
-                cd.RandomSource = new Random();
+                cd.RandomSource = MersenneTwister.Default;
             }
         }
 
diff --git a/src/UnitTests/DistributionTests/Continuous/BetaTests.cs b/src/UnitTests/DistributionTests/Continuous/BetaTests.cs
index c245d636..08bff7ac 100644
--- a/src/UnitTests/DistributionTests/Continuous/BetaTests.cs
+++ b/src/UnitTests/DistributionTests/Continuous/BetaTests.cs
@@ -274,7 +274,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleStatic()
         {
-            Beta.Sample(new Random(), 2.0, 3.0);
+            Beta.Sample(new Random(0), 2.0, 3.0);
         }
 
         /// <summary>
@@ -283,7 +283,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = Beta.Samples(new Random(), 2.0, 3.0);
+            var ied = Beta.Samples(new Random(0), 2.0, 3.0);
             ied.Take(5).ToArray();
         }
 
@@ -293,7 +293,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Beta.Sample(new Random(), 1.0, -1.0));
+            Assert.Throws<ArgumentOutOfRangeException>(() => Beta.Sample(new Random(0), 1.0, -1.0));
         }
 
         /// <summary>
@@ -302,7 +302,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Beta.Samples(new Random(), 1.0, -1.0).First());
+            Assert.Throws<ArgumentOutOfRangeException>(() => Beta.Samples(new Random(0), 1.0, -1.0).First());
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Continuous/ChiSquareTests.cs b/src/UnitTests/DistributionTests/Continuous/ChiSquareTests.cs
index a0d9c8de..6135cd1d 100644
--- a/src/UnitTests/DistributionTests/Continuous/ChiSquareTests.cs
+++ b/src/UnitTests/DistributionTests/Continuous/ChiSquareTests.cs
@@ -289,7 +289,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleStatic()
         {
-            ChiSquared.Sample(new Random(), 2.0);
+            ChiSquared.Sample(new Random(0), 2.0);
         }
 
         /// <summary>
@@ -298,7 +298,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => ChiSquared.Sample(new Random(), -1.0));
+            Assert.Throws<ArgumentOutOfRangeException>(() => ChiSquared.Sample(new Random(0), -1.0));
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Continuous/ContinuousUniformTests.cs b/src/UnitTests/DistributionTests/Continuous/ContinuousUniformTests.cs
index 1b1bac88..31defe40 100644
--- a/src/UnitTests/DistributionTests/Continuous/ContinuousUniformTests.cs
+++ b/src/UnitTests/DistributionTests/Continuous/ContinuousUniformTests.cs
@@ -322,7 +322,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleStatic()
         {
-            ContinuousUniform.Sample(new Random(), 0.0, 1.0);
+            ContinuousUniform.Sample(new Random(0), 0.0, 1.0);
         }
 
         /// <summary>
@@ -331,7 +331,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = ContinuousUniform.Samples(new Random(), 0.0, 1.0);
+            var ied = ContinuousUniform.Samples(new Random(0), 0.0, 1.0);
             ied.Take(5).ToArray();
         }
 
@@ -341,7 +341,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => ContinuousUniform.Sample(new Random(), 0.0, -1.0));
+            Assert.Throws<ArgumentOutOfRangeException>(() => ContinuousUniform.Sample(new Random(0), 0.0, -1.0));
         }
 
         /// <summary>
@@ -350,7 +350,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => ContinuousUniform.Samples(new Random(), 0.0, -1.0).First());
+            Assert.Throws<ArgumentOutOfRangeException>(() => ContinuousUniform.Samples(new Random(0), 0.0, -1.0).First());
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Continuous/GammaTests.cs b/src/UnitTests/DistributionTests/Continuous/GammaTests.cs
index 98e2e07f..d7ecf57e 100644
--- a/src/UnitTests/DistributionTests/Continuous/GammaTests.cs
+++ b/src/UnitTests/DistributionTests/Continuous/GammaTests.cs
@@ -419,7 +419,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleStatic()
         {
-            Gamma.Sample(new Random(), 1.0, 1.0);
+            Gamma.Sample(new Random(0), 1.0, 1.0);
         }
 
         /// <summary>
@@ -428,7 +428,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = Gamma.Samples(new Random(), 1.0, 1.0);
+            var ied = Gamma.Samples(new Random(0), 1.0, 1.0);
             ied.Take(5).ToArray();
         }
 
@@ -438,7 +438,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Normal.Sample(new Random(), 1.0, -1.0));
+            Assert.Throws<ArgumentOutOfRangeException>(() => Normal.Sample(new Random(0), 1.0, -1.0));
         }
 
         /// <summary>
@@ -447,7 +447,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Normal.Samples(new Random(), 1.0, -1.0).First());
+            Assert.Throws<ArgumentOutOfRangeException>(() => Normal.Samples(new Random(0), 1.0, -1.0).First());
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Continuous/LogNormalTests.cs b/src/UnitTests/DistributionTests/Continuous/LogNormalTests.cs
index d348eee3..26462309 100644
--- a/src/UnitTests/DistributionTests/Continuous/LogNormalTests.cs
+++ b/src/UnitTests/DistributionTests/Continuous/LogNormalTests.cs
@@ -429,7 +429,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleStatic()
         {
-            LogNormal.Sample(new Random(), 0.0, 1.0);
+            LogNormal.Sample(new Random(0), 0.0, 1.0);
         }
 
         /// <summary>
@@ -438,7 +438,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = LogNormal.Samples(new Random(), 0.0, 1.0);
+            var ied = LogNormal.Samples(new Random(0), 0.0, 1.0);
             ied.Take(5).ToArray();
         }
 
@@ -448,7 +448,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => { var d = LogNormal.Sample(new Random(), 0.0, -1.0); });
+            Assert.Throws<ArgumentOutOfRangeException>(() => { var d = LogNormal.Sample(new Random(0), 0.0, -1.0); });
         }
 
         /// <summary>
@@ -457,7 +457,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => { var ied = LogNormal.Samples(new Random(), 0.0, -1.0).First(); });
+            Assert.Throws<ArgumentOutOfRangeException>(() => { var ied = LogNormal.Samples(new Random(0), 0.0, -1.0).First(); });
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Continuous/NormalTests.cs b/src/UnitTests/DistributionTests/Continuous/NormalTests.cs
index c44cef16..dfa427f6 100644
--- a/src/UnitTests/DistributionTests/Continuous/NormalTests.cs
+++ b/src/UnitTests/DistributionTests/Continuous/NormalTests.cs
@@ -402,7 +402,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleStatic()
         {
-            Normal.Sample(new Random(), 0.0, 1.0);
+            Normal.Sample(new Random(0), 0.0, 1.0);
         }
 
         /// <summary>
@@ -411,7 +411,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = Normal.Samples(new Random(), 0.0, 1.0);
+            var ied = Normal.Samples(new Random(0), 0.0, 1.0);
             ied.Take(5).ToArray();
         }
 
@@ -421,7 +421,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => { var d = Normal.Sample(new Random(), 0.0, -1.0); });
+            Assert.Throws<ArgumentOutOfRangeException>(() => { var d = Normal.Sample(new Random(0), 0.0, -1.0); });
         }
 
         /// <summary>
@@ -430,7 +430,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => { var ied = Normal.Samples(new Random(), 0.0, -1.0).First(); });
+            Assert.Throws<ArgumentOutOfRangeException>(() => { var ied = Normal.Samples(new Random(0), 0.0, -1.0).First(); });
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Continuous/StudentTTests.cs b/src/UnitTests/DistributionTests/Continuous/StudentTTests.cs
index c75ddcf4..7e0aa958 100644
--- a/src/UnitTests/DistributionTests/Continuous/StudentTTests.cs
+++ b/src/UnitTests/DistributionTests/Continuous/StudentTTests.cs
@@ -377,7 +377,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleStatic()
         {
-            StudentT.Sample(new Random(), 0.0, 1.0, 3.0);
+            StudentT.Sample(new Random(0), 0.0, 1.0, 3.0);
         }
 
         /// <summary>
@@ -386,7 +386,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = StudentT.Samples(new Random(), 0.0, 1.0, 3.0);
+            var ied = StudentT.Samples(new Random(0), 0.0, 1.0, 3.0);
             ied.Take(5).ToArray();
         }
 
@@ -396,7 +396,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => StudentT.Sample(new Random(), Double.NaN, 1.0, Double.NaN));
+            Assert.Throws<ArgumentOutOfRangeException>(() => StudentT.Sample(new Random(0), Double.NaN, 1.0, Double.NaN));
         }
 
         /// <summary>
@@ -405,7 +405,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleSequenceStatic()
         {
-            var ied = StudentT.Samples(new Random(), 0.0, 1.0, Double.NaN);
+            var ied = StudentT.Samples(new Random(0), 0.0, 1.0, Double.NaN);
             Assert.Throws<ArgumentOutOfRangeException>(() => ied.Take(5).ToArray());
         }
 
diff --git a/src/UnitTests/DistributionTests/Continuous/WeibullTests.cs b/src/UnitTests/DistributionTests/Continuous/WeibullTests.cs
index bb76f330..820c0d2f 100644
--- a/src/UnitTests/DistributionTests/Continuous/WeibullTests.cs
+++ b/src/UnitTests/DistributionTests/Continuous/WeibullTests.cs
@@ -323,7 +323,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleStatic()
         {
-            Weibull.Sample(new Random(), 1.0, 1.0);
+            Weibull.Sample(new Random(0), 1.0, 1.0);
         }
 
         /// <summary>
@@ -332,7 +332,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = Weibull.Samples(new Random(), 1.0, 1.0);
+            var ied = Weibull.Samples(new Random(0), 1.0, 1.0);
             ied.Take(5).ToArray();
         }
 
@@ -342,7 +342,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Normal.Sample(new Random(), 1.0, -1.0));
+            Assert.Throws<ArgumentOutOfRangeException>(() => Normal.Sample(new Random(0), 1.0, -1.0));
         }
 
         /// <summary>
@@ -351,7 +351,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Continuous
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Normal.Samples(new Random(), 1.0, -1.0).First());
+            Assert.Throws<ArgumentOutOfRangeException>(() => Normal.Samples(new Random(0), 1.0, -1.0).First());
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Discrete/BernoulliTests.cs b/src/UnitTests/DistributionTests/Discrete/BernoulliTests.cs
index 2c2b260b..68ac7408 100644
--- a/src/UnitTests/DistributionTests/Discrete/BernoulliTests.cs
+++ b/src/UnitTests/DistributionTests/Discrete/BernoulliTests.cs
@@ -233,7 +233,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void CanSampleStatic()
         {
-            Bernoulli.Sample(new Random(), 0.3);
+            Bernoulli.Sample(new Random(0), 0.3);
         }
 
         /// <summary>
@@ -242,7 +242,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = Bernoulli.Samples(new Random(), 0.3);
+            var ied = Bernoulli.Samples(new Random(0), 0.3);
             ied.Take(5).ToArray();
         }
 
@@ -252,7 +252,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Bernoulli.Sample(new Random(), -1.0));
+            Assert.Throws<ArgumentOutOfRangeException>(() => Bernoulli.Sample(new Random(0), -1.0));
         }
 
         /// <summary>
@@ -261,7 +261,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Bernoulli.Samples(new Random(), -1.0).First());
+            Assert.Throws<ArgumentOutOfRangeException>(() => Bernoulli.Samples(new Random(0), -1.0).First());
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Discrete/BinomialTests.cs b/src/UnitTests/DistributionTests/Discrete/BinomialTests.cs
index 98cc03b1..e599f099 100644
--- a/src/UnitTests/DistributionTests/Discrete/BinomialTests.cs
+++ b/src/UnitTests/DistributionTests/Discrete/BinomialTests.cs
@@ -257,7 +257,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void CanSampleStatic()
         {
-            Binomial.Sample(new Random(), 0.3, 5);
+            Binomial.Sample(new Random(0), 0.3, 5);
         }
 
         /// <summary>
@@ -266,7 +266,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = Binomial.Samples(new Random(), 0.3, 5);
+            var ied = Binomial.Samples(new Random(0), 0.3, 5);
             ied.Take(5).ToArray();
         }
 
@@ -276,7 +276,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Binomial.Sample(new Random(), -1.0, 5));
+            Assert.Throws<ArgumentOutOfRangeException>(() => Binomial.Sample(new Random(0), -1.0, 5));
         }
 
         /// <summary>
@@ -285,7 +285,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Binomial.Samples(new Random(), -1.0, 5).First());
+            Assert.Throws<ArgumentOutOfRangeException>(() => Binomial.Samples(new Random(0), -1.0, 5).First());
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Discrete/CategoricalTests.cs b/src/UnitTests/DistributionTests/Discrete/CategoricalTests.cs
index 8c2037f0..5c56956f 100644
--- a/src/UnitTests/DistributionTests/Discrete/CategoricalTests.cs
+++ b/src/UnitTests/DistributionTests/Discrete/CategoricalTests.cs
@@ -161,7 +161,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void CanSampleStatic()
         {
-            Categorical.SampleWithProbabilityMass(new Random(), _largeP);
+            Categorical.SampleWithProbabilityMass(new Random(0), _largeP);
         }
 
         /// <summary>
@@ -170,7 +170,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Categorical.SampleWithProbabilityMass(new Random(), _badP));
+            Assert.Throws<ArgumentOutOfRangeException>(() => Categorical.SampleWithProbabilityMass(new Random(0), _badP));
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Discrete/DiscreteUniformTests.cs b/src/UnitTests/DistributionTests/Discrete/DiscreteUniformTests.cs
index 16f2ef56..1580af3b 100644
--- a/src/UnitTests/DistributionTests/Discrete/DiscreteUniformTests.cs
+++ b/src/UnitTests/DistributionTests/Discrete/DiscreteUniformTests.cs
@@ -279,7 +279,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void CanSampleStatic()
         {
-            DiscreteUniform.Sample(new Random(), 0, 10);
+            DiscreteUniform.Sample(new Random(0), 0, 10);
         }
 
         /// <summary>
@@ -288,7 +288,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void CanSampleSequenceStatic()
         {
-            var ied = DiscreteUniform.Samples(new Random(), 0, 10);
+            var ied = DiscreteUniform.Samples(new Random(0), 0, 10);
             ied.Take(5).ToArray();
         }
 
@@ -298,7 +298,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => DiscreteUniform.Sample(new Random(), 20, 10));
+            Assert.Throws<ArgumentOutOfRangeException>(() => DiscreteUniform.Sample(new Random(0), 20, 10));
         }
 
         /// <summary>
@@ -307,7 +307,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Discrete
         [Test]
         public void FailSampleSequenceStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => DiscreteUniform.Samples(new Random(), 20, 10).First());
+            Assert.Throws<ArgumentOutOfRangeException>(() => DiscreteUniform.Samples(new Random(0), 20, 10).First());
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Multivariate/DirichletTests.cs b/src/UnitTests/DistributionTests/Multivariate/DirichletTests.cs
index d8b08bbc..cf74eba4 100644
--- a/src/UnitTests/DistributionTests/Multivariate/DirichletTests.cs
+++ b/src/UnitTests/DistributionTests/Multivariate/DirichletTests.cs
@@ -110,7 +110,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         {
             new Dirichlet(0.3, 5)
             {
-                RandomSource = new Random()
+                RandomSource = new Random(0)
             };
         }
 
diff --git a/src/UnitTests/DistributionTests/Multivariate/InverseWishartTests.cs b/src/UnitTests/DistributionTests/Multivariate/InverseWishartTests.cs
index 4d50814f..9bc60908 100644
--- a/src/UnitTests/DistributionTests/Multivariate/InverseWishartTests.cs
+++ b/src/UnitTests/DistributionTests/Multivariate/InverseWishartTests.cs
@@ -117,7 +117,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         {
             new InverseWishart(1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2))
             {
-                RandomSource = new Random()
+                RandomSource = new Random(0)
             };
         }
 
@@ -303,7 +303,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         [Test]
         public void CanSampleStatic()
         {
-            InverseWishart.Sample(new Random(), 1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2));
+            InverseWishart.Sample(new Random(0), 1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2));
         }
 
         /// <summary>
@@ -312,7 +312,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => InverseWishart.Sample(new Random(), -1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2)));
+            Assert.Throws<ArgumentOutOfRangeException>(() => InverseWishart.Sample(new Random(0), -1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2)));
         }
     }
 }
diff --git a/src/UnitTests/DistributionTests/Multivariate/MatrixNormalTests.cs b/src/UnitTests/DistributionTests/Multivariate/MatrixNormalTests.cs
index 8275162b..2c835e92 100644
--- a/src/UnitTests/DistributionTests/Multivariate/MatrixNormalTests.cs
+++ b/src/UnitTests/DistributionTests/Multivariate/MatrixNormalTests.cs
@@ -126,7 +126,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
             const int P = 3;
             new MatrixNormal(MatrixLoader.GenerateRandomDenseMatrix(N, P), MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(N), MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(P))
             {
-                RandomSource = new Random()
+                RandomSource = new Random(0)
             };
         }
 
@@ -321,7 +321,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         [TestCase(10, 10)]
         public void CanSampleStatic(int n, int p)
         {
-            MatrixNormal.Sample(new Random(), MatrixLoader.GenerateRandomDenseMatrix(n, p), MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(n), MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(p));
+            MatrixNormal.Sample(new Random(0), MatrixLoader.GenerateRandomDenseMatrix(n, p), MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(n), MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(p));
         }
 
         /// <summary>
@@ -343,7 +343,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         [TestCase(5, 2, 5, 5, 2, 3)]
         public void FailSampleStatic(int rowsOfM, int columnsOfM, int rowsOfV, int columnsOfV, int rowsOfK, int columnsOfK)
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => MatrixNormal.Sample(new Random(), MatrixLoader.GenerateRandomDenseMatrix(rowsOfM, columnsOfM), MatrixLoader.GenerateRandomDenseMatrix(rowsOfV, columnsOfV), MatrixLoader.GenerateRandomDenseMatrix(rowsOfK, columnsOfK)));
+            Assert.Throws<ArgumentOutOfRangeException>(() => MatrixNormal.Sample(new Random(0), MatrixLoader.GenerateRandomDenseMatrix(rowsOfM, columnsOfM), MatrixLoader.GenerateRandomDenseMatrix(rowsOfV, columnsOfV), MatrixLoader.GenerateRandomDenseMatrix(rowsOfK, columnsOfK)));
         }
     }
 }
diff --git a/src/UnitTests/DistributionTests/Multivariate/MultinomialTests.cs b/src/UnitTests/DistributionTests/Multivariate/MultinomialTests.cs
index 26df25d9..2d7d3775 100644
--- a/src/UnitTests/DistributionTests/Multivariate/MultinomialTests.cs
+++ b/src/UnitTests/DistributionTests/Multivariate/MultinomialTests.cs
@@ -245,7 +245,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         [Test]
         public void CanSampleStatic()
         {
-            Multinomial.Sample(new Random(), _largeP, 4);
+            Multinomial.Sample(new Random(0), _largeP, 4);
         }
 
         /// <summary>
@@ -254,7 +254,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Multinomial.Sample(new Random(), _badP, 4));
+            Assert.Throws<ArgumentOutOfRangeException>(() => Multinomial.Sample(new Random(0), _badP, 4));
         }
 
         /// <summary>
diff --git a/src/UnitTests/DistributionTests/Multivariate/NormalGammaTests.cs b/src/UnitTests/DistributionTests/Multivariate/NormalGammaTests.cs
index a04c28b6..8df7afec 100644
--- a/src/UnitTests/DistributionTests/Multivariate/NormalGammaTests.cs
+++ b/src/UnitTests/DistributionTests/Multivariate/NormalGammaTests.cs
@@ -304,7 +304,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         {
             new NormalGamma(0.0, 1.0, 1.0, 1.0)
             {
-                RandomSource = new Random()
+                RandomSource = new Random(0)
             };
         }
 
diff --git a/src/UnitTests/DistributionTests/Multivariate/WishartTests.cs b/src/UnitTests/DistributionTests/Multivariate/WishartTests.cs
index 4c6dd717..cf22e27e 100644
--- a/src/UnitTests/DistributionTests/Multivariate/WishartTests.cs
+++ b/src/UnitTests/DistributionTests/Multivariate/WishartTests.cs
@@ -119,7 +119,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         {
             new Wishart(1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2))
             {
-                RandomSource = new Random()
+                RandomSource = new Random(0)
             };
         }
 
@@ -305,7 +305,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         [Test]
         public void CanSampleStatic()
         {
-            Wishart.Sample(new Random(), 1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2));
+            Wishart.Sample(new Random(0), 1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2));
         }
 
         /// <summary>
@@ -314,7 +314,7 @@ namespace MathNet.Numerics.UnitTests.DistributionTests.Multivariate
         [Test]
         public void FailSampleStatic()
         {
-            Assert.Throws<ArgumentOutOfRangeException>(() => Wishart.Sample(new Random(), -1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2)));
+            Assert.Throws<ArgumentOutOfRangeException>(() => Wishart.Sample(new Random(0), -1.0, MatrixLoader.GenerateRandomPositiveDefiniteDenseMatrix(2)));
         }
     }
 }
diff --git a/src/UnitTests/LinearAlgebraTests/Complex/SparseMatrixTests.cs b/src/UnitTests/LinearAlgebraTests/Complex/SparseMatrixTests.cs
index 25a2cfd5..555f5475 100644
--- a/src/UnitTests/LinearAlgebraTests/Complex/SparseMatrixTests.cs
+++ b/src/UnitTests/LinearAlgebraTests/Complex/SparseMatrixTests.cs
@@ -191,7 +191,7 @@ namespace MathNet.Numerics.UnitTests.LinearAlgebraTests.Complex
         {
             var matrix = new SparseMatrix(500, 1000);
             var nonzero = 0;
-            var rnd = new System.Random();
+            var rnd = new System.Random(0);
 
             for (var i = 0; i < matrix.RowCount; i++)
             {
diff --git a/src/UnitTests/LinearAlgebraTests/Complex32/SparseMatrixTests.cs b/src/UnitTests/LinearAlgebraTests/Complex32/SparseMatrixTests.cs
index 1166b1cb..e8e904a7 100644
--- a/src/UnitTests/LinearAlgebraTests/Complex32/SparseMatrixTests.cs
+++ b/src/UnitTests/LinearAlgebraTests/Complex32/SparseMatrixTests.cs
@@ -187,7 +187,7 @@ namespace MathNet.Numerics.UnitTests.LinearAlgebraTests.Complex32
         {
             var matrix = new SparseMatrix(500, 1000);
             var nonzero = 0;
-            var rnd = new System.Random();
+            var rnd = new System.Random(0);
 
             for (var i = 0; i < matrix.RowCount; i++)
             {
diff --git a/src/UnitTests/LinearAlgebraTests/Double/SparseMatrixTests.cs b/src/UnitTests/LinearAlgebraTests/Double/SparseMatrixTests.cs
index 5941df9d..078e8e3b 100644
--- a/src/UnitTests/LinearAlgebraTests/Double/SparseMatrixTests.cs
+++ b/src/UnitTests/LinearAlgebraTests/Double/SparseMatrixTests.cs
@@ -185,7 +185,7 @@ namespace MathNet.Numerics.UnitTests.LinearAlgebraTests.Double
         {
             var matrix = new SparseMatrix(500, 1000);
             var nonzero = 0;
-            var rnd = new System.Random();
+            var rnd = new System.Random(0);
 
             for (var i = 0; i < matrix.RowCount; i++)
             {
diff --git a/src/UnitTests/LinearAlgebraTests/Single/SparseMatrixTests.cs b/src/UnitTests/LinearAlgebraTests/Single/SparseMatrixTests.cs
index 45182848..8785b0d0 100644
--- a/src/UnitTests/LinearAlgebraTests/Single/SparseMatrixTests.cs
+++ b/src/UnitTests/LinearAlgebraTests/Single/SparseMatrixTests.cs
@@ -185,7 +185,7 @@ namespace MathNet.Numerics.UnitTests.LinearAlgebraTests.Single
         {
             var matrix = new SparseMatrix(500, 1000);
             var nonzero = 0;
-            var rnd = new System.Random();
+            var rnd = new System.Random(0);
 
             for (var i = 0; i < matrix.RowCount; i++)
             {
diff --git a/src/UnitTests/Random/SystemRandomExtensionTests.cs b/src/UnitTests/Random/SystemRandomExtensionTests.cs
index 7e9bf07d..9c2e8bd6 100644
--- a/src/UnitTests/Random/SystemRandomExtensionTests.cs
+++ b/src/UnitTests/Random/SystemRandomExtensionTests.cs
@@ -42,7 +42,7 @@ namespace MathNet.Numerics.UnitTests.Random
         [Test]
         public void CanSampleInt64()
         {
-            var rnd = new Random();
+            var rnd = new Random(0);
             rnd.NextInt64();
         }
 
@@ -52,7 +52,7 @@ namespace MathNet.Numerics.UnitTests.Random
         [Test]
         public void CanSampleFullRangeInt32()
         {
-            var rnd = new Random();
+            var rnd = new Random(0);
             rnd.NextFullRangeInt32();
         }
 
@@ -62,7 +62,7 @@ namespace MathNet.Numerics.UnitTests.Random
         [Test]
         public void CanSampleFullRangeInt64()
         {
-            var rnd = new Random();
+            var rnd = new Random(0);
             rnd.NextFullRangeInt64();
         }
 
@@ -72,7 +72,7 @@ namespace MathNet.Numerics.UnitTests.Random
         [Test]
         public void CanSampleDecimal()
         {
-            var rnd = new Random();
+            var rnd = new Random(0);
             rnd.NextDecimal();
         }
     }
diff --git a/src/UnitTests/SortingTests.cs b/src/UnitTests/SortingTests.cs
index 7fe5133d..e42c69cc 100644
--- a/src/UnitTests/SortingTests.cs
+++ b/src/UnitTests/SortingTests.cs
@@ -43,7 +43,7 @@ namespace MathNet.Numerics.UnitTests
         public void TestRandomTupleArraySorting()
         {
             const int Len = 0x1 << 10;
-            var random = new System.Random();
+            var random = new System.Random(0);
 
             var keys = new int[Len];
             var items = new int[Len];
@@ -77,7 +77,7 @@ namespace MathNet.Numerics.UnitTests
         public void TestRandomTupleListSorting()
         {
             const int Len = 0x1 << 10;
-            var random = new System.Random();
+            var random = new System.Random(0);
 
             var keys = new List<int>(Len);
             var items = new List<int>(Len);
@@ -112,7 +112,7 @@ namespace MathNet.Numerics.UnitTests
         public void TestRandomTripleArraySorting()
         {
             const int Len = 0x1 << 10;
-            var random = new System.Random();
+            var random = new System.Random(0);
 
             var keys = new int[Len];
             var items1 = new int[Len];
@@ -149,7 +149,7 @@ namespace MathNet.Numerics.UnitTests
         public void TestAppliedListSorting()
         {
             const int Len = 0x1 << 10;
-            var random = new System.Random();
+            var random = new System.Random(0);
 
             var list = new List<int>();
 
diff --git a/src/UnitTests/StatisticsTests/MCMCTests/HybridMCTest.cs b/src/UnitTests/StatisticsTests/MCMCTests/HybridMCTest.cs
index dd55ebe3..6bb0ff08 100644
--- a/src/UnitTests/StatisticsTests/MCMCTests/HybridMCTest.cs
+++ b/src/UnitTests/StatisticsTests/MCMCTests/HybridMCTest.cs
@@ -55,7 +55,7 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
             var hybrid = new HybridMC(new double[] { 0 }, x => _normal.DensityLn(x[0]), 10, 0.1);
             Assert.IsNotNull(hybrid.RandomSource);
 
-            hybrid.RandomSource = new System.Random();
+            hybrid.RandomSource = new System.Random(0);
             Assert.IsNotNull(hybrid.RandomSource);
         }
 
diff --git a/src/UnitTests/StatisticsTests/MCMCTests/MCMCDiagnosticsTest.cs b/src/UnitTests/StatisticsTests/MCMCTests/MCMCDiagnosticsTest.cs
index 99075ee3..80f93435 100644
--- a/src/UnitTests/StatisticsTests/MCMCTests/MCMCDiagnosticsTest.cs
+++ b/src/UnitTests/StatisticsTests/MCMCTests/MCMCDiagnosticsTest.cs
@@ -45,7 +45,7 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
         /// <summary>
         /// For generation of a random series to test the methods.
         /// </summary>
-        private readonly System.Random _rnd = new System.Random();
+        private readonly System.Random _rnd = new System.Random(0);
         /// <summary>
         /// Distribution to sample the entries of the random series from.
         /// </summary>
diff --git a/src/UnitTests/StatisticsTests/MCMCTests/MetropolisHastingsSamplerTests.cs b/src/UnitTests/StatisticsTests/MCMCTests/MetropolisHastingsSamplerTests.cs
index 089fdd36..4a490a52 100644
--- a/src/UnitTests/StatisticsTests/MCMCTests/MetropolisHastingsSamplerTests.cs
+++ b/src/UnitTests/StatisticsTests/MCMCTests/MetropolisHastingsSamplerTests.cs
@@ -36,8 +36,6 @@ using NUnit.Framework;
 
 namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
 {
-    using Random = System.Random;
-
     /// <summary>
     /// Metropolis hastings sampler tests.
     /// </summary>
@@ -53,13 +51,13 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
             var normal = new Normal(0.0, 1.0);
             var rnd = new MersenneTwister();
 
-            var ms = new MetropolisHastingsSampler<double>(0.2, normal.Density, (x, y) => (new Normal(x, 0.1)).Density(y), x => Normal.Sample(rnd, x, 0.1), 10)
+            var ms = new MetropolisHastingsSampler<double>(0.2, normal.Density, (x, y) => Normal.PDF(x, 0.1, y), x => Normal.Sample(rnd, x, 0.1), 10)
                 {
                     RandomSource = rnd
                 };
             Assert.IsNotNull(ms.RandomSource);
 
-            ms.RandomSource = new System.Random();
+            ms.RandomSource = new System.Random(0);
             Assert.IsNotNull(ms.RandomSource);
         }
 
@@ -72,7 +70,7 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
             var normal = new Normal(0.0, 1.0);
             var rnd = new MersenneTwister();
 
-            var ms = new MetropolisHastingsSampler<double>(0.2, normal.Density, (x, y) => (new Normal(x, 0.1)).Density(y), x => Normal.Sample(rnd, x, 0.1), 10)
+            var ms = new MetropolisHastingsSampler<double>(0.2, normal.Density, (x, y) => Normal.PDF(x, 0.1, y), x => Normal.Sample(rnd, x, 0.1), 10)
                 {
                     RandomSource = rnd
                 };
@@ -89,7 +87,7 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
             var normal = new Normal(0.0, 1.0);
             var rnd = new MersenneTwister();
 
-            var ms = new MetropolisHastingsSampler<double>(0.2, normal.Density, (x, y) => (new Normal(x, 0.1)).Density(y), x => Normal.Sample(rnd, x, 0.1), 10)
+            var ms = new MetropolisHastingsSampler<double>(0.2, normal.Density, (x, y) => Normal.PDF(x, 0.1, y), x => Normal.Sample(rnd, x, 0.1), 10)
                 {
                     RandomSource = rnd
                 };
@@ -103,8 +101,9 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
         [Test]
         public void NullRandomNumberGenerator()
         {
-            var normal = new Normal(0.0, 1.0);
-            var ms = new MetropolisHastingsSampler<double>(0.2, normal.Density, (x, y) => (new Normal(x, 0.1)).Density(y), x => Normal.Sample(new Random(), x, 0.1), 10);
+            var random = MersenneTwister.Default;
+            var normal = new Normal(0.0, 1.0, random);
+            var ms = new MetropolisHastingsSampler<double>(0.2, normal.Density, (x, y) => Normal.PDF(x, 0.1, y), x => Normal.Sample(random, x, 0.1), 10);
             Assert.Throws<ArgumentNullException>(() => ms.RandomSource = null);
         }
     }
diff --git a/src/UnitTests/StatisticsTests/MCMCTests/MetropolisSamplerTests.cs b/src/UnitTests/StatisticsTests/MCMCTests/MetropolisSamplerTests.cs
index e7112393..53462caf 100644
--- a/src/UnitTests/StatisticsTests/MCMCTests/MetropolisSamplerTests.cs
+++ b/src/UnitTests/StatisticsTests/MCMCTests/MetropolisSamplerTests.cs
@@ -101,7 +101,7 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
         public void NullRandomNumberGenerator()
         {
             var normal = new Normal(0.0, 1.0);
-            var ms = new MetropolisSampler<double>(0.2, normal.Density, x => Normal.Sample(new Random(), x, 0.1), 10);
+            var ms = new MetropolisSampler<double>(0.2, normal.Density, x => Normal.Sample(new Random(0), x, 0.1), 10);
             Assert.Throws<ArgumentNullException>(() => ms.RandomSource = null);
         }
     }
diff --git a/src/UnitTests/StatisticsTests/MCMCTests/UnivariateHybridMCTest.cs b/src/UnitTests/StatisticsTests/MCMCTests/UnivariateHybridMCTest.cs
index fe750892..5e97a94e 100644
--- a/src/UnitTests/StatisticsTests/MCMCTests/UnivariateHybridMCTest.cs
+++ b/src/UnitTests/StatisticsTests/MCMCTests/UnivariateHybridMCTest.cs
@@ -58,7 +58,7 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
             var hybrid = new UnivariateHybridMC(0, _normal.DensityLn, 10, 0.1);
             Assert.IsNotNull(hybrid.RandomSource);
 
-            hybrid.RandomSource = new System.Random();
+            hybrid.RandomSource = new System.Random(0);
             Assert.IsNotNull(hybrid.RandomSource);
         }
 
diff --git a/src/UnitTests/StatisticsTests/MCMCTests/UnivariateSliceSamplerTests.cs b/src/UnitTests/StatisticsTests/MCMCTests/UnivariateSliceSamplerTests.cs
index b23b7876..24e7c672 100644
--- a/src/UnitTests/StatisticsTests/MCMCTests/UnivariateSliceSamplerTests.cs
+++ b/src/UnitTests/StatisticsTests/MCMCTests/UnivariateSliceSamplerTests.cs
@@ -78,7 +78,7 @@ namespace MathNet.Numerics.UnitTests.StatisticsTests.McmcTests
             var ss = new UnivariateSliceSampler(0.1, x => -0.5*x*x, 5, 1.0);
 
             Assert.IsNotNull(ss.RandomSource);
-            ss.RandomSource = new Random();
+            ss.RandomSource = new Random(0);
             Assert.IsNotNull(ss.RandomSource);
         }