CouenneDomain.hpp

Go to the documentation of this file.

/* $Id: CouenneDomain.hpp 927 2012-11-28 15:25:47Z stefan $
 *
 * Name:    domain.hpp
 * Author:  Pietro Belotti
 * Purpose: class for point and bounding box
 *
 * (C) Carnegie-Mellon University, 2008-09.
 * This file is licensed under the Eclipse Public License (EPL)
 */

#ifndef COUENNE_DOMAIN_HPP
#define COUENNE_DOMAIN_HPP

#include <stdlib.h>
#include <stack>

#include "CouenneTypes.hpp"

namespace Osi {

  class OsiSolverInterface;
  class OsiCuts;
}


namespace Couenne {


class DomainPoint {

  friend class Domain;

protected:

  int dimension_; 

  CouNumber *x_;  
  CouNumber *lb_; 
  CouNumber *ub_; 

  bool copied_;   

  bool isNlp_;    
public:

  DomainPoint (int dim, 
               CouNumber *x, 
               CouNumber *lb,
               CouNumber *ub,
               bool copy = true);

  DomainPoint (int dim = 0, 
               const CouNumber *x   = NULL, 
               const CouNumber *lb  = NULL, 
               const CouNumber *ub  = NULL,
               bool copy = true);

  ~DomainPoint () {
    if (copied_) {
      if (x_)  free (x_);
      if (lb_) free (lb_);
      if (ub_) free (ub_);
    }
  }

  DomainPoint (const DomainPoint &src);

  void resize (int newdim);

  int size () const {return dimension_;}

  inline int Dimension () {return dimension_;}

  inline CouNumber &x  (register int index) {return x_  [index];} 
  inline CouNumber &lb (register int index) {return lb_ [index];} 
  inline CouNumber &ub (register int index) {return ub_ [index];} 

  inline CouNumber *x  () {return x_ ;} 
  inline CouNumber *lb () {return lb_;} 
  inline CouNumber *ub () {return ub_;} 

  DomainPoint &operator= (const DomainPoint &src);

  bool &isNlp () 
  {return isNlp_;}
};



class Domain {

protected:

  DomainPoint *point_;                  
  std::stack <DomainPoint *> domStack_; 

public:

  Domain (): point_ (NULL) {}

  Domain (const Domain &src) {
    point_ = new DomainPoint (*(src.point_));
    // TODO -- not important, discard previous points when copying problem
    /*for (std::stack <DomainPoint *>::iterator i = src.domStack_.begin ();
         i != src.domStack_.end (); ++i)
         domStack_.push (new DomainPoint (**i));*/
  } 

  ~Domain ();

  void push (int dim, 
             CouNumber *x, 
             CouNumber *lb, 
             CouNumber *ub, 
             bool copy = true);

  void push (int dim, 
             const CouNumber *x, 
             const CouNumber *lb, 
             const CouNumber *ub,
             bool copy = true);

  void push (const OsiSolverInterface *si,
             OsiCuts *cs = NULL, 
             bool copy = true);

  void push (const DomainPoint &dp, bool copy = true);

  void pop ();

  inline DomainPoint *current ()   {return point_;}                          

  inline CouNumber &x  (register int index) {return point_ -> x  (index);}   
  inline CouNumber &lb (register int index) {return point_ -> lb (index);}   
  inline CouNumber &ub (register int index) {return point_ -> ub (index);}   

  inline CouNumber *x  () {return point_ -> x  ();}   
  inline CouNumber *lb () {return point_ -> lb ();}   
  inline CouNumber *ub () {return point_ -> ub ();}   
};

}

#endif