...

tests/plot_stream.cpp

@@ -39,7 +39,7 @@
 //set a global seed
 static const uint64_t default_seed= 0x548c9decbce65295  ;
 
-#define NCHECKS  1024
+#define NSAMPLES  1024
 
 
 
@@ -48,14 +48,36 @@ int main(int argv, char** argc)
 
 	    random_iterator::philox RNG(seed);
 
-		std::uniform_real_distribution<double> udist(0.0, 10.0);
-		std::normal_distribution<double> ndist(0.0, 1.0);
+		std::uniform_real_distribution<double> udist(0.0, 1.0);
+		std::normal_distribution<double>       ndist(0.0, 1.0);
+		std::exponential_distribution<double>  edist(1.0);
 
+		auto urng_stream = random_iterator::make_stream( udist, RNG, 0);
+		auto nrng_stream = random_iterator::make_stream( ndist, RNG, 1);
+		auto erng_stream = random_iterator::make_stream( edist, RNG, 1);
 
+		std::vector<double> data_normal_distribution;
+		std::vector<double> data_uniform_distribution;
+		std::vector<double> data_exponential_distribution;
 
-		auto rng_stream = random_iterator::make_stream( dist, RNG, stream);
 
-		//histogram as a simple vector
-		std::vector<size_t> histogram(10);
+		for(size_t i =0; i<NSAMPLES; ++i){
+
+			     data_normal_distribution.push_back(nrng_stream[i]);
+		    	data_uniform_distribution.push_back(urng_stream[i]);
+			data_exponential_distribution.push_back(erng_stream[i]);
+
+		}
+
+		subplot(1, 3, 0);
+		auto data_normal_distribution  = hist(data_normal_distribution);
+		subplot(1, 3, 1);
+		auto data_uniform_distribution = hist(data_uniform_distribution);
+		subplot(1, 3, 2);
+
+
+		show();
+
+		return 0;
 
 }