CouenneExprSum.hpp

Go to the documentation of this file.

/* $Id: CouenneExprSum.hpp 615 2011-06-08 20:36:24Z pbelotti $
 *
 * Name:    exprSum.hpp
 * Author:  Pietro Belotti
 * Purpose: definition of sum expressions
 *
 * (C) Carnegie-Mellon University, 2006-10.
 * This file is licensed under the Eclipse Public License (EPL)
 */

#ifndef COUENNE_EXPRSUM_H
#define COUENNE_EXPRSUM_H

#include <vector>

#include "CouenneExprOp.hpp"

namespace Couenne {


class exprSum: public exprOp {

 public:

  exprSum  (expression ** = NULL, int = 0);

  exprSum (expression *, expression *);

  virtual ~exprSum () {}
 
  virtual expression *clone (Domain *d = NULL) const
    {return new exprSum (clonearglist (d), nargs_);}

  std::string printOp () const
    {return "+";}

  virtual CouNumber operator () ();

  virtual expression *differentiate (int index); 

  virtual expression *simplify ();

  virtual int Linearity ();

  virtual void getBounds (expression *&, expression *&);

  virtual void getBounds (CouNumber &, CouNumber &);

  virtual exprAux *standardize (CouenneProblem *p, bool addAux = true);

  virtual void generateCuts (expression *, //const OsiSolverInterface &, 
                             OsiCuts &, const CouenneCutGenerator *,
                             t_chg_bounds * = NULL, int = -1, 
                             CouNumber = -COUENNE_INFINITY, 
                             CouNumber =  COUENNE_INFINITY);

  virtual enum expr_type code () 
    {return COU_EXPRSUM;}

  virtual bool impliedBound (int, CouNumber *, CouNumber *, t_chg_bounds *, enum auxSign = expression::AUX_EQ);

  exprAux *createQuadratic (CouenneProblem *);

protected:

  int impliedBoundSum (CouNumber wl, 
                       CouNumber wu, 
                       std::vector <CouNumber> &xl,
                       std::vector <CouNumber> &xu,
                       std::vector <std::pair <int, CouNumber> > &nl,
                       std::vector <std::pair <int, CouNumber> > &nu);
};



inline CouNumber exprSum::operator () () {

  CouNumber ret = 0;

  expression **al = arglist_;

  for (int n = nargs_; n--;)
    ret += (**al++) ();

  return ret;
}

}

#endif