Changeset 2213:2c094dfa176d in lemon-0.x for lemon/mip_glpk.cc

Timestamp:

09/08/06 17:51:32 (18 years ago)

Author:

athos

Branch:

default

Phase:

public

Convert:

svn:c9d7d8f5-90d6-0310-b91f-818b3a526b0e/lemon/trunk@2941

Message:

Some corrections.

File:

• lemon/mip_glpk.cc (modified) (2 diffs)

lemon/mip_glpk.cc

@@ -58 +58 @@
   LpGlpk::SolveExitStatus MipGlpk::_solve(){
     int result = lpx_simplex(lp);
-    result = lpx_integer(lp);
-    switch (result){
-      case LPX_E_OBJLL:
-      case LPX_E_OBJUL:
-      case LPX_E_ITLIM:
-      case LPX_E_TMLIM:
+    //
+    if (lpx_get_status(lp)==LPX_OPT){
+      //Maybe we could try the routine lpx_intopt(lp), a revised
+      //version of lpx_integer
+      result = lpx_integer(lp);
+      switch (result){
       case LPX_E_OK:
-        return SOLVED;
+        return SOLVED;
       default:
-        return UNSOLVED;
+        return UNSOLVED;
+      }
+      
     }
+    return UNSOLVED;
   }
 
@@ -74 +77 @@
   LpGlpk::SolutionStatus MipGlpk::_getMipStatus(){
 
-    //Meg kell nezni: ha az LP is infinite, akkor ez is, ha az is
-    //infeasible, akkor ez is, de ez lehet maskepp is infeasible.
-    int stat=  lpx_mip_status(lp);
-    switch (stat) {
-    case LPX_I_UNDEF://Undefined (no solve has been run yet)
-      return UNDEFINED;
-   case LPX_I_NOFEAS://There is no feasible integral solution (primal, I guess)
-      return INFEASIBLE;
-//     case LPX_UNBND://Unbounded
-//       return INFINITE;
-    case LPX_I_FEAS://Feasible
-      return FEASIBLE;
-    case LPX_I_OPT://Feasible
-      return OPTIMAL;
-    default:
+    if (lpx_get_status(lp)==LPX_OPT){
+      //Meg kell nezni: ha az LP is infinite, akkor ez is, ha az is
+      //infeasible, akkor ez is, de ez lehet maskepp is infeasible.
+      int stat=  lpx_mip_status(lp);
+      
+      switch (stat) {
+      case LPX_I_UNDEF://Undefined (no solve has been run yet)
+        return UNDEFINED;
+      case LPX_I_NOFEAS://There is no feasible integral solution
+        return INFEASIBLE;
+        //     case LPX_UNBND://Unbounded
+        //       return INFINITE;
+      case LPX_I_FEAS://Feasible
+        return FEASIBLE;
+      case LPX_I_OPT://Feasible
+        return OPTIMAL;
+      default:
       return UNDEFINED; //to avoid gcc warning
       //FIXME error
+      }
     }
+    else 
+      return UNDEFINED; //Maybe we could refine this: what does the LP
+                        //relaxation look like
+      
   }