probdist.h

   1 #ifndef ZZGO_PROBDIST_H
   2 #define ZZGO_PROBDIST_H
   3
   4 /* Tools for picking an item according to a probability distribution. */
   5
   6 /* The probability distribution structure is designed to be once
   7  * initialized, then random items assigned a value repeatedly and
   8  * random items picked repeatedly as well. */
   9
  10 #include "move.h"
  11 #include "util.h"
  12
  13 /* The interface looks a bit funny-wrapped since we used to switch
  14  * between different probdist representations. */
  15
  16 struct probdist {
  17         int n;
  18         double *items; // [n], items[i] = P(pick==i)
  19         double total;
  20 };
  21 #define probdist_total(pd) ((pd)->total)
  22 #define probdist_one(pd, i) ((pd)->items[i])
  23 /* Probability so small that it's same as zero; used to compensate
  24  * for probdist.total inaccuracies. */
  25 #define PROBDIST_EPSILON 0.05
  26
  27 static void probdist_set(struct probdist *pd, int i, double val);
  28
  29 /* Pick a random item. ignore is a zero-terminated sorted array of items
  30  * that are not to be considered (and whose values are not in @total). */
  31 int probdist_pick(struct probdist *pd, int *ignore);
  32
  33
  34 /* We disable the assertions here since this is quite time-critical
  35  * part of code, and also the compiler is reluctant to inline the
  36  * functions otherwise. */
  37 static inline void
  38 probdist_set(struct probdist *pd, int i, double val)
  39 {
  40 #if 0
  41         assert(i >= 0 && i < pd->n);
  42         assert(val >= 0);
  43 #endif
  44         pd->total += val - pd->items[i];
  45         pd->items[i] = val;
  46 }
  47
  48 #endif