Changeset 3874

Timestamp:

07/08/05 12:48:51 (4 years ago)

Author:

josyoun

Message:

( Links to Acro-related mail archives at http://software.sandia.gov/Acro/ )

( Links to Acro-related mail archives at http://software.sandia.gov/Acro/ )

Initial version of the GLPK derived class. This will not work with the current version of COIN since the function getModelPtr is public, noncost and getMutablePtr is private.

Files:

• pico/trunk/src/lp/pico/PicoGlpkLP.h (modified) (1 diff)

pico/trunk/src/lp/pico/PicoGlpkLP.h

@@ -20 +20 @@
 #include <pico/PicoLPSubclass.h>
 
-namespace pico {
-
-
-template <>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::initialize()
+namespace pico
 {
-}
+
+
+#if 0
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::initialize ()
+    {
+    }
+#endif
 
 //
 //
-template<>
-bool PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::isTimeLimitReached() const
-{
+    template <>
+        bool PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::isTimeLimitReached () const
+    {
+        LPX *const model = getModelPtr ();
+        double tmlim = lpx_get_real_parm (model, LPX_K_TMLIM);
+
+        if (fabs (tmlim) < 1e-8)
+              return true;
+        else
+              return false;
+
+    }
+
+
+    template <>
+        bool PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::isProvenPrimalUnbounded () const
+    {
+        LPX *model = getModelPtr ();
+        int prim = lpx_get_prim_stat (model);
+        int dual = lpx_get_dual_stat (model);
+
+        if (prim == LPX_P_FEAS && dual == LPX_D_NOFEAS)
+            return true;
+        else
+              return false;
+    }
+
+
+    template <>
+        bool PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::isProvenDualUnbounded () const
+    {
+        LPX *model = getModelPtr ();
+        int prim = lpx_get_prim_stat (model);
+        int dual = lpx_get_dual_stat (model);
+
+        if (prim == LPX_P_NOFEAS && dual == LPX_D_FEAS)
+            return true;
+        else
+              return false;
+    }
+
+
 #if 0
-int stat = CPXgetstat( const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(), getConstLpPtr() );
-#if (CPX_VERSION >= 800)
-return stat == CPXMIP_TIME_LIM_INFEAS || stat == CPXMIP_TIME_LIM_FEAS;
-#else
-return stat == CPX_TIME_LIM_INFEAS || stat == CPX_TIME_LIM_FEAS;
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getColLower (int whichVariable) const
+    {
+        return getColLower ()[whichVariable];
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getColUpper (int whichVariable) const
+    {
+        return getColUpper ()[whichVariable];
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getRowLower (int whichVariable) const
+    {
+        return getRowLower ()[whichVariable];
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getRowUpper (int whichVariable) const
+    {
+        return getRowUpper ()[whichVariable];
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getRightHandSide (int whichConstraint) const
+    {
+        return getRightHandSide ()[whichConstraint];
+    }
+
+
+    template <>
+        PicoLPInterface::constraintSense PicoLPSubclass <
+        OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getRowSense (int whichConstraint) const
+    {
+        return getRowSense ()[whichConstraint];
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getObjCoefficient (int whichVariable) const
+    {
+        return getObjCoefficients ()[whichVariable];
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getColSolution (const int whichCol) const
+    {
+        return getColSolution ()[whichCol];
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getRowPrice (const int whichRow) const
+    {
+        return getRowPrice ()[whichRow];
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getReducedCost (const int whichCol) const
+    {
+        return getReducedCost ()[whichCol];
+    }
 #endif
+
+
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getBasis (BasisArray & basisBuffer) const
+    {
+        LPX *model = getModelPtr ();
+
+        /* Get the basis status for all the structural variables */
+        for (int i = 0; i < getNumCols (); i++)
+            switch (lpx_get_col_stat (model, i + 1))
+              {
+              case LPX_BS:
+                  basisBuffer.put (i, basic);
+                  break;
+                  case LPX_NL:basisBuffer.put (i, atLower);
+                  break;
+                  case LPX_NU:basisBuffer.put (i, atUpper);
+                  break;
+                  default:basisBuffer.put (i, other);
+                  break;
+              }
+
+        /* Get the basis status for all the auxillary variables */
+        for (int j = 0; j < getNumRows (); j++)
+            switch (lpx_get_row_stat (model, j + 1))
+              {
+              case LPX_BS:
+                  basisBuffer.put (getNumCols () + j, basic);
+                  break;
+              case LPX_NL:
+                  basisBuffer.put (getNumCols () + j, atLower);
+                  break;
+              case LPX_NU:
+                  basisBuffer.put (getNumCols () + j, atUpper);
+                  break;
+              default:
+                  basisBuffer.put (getNumCols () + j, other);
+                  break;
+              }
+    }
+
+    template <>
+        basisState PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getBasis (int whichVariable) const
+    {
+        LPX *model = getModelPtr ();
+        int var;
+
+        // If we need an structural variable */
+        if (whichVariable < getNumCols ())
+              var = lpx_get_col_stat (model, whichVariable + 1);
+        else
+              var =
+                lpx_get_row_stat (model, whichVariable - getNumCols () + 1);
+        switch (var)
+          {
+          case LPX_BS:
+              return basic;
+              case LPX_NL:return atLower;
+              case LPX_NU:return atUpper;
+              default:return other;
+          }
+    }
+
+
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::setBasis (BasisArray & basisBuffer)
+    {
+        LPX *model = getModelPtr ();
+
+        /* Set the basis status of all the structural variables */
+        for (int i = 0; i < getNumCols (); i++)
+            switch (basisBuffer (i))
+              {
+              case basic:
+                  lpx_set_col_stat (model, i + 1, LPX_BS);
+                  break;
+              case atLower:
+                  lpx_set_col_stat (model, i + 1, LPX_NL);
+                  break;
+              case atUpper:
+                  lpx_set_col_stat (model, i + 1, LPX_NU);
+                  break;
+              }
+
+        /* Set the basis status for all the auxillary variables */
+        for (int j = 0; j < getNumRows (); j++)
+            switch (basisBuffer (getNumCols () + j))
+              {
+              case basic:
+                  lpx_set_col_stat (model, j + 1, LPX_BS);
+                  break;
+              case atLower:
+                  lpx_set_col_stat (model, j + 1, LPX_NL);
+                  break;
+              case atUpper:
+                  lpx_set_col_stat (model, j + 1, LPX_NU);
+                  break;
+              }
+    }
+
+
+    template <>
+        priceRule PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getPricingRule ()const
+    {
+        LPX *model = getModelPtr ();
+
+        switch (lpx_get_int_parm (model, LPX_K_PRICE))
+          {
+              /* Textbook pricing */
+          case 0:
+              return defaultPricing;
+              /* Steepest edge pricing */
+              case 1:return steepestEdgePricing;
+          }
+    }
+
+
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::setPricingRule (priceRule rule)
+    {
+        LPX *model = getModelPtr ();
+        switch (rule)
+          {
+              /* Textbook pricing */
+          case defaultPricing:
+              lpx_set_int_parm (model, LPX_K_PRICE, 0);
+              break;
+              /* Steepest edge pricing */
+          case steepestEdgePricing:
+              lpx_set_int_parm (model, LPX_K_PRICE, 1);
+              break;
+          }
+    }
+
+
+    template <>
+        double PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getTimeLimit () const
+    {
+        LPX *model = getModelPtr ();
+          return lpx_get_real_parm (model, LPX_K_TMLIM);
+    }
+
+
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::setTimeLimit (double time)
+    {
+        LPX *model = getModelPtr ();
+        lpx_set_real_parm (model, LPX_K_TMLIM, time);
+    }
+
+
+    template <>
+        clockType PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::getClockType ()const
+    {
+        return wallTime;
+    }
+
+
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::setClockType (clockType cType)
+    {
+        /* GLPK only allows a time limit to be set using the wall clock */
+    }
+
+
+
+    template <>
+        PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::PicoLPSubclass ()
+    {
+    }
+
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::write (char *filename, ilp_format format)
+    {
+        LPX *model = getModelPtr ();
+
+        switch (format)
+          {
+          case mps_format:
+          case glpk_mps_format:
+              lpx_write_mps (model, filename);
+              break;
+          case cplex_mps_format:
+              lpx_write_cpxlp (model, filename);
+              break;
+          }
+
+    }
+
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::printBrokenInfo ()
+    {
+    }
+
+    template <>
+        void PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt >::setMatrixCoeff (int row, int col, double newValue)
+    {
+    }
+
+    typedef PicoLPSubclass < OsiGlpkSolverInterface,
+        PicoLPCutMngmnt > PicoGlpkLP;
+}                               // namespace pico
+
 #endif
-return false;
-}
-
-
-template<>
-bool PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::isProvenPrimalUnbounded() const
-{
-#if 0
-int stat = CPXgetstat( const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(), getConstLpPtr() );
-#if (CPX_VERSION >= 800)
-return stat == CPX_STAT_UNBOUNDED;
-#else
-int method = CPXgetmethod(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(), getConstLpPtr() );
-// stat for version < 8 is relative to solver. If running dual simplex and stat is
-// infeasible, we do not necessarily know the primal is unbounded (could be infeasible).
- return ((method == CPX_ALG_PRIMAL && stat == CPX_UNBOUNDED));
+
 #endif
-#endif
-return false;
-}
-
-
-template<>
-bool PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::isProvenDualUnbounded() const
-{
-#if 0
-  // Version 8.0 and up returns status wrt the primal only.  There is no primal status that
-  // proves dual unboundedness.  CAP.
-#if (CPX_VERSION < 800)
-int stat = 
-#endif
-CPXgetstat( const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(), getConstLpPtr() );
-#if (CPX_VERSION >= 800)
-return false;
-#else
-int method = CPXgetmethod(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(), getConstLpPtr() );
-// stat for version < 8 is relative to solver. If running primal simplex and stat is
-// infeasible, we do not necessarily know the dual is unbounded (could be infeasible).
- return ((method == CPX_ALG_DUAL && stat == CPX_UNBOUNDED));
-#endif
-#endif
-return false;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getColLower(int whichVariable) const
-{
-#if 0
-double ans;
-CPXgetlb(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&ans,whichVariable,whichVariable);
-return ans;
-#endif
-return -999.0;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getColUpper(int whichVariable) const
-{
-#if 0
-double ans;
-CPXgetub(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&ans,whichVariable,whichVariable);
-return ans;
-#endif
-return -999.0;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getRowLower(int whichVariable) const
-{
-#if 0
-char sense;
-CPXgetsense(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&sense,whichVariable,whichVariable);
-double rhs;
-CPXgetrhs(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&rhs,whichVariable,whichVariable);
-double range;
-CPXgetrngval(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&range,whichVariable,whichVariable);
-double lower, upper;
-convertSenseToBound(sense, rhs, range, lower, upper);
-return lower;
-#endif
-return -999.0;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getRowUpper(int whichVariable) const
-{
-#if 0
-char sense;
-CPXgetsense(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&sense,whichVariable,whichVariable);
-double rhs;
-CPXgetrhs(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&rhs,whichVariable,whichVariable);
-double range;
-CPXgetrngval(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&range,whichVariable,whichVariable);
-double lower, upper;
-convertSenseToBound(sense, rhs, range, lower, upper);
-return upper;
-#endif
-return -999.0;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getRightHandSide(int whichConstraint) const
-{
-#if 0
-double rhs;
-CPXgetrhs(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&rhs,whichConstraint,whichConstraint);
-return rhs;
-#endif
-return -999.0;
-}
-
-
-template<>
-PicoLPInterface::constraintSense PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getRowSense(int whichConstraint) const
-{
-#if 0
-char csense;
-CPXgetsense( const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(), &csense, whichConstraint, whichConstraint);
-
-PicoLPInterface::constraintSense sense = notRestricted;
-switch (csense) {
-  case 'L':     sense = lessOrEqual; break;
-  case 'G':     sense = greaterOrEqual; break;
-  case 'E':     sense = equal; break;
-  case 'R':     sense = ranged; break;
-  case 'N':     sense = notRestricted; break;
-  default:
-        EXCEPTION_MNGR(runtime_error,"PicoLPSubclass::getRowSense - bad sense type=" << csense);
-  };
-return sense;
-#endif
-return notRestricted;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getObjCoefficient(int whichVariable) const
-{
-#if 0
-double ans;
-int err = CPXgetobj(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&ans,whichVariable,whichVariable);
-checkCPXerror( err, "CPXgetobj", "getObjCoefficient" );
-return ans;
-#endif 
-return -999.0;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getColSolution(const int whichCol) const
-{
-#if 0
-double ans = 0.0;
-int probType = CPXgetprobtype(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr());
-#if (CPX_VERSION >= 800)
-if ( probType == CPXPROB_MILP ) {
-#else
-if ( probType == CPXPROB_MIP ) {
-#endif
-   int err = CPXgetmipx( const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(), &ans, whichCol, whichCol);
-   if ( err == CPXERR_NO_INT_SOLN )
-      ans = 0.0;
-   else
-      checkCPXerror( err, "CPXgetmipx", "getColSolution" );
-   }
-else {
-   int err = CPXgetx( const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(), &ans, whichCol, whichCol );
-   if ( err == CPXERR_NO_SOLN )
-      ans = 0.0;
-   else
-      checkCPXerror( err, "CPXgetx", "getColSolution" );
-   }
-return ans;
-#endif
-return -999.0;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getRowPrice(const int whichRow) const
-{
-#if 0
-double ans;
-int err = CPXgetpi(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&ans,whichRow,whichRow);
-if (err == CPXERR_NO_SOLN)
-   ans = 0.0;
-else
-   checkCPXerror( err, "CPXgetpi", "getRowPrice" );
-return ans;
-#endif
-return -999.0;
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getReducedCost(const int whichCol) const
-{
-#if 0
-double ans;
-int err = CPXgetdj(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),&ans,whichCol,whichCol);
-if (err == CPXERR_NO_SOLN)
-   ans = 0.0;
-else
-   checkCPXerror( err, "CPXgetdj", "getReducedCost" );
-return ans;
-#endif
-return -999.0;
-}
-
-
-template <>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getBasis(BasisArray& basisBuffer) const
-{
-#if 0
-int numCols = getNumCols();
-int basisSize = numCols + getNumRows();
-utilib::IntVector tempBasis(basisSize);
-int tmp = CPXgetbase(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),(int*)(tempBasis.data()),(int*)(&tempBasis[numCols]));
-if (tmp)
-   EXCEPTION_MNGR(runtime_error,"CPLEX would not retreive basis");
-for (int i = 0; i < basisSize; i++) {
-  int val = tempBasis[i];
-  if (val == 0) basisBuffer.put(i, atLower);
-  else if (val == 3) basisBuffer.put(i, other);
-  else basisBuffer.put(i, (basisState)val);
-  }
-#endif
-}
-
-template <>
-basisState PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getBasis(int whichVariable) const
-{
-#if 0
-int numCols = getNumCols();
-utilib::IntVector tempBasis(numCols + getNumRows());
-int tmp = CPXgetbase(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),
-                (int*)(tempBasis.data()),(int*)(&tempBasis[numCols]));
-if (tmp)
-  EXCEPTION_MNGR(runtime_error, "PicoLPSubclass<T,V>::getBasis -- "
-                 "Couldn't retreive basis");
-int val = tempBasis[whichVariable];
-//
-// Convert to Osi standard representation for the basis
-//
-if (val == 0) return atLower;
-if (val == 3) return other;
-return (basisState)val;
-#endif
-return other;
-}
-
-
-template <>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::setBasis(BasisArray& basisBuffer)
-{
-#if 0
-static BasicArray<int> tempBasis;
-
-int basisSize = getNumRows() + getNumCols();
-tempBasis.resize(basisSize);
-for (int i = 0; i < basisSize; i++)
-  tempBasis[i] = (int)basisBuffer(i);
-status = CPXcopybase(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),getConstLpPtr(),(int*)(tempBasis.data()),
-                                        (int*)(&tempBasis[getNumCols() ]));
-if (status && lpWarnings)
-   ucout << "Warning: cplex would not load basis, status = " << status;
-else
-  {
-  // We may want to rethink at some point.  I don't think the MIP uses this 
-  // and I don't think we suppress basis copying if this is false.
-  have_basis = true;
-  // If you change the basis, presumably there needs to be some
-  // computation to make it a solution to the problem.
-  myState= unsolved;
-  }
-#endif
-}
-
-
-template<>
-priceRule PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getPricingRule() const
-{
-/// TODO
-return defaultPricing;
-}
-
-
-template<>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::setPricingRule(priceRule rule)
-{
-/// TODO
-}
-
-
-template<>
-double PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getTimeLimit() const
-{
-#if 0
-double ans;
-int tmp = CPXgetdblparam(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),CPX_PARAM_TILIM,&ans);
-if (tmp != 0)
-  {
-    EXCEPTION_MNGR(runtime_error, "Unable to get time limit");
-  }
-return ans;
-#endif
-return -999.0;
-}
-
-
-template<>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::setTimeLimit(double time)
-{
-#if 0
-int tmp = CPXsetdblparam(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),CPX_PARAM_TILIM,time);
-if (tmp)
-  EXCEPTION_MNGR(runtime_error,"Unable to set time limit");
-#endif
-}
-
-
-template<>
-clockType PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::getClockType() const
-{
-#if 0
-int ans;
-int tmp = CPXgetintparam(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),CPX_PARAM_CLOCKTYPE,&ans);
-if (tmp != 0)
-  {
-    EXCEPTION_MNGR(runtime_error,"Unable to get clock type");
-  }
-if ((clockType) ans != cpuTime && (clockType)ans != wallTime)
-  EXCEPTION_MNGR(runtime_error, "Unknown clock type " <<  ans)
-return (clockType)ans;
-#endif
-return (clockType)0.0;
-}
-
-
-template<>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::setClockType(clockType cType)
-{
-#if 0
-if (cType != cpuTime && cType != wallTime)
-  EXCEPTION_MNGR(runtime_error, "Improper clock type " << cType)
-int tmp = CPXsetintparam(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(),CPX_PARAM_CLOCKTYPE, (int)cType);
-if (tmp)
-  EXCEPTION_MNGR(runtime_error, "Unable to set clock type")
-#endif
-}
-
-
-
-template<>
-PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::PicoLPSubclass()
-{
-#if 0
-solverMethod = primalSimplex;
-// TODO
-// Turn cplex preprocessing off till we handle presolve infeasibility
-// errors more gracefully
-CPXsetintparam(const_cast<PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>*>(this)->OsiGlpkSolverInterface::getEnvironmentPtr(), CPX_PARAM_PREIND, 0);
-#endif
-}
-
-template<>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::write(char *filename, ilp_format format)
-{
-#if 0
-  // actually just use cplex mps and lp formats for now.
-int writeStatus;
-if (format == mps_format)
-  {
-    // This will write with original names, but since PICO currently
-    // doesn't give the real variable names (etc) to the cplex solver,
-    // it will be equivalent to the "REW" filetype = generic names
-  writeStatus = CPXwriteprob(getEnvironmentPtr(), getConstLpPtr(), filename, "MPS");
-  if (writeStatus)
-    ucout << "Error writing lp to file\n";
-  return;
-  }
-if (format == lp_format)
-  {
-  writeStatus = CPXwriteprob(getEnvironmentPtr(), getConstLpPtr(), filename, "LP");
-  if (status)
-    ucout << "Error writing lp to file\n";
-  return;
-  }
-// default for the other formats.
- PicoLPInterface::write(filename, format);
-#endif
-}
-
-template<>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::printBrokenInfo()
-{
-#if 0
-    int stat = CPXgetstat( getEnvironmentPtr(), getConstLpPtr() );
-    cout << "Broken LP has cplex status " << stat << "\n";
-#endif
-}
-
-template<>
-void PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt>::setMatrixCoeff(int row, int col, double newValue)
-{
-#if 0
-    int stat = CPXchgcoef(getEnvironmentPtr(), getConstLpPtr(), row, col, newValue);
-    if (stat != 0) EXCEPTION_MNGR(runtime_error, "Cplex error setting matrix coefficient for (" << row << ", " << col << ")")
-#endif
-}
-
-typedef PicoLPSubclass<OsiGlpkSolverInterface,PicoLPCutMngmnt> PicoGlpkLP;
-
-} // namespace pico
-
-#endif
-
-#endif