[GNUnet-SVN] r19111 - gnunet/src/ats

[gnunet]

Author: wachs
Date: 2012-01-12 17:27:58 +0100 (Thu, 12 Jan 2012)
New Revision: 19111

Modified:
   gnunet/src/ats/gnunet-service-ats_addresses_mlp.c
   gnunet/src/ats/gnunet-service-ats_addresses_mlp.h
Log:
- more mlp


Modified: gnunet/src/ats/gnunet-service-ats_addresses_mlp.c
===================================================================
--- gnunet/src/ats/gnunet-service-ats_addresses_mlp.c   2012-01-12 14:52:08 UTC 
(rev 19110)
+++ gnunet/src/ats/gnunet-service-ats_addresses_mlp.c   2012-01-12 16:27:58 UTC 
(rev 19111)
@@ -39,6 +39,85 @@
 
 
 /**
+ * Solves the MLP problem
+ * @return GNUNET_OK if could be solved, GNUNET_SYSERR on failure
+ */
+int
+mlp_solve_lp_problem (struct GAS_MLP_Handle *mlp)
+{
+  int res;
+
+  /* LP presolver?
+   * Presolver is required if the problem was modified and an existing
+   * valid basis is now invalid */
+  if (mlp->presolver_required == GNUNET_YES)
+    mlp->control_param_lp.presolve = GLP_ON;
+  else
+    mlp->control_param_lp.presolve = GLP_OFF;
+
+
+  /* Solve LP problem to have initial valid solution */
+lp_solv:
+  res = glp_simplex(mlp->prob, &mlp->control_param_lp);
+  if (res == 0)
+  {
+    /* The LP problem instance has been successfully solved. */
+  }
+  else if (res == GLP_EITLIM)
+  {
+    /* simplex iteration limit has been exceeded. */
+    // TODO Increase iteration limit?
+  }
+  else if (res == GLP_ETMLIM)
+  {
+    /* Time limit has been exceeded.  */
+    // TODO Increase time limit?
+  }
+  else
+  {
+    /* Problem was ill-defined, retry with presolver */
+    if (mlp->presolver_required == GNUNET_NO)
+    {
+      mlp->presolver_required = GNUNET_YES;
+      goto lp_solv;
+    }
+    else
+    {
+      /* Problem was ill-defined, no way to handle that */
+      GNUNET_log_from (GNUNET_ERROR_TYPE_ERROR,
+          "ats-mlp",
+          "Solving LP problem failed: glp_simplex error 0x%X", res);
+      return GNUNET_SYSERR;
+    }
+  }
+
+  /* Analyze problem status  */
+  res = glp_get_status (mlp->prob);
+  switch (res) {
+    /* solution is optimal */
+    case GLP_OPT:
+    /* solution is feasible */
+    case GLP_FEAS:
+      break;
+
+    /* Problem was ill-defined, no way to handle that */
+    default:
+      GNUNET_log_from (GNUNET_ERROR_TYPE_ERROR,
+          "ats-mlp",
+          "Solving LP problem failed, no solution: glp_get_status 0x%X", res);
+      return GNUNET_SYSERR;
+      break;
+  }
+
+  /* solved sucessfully, no presolver required next time */
+  mlp->presolver_required = GNUNET_NO;
+
+  return GNUNET_OK;
+}
+
+
+
+/**
  * Init the MLP problem solving component
  *
  * @param max_duration maximum numbers of iterations for the LP/MLP Solver
@@ -62,10 +141,21 @@
 
   /* Init LP solving parameters */
   glp_init_smcp(&GAS_mlp->control_param_lp);
+#if DEBUG_MLP
+  GAS_mlp->control_param_lp.msg_lev = GLP_MSG_ALL;
+#else
+  GAS_mlp->control_param_lp.msg_lev = GLP_MSG_OFF;
+#endif
   GAS_mlp->control_param_lp.it_lim = max_iterations;
   GAS_mlp->control_param_lp.tm_lim = max_duration.rel_value;
+
   /* Init MLP solving parameters */
   glp_init_iocp(&GAS_mlp->control_param_mlp);
+#if DEBUG_MLP
+  GAS_mlp->control_param_mlp.msg_lev = GLP_MSG_ALL;
+#else
+  GAS_mlp->control_param_mlp.msg_lev = GLP_MSG_OFF;
+#endif
   GAS_mlp->control_param_mlp.tm_lim = max_duration.rel_value;
 
   return GNUNET_OK;

Modified: gnunet/src/ats/gnunet-service-ats_addresses_mlp.h
===================================================================
--- gnunet/src/ats/gnunet-service-ats_addresses_mlp.h   2012-01-12 14:52:08 UTC 
(rev 19110)
+++ gnunet/src/ats/gnunet-service-ats_addresses_mlp.h   2012-01-12 16:27:58 UTC 
(rev 19111)
@@ -33,6 +33,8 @@
 #ifndef GNUNET_SERVICE_ATS_ADDRESSES_MLP_H
 #define GNUNET_SERVICE_ATS_ADDRESSES_MLP_H
 
+#define VERBOSE GNUNET_EXTRA_LOGGING
+#define DEBUG_MLP GNUNET_EXTRA_LOGGING
 
 #define MLP_MAX_EXEC_DURATION   
GNUNET_TIME_relative_multiply(GNUNET_TIME_UNIT_SECONDS, 3)
 #define MLP_MAX_ITERATIONS      INT_MAX
@@ -68,6 +70,20 @@
    */
   unsigned int max_iterations;
 
+
+  /* state information */
+
+  /**
+   * Do we need to use the LP presolver?
+   *
+   * If the problem addresses were added or removed and the last basis was we
+   * need to use the presolver.
+   * presolver_required == GNUNET_YES
+   *
+   * If values were modified, we can reuse a valid basis
+   * presolver_required == GNUNET_NO
+   */
+  int presolver_required;
 };
 
 /**