Changes in src/LinearAlgebra/RealSpaceMatrix.cpp [66fd49:a5028f]

File:

• src/LinearAlgebra/RealSpaceMatrix.cpp (modified) (2 diffs)

src/LinearAlgebra/RealSpaceMatrix.cpp

@@ -28 +28 @@
 #include "LinearAlgebra/Vector.hpp"
 #include "LinearAlgebra/VectorContent.hpp"
+#include "RandomNumbers/RandomNumberGeneratorFactory.hpp"
+#include "RandomNumbers/RandomNumberGenerator.hpp"
 
 #include <gsl/gsl_blas.h>
@@ -125 +127 @@
 {
   double phi[NDIM];
+  RandomNumberGenerator &random = RandomNumberGeneratorFactory::getInstance().makeRandomNumberGenerator();
+  const double rng_min = random.min();
+  const double rng_max = random.max();
+
 
   for (int i=0;i<NDIM;i++) {
-    phi[i] = rand()/(RAND_MAX/(2.*M_PI));
+    phi[i] = (random()/(rng_max-rng_min))*(2.*M_PI);
+    std::cout << "Random angle is " << phi[i] << std::endl;
   }