OSDecompSolver.h

Go to the documentation of this file.

/* $Id: OSDecompSolver.h 3038 2009-11-07 11:43:44Z Gassmann $ */
#ifndef OSDECOMPSOLVER_H
#define OSDECOMPSOLVER_H

// --------------------------------------------------------------------- //
#include "OSInstance.h"
#include "OSResult.h"
#include "OSDataStructures.h"
#include "OSErrorClass.h"
#include "OSOption.h"
#include "OSCoinSolver.h"
#include "OSDecompParam.h"
#include <vector>
#include <string>
#include <map>




// --------------------------------------------------------------------- //
// --------------------------------------------------------------------- //
class OSDecompSolver{
public:
        

        OSInstance *m_osinstanceMaster;
        
        int m_multiCommodCutLimit;
        int m_numMultCuts;
        
        
        OSDecompParam m_osDecompParam;
        
        double m_bestIPValue;
        double m_bestLPValue;
        double m_rootLPValue;
        
        
        
        //the transformation matrix
        int* m_thetaPnt;
        int* m_thetaIndex;
        int m_numThetaVar;
        int m_numThetaNonz;
        
        
        //arrays for the added constraints
        //for now the added constraints are 
        //tour breaking and variable branching
        //constraints
        //
        int* m_pntBmatrix;
        //m_Bmatrix is the matrix of Xij indexes
        int* m_BmatrixIdx;  
        //m_Bmatrix is the matrix of Xij coefficients
        double* m_BmatrixVal;  
        
        
        std::set<std::pair<int, double> > intVarSet;
        
        int m_numHubs;
        
        int m_numNodes;
        
        
        
        //arrays for the coupling constraint matrix
        //we store indexes since values are 1.0
        //also the RHS is 1.0
        //the matrix has m_numNodes - m_numHubs rows
        int* m_pntAmatrix;
        //m_Amatrix holds the column indexes for each row
        int* m_Amatrix;
        //
        
        int m_numBmatrixCon;   
        int m_numBmatrixNonz;
        
        
        int m_maxBmatrixNonz;
        //
        //end arrays for added constaints
        
        
        int m_maxBmatrixCon;
        
        
        int m_maxMasterColumns;
        
        int m_maxMasterRows;
        
        
        
        //these variable names are in the original variable space
        std::string* m_variableNames;
        
        
        OSOption *m_osoption;
        

        virtual OSInstance* getInitialRestrictedMaster( ) = 0;
        
        
        virtual void getCutsTheta(const  double* thetaVar, const int numThetaVar,
                                          int &numNewRows, int*  &numNonz, int** &colIdx,
                                          double** &values, double* &rowLB, double* &rowUB) = 0 ;       
        

        virtual void getCutsMultiCommod(const  double* thetaVar, const int numThetaVar,
                        int &numNewRows, int*  &numNonz, int** &colIdx,
                        double** &values, double* &rowLB, double* &rowUB) = 0;  
        
        
        
        
        virtual void getColumns(const  double* yA, const int numARows,
                                        const  double* yB, const int numBRows,
                                        int &numNewColumns, int* &numNonz, double* &cost, 
                                        int** &rowIdx, double** &values, double &lowerBound)  = 0;
        
        
        
        virtual void getBranchingCut(const double* thetaVar, const int numThetaVar,
                                const std::map<int, int> &varConMap, int &varIdx,  int &numNonz, 
                                int* &indexes,  double* &values) = 0 ;  
        
        
        
        virtual void getBranchingCut(const int* thetaIdx, const double* theta, 
                                const int numThetaVar, const std::map<int, int> &varConMap, 
                                int &varIdx,  int &numNonz, int* &indexes, double* &values) = 0 ;
        

        
        
        
        //this method gets called when we are done
        virtual void pauHana(std::vector<int> &m_zOptIndexes, 
                        std::vector<double> &m_zRootLPx_vals,
                        int numNodes, int numColsGen, std::string message) = 0;
        
        

        virtual void resetMaster(std::map<int, int> &inVars, OsiSolverInterface *si ) = 0;
        
        
        
        virtual void initializeDataStructures() = 0;
        
        
        OSDecompSolver();
        
        OSDecompSolver( OSOption* osoption);
                
        
        virtual ~OSDecompSolver() = 0;  
        

        
        //
};//end class OSDecompSolver



#endif