makefile + gaissoa

Makefile

+CC=g++
+LFLAGS=-lm
+EXEC= r_gaussian \
+			r_monte_carlo \
+			r_PM
+
+all: $(EXEC)
+
+r_%: %.cpp
+	$(CC) $(LFLAGS) $< -o $@

fcts.h

@@ -20,8 +20,8 @@ long uniform() {
    return seed ;
 }
 void gaussian(double U0, double U1, double *res){
-  res[0] = sqrt(-2 * log(U0));// * cos(2 * M_PI * U1);
-  res[1] = sqrt(-2 * log(U0));// * sin(2 * M_PI * U1);
+  res[0] = sqrt(-2 * log(U0)) * cos(2 * M_PI * U1);
+  res[1] = sqrt(-2 * log(U0)) * sin(2 * M_PI * U1);
 }
 
 void store_freqs(int freqs[], int nCuts){

gaussian.cpp

@@ -22,12 +22,12 @@ int main(int argc, char **argv) {
      freqs[i] = 0;
 
   for (int i = 0; i < nRNG; i++){
-    U0 = uniform();
-    U1 = uniform();
+    U0 = uniform()/(double)0x7FFFFFFF;
+    U1 = uniform()/(double)0x7FFFFFFF;
     gaussian(U0, U1, Y);
-    printf("%lf, %lf\n", Y[0], sqrt(2 * log(U1)) * cos(2 * M_PI * U1));
+    printf("%lf, %lf\n", Y[0], Y[1]);
   }
-  // store_res_double(Y)
+  store_res_double(Y)
   rescsv.close();
   return 0;
 }