Changeset e59735

Timestamp:

May 9, 2017, 11:06:46 AM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_goodresults, ForceAnnealing_tocheck

Children:

d2a077a

Parents:

ba7a54

Message:

tempcommit: Fixing RandomPerturbationAction.

File:

• src/Actions/AtomAction/RandomPerturbationAction.cpp (modified) (2 diffs)

src/Actions/AtomAction/RandomPerturbationAction.cpp

@@ -78 +78 @@
   const double rng_max = random.max();
   const double range = (rng_max - rng_min);
-  const double offset = (rng_max - rng_min);
+  const double offset = rng_min;
 
   std::vector<AtomicInfo> NewWalkers;
@@ -85 +85 @@
     double noisecomponents[NDIM];
     for (size_t i=0;i<NDIM;++i)
-      noisecomponents[i] = params.noise.get()*(random()*2./range - offset);
+      noisecomponents[i] = params.noise.get()*((random()-offset)*2./range - 1.);
     const Vector noisevector(noisecomponents);
     LOG(3, "DEBUG: Changing atom " << *first << "'s position by " << noisevector << ".");