diff --git a/freqs1M.odc b/freqs1M.odc deleted file mode 100644 index 0356a12ab9089bd159ab9b48dd75251ccabf3a91..0000000000000000000000000000000000000000 Binary files a/freqs1M.odc and /dev/null differ diff --git a/gva.cpp b/gva.cpp deleted file mode 100644 index ebd039217cd4876b8cfc83229d22e7009fb375a0..0000000000000000000000000000000000000000 --- a/gva.cpp +++ /dev/null @@ -1,60 +0,0 @@ -#include <math.h> -#include <stdio.h> -#include <iostream> -#include <fstream> - -int x0PM = 110425; -unsigned int GVA(int n,int a,int b,long int m){ - std::ofstream myfile, freqcsv; - myfile.open ("gva.csv"); - freqcsv.open ("freqs.csv"); - myfile << n << ", " << x0PM << ", ParkMiller \n"; - unsigned int i, res=x0PM; - unsigned int freqs[10]; - for (i =0; i < 10; i ++){ - freqs[i] = 0; - } - - for (i=0; i< n; i ++){ - res=(a*res + b)% (m); - // res=16807*res % ((long int) (pow(2,31)-1)); - freqs[(unsigned int) ((res/(pow(2,31)-1))*10)]+=1; - myfile << ((res/(pow(2,31)-1))); - myfile << "\n"; - } - for (i =0; i < 10; i ++){ - freqcsv << "0."; - freqcsv << i; - freqcsv << ", "; - freqcsv << freqs[i]; - freqcsv << "\n"; - } - myfile.close(); - freqcsv.close(); - return res; -} - -int main(int argc, char ** argv) -{ - int n =10000000; - int a = 16807; - int b = 0; - unsigned int m = (pow(2,31)-1); - - GVA(n,a,b,m); - - /* int nIter; - if (argc < 2 ){ - nIter=100; - } - else{ - nIter = atoi(argv[1]); - } - int data[nIter], i; - for (i = 0; i < nIter; i ++) - data[i] = ParkMiller(i); - //printf("%d \n", ParkMiller(i));*/ - return 0; -} - - diff --git a/gva.csv b/gva.csv deleted file mode 100644 index adaae02571be14045d020d4e466369b3f040d8e0..0000000000000000000000000000000000000000 Binary files a/gva.csv and /dev/null differ diff --git a/gva.odc b/gva.odc deleted file mode 100644 index cc1fd29ff5b0262a40f7714f398aa116b554109a..0000000000000000000000000000000000000000 Binary files a/gva.odc and /dev/null differ diff --git a/gva1M.csv b/gva1M.csv deleted file mode 100644 index e59c504fd79c59aa9eb1867456b018f5032bc280..0000000000000000000000000000000000000000 Binary files a/gva1M.csv and /dev/null differ diff --git a/gva50.csv b/gva50.csv deleted file mode 100644 index 87af24b62cd6a089fa3114d59d7e60382c7d2a64..0000000000000000000000000000000000000000 --- a/gva50.csv +++ /dev/null @@ -1,60 +0,0 @@ -7.82637e-06 -0.131538 -0.755605 -0.45865 -0.532767 -0.218959 -0.0470446 -0.678865 -0.679296 -0.934693 -0.383502 -0.519416 -0.830965 -0.0345721 -0.0534616 -0.5297 -0.671149 -0.00769819 -0.383416 -0.0668422 -0.417486 -0.686773 -0.588977 -0.930436 -0.846167 -0.526929 -0.0919649 -0.653919 -0.415999 -0.701191 -0.910321 -0.762198 -0.262453 -0.0474645 -0.736082 -0.328234 -0.632639 -0.75641 -0.991037 -0.365339 -0.247039 -0.98255 -0.72266 -0.753356 -0.651519 -0.0726859 -0.631635 -0.884707 -0.27271 -0.436411 -9 -1 -4 -4 -4 -5 -8 -7 -3 -5 diff --git a/gvg.cpp b/gvg.cpp deleted file mode 100644 index 58c1b2838f9018f0d3fded7df858c34d7f790a48..0000000000000000000000000000000000000000 --- a/gvg.cpp +++ /dev/null @@ -1,29 +0,0 @@ -#include <cstdlib> -#include <cmath> -#include <limits> -double generateGaussianNoise(double mu, double sigma) -{ - static const double epsilon = std::numeric_limits<double>::min(); - static const double two_pi = 2.0*3.14159265358979323846; - - thread_local double z1; - thread_local bool generate; - generate = !generate; - - if (!generate) - return z1 * sigma + mu; - - double u1, u2; - do - { - u1 = rand() * (1.0 / RAND_MAX); - u2 = rand() * (1.0 / RAND_MAX); - } - while ( u1 <= epsilon ); - - double z0; - z0 = sqrt(-2.0 * log(u1)) * cos(two_pi * u2); - z1 = sqrt(-2.0 * log(u1)) * sin(two_pi * u2); - return z0 * sigma + mu; -} -