eoRealAtomXover.h

//-----------------------------------------------------------------------------

#ifndef _eoRealAtomXover_H
#define _eoRealAtomXover_H

#include <utils/eoRNG.h>

#include <eoOp.h>

class eoDoubleExchange: public eoBinOp<double>
{
public:
  eoDoubleExchange() {}

  virtual std::string className() const { return "eoDoubleExchange"; }

  bool operator()(double& r1, const double& r2)
  {
    if (eo::rng.flip())
      if (r1 != r2)        // if r1 == r2 you must return false
        {
          r1 = r2;
          return true;
        }
    return false;
  }

};

class eoDoubleIntermediate: public eoBinOp<double>
{
public:
  eoDoubleIntermediate() {}

  virtual std::string className() const { return "eoDoubleIntermediate"; }

  bool operator()(double& r1, const double& r2)
  {
    double alpha = eo::rng.uniform();
    r1 = alpha * r2 + (1-alpha) * r1;
    return true;
  }

};

#endif