OSColGenApp.h

Go to the documentation of this file.

/* $Id: OSColGenApp.h 3038 2009-11-07 11:43:44Z kmartin $ */
#ifndef OSCOLGENAPP_H
#define OSCOLGENAPP_H

// --------------------------------------------------------------------- //
#include "OSDataStructures.h"
#include "OSDefaultSolver.h" 
#include "OSrLReader.h"
#include "OSOption.h"
#include "OSInstance.h"
#include "OSoLReader.h"
#include "OSCoinSolver.h"
#include "OSBearcatSolverXij.h"
#include "OSDecompParam.h"
#include "OSDecompFactoryInitializer.h"

#include "OSNode.h"

#include <vector>

// --------------------------------------------------------------------- //
// --------------------------------------------------------------------- //
class OSColGenApp {
public:
        
        OSDecompFactoryInitializer *m_factoryInit;
        

        OSDecompSolver *m_osrouteSolver;
        
        //OSBearcatSolverXij *m_osrouteSolver;
        //OSDecompSolver *m_osrouteSolver;
        
        std::ostringstream outStr;
        
        OSInstance *m_osinstanceMaster;
        OSOption *m_osoption;
        OSResult *m_osresult;
        
        std::vector<double> m_zRootLPx_vals;
        
        
        std::vector<int> m_zRootLPx_idx;
        

        
        bool m_calledBranchAndBound;
        
        std::map<int, int> m_rowIdxVarMap;

        
        //the OS solver wrapper
        CoinSolver *m_solver;

        OsiSolverInterface *m_si;

        
        OSDecompParam m_osDecompParam;
        
        int m_numNodesGenerated;
        
        int m_numColumnsGenerated;
        
        int m_numColumnsOld;
        
        std::string m_message;
        
        double m_zUB;
        
        double m_zLB;
        
        double m_zRootLP;
        
        std::map<int, OSNode*> m_nodeMap;
        
        std::vector<int> m_zOptIndexes;
        
        std::vector<int> m_zOptRootLP;
        
        
        std::map<int, int> inVars;
        

        
        
        double *m_yA;
        
        double *m_yB;
        
        double* m_theta;
        
        int m_maxCols;
        
        int m_maxRows;
        
        void solveRestrictedMasterRelaxation( );

        
        void getOptions( OSOption *osoption);
        
        bool branchAndBound();
        
        void solve();
        
        /***  
         * print the information about the current tree
         */
        void printTreeInfo();
        
        
        void getInitialRestrictedMaster( );
        
        bool isInteger( const double *thetaVar, const int numThetaVar, 
                        const double tol);
        

        void printDebugInfo( );
        
        
        void checkNodeConsistency( const int rowIdx, const OSNode *osnode);
        
        void getCuts(const  double* thetaVar, const int numThetaVar,
                        int &numNewRows, int*  &numNonz, int** &colIdx,
                        double** &values, double* &rowLB, double* &rowUB) ;     
        
        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) ;
        
        
        OSNode* createChild(const OSNode *osnode, std::map<int, int> &varConMap,
                        const int rowIdx, const double rowLB, const double rowUB);
        
        
        
        void createBranchingCut(const int* thetaIdx, const double* theta, 
                        const int numThetaVar, std::map<int, int> &varConMap, int &rowIdx);
        
        
        void createBranchingCut(const double* theta, const int numThetaVar, 
                        std::map<int, int> &varConMap, int &rowIdx);
        
        
        void resetMaster( );
        
        OSColGenApp();
        
        OSColGenApp( OSOption *osption);
        

        ~OSColGenApp(); 

        //
};//end class OSColGenApp

#endif