Changes in src/Thermostats/Langevin.cpp [ad011c:a5028f]

File:

• src/Thermostats/Langevin.cpp (modified) (3 diffs)

src/Thermostats/Langevin.cpp

@@ -19 +19 @@
 #include "CodePatterns/Log.hpp"
 #include "ThermoStatContainer.hpp"
+
+#include "RandomNumbers/RandomNumberGeneratorFactory.hpp"
+#include "RandomNumbers/RandomNumberGenerator.hpp"
 
 Langevin::Langevin(double _TempFrequency,double _alpha) :
@@ -73 +76 @@
 double Langevin::doScaleAtoms(unsigned int step,double ActualTemp,ForwardIterator begin, ForwardIterator end){
   DoLog(2) && (Log() << Verbose(2) <<  "Applying Langevin thermostat..." << endl);
+  RandomNumberGenerator &random = RandomNumberGeneratorFactory::getInstance().makeRandomNumberGenerator();
+  const double rng_min = random.min();
+  const double rng_max = random.max();
   double ekin=0;
   for(ForwardIterator iter=begin;iter!=end;++iter){
@@ -79 +85 @@
     if ((*iter)->FixedIon == 0) { // even FixedIon moves, only not by other's forces
       // throw a dice to determine whether it gets hit by a heat bath particle
-      if (((((rand()/(double)RAND_MAX))*TempFrequency) < 1.)) {
+      if (((((random()/(rng_max-rng_min)))*TempFrequency) < 1.)) {
         DoLog(3) && (Log() << Verbose(3) << "Particle " << (**iter) << " was hit (sigma " << sigma << "): " << U.Norm() << " -> ");
         // pick three random numbers from a Boltzmann distribution around the desired temperature T for each momenta axis