eoIntBounds.h

// -*- mode: c++; c-indent-level: 4; c++-member-init-indent: 8; comment-column: 35; -*-

//-----------------------------------------------------------------------------
// eoIntBounds.h
// (c) Marc Schoenauer 2001, Maarten Keijzer 2000, GeNeura Team, 1998
/*
    This library is free software; you can redistribute it and/or
    modify it under the terms of the GNU Lesser General Public
    License as published by the Free Software Foundation; either
    version 2 of the License, or (at your option) any later version.

    This library is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public
    License along with this library; if not, write to the Free Software
    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

    Contact: todos@geneura.ugr.es, http://geneura.ugr.es
             Marc.Schoenauer@polytechnique.fr
             mak@dhi.dk
 */
//-----------------------------------------------------------------------------

#ifndef _eoIntBounds_h
#define _eoIntBounds_h

#include <stdexcept>               // std::exceptions!
#include <utils/eoRNG.h>

class eoIntBounds : public eoPersistent
{
public:
  virtual ~eoIntBounds(){}

  virtual bool isBounded(void) const  = 0;

  virtual bool hasNoBoundAtAll(void) const  = 0;

  virtual bool isMinBounded(void) const = 0;

  virtual bool isMaxBounded(void) const = 0;

  virtual bool isInBounds(double)  const = 0;

  virtual void foldsInBounds(double &)  const = 0;

  virtual void foldsInBounds(long int & i)  const
  {
    double r = double(i);
    foldsInBounds(r);
    i = (long int)(r);
  }

  virtual void truncate(double &)  const = 0;

  virtual void truncate(long int & i)  const
  {
    double r = double(i);
    truncate(r);
    i = (long int)(r);
  }

  virtual long int minimum()  const = 0 ;
  virtual long int maximum()  const = 0 ;
  virtual long int range()  const = 0;

  virtual double uniform(eoRng & _rng = eo::rng)  const = 0;
  virtual long int random(eoRng & _rng = eo::rng)  const = 0;

  virtual eoIntBounds * dup()  const = 0;
};

class eoIntNoBounds : public eoIntBounds
{
public:
  virtual ~eoIntNoBounds(){}

  virtual bool isBounded(void)  const {return false;}
  virtual bool hasNoBoundAtAll(void) const  {return true;}
  virtual bool isMinBounded(void)  const {return false;}
  virtual bool isMaxBounded(void)  const {return false;}
  virtual void foldsInBounds(double &)  const {return;}
  virtual void truncate(double &)  const {return;}
  virtual bool isInBounds(double)  const {return true;}

  virtual long int minimum() const
  {
    throw std::logic_error("Trying to get minimum of unbounded eoIntBounds");
  }
  virtual long int maximum() const
  {
    throw std::logic_error("Trying to get maximum of unbounded eoIntBounds");
  }
  virtual long int range() const
  {
    throw std::logic_error("Trying to get range of unbounded eoIntBounds");
  }

  virtual double uniform(eoRng & _rng = eo::rng) const
  {
    (void)_rng;

    throw std::logic_error("Trying to generate uniform values in unbounded eoIntBounds");
  }

  virtual long int random(eoRng & _rng = eo::rng) const
  {
    (void)_rng;

    throw std::logic_error("Trying to generate uniform values in unbounded eoIntBounds");
  }

  // methods from eoPersistent
  virtual void readFrom(std::istream& _is)
  {
    (void)_is;

    throw std::runtime_error("Should not use eoIntBounds::readFrom");
  }

  virtual void printOn(std::ostream& _os) const
  {
    _os << "[-inf,+inf]";
  }

  virtual eoIntBounds * dup() const
  {
    return new eoIntNoBounds(*this);
  }

};

extern eoIntNoBounds eoDummyIntNoBounds;

class eoIntInterval : public eoIntBounds
{
public :
  virtual ~eoIntInterval(){}

  eoIntInterval(long int _min=0, long int _max=1) :
    repMinimum(_min), repMaximum(_max), repRange(_max-_min)
  {
    if (repRange<=0)
      throw std::logic_error("Void range in eoIntBounds");
  }

  // accessors
  virtual long int minimum() const { return repMinimum; }
  virtual long int maximum() const { return repMaximum; }
  virtual long int range()  const { return repRange; }

  // description
  virtual bool isBounded(void)  const {return true;}
  virtual bool hasNoBoundAtAll(void)  const {return false;}
  virtual bool isMinBounded(void)  const {return true;}
  virtual bool isMaxBounded(void)  const {return true;}

  virtual double uniform(eoRng & _rng = eo::rng) const
  {
    return repMinimum + _rng.uniform(repRange);
  }

  virtual long int random(eoRng & _rng = eo::rng) const
  {
    return repMinimum + _rng.random(repRange);
  }

  // says if a given double is within the bounds
  virtual bool isInBounds(double _r) const
  {
    if (_r < repMinimum)
      return false;
    if (_r > repMaximum)
      return false;
    return true;
  }

  // folds a value into bounds
  virtual void foldsInBounds(double &  _r) const
  {
    long iloc;
    double dlargloc = 2 * range() ;

    if (fabs(_r) > 1.0E9)               // iloc too large!
      {
        _r = uniform();
        return;
      }

    if ( (_r > maximum()) )
      {
        iloc = (long) ( (_r-minimum()) / dlargloc ) ;
        _r -= dlargloc * iloc ;
        if ( _r > maximum() )
          _r = 2*maximum() - _r ;
      }

    if (_r < minimum())
      {
        iloc = (long) ( (maximum()-_r) / dlargloc ) ;
        _r += dlargloc * iloc ;
        if (_r < minimum())
          _r = 2*minimum() - _r ;
      }
  }

  // truncates to the bounds
  virtual void truncate(double & _r) const
  {
    if (_r < repMinimum)
      _r = repMinimum;
    else if (_r > repMaximum)
      _r = repMaximum;
    return;
  }

  // methods from eoPersistent
  virtual void readFrom(std::istream& _is)
  {
    (void)_is;

    throw std::runtime_error("Should not use eoIntInterval::readFrom");
  }

  virtual void printOn(std::ostream& _os) const
  {
    _os << "[" << repMinimum << "," << repMaximum << "]";
  }

  virtual eoIntBounds * dup() const
  {
    return new eoIntInterval(*this);
  }

private :
  long int repMinimum;
  long int repMaximum;
  long int repRange;                       // to minimize operations ???
};

class eoIntBelowBound : public eoIntBounds
{
public :
  virtual ~eoIntBelowBound(){}
  eoIntBelowBound(long int _min=0) :
    repMinimum(_min)
  {}

  // accessors
  virtual long int minimum() const { return repMinimum; }

  virtual long int maximum() const
  {
    throw std::logic_error("Trying to get maximum of eoIntBelowBound");
  }
  virtual long int range() const
  {
    throw std::logic_error("Trying to get range of eoIntBelowBound");
  }

  virtual double uniform(eoRng & _rng = eo::rng) const
  {
    (void)_rng;

    throw std::logic_error("Trying to generate uniform values in eoIntBelowBound");
  }

  virtual long int random(eoRng & _rng = eo::rng) const
  {
    (void)_rng;

    throw std::logic_error("Trying to generate uniform values in eoIntBelowBound");
  }

  // description
  virtual bool isBounded(void)  const {return false;}
  virtual bool hasNoBoundAtAll(void)  const {return false;}
  virtual bool isMinBounded(void)  const {return true;}
  virtual bool isMaxBounded(void) const  {return false;}

  // says if a given double is within the bounds
  virtual bool isInBounds(double _r) const
  {
    if (_r < repMinimum)
      return false;
    return true;
  }

  // folds a value into bounds
  virtual void foldsInBounds(double &  _r) const
  {
    // easy as a pie: symmetry w.r.t. minimum
    if (_r < repMinimum)           // nothing to do otherwise
      _r = 2*repMinimum - _r;
    return ;
  }

  // truncates to the bounds
  virtual void truncate(double & _r) const
  {
    if (_r < repMinimum)
      _r = repMinimum;
    return;
  }

  // methods from eoPersistent
  virtual void readFrom(std::istream& _is)
  {
    (void)_is;

    throw std::runtime_error("Should not use eoIntBelowBound::readFrom");
  }

  virtual void printOn(std::ostream& _os) const
  {
    _os << "[" << repMinimum << ",+inf]";
  }

  virtual eoIntBounds * dup() const
  {
    return new eoIntBelowBound(*this);
  }

private :
  long int repMinimum;
};

class eoIntAboveBound : public eoIntBounds
{
public :
  virtual ~eoIntAboveBound(){}

  eoIntAboveBound(long int _max=0) :
    repMaximum(_max)
  {}

  // accessors
  virtual long int maximum() const  { return repMaximum; }

  virtual long int minimum() const
  {
    throw std::logic_error("Trying to get minimum of eoIntAboveBound");
  }
  virtual long int range() const
  {
    throw std::logic_error("Trying to get range of eoIntAboveBound");
  }

  virtual double uniform(eoRng & _rng = eo::rng) const
  {
    (void)_rng;

    throw std::logic_error("Trying to generate uniform values in eoIntAboveBound");
  }

  virtual long int random(eoRng & _rng = eo::rng) const
  {
    (void)_rng;

    throw std::logic_error("Trying to generate uniform values in eoIntAboveBound");
  }

  // description
  virtual bool isBounded(void)  const {return false;}
  virtual bool hasNoBoundAtAll(void)  const {return false;}
  virtual bool isMinBounded(void)  const {return false;}
  virtual bool isMaxBounded(void)  const {return true;}

  // says if a given double is within the bounds
  virtual bool isInBounds(double _r) const
  {
    if (_r > repMaximum)
      return false;
    return true;
  }

  // folds a value into bounds
  virtual void foldsInBounds(double &  _r) const
  {
    // easy as a pie: symmetry w.r.t. maximum
    if (_r > repMaximum)           // nothing to do otherwise
      _r = 2*repMaximum - _r;
    return ;
  }

  // truncates to the bounds
  virtual void truncate(double & _r) const
  {
    if (_r > repMaximum)
      _r = repMaximum;
    return;
  }

  // methods from eoPersistent
  virtual void readFrom(std::istream& _is)
  {
    (void)_is;

    throw std::runtime_error("Should not use eoIntAboveBound::readFrom");
  }

  virtual void printOn(std::ostream& _os) const
  {
    _os << "[-inf," << repMaximum << "]";
  }

  virtual eoIntBounds * dup() const
  {
    return new eoIntAboveBound(*this);
  }

private :
  long int repMaximum;
};


class eoGeneralIntBounds : public eoIntBounds
{
public:
  eoGeneralIntBounds(std::string _s = "[-infinity,+infinity]")
  {
    repBound = getBoundsFromString(_s);
  }

  eoGeneralIntBounds(const eoGeneralIntBounds & _b) : eoIntBounds(_b)
  {
    // replicate the embedded bound (I'm pretty sure there is another
    // way to do that !!!

    bool minBounded = _b.isMinBounded();
    bool maxBounded = _b.isMaxBounded();
    long int minimum, maximum;
    const eoIntBounds & bb = _b.theBounds();
    if (minBounded) minimum = bb.minimum();
    if (maxBounded) maximum = bb.maximum();

      if (minBounded && maxBounded)
        repBound = new eoIntInterval(minimum, maximum);
      else if (!minBounded && !maxBounded)      // no bound at all
        repBound = new eoIntNoBounds;
      else if (!minBounded && maxBounded)
        repBound = new eoIntAboveBound(maximum);
      else if (minBounded && !maxBounded)
        repBound = new eoIntBelowBound(minimum);
  }

  eoGeneralIntBounds& operator=(const eoGeneralIntBounds& _b)
  {
    // replicate the embedded bound (I'm pretty sure there is another
    // way to do that !!!

    bool minBounded = _b.isMinBounded();
    bool maxBounded = _b.isMaxBounded();
    long int minimum, maximum;
    const eoIntBounds & bb = _b.theBounds();
    if (minBounded) minimum = bb.minimum();
    if (maxBounded) maximum = bb.maximum();

    // first delete the embedded bounds if necessary
    if (repBound)
      delete repBound;
    // now reallocate
      if (minBounded && maxBounded)
        repBound = new eoIntInterval(minimum, maximum);
      else if (!minBounded && !maxBounded)      // no bound at all
        repBound = new eoIntNoBounds;
      else if (!minBounded && maxBounded)
        repBound = new eoIntAboveBound(maximum);
      else if (minBounded && !maxBounded)
        repBound = new eoIntBelowBound(minimum);
      return (*this);
  }


  ~eoGeneralIntBounds()
  {
    delete repBound;
  }


  virtual bool isBounded(void)  const {return repBound->isBounded();}

  virtual bool hasNoBoundAtAll(void)  const {return repBound->hasNoBoundAtAll();}

  virtual bool isMinBounded(void)  const {return repBound->isMinBounded();}

  virtual bool isMaxBounded(void) const {return repBound->isMaxBounded();}

  virtual bool isInBounds(double _x)  const {return repBound->isInBounds(_x);}

  virtual void foldsInBounds(double & _x) const {return repBound->foldsInBounds(_x);}

  virtual void truncate(double & _x)  const {return repBound->truncate(_x);}

  virtual long int minimum()  const {return repBound->minimum();}
  virtual long int maximum() const {return repBound->maximum();}
  virtual long int range()  const {return repBound->range();}

  virtual double uniform(eoRng & _rng = eo::rng)  const {(void)_rng; return repBound->uniform();}

  virtual long int random(eoRng & _rng = eo::rng)  const {(void)_rng; return repBound->random();}

  virtual eoIntBounds * dup() const  {return repBound->dup();}

  const eoIntBounds & theBounds()  const { return *repBound;}

  virtual void printOn(std::ostream& _os) const
  {
    repBound->printOn(_os);
  }

  virtual void readFrom(std::istream& _is)
  {
    std::string s;
    _is >> s;
    if (repBound)
      delete repBound;
    repBound = getBoundsFromString(s);
  }

private:
  // reading from a string
  eoIntBounds * getBoundsFromString(std::string);

  eoIntBounds * repBound;
};


#endif