nonlinear_inst.h

#ifndef NONLINEAR_INST_H
#define NONLINEAR_INST_H

#include<scil/scil.h>
#include<scil/constraints/cut.h>
#include<scil/core/monomial.h>
#include<scil/core/qelement.h>
#include<scil/core/qmonomial.h>
#include<scil/core/qboolfunction.h>
#include<scil/variable.h>
#include<boost/graph/adjacency_list.hpp>
#include<boost/graph/graph_traits.hpp>
#include<boost/graph/iteration_macros.hpp>

#include<vector>
#include<list>
#include<map>

namespace SCIL{

  class nonlinear_inst : public sym_constraint {
    //  private:
     protected:
         typedef adjacency_list<vecS, vecS, undirectedS> Graph;
         typedef graph_traits<Graph> GraphTraits;
         typedef GraphTraits::vertex_descriptor vertex_descriptor;
         typedef GraphTraits::edge_descriptor edge_descriptor;

         std::list<qelement<Graph>*> qm; //the elementss comprising the quadratic instance
         Graph G; //the separation graph
         row_map<edge_descriptor> VM; //the mapping between edges and variables
         //assign nodes and edges to qelementss
         virtual void construct_separation_graph(ILP_Problem& IP);
         //output the edge of the separation graph defined by a and b
         //(if QUADINST_DEBUG)
         void printEdge(qelement<Graph>* a, qelement<Graph>* b);
         //output the quadratic decomposition (if QUADINST_DEBUG)
         void printDecomposition();

   //   protected:
         SCIL::CUT<Graph>* scc; //the symbolic cut constraint for separation
         //check if the values defined by mapping the fractional solution on the edges are valid
         bool checkBounds(subproblem& S);

      public:
         //trivial default constructor
         nonlinear_inst()
            : scc(NULL) {}

         //destructor
         ~nonlinear_inst();

         //initialization, derived from sym_constraint
         virtual void init(subproblem& S);

         //separation routine, derived from sym_constraint
         virtual status standard_separation(subproblem& S);

         //feasibility test, derived from sym_constraint
         virtual status feasible(solution& S){};
   };
}
#endif

---