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

[gnunet]

Author: wachs
Date: 2012-01-19 13:45:08 +0100 (Thu, 19 Jan 2012)
New Revision: 19254

Modified:
   gnunet/src/ats/gnunet-service-ats_addresses_mlp.c
   gnunet/src/ats/gnunet-service-ats_addresses_mlp.h
Log:
- disabling code : debugging was required


Modified: gnunet/src/ats/gnunet-service-ats_addresses_mlp.c
===================================================================
--- gnunet/src/ats/gnunet-service-ats_addresses_mlp.c   2012-01-19 11:06:37 UTC 
(rev 19253)
+++ gnunet/src/ats/gnunet-service-ats_addresses_mlp.c   2012-01-19 12:45:08 UTC 
(rev 19254)
@@ -29,11 +29,10 @@
 #include "gnunet-service-ats_addresses.h"
 #include "gnunet-service-ats_addresses_mlp.h"
 #include "gnunet_statistics_service.h"
-#if HAVE_LIBGLPK
 #include "glpk.h"
-#endif
 #include "float.h"
 
+#define WRITE_MLP GNUNET_NO
 #define DEBUG_ATS GNUNET_YES
 /* A very big value */
 #define M DBL_MAX
@@ -199,7 +198,7 @@
     if (mlp->ia != NULL)
     {
       GNUNET_free (mlp->ia);
-      mlp->ja = NULL;
+      mlp->ia = NULL;
     }
 
     /* delete column index */
@@ -255,10 +254,10 @@
   mlp->ci++;
 
   mlp->ia[mlp->ci] = row_index;
-  mlp->ja[mlp->ci] = mlpi->c_b;
+  mlp->ja[mlp->ci] = mlpi->c_n;
   mlp->ar[mlp->ci] = -M;
   mlp->ci++;
-
+#if 0
   /* c 3) minimum bandwidth
    * b_t + (-n_t * b_min) >= 0
    */
@@ -274,7 +273,7 @@
   mlp->ci++;
 
   mlp->ia[mlp->ci] = row_index;
-  mlp->ja[mlp->ci] = mlpi->c_b;
+  mlp->ja[mlp->ci] = mlpi->c_n;
   mlp->ar[mlp->ci] = -mlp->b_min;
   mlp->ci++;
 
@@ -286,6 +285,15 @@
   mlp->ar[mlp->ci] = 1;
   mlp->ci++;
 
+  /* c 6) maximize diversity
+   * (1)*n_1 + ... + (1)*n_m - d == 0
+   */
+  mlp->ia[mlp->ci] = mlp->r_c6;
+  mlp->ja[mlp->ci] = mlpi->c_n;
+  mlp->ar[mlp->ci] = 1;
+  mlp->ci++;
+#endif
+
   return GNUNET_OK;
 }
 
@@ -302,6 +310,7 @@
 {
   unsigned int n_addresses;
   int row_index;
+  //int c;
 
   /* Problem matrix*/
   n_addresses = GNUNET_CONTAINER_multihashmap_size(addresses);
@@ -335,10 +344,17 @@
    * #indices: 7 * |n_addresses|
    *
    * optimality constraints:
-   * tbc
+   *
+   * c 6) diversity
+   * #rows: 1
+   * #indices: |n_addresses| + 1
+   *
+   * c 7) quality
+   * #rows: |quality properties|
+   * #indices:|quality properties| + |n_addresses|
    * */
 
-  int pi = (7 * n_addresses);
+  int pi = ((7 * n_addresses) /*+ (2 * n_addresses +  mlp->m_q + 1)*/);
   mlp->cm_size = pi;
   mlp->ci = 0;
 
@@ -361,10 +377,27 @@
    * c 1) bandwidth capping
    * c 3) minimum bandwidth
    * c 4) minimum number of connections
+   * c 6) maximize diversity
    */
+
+  int min = mlp->n_min;
+  if (mlp->n_min > mlp->c_p)
+    min = mlp->c_p;
+/*
   mlp->r_c4 = glp_add_rows (mlp->prob, 1);
-  glp_set_row_bnds (mlp->prob, mlp->r_c4, GLP_LO, mlp->n_min, 0.0);
-
+  glp_set_row_bnds (mlp->prob, mlp->r_c4, GLP_FX, min, min);
+*/
+  /* Add row for c6) */
+#if 0
+  mlp->r_c6 = glp_add_rows (mlp->prob, 1);
+  /* Set type type to fix */
+  glp_set_row_bnds (mlp->prob, mlp->r_c6, GLP_FX, 0.0, 0.0);
+  /* Setting -D */
+  ia[mlp->ci] = mlp->r_c6 ;
+  ja[mlp->ci] = mlp->c_d;
+  ar[mlp->ci] = -1;
+  mlp->ci++;
+#endif
   GNUNET_CONTAINER_multihashmap_iterate (addresses, create_constraint_it, mlp);
 
   /* Adding constraint rows
@@ -374,7 +407,7 @@
    * c 2) 1 address per peer
    * sum (n_p1_1 + ... + n_p1_n) = 1
    */
-
+return;
   /* Adding rows for c 2) */
   row_index = glp_add_rows (mlp->prob, mlp->c_p);
 
@@ -396,9 +429,26 @@
 
       addr = addr->next;
     }
-
     peer = peer->next;
   }
+
+
+  /* For all quality metrics */
+/*
+  for (c = 0; c < mlp->m_q; c++)
+  {
+    struct ATS_Peer *p = mlp->peer_head;
+    while (p != NULL)
+    {
+      ia[mlp->ci] = row_index;
+      ja[mlp->ci] = mlp->c_q[c];
+      ar[mlp->ci] = mlp->;
+      mlp->ci++;
+
+      p = p->next;
+    }
+  }
+*/
 }
 
 
@@ -428,6 +478,10 @@
   mlpi->c_b = col;
   mlpi->c_n = col + 1;
 
+  GNUNET_log_from (GNUNET_ERROR_TYPE_ERROR,
+      "ats-mlp",
+      "Culoumn %i %i\n", mlpi->c_b, mlpi->c_n);
+
   GNUNET_asprintf (&name, "b_%s_%s", GNUNET_i2s (&address->peer), 
address->plugin);
   glp_set_col_name (mlp->prob, mlpi->c_b , name);
   GNUNET_free (name);
@@ -438,6 +492,7 @@
   /* Objective function coefficient == 0 */
   glp_set_obj_coef (mlp->prob, mlpi->c_b , 0);
 
+
   /* Add usage column */
   GNUNET_asprintf (&name, "n_%s_%s", GNUNET_i2s (&address->peer), 
address->plugin);
   glp_set_col_name (mlp->prob, mlpi->c_n, name);
@@ -447,7 +502,7 @@
   /* Integer value*/
   glp_set_col_kind (mlp->prob, mlpi->c_n, GLP_IV);
   /* Objective function coefficient == 0 */
-  glp_set_obj_coef (mlp->prob, mlpi->c_n, 0);
+  glp_set_obj_coef (mlp->prob, mlpi->c_n, 1);
 
   return GNUNET_OK;
 }
@@ -464,10 +519,15 @@
 mlp_create_problem (struct GAS_MLP_Handle *mlp, struct 
GNUNET_CONTAINER_MultiHashMap * addresses)
 {
   int res = GNUNET_OK;
+/*
   int col;
   int c;
   char *name;
+*/
+  GNUNET_assert (mlp->prob == NULL);
 
+  /* create the glpk problem */
+  mlp->prob = glp_create_prob ();
 
   /* Set a problem name */
   glp_set_prob_name (mlp->prob, "gnunet ats bandwidth distribution");
@@ -478,7 +538,7 @@
   /* Adding invariant columns */
 
   /* Diversity d column  */
-
+#if 0
   col = glp_add_cols (mlp->prob, 1);
   mlp->c_d = col;
   /* Column name */
@@ -522,13 +582,16 @@
     /* Coefficient == Qm */
     glp_set_obj_coef (mlp->prob, col + c, mlp->co_Q[c]);
   }
-
+#endif
   /* Add columns for addresses */
   GNUNET_CONTAINER_multihashmap_iterate (addresses, create_columns_it, mlp);
 
   /* Add constraints */
   mlp_add_constraints_all_addresses (mlp, addresses);
 
+  /* Load the matrix */
+  glp_load_matrix(mlp->prob, (mlp->ci - 1), mlp->ia, mlp->ja, mlp->ar);
+
   return res;
 }
 
@@ -584,7 +647,7 @@
       /* Problem was ill-defined, no way to handle that */
       GNUNET_log_from (GNUNET_ERROR_TYPE_ERROR,
           "ats-mlp",
-          "Solving LP problem failed:  %s\n", mlp_solve_to_string(res));
+          "Solving LP problem failed: %i %s\n", res, mlp_solve_to_string(res));
       return GNUNET_SYSERR;
     }
   }
@@ -727,11 +790,66 @@
 {
   int res;
   mlp->last_execution = GNUNET_TIME_absolute_get ();
+#if DEBUG_ATS
+  GNUNET_log (GNUNET_ERROR_TYPE_DEBUG, "Problem solving\n");
+#endif
+
+#if WRITE_MLP
+  char * name;
+  static int i;
+  i++;
+  GNUNET_asprintf(&name, "problem_%i", i);
+  glp_write_lp (mlp->prob, 0, name);
+  GNUNET_free (name);
+# endif
+
   res = mlp_solve_lp_problem (mlp);
-  if (res == GNUNET_OK)
-    res = mlp_solve_mlp_problem (mlp);
 
+#if WRITE_MLP
+  GNUNET_asprintf(&name, "problem_%i_lp_solution", i);
+  glp_print_sol (mlp->prob,  name);
+  GNUNET_free (name);
+# endif
 
+  if (res != GNUNET_OK)
+  {
+#if DEBUG_ATS
+  GNUNET_log (GNUNET_ERROR_TYPE_DEBUG, "LP Problem solving failed\n");
+#endif
+    return GNUNET_SYSERR;
+  }
+
+  res = mlp_solve_mlp_problem (mlp);
+
+#if WRITE_MLP
+  GNUNET_asprintf(&name, "problem_%i_mlp_solution", i);
+  glp_print_mip (mlp->prob, name);
+  GNUNET_free (name);
+# endif
+
+
+  if (res != GNUNET_OK)
+  {
+#if DEBUG_ATS
+  GNUNET_log (GNUNET_ERROR_TYPE_DEBUG, "MLP Problem solving failed\n");
+#endif
+    return GNUNET_SYSERR;
+  }
+
+  res = glp_mip_status(mlp->prob);
+
+  if (res != GNUNET_OK)
+  {
+#if DEBUG_ATS
+  GNUNET_log (GNUNET_ERROR_TYPE_DEBUG, "LP Problem solving failed\n");
+#endif
+    return GNUNET_SYSERR;
+  }
+
+#if DEBUG_ATS
+  GNUNET_log (GNUNET_ERROR_TYPE_DEBUG, "Problem solved: %i %s\n", res, 
mlp_status_to_string(res));
+#endif
+
   /* Process result */
 
   if (mlp->mlp_task != GNUNET_SCHEDULER_NO_TASK)
@@ -863,7 +981,7 @@
 
   /* Init LP solving parameters */
   glp_init_smcp(&mlp->control_param_lp);
-#if DEBUG_MLP
+#if DEBUG_ATS
   mlp->control_param_lp.msg_lev = GLP_MSG_ALL;
 #else
   mlp->control_param_lp.msg_lev = GLP_MSG_OFF;
@@ -873,7 +991,7 @@
 
   /* Init MLP solving parameters */
   glp_init_iocp(&mlp->control_param_mlp);
-#if DEBUG_MLP
+#if DEBUG_ATS
   mlp->control_param_mlp.msg_lev = GLP_MSG_ALL;
 #else
   mlp->control_param_mlp.msg_lev = GLP_MSG_OFF;
@@ -968,7 +1086,11 @@
 
   /* Recalculate */
   if (new == GNUNET_YES)
+  {
+    mlp_delete_problem (mlp);
+    mlp_create_problem (mlp, addresses);
     mlp->presolver_required = GNUNET_YES;
+  }
   mlp_solve_problem (mlp);
 }
 

Modified: gnunet/src/ats/gnunet-service-ats_addresses_mlp.h
===================================================================
--- gnunet/src/ats/gnunet-service-ats_addresses_mlp.h   2012-01-19 11:06:37 UTC 
(rev 19253)
+++ gnunet/src/ats/gnunet-service-ats_addresses_mlp.h   2012-01-19 12:45:08 UTC 
(rev 19254)
@@ -34,7 +34,6 @@
 #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)
@@ -168,6 +167,8 @@
 
   /* Row index constraint 4: minimum connections */
   unsigned int r_c4;
+  /* Row index constraint 6: maximize diversity */
+  unsigned int r_c6;
 
   /* column index Diversity (D) column */
   int c_d;