src/util/load.c

   1 /*
   2      This file is part of GNUnet.
   3      (C) 2010 Christian Grothoff (and other contributing authors)
   4
   5      GNUnet is free software; you can redistribute it and/or modify
   6      it under the terms of the GNU General Public License as published
   7      by the Free Software Foundation; either version 2, or (at your
   8      option) any later version.
   9
  10      GNUnet is distributed in the hope that it will be useful, but
  11      WITHOUT ANY WARRANTY; without even the implied warranty of
  12      MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13      General Public License for more details.
  14
  15      You should have received a copy of the GNU General Public License
  16      along with GNUnet; see the file COPYING.  If not, write to the
  17      Free Software Foundation, Inc., 59 Temple Place - Suite 330,
  18      Boston, MA 02111-1307, USA.
  19 */
  20
  21 /**
  22  * @file util/load.c
  23  * @brief functions related to load calculations
  24  * @author Christian Grothoff
  25  */
  26 #include "platform.h"
  27 #include "gnunet_load_lib.h"
  28
  29 #define DEBUG_LOAD GNUNET_NO
  30
  31 /**
  32  * Values we track for load calculations.
  33  */
  34 struct GNUNET_LOAD_Value
  35 {
  36
  37   /**
  38    * How fast should the load decline if no values are added?
  39    */
  40   struct GNUNET_TIME_Relative autodecline;
  41
  42   /**
  43    * Last time this load value was updated by an event.
  44    */
  45   struct GNUNET_TIME_Absolute last_update;
  46
  47   /**
  48    * Sum of all datastore delays ever observed (in ms).  Note that
  49    * delays above 64k ms are excluded (to avoid overflow within
  50    * first 4 billion requests).
  51    */
  52   uint64_t cummulative_delay;
  53
  54   /**
  55    * Sum of squares of all datastore delays ever observed (in ms).   Note that
  56    * delays above 64k ms are excluded (to avoid overflow within
  57    * first 4 billion requests).
  58    */
  59   uint64_t cummulative_squared_delay;
  60
  61   /**
  62    * Total number of requests included in the cummulative datastore delay values.
  63    */
  64   uint64_t cummulative_request_count;
  65
  66   /**
  67    * Current running average datastore delay.  Its relation to the
  68    * average datastore delay and it std. dev. (as calcualted from the
  69    * cummulative values) tells us our current load.
  70    */
  71   double runavg_delay;
  72
  73   /**
  74    * How high is the load?  0 for below average, otherwise
  75    * the number of std. devs we are above average, or 100 if the
  76    * load is so high that we currently cannot calculate it.
  77    */
  78   double load;
  79
  80 };
  81
  82
  83 static void
  84 internal_update (struct GNUNET_LOAD_Value *load)
  85 {
  86   struct GNUNET_TIME_Relative delta;
  87   unsigned int n;
  88
  89   if (load->autodecline.rel_value == GNUNET_TIME_UNIT_FOREVER_REL.rel_value)
  90     return;
  91   delta = GNUNET_TIME_absolute_get_duration (load->last_update);
  92   if (delta.rel_value < load->autodecline.rel_value)
  93     return;
  94   if (load->autodecline.rel_value == 0)
  95   {
  96     load->runavg_delay = 0.0;
  97     load->load = 0;
  98     return;
  99   }
 100   n = delta.rel_value / load->autodecline.rel_value;
 101   if (n > 16)
 102   {
 103     load->runavg_delay = 0.0;
 104     load->load = 0;
 105     return;
 106   }
 107   while (n > 0)
 108   {
 109     n--;
 110     load->runavg_delay = (load->runavg_delay * 7.0) / 8.0;
 111   }
 112 }
 113
 114
 115 /**
 116  * Create a new load value.
 117  *
 118  * @param autodecline speed at which this value should automatically
 119  *        decline in the absence of external events; at the given
 120  *        frequency, 0-load values will be added to the load
 121  * @return the new load value
 122  */
 123 struct GNUNET_LOAD_Value *
 124 GNUNET_LOAD_value_init (struct GNUNET_TIME_Relative autodecline)
 125 {
 126   struct GNUNET_LOAD_Value *ret;
 127
 128   ret = GNUNET_malloc (sizeof (struct GNUNET_LOAD_Value));
 129   ret->autodecline = autodecline;
 130   ret->last_update = GNUNET_TIME_absolute_get ();
 131   return ret;
 132 }
 133
 134
 135 /**
 136  * Change the value by which the load automatically declines.
 137  *
 138  * @param load load to update
 139  * @param autodecline frequency of load decline
 140  */
 141 void
 142 GNUNET_LOAD_value_set_decline (struct GNUNET_LOAD_Value *load,
 143                                struct GNUNET_TIME_Relative autodecline)
 144 {
 145   internal_update (load);
 146   load->autodecline = autodecline;
 147 }
 148
 149
 150 /**
 151  * Recalculate our load value.
 152  *
 153  * @param load load to update
 154  */
 155 static void
 156 calculate_load (struct GNUNET_LOAD_Value *load)
 157 {
 158   double stddev;
 159   double avgdel;
 160   double sum_val_i;
 161   double n;
 162   double nm1;
 163
 164   if (load->cummulative_request_count <= 1)
 165     return;
 166   /* calcuate std dev of latency; we have for n values of "i" that:
 167    *
 168    * avg = (sum val_i) / n
 169    * stddev = (sum (val_i - avg)^2) / (n-1)
 170    * = (sum (val_i^2 - 2 avg val_i + avg^2) / (n-1)
 171    * = (sum (val_i^2) - 2 avg sum (val_i) + n * avg^2) / (n-1)
 172    */
 173   sum_val_i = (double) load->cummulative_delay;
 174   n = ((double) load->cummulative_request_count);
 175   nm1 = n - 1.0;
 176   avgdel = sum_val_i / n;
 177   stddev =
 178       (((double) load->cummulative_squared_delay) - 2.0 * avgdel * sum_val_i +
 179        n * avgdel * avgdel) / nm1;
 180   if (stddev <= 0)
 181     stddev = 0.01;              /* must have been rounding error or zero; prevent division by zero */
 182   /* now calculate load based on how far out we are from
 183    * std dev; or if we are below average, simply assume load zero */
 184   if (load->runavg_delay < avgdel)
 185     load->load = 0.0;
 186   else
 187     load->load = (load->runavg_delay - avgdel) / stddev;
 188 }
 189
 190
 191 /**
 192  * Get the current load.
 193  *
 194  * @param load load handle
 195  * @return zero for below-average load, otherwise
 196  *         number of std. devs we are above average;
 197  *         100 if the latest updates were so large
 198  *         that we could not do proper calculations
 199  */
 200 double
 201 GNUNET_LOAD_get_load (struct GNUNET_LOAD_Value *load)
 202 {
 203   internal_update (load);
 204   calculate_load (load);
 205   return load->load;
 206 }
 207
 208
 209 /**
 210  * Get the average value given to update so far.
 211  *
 212  * @param load load handle
 213  * @return zero if update was never called
 214  */
 215 double
 216 GNUNET_LOAD_get_average (struct GNUNET_LOAD_Value *load)
 217 {
 218   double n;
 219   double sum_val_i;
 220
 221   internal_update (load);
 222   if (load->cummulative_request_count == 0)
 223     return 0.0;
 224   n = ((double) load->cummulative_request_count);
 225   sum_val_i = (double) load->cummulative_delay;
 226   return sum_val_i / n;
 227 }
 228
 229
 230 /**
 231  * Update the current load.
 232  *
 233  * @param load to update
 234  * @param data latest measurement value (for example, delay)
 235  */
 236 void
 237 GNUNET_LOAD_update (struct GNUNET_LOAD_Value *load, uint64_t data)
 238 {
 239   uint32_t dv;
 240
 241   internal_update (load);
 242   load->last_update = GNUNET_TIME_absolute_get ();
 243   if (data > 64 * 1024)
 244   {
 245     /* very large */
 246     load->load = 100.0;
 247     return;
 248   }
 249   dv = (uint32_t) data;
 250   load->cummulative_delay += dv;
 251   load->cummulative_squared_delay += dv * dv;
 252   load->cummulative_request_count++;
 253   load->runavg_delay = ((load->runavg_delay * 7.0) + dv) / 8.0;
 254 }
 255
 256
 257
 258 /* end of load.c */