gva rien a dire

13864eb7 · Hugo TRACHINO · 03da37fb · 13864eb7 · 13864eb7 · 13864eb7
--- a/PM.cpp
+++ b/PM.cpp
-#include <iostream>
-#include <fstream>
-#include <stdio.h>
-#include <math.h>
-using namespace std;
-const unsigned int n = 0xFFFF;
-const unsigned int m = 16807;
-const unsigned int b = 0;
-static long  seed = 0xAAAA;
-long uniform() {
-   long unsigned int hi = m * (seed & n);
-   long unsigned int lo = m * (seed >> 16);
-
-   lo+= (hi & 0x7FFF) << 16;
-   lo+= hi >> 15;   
-   
-   if (lo > 0x7FFFFFFFF)
-      lo -= 0x7FFFFFFF;
-   seed=(long) lo;
-   return seed %10;
-}
-void store_res(long * res, int freqs[10], int size){
-	std::ofstream myfile, freqcsv;
-	myfile.open("gva_res.csv");
-	freqcsv.open("gva_freqs.csv");
-	int i;
-	for (i = 0; i < size; i ++){
-		freqs[res[i]] ++;
-		myfile << res[i];
-		myfile << "\n";
-	}
-	for (i=0; i < 10; i++){
-		freqcsv << "0.";
-		freqcsv << i;
-		freqcsv << ", ";
-		freqcsv << freqs[i];
-		freqcsv << "\n";
-	}
-	myfile.close();
-	freqcsv.close();
-}
+#include "fcts.h"

 int main(int argc, char **argv) {
-   if(argc <2)
+   int nRNG, nCuts;
+   if(argc <2){
+	   printf("Please give the number of random generated number");
 	   return -1;
-   long unif[atoi(argv[1])];
+   }
+   nRNG = atoi(argv[1]);
+   if(argc <3){
+   	   nCuts = 10;
+   }
+   else
+	   nCuts = atoi(argv[2]);
+
+   long unif[nRNG];
   // create the outputs
-   for (int i = 0; i < atoi(argv[1]); i++)
+   for (int i = 0; i < nRNG; i++)
      unif[i] = uniform();
   
   // store hist freqs
-   int freqs[10];
-   for(int i=0; i < 10; i ++)
+   int freqs[nCuts];
+   for(int i=0; i < nCuts; i ++)
 	   freqs[i] = 0;
-   store_res(unif, freqs, atoi(argv[1]));
-   for (int i = 0; i < atoi(argv[1]); i++)
-      cout << unif[i] << endl;
+   printf("%d \n", nRNG);
+   store_res(unif, freqs, nRNG, nCuts);
 }
--- a/fcts.h
+++ b/fcts.h
+#include <iostream>
+#include <fstream>
+#include <stdio.h>
+#include <math.h>
+using namespace std;
+const unsigned int n = 0xFFFF;
+const unsigned int m = 16807;
+const unsigned int b = 0;
+static long  seed = 0xAAAA;
+long uniform() {
+   long unsigned int hi = m * (seed & n);
+   long unsigned int lo = m * (seed >> 16);
+   
+   lo+= (hi & 0x7FFF) << 16;
+   lo+= hi >> 15;
+   
+   if (lo > 0x7FFFFFFFF)
+      lo -= 0x7FFFFFFF;
+   seed=(long) lo;
+   return seed ;
+}
+//void gaussian(){
+//  Y0 = sqrt(-2* log(U0))*cos(2* pi * U1);
+//  Y1 = sqrt(-2* log(U0))*sin(2* pi * U1);
+
+//}
+
+void store_res(long * res, int freqs[], int nRNG, int nCuts){
+	std::ofstream myfile, freqcsv;
+	myfile.open("gva_res.csv");
+	freqcsv.open("gva_freqs.csv");
+	int i;
+	for (i = 0; i < nRNG; i ++){
+		printf("i = %d, nRNG= %d\n ", i, nRNG);
+		freqs[res[i]%nCuts] ++;
+		myfile << res[i];
+		myfile << "\n";
+	}
+	for (i=0; i < nCuts; i++){
+		freqcsv << "0.";
+		freqcsv << i;
+		freqcsv << ", ";
+		freqcsv << freqs[i];
+		freqcsv << "\n";
+	}
+	freqcsv << "1.0, 0";
+	myfile.close();
+	freqcsv.close();
+}
--- a/gva_freqs.csv
+++ b/gva_freqs.csv
+0.0, 42
+0.1, 49
+0.2, 55
+0.3, 45
+0.4, 50
+0.5, 54
+0.6, 52
+0.7, 47
+0.8, 48
+0.9, 55
+0.10, 47
+0.11, 41
+0.12, 57
+0.13, 49
+0.14, 50
+0.15, 60
+0.16, 61
+0.17, 50
+0.18, 47
+0.19, 41
+1.0, 0
\ No newline at end of file
--- a/gva_res.csv
+++ b/gva_res.csv
+2129024904
+1111710447
+300378969
+614894528
+1940280823
+1758053217
+807848021
+1355576714
+1270791517
+941099026
+1679864779
+1847250199
+1027462695
+707124642
+1098732752
+880517387
+1721944534
+451706164
+1741929515
+1169010216
+1442232878
+940435704
+2150126840
+2259021228
+1505506442
+1728667798
+843493737
+1902374056
+648844279
+1191981592
+1474330245
+1682732036
+1280636181
+1229364352
+709550831
+1136776363
+938320123
+1293220652
+427351264
+2091419281
+965688989
+1143864166
+1996908364
+1293574117
+1193875237
+2422399133
+2289206562
+2653698214
+1121613762
+1924346786
+1847091264
+637278982
+816304448
+1295683396
+1294623825
+362886106
+851580831
+1421699755
+373848833
+2220757006
+689192846
+237699237
+1413818721
+887329596
+1437888153
+2496255269
+2521531761
+1682070110
+1795708775
+1185154828
+921812121
+1860593906
+2219766749
+1410963958
+2064352326
+792220421
+2036281372
+2640602277
+889209819
+1819588045
+1099665643
+1159485185
+341859812
+687200968
+184114442
+1045860647
+661506701
+2031602942
+1455439397
+2237822751
+913187335
+2273739061
+1632940207
+1243556971
+936333313
+2146902033
+1558725011
+834184442
+634813350
+1493058565
+991276693
+1197095230
+666796818
+485453356
+304072555
+1181947777
+674771987
+1437962045
+2210801546
+1993429105
+1227865241
+1939082596
+1038090346
+1057242454
+1766014737
+2417348109
+1550761058
+1570937441
+1850377692
+1254949102
+477157685
+769551480
+1175159156
+2023931538
+2184066412
+1725347176
+1496828069
+1602523338
+1447370079
+1578687251
+1418127877
+2213585265
+596464734
+1886398784
+1167449656
+1080077101
+668841580
+1991080891
+1349293905
+1534537744
+2179267143
+1822498842
+651284094
+1009175872
+875383824
+1003593054
+363423502
+1303396122
+1474461229
+1390737189
+862287231
+1913088038
+872963836
+2143037935
+715859238
+1824215962
+1062769630
+1678833215
+1002624189
+1402907307
+1442775451
+2099960107
+908509129
+855792424
+741144335
+2233301867
+1341157944
+1325944869
+1097318227
+1928016722
+1777777823
+2246440504
+2027879340
+2016655335
+1032108696
+1626277072
+2357980512
+1571517293
+2434842816
+2725785276
+1649314880
+1458976465
+1772375191
+1754838340
+506408993
+1313126488
+763015401
+480837459
+696179853
+933199953
+1411049865
+921759422
+2315862275
+1331394103
+691739647
+216788697
+1809814222
+772556997
+2303867686
+1593935060
+2460581575
+1617811013
+1690173953
+696055545
+1611860624
+2358478031
+1444301571
+1208542419
+1797159573
+1739308401
+1615686839
+2064881405
+1536911830
+2014860133
+716151406
+917804638
+2203693195
+2539430555
+2498653290
+1230490845
+1408915412
+886402850
+783992051
+964762202
+427119019
+940866102
+2107369402
+2223817207
+1629447888
+2056810228
+1956432400
+2119701863
+1033358946
+330952284
+1595101424
+1425146344
+1008786648
+2018757157
+1644099795
+2110188113
+2350434904
+2069226352
+1840348698
+2031599652
+2592749453
+832370308
+1624602420
+2411822358
+1412958223
+643968170
+1455293142
+1171838310
+762565980
+380388980
+1870436950
+548504845
+551389365
+729862985
+1757891806
+498141246
+638538673
+322555318
+1822977201
+1951711575
+1419155958
+2066453201
+1995280218
+2083916577
+1130480367
+1781594434
+2001335592
+867161949
+669815967
+1559642926
+416890095
+1967680975
+1327301194
+1552939543
+2478715022
+1619376089
+542120463
+1460826509
+1984927995
+622103296
+914525015
+2327566405
+1033333226
+292666068
+952463754
+697649087
+1464015359
+1253702922
+783300984
+2168537673
+1740960870
+1864823350
+1404806149
+1505521181
+1223589406
+1984315889
+1401687932
+726733893
+1226759484
+551862795
+1755102274
+736246355
+711336024
+2085068576
+2378132206
+1402744177
+1797799442
+725167854
+374850369
+1713507814
+2436581710
+1159781910
+2283055560
+744370209
+1374152705
+782797744
+1539792701
+894736042
+2157144850
+2344053199
+954061222
+2011144370
+618406885
+416741869
+1023019466
+1013532187
+639317364
+1550963207
+1632004720
+2013425985
+2167276900
+2304570821
+2076013842
+2423910742
+1780970423
+2392482250
+962089630
+298248904
+1963415189
+1849051897
+1035125166
+374658165
+1393381828
+807969899
+2032592492
+1535517689
+1206373597
+1402730436
+2051892273
+1817479326
+1490699486
+1637577608
+2126529118
+1238602058
+607454762
+950998716
+523072722
+251461446
+1350714706
+522431068
+2248182325
+1072738531
+393936970
+2152434786
+2530542573
+1665760246
+1827704115
+724912380
+595118594
+1382482660
+1541269932
+851673605
+591913701
+829390513
+589606756
+1463773809
+924262321
+731371330
+960247017
+1001172295
+530510356
+2267560507
+2586210549
+2354802361
+1471013671
+736990045
+1240408112
+1111882347
+399261812
+2244375330
+1602012471
+912141067
+1830710117
+1340015527
+1567406936
+643677579
+2055472354
+819557021
+770847064
+338723885
+730373915
+818364243
+1856469323
+1196282919
+1589280011
+1970816376
+1734884585
+2139356181
+1734800064
+902092300
+1000228005
+1005397319
+1739544505
+1296449319
+759320334
+733835418
+430817036
+862588632
+337083331
+1146892667
+331681125
+1727328609
+833526010
+1372841208
+1791256853
+2098503617
+1476066812
+1106534491
+2269588375
+1764908204
+992903462
+1962380095
+1930972044
+1910002535
+744703663
+378223574
+1684151901
+2238279040
+2427906216
+2461340686
+1220653916
+1873597157
+1426518367
+700938423
+1092093371
+1219400088
+617341487
+2233007983
+1921740203
+2465020908
+681714775
+1445770833
+2515583030
+1068376614
+2099169797
+777626672
+1462970693
+2271210095
+1075897911
+894659316
+1188631296
+1999329869
+2479543901
+2576490680
+929714438
+1864382900
+1961606008
+1631866996
+2292048031
+2311220896
+1840535627
+1074753973
+847303451
+512810439
+1000862128
+1763495143
+1822824028
+2380337601
+977909730
+1700850927
+1491663654
+2459396136
+1068510445
+2347782874
+2316919822
+1988928160
+1422341309
+2047406045
+896998975
+1204758308
+1055292188
+2070261038
+1185381558
+584629736
+1397224527
+2212812703
+697064532
+427021743
+1521498782
+1918105312
+2305959364
+2317070012
+1544506428
+2394948578
+2399789373
+2423917413
+955268974
+1566323992
+1385782847
+2403855141
+1860302572
+2623796332
+2257560818
+778070755
+1761204082
+2554995183
+2029466541
+1105246882
+1922250969
+1274097165
+1794429653
+962667855
+1950360169
+961488742
+675055038
+2269364156
+1918149132
+1596341225
+1449514628
+2319740649
+645181247
+1526064014
+2532704137
+2719363555
+1470524400
+1879739046
+889579427
+2347651601
+2113527897
+2472401280
+2286628409
+1126370210
+1239790796
+1342682654
+2149079289
+1682494794
+2264801373
+1548328523
+1704463823
+722920765
+316155630
+1477916383
+2165336613
+1144813537
+1715928224
+2258304828
+1588161375
+1111482449
+1993648422
+2621681421
+2794343549
+1296161757
+503011116
+291796595
+2088433916
+2286988618
+2017170705
+770487794
+782581318
+1906340499
+772340260
+1565167492
+990441186
+1183974590
+1745558956
+568506725
+513006469
+2157622652
+2162121611
+599383917
+1984215367
+2191903075
+605455854
+1434674443
+1260076872
+591067666
+1019688217
+2347059212
+800361559
+805116519
+1820301461
+2047235099
+1474841357
+1493591766
+2188299688
+1913351485
+1560963124
+1170764024
+1001228137
+2379030606
+2235539020
+834681930
+721965868
+2025669856
+1830216606
+906769433
+1663789196
+1489162668
+2314703896
+1460273942
+799822619
+500625364
+1324566428
+1728272447
+713302840
+1013928492
+1446233486
+2210239384
+2180069630
+1292178034
+2384366776
+679123903
+586854776
+1178641948
+2320001266
+2337583819
+1361613696
+1297114303
+1265431260
+1557928964
+1651277282
+430694661
+366517495
+565930090
+265073035
+549090596
+383831915
+2141927042
+831246210
+1820526363
+2574323737
+2158547448
+2445738313
+2231476345
+2315998965
+2453275846
+2292850760
+2047107444
+636931736
+1726260080
+678647440
+307217003
+2206164032
+1199127108
+1487959761
+706070465
+1521626919
+1085502541
+2077418347
+630099107
+720819652
+1105259038
+1752650036
+1203144263
+833962649
+1468957755
+1408331636
+1043541089
+607821064
+2172824672
+1943459383
+1385583737
+1763516717
+743440367
+793147854
+1333910307
+1798901388
+617837221
+940888946
+1874859391
+1635889478
+2209068819
+1646898547
+1915865781
+1314674415
+2030151322
+1736345114
+796968442
+742711819
+1267177548
+955439662
+278726598
+1550628711
+1659180165
+1092596457
+1873982972
+1342794869
+2312746623
+1865638697
+893628449
+1412439051
+734276018
+1884775998
+2369885233
+1093726825
+1547121317
+2252945626
+2460368839
+1198321920
+1230077712
+2252187138
+1874551792
+1043829055
+1326307168
+1944472818
+1100425345
+1039742279
+1479918018
+1043161844
+287190812
+1537728699
+746497077
+519854187
+1623154621
+2317012147
+2670729715
+1230509471
+2357738366
+1899778875
+1669809364
+1406293854
+1909533901
+1676371793
+2423721546
+1095924725
+580103698
+681327149
+1791309722
+2060961799
+1662132185
+2029470483
+862634632
+2078791660
+1723525488
+1382593809
+2218362550
+950206279
+483890845
+617657473
+177745989
+79337962
+1390190151
+1553021618
+2480624049
+1633827011
+590251531
+1160963522
+894258038
+1034254575
+582506922
+818777093
+684141196
+1472804515
+500721957
+2211087829
+2227144554
+1647402253
+430550705
+1829502936
+1507815526
+1335261471
+1990344441
+1809556986
+1539263486
+674856851
+372635550
+734294164
+264604314
+142710170
+866425314
+1370268393
+2213627586
+2355390337
+1781018194
+616005594
+1252043764
+860858002
+1214700792
+811176416
+361131489
+1414285509
+2468441830
+1673892968
+1366189729
+1760117093
+1758116489
+2470380914
+2770407395
+711860653
+1874636039
+1127807075
+1020432622
+2296072531
+2336098168
+2063439686
+2486486433
+1929987913
+656373871
+1232183806
+1032318157
+377676678
+2016421248
+2505222795
+838373879
+2079447973
+1298955231
+2169777404
+1972291806
+2230854824
+867289237
+997527640
+279414453
+1498974472
+2200037040
+1316053465
+2175602406
+1061915775
+1981062396
+2192337808
+1551049055
+2075038498
+1021716473
+1544369290
+1738634349
+984528574
+2197738300
+700210698
+1832527327
+1870401090
+732249396
+740130055
+1642161163
+658976279
+1913896436
+2070781650
+2359657104
+1610737562
+1209358449
+1966311439
+785886478
+2116386307
+2538532804
+1403649244
+1694039536
+829773831
+1077945721
+996752632
+1292143044
+916228165
+2013561256
+1532989724
+1555235942
+1548742150
+1670824094
+1587310501
+1374063879
+1536301919
+2440354887
+1587464103
+1026899321
+1654751770
+2453764801
+979970141
+2158352434
+1536966534
+2380824408
+2395782056
+926372692
+1039476382
+1593160609
+660701941
+753702002
+1407419927
+1703190747
+602279837
+362797324
+677374220
+489614465
+279130192
+1484019228
+1593005014
+519294597
+2243120497
+1476448893
+589091830
+1148940283
+1498226860
+2468011772
+2552103508
+2144262851
+1492917654
+2041321360
+2384729929
+1578297709
+2335926135
+1498809230
+1988828908
+911912841
+1498392766
+1759094805
+2509690052
+909711355
+1839522939
+1666948908
+456080641
+329239027
+1972212557
+1113638792
+1119752951
+658428332
+2101658364
+1317298014
+678751494
+793399614
+1234357489
+1964341595
+1667494897
+1823765939
+1453709770
+1929304663
+591327663
+557004648
+358332336
+155871124
+669912678
+1422360996
+1270486863
+1458875466
+1150485898
+949752252
+1546162908
+1093046635
+411201632
+887699569
+273205746
+856379020
+1449869174
+1914940249
+2102702007
+931491820
+2100377130
+2676045758
+1138354878
+1901646837
+2531570974
+1799663521
+2391381872
+1587417611
+678486099
+1938048030
+724697612
+1827135658
+1423189228
+397753660
+876120955
+1168216082
+530159340
+1745797943
+1787876162
+761430091
+1200668347
+2035785923
+2169737525
+1337834143
+824992003
+412184950
+2152141256
+2187181679
+2329413975
+1381789618
+607997774
+321735032
+1781732617
+488326415
+1665137416
+1672232410
+1288032493
+2336982187
+2588290296
+1968734394
+1366173294
+2493784183
+2327817122
+1011046817
+1752924333
+2024187196
--- a/gva_res_100000.csv
+++ b/gva_res_100000.csv