CouenneExprSin.hpp

Go to the documentation of this file.

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

#ifndef COUENNE_EXPRSIN_HPP
#define COUENNE_EXPRSIN_HPP

#include <math.h>
#include <assert.h>

#include "CouenneExprUnary.hpp"
#include "CouenneExprConst.hpp"

namespace Couenne {

enum cou_trig {COU_SINE, COU_COSINE};


inline CouNumber modulo (register CouNumber a, register CouNumber b)
  {return a - b * floor (a/b);}


CouNumber trigSelBranch (const CouenneObject *obj, 
                         const OsiBranchingInformation *info,
                         expression * &var,
                         double * &brpts,
                         double * &brDist, // distance of current LP
                                           // point to new convexifications
                         int &way,
                         enum cou_trig type);


bool trigImpliedBound (enum cou_trig, int, int, CouNumber *, CouNumber *, t_chg_bounds *);


class exprSin: public exprUnary {

 public:

  exprSin (expression *al): 
    exprUnary (al) {} //< non-leaf expression, with argument list

  expression *clone (Domain *d = NULL) const
  {return new exprSin (argument_ -> clone (d));}

  inline unary_function F () 
  {return sin;}

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

  inline CouNumber gradientNorm (const double *x) {
    return (argument_ -> Index () < 0) ? 
      0. : fabs (cos (x [argument_ -> Index ()]));
  }

  expression *differentiate (int index); 

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

  void getBounds (CouNumber &lb, CouNumber &ub);

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

  virtual enum expr_type code () 
  {return COU_EXPRSIN;}

  bool impliedBound (int index, CouNumber *l, CouNumber *u, t_chg_bounds *chg, enum auxSign = expression::AUX_EQ) {

    bool impl = trigImpliedBound (COU_SINE, index, argument_ -> Index (), l, u, chg);

    if (impl && argument_ -> isInteger ()) {

      int ind = argument_ -> Index ();
      assert (ind >= 0);
      l [ind] = ceil  (l [ind] - COUENNE_EPS);
      u [ind] = floor (u [ind] + COUENNE_EPS);
    }

    return impl;
  }

  virtual CouNumber selectBranch (const CouenneObject *obj, 
                                  const OsiBranchingInformation *info,
                                  expression * &var, 
                                  double * &brpts, 
                                  double * &brDist, // distance of current LP
                                                    // point to new convexifications
                                  int &way)

  {return trigSelBranch (obj, info, var, brpts, brDist, way, COU_SINE);}

  virtual void closestFeasible (expression *varind, expression *vardep,
                                CouNumber& left, CouNumber& right) const;

  virtual bool isCuttable (CouenneProblem *problem, int index) const
  {return false;}
};

}

#endif