/Users/kmartin/Documents/files/code/cpp/OScpp/COIN-OS/OS/src/OSModelInterfaces/OSmps2osil.cpp

Go to the documentation of this file.

#include "OSmps2osil.h"
#include <iostream> 



using std::cout;
using std::endl;

OSmps2osil::OSmps2osil( std::string mpsfilename)
{
        m_MpsData = new CoinMpsIO();
        m_MpsData->readMps( &mpsfilename[ 0] ); 
        m_CoinPackedMatrix  =  new CoinPackedMatrix( *(m_MpsData->getMatrixByCol()));
}// end constructor

OSmps2osil::~OSmps2osil(){
        delete m_MpsData;
        m_MpsData = NULL;
        delete m_CoinPackedMatrix;
        m_CoinPackedMatrix = NULL;
        osinstance->instanceData->linearConstraintCoefficients->value->el = NULL;
        osinstance->instanceData->linearConstraintCoefficients->start->el = NULL;
        osinstance->instanceData->linearConstraintCoefficients->colIdx->el = NULL;
        osinstance->instanceData->linearConstraintCoefficients->rowIdx->el = NULL;
        delete osinstance;     
        osinstance = NULL;
        #ifdef DEBUG
        cout << "now delete m_MpsData" << endl;
        #endif
        delete m_MpsData;

} // end destructor

// Kipp-- put in integer programming
bool OSmps2osil::createOSInstance( )
{
        osinstance = new OSInstance();
        int i;
        int numvar =  m_MpsData->getNumCols();
        int numrows =  m_MpsData->getNumRows();
        int numnonz = m_MpsData->getNumElements();
        int numberObj = 1;
        int objIndex = -1;
        double objWeight = 1.0;
        osinstance->setInstanceName( const_cast<char*>(m_MpsData->getProblemName()));
        //
        // get the variable information
        //
        std::string initString;
        double init = OSNAN;
        osinstance->setVariableNumber( numvar);
        for(i = 0; i < numvar; i++)
        {
                osinstance->addVariable(i, m_MpsData->columnName( i), 
                        m_MpsData->getColLower()[ i], m_MpsData->getColUpper()[ i], 
                        m_MpsData->isInteger( i)? 'I':'C', init, initString);
        }
        //
        // now create the objective function
        //      
        SparseVector* objectiveCoefficients = NULL;
        objectiveCoefficients = new SparseVector( numvar);
        double *p = const_cast<double*>(m_MpsData->getObjCoefficients());
        for(i = 0; i < numvar; i++){
                objectiveCoefficients->indexes[i] = i;
                objectiveCoefficients->values[i] = *(p++);
        }
        // kipp -- check on this -- does MPS assume a min???
        osinstance->setObjectiveNumber( numberObj) ;
        osinstance->addObjective(objIndex, m_MpsData->getObjectiveName(), 
        "min", m_MpsData->objectiveOffset(),  objWeight, objectiveCoefficients) ;
        delete objectiveCoefficients;
        objectiveCoefficients = NULL;
        //
        // now fill in row information
        //
        osinstance->setConstraintNumber( numrows);
        double constant = 0.0;
        for(i = 0; i < numrows; i++){
                osinstance->addConstraint(i, m_MpsData->rowName( i), m_MpsData->getRowLower()[i], 
                m_MpsData->getRowUpper()[i], constant);
        }
        //
        // now fill in the lp arrays
        //
        int valuesBegin = 0;
        int valuesEnd = numnonz - 1;
        int startsBegin = 0;
        int indexesBegin = 0;
        int indexesEnd = numnonz - 1;
        int startsEnd = m_CoinPackedMatrix->isColOrdered()?numvar:numrows;
        osinstance->setLinearConstraintCoefficients(numnonz,  m_CoinPackedMatrix->isColOrdered(), 
                const_cast<double*>(m_CoinPackedMatrix->getElements()), valuesBegin,  valuesEnd, 
                const_cast<int*>(m_CoinPackedMatrix->getIndices()),  indexesBegin,  indexesEnd,                         
                const_cast<int*>(m_CoinPackedMatrix->getVectorStarts()),  startsBegin,  startsEnd);
        return true;
}

---