33 #ifndef CBC_SYMMETRY_HPP
34 #define CBC_SYMMETRY_HPP
37 #include "nausparse.h"
52 #define NTY_BAD_DEPTH 10
58 #define COUENNE_HACKED_EPS 1.e-07
59 #define COUENNE_HACKED_EPS_SYMM 1e-8
60 #define COUENNE_HACKED_EXPRGROUP 8
79 void node(
int,
double,
double,
double,
int,
int);
88 inline void bounds(
register double a,
register double b){
lb = a;
ub = b;}
115 inline bool operator() (
register const char *a,
register const char *b)
const
116 {
return strcmp (a, b) < 0;}
156 int largestOrbit(
const double * lower,
const double * upper)
const;
157 void ChangeBounds (
const double * lower,
const double * upper,
158 int numberColumns,
bool justFixedAtOne)
const;
159 inline bool compare (
register Node &a,
register Node &b)
const;
189 CbcNauty (
int n,
const size_t * v,
const int * d,
const int * e);
218 std::vector<std::vector<int> > *
getOrbits()
const;
224 {
return stats_->errstatus;}
269 std::multimap<int,int>::iterator
it;
271 std::pair<std::multimap<int,int>::iterator,
272 std::multimap<int,int>::iterator>
ret;
295 int numberExtra,
const int * extraToZero);
315 double * lower,
double * upper,
316 int branchState)
const ;
328 virtual void print();
void setupSymmetry(const OsiSolverInterface &solver)
empty if no NTY, symmetry data structure setup otherwise
void color_node(int ix, int color)
virtual ~CbcOrbitalBranchingObject()
void addElement(int ix, int jx)
virtual CbcBranchingObject * clone() const
Clone.
bool operator()(register const char *a, register const char *b) const
std::multimap< int, int > constr_rhs
void ChangeBounds(const double *lower, const double *upper, int numberColumns, bool justFixedAtOne) const
#define COUENNE_HACKED_EPS_SYMM
void color_vertex(register int k)
void deleteElement(int ix, int jx)
std::vector< std::vector< int > > * getOrbits() const
Returns the orbits in a "convenient" form.
std::multimap< int, int >::iterator it
virtual double branch()
Does next branch and updates state.
Branching object for Orbital branching.
~CbcSymmetry()
Destructor.
Abstract Base Class for describing an interface to a solver.
void Compute_Symmetry() const
bool compare(register Node &a, register Node &b) const
Class to deal with symmetry.
double getGroupSize() const
int * fixToZero_
Fix to zero.
void setWriteAutoms(const std::string &afilename)
Methods to classify orbits.
int way() const
Get the state of the branching object.
Abstract branching object base class Now just difference with OsiBranchingObject. ...
void bounds(register double a, register double b)
int numberUsefulOrbits() const
virtual void previousBranch()
Reset every information so that the branching object appears to point to the previous child...
void getVstat(double *v, int nv)
virtual CbcRangeCompare compareBranchingObject(const CbcBranchingObject *brObj, const bool replaceIfOverlap=false)
Compare the this with brObj.
void node(int, double, double, double, int, int)
int column_
Column to go to 1.
void Print_Orbits() const
CbcOrbitalBranchingObject & operator=(const CbcOrbitalBranchingObject &rhs)
CbcModel * model() const
Return model.
virtual double branch()=0
Execute the actions required to branch, as specified by the current state of the branching object...
virtual CbcBranchObjType type() const
Return the type (an integer identifier) of this.
std::vector< Node > node_info_
virtual void print() const
Print something about branch - only if log level high.
void insertRHS(int rhs, int cons)
std::pair< std::multimap< int, int >::iterator, std::multimap< int, int >::iterator > ret
int numberOther_
Number (without column) going to zero on down branch.
int numberExtra_
Number extra.
int statsOrbits(CbcModel *model, int type) const
CbcNauty & operator=(const CbcNauty &rhs)
Assignment operator.
bool operator()(register const Node &a, register const Node &b)
virtual int compareOriginalObject(const CbcBranchingObject *brObj) const
Compare the original object of this with the original object of brObj.
CbcOrbitalBranchingObject()
virtual void print()
Print something about branch - only if log level high.
virtual void previousBranch()
Reset every information so that the branching object appears to point to the previous child...
virtual void fix(OsiSolverInterface *solver, double *lower, double *upper, int branchState) const
Update bounds in solver as in 'branch' and update given bounds.
CbcSymmetry & operator=(const CbcSymmetry &rhs)
Assignment operator.
int orbitalFixing(OsiSolverInterface *solver)
Fixes variables using orbits (returns number fixed)
Abstract base class for `objects'.
int largestOrbit(const double *lower, const double *upper) const
Simple Branch and bound class.
int getNumGenerators() const
std::vector< int > * Find_Orbit(int) const
CbcSymmetry()
Default constructor.
CbcNauty()
Default constructor.
int numberUsefulObjects() const
bool operator()(register const Node &a, register const Node &b)