src/cadet/desirability_table.c

   1 /* This file is in the public domain. */
   2
   3 /**
   4  * @brief Program to simulate results from #GCP_get_desirability_of_path()
   5  * for various plausible inputs.
   6  * @author Christian Grothoff
   7  */
   8 #include <stdio.h>
   9
  10 int
  11 main ()
  12 {
  13   for (unsigned int num_alts=1; num_alts<10; num_alts++)
  14     for (unsigned int off=0; off<10; off++)
  15       for (double delta=-(int) off;delta<=5;delta += 0.25)
  16       {
  17         double weight_alts;
  18
  19         if (delta <= - 1.0)
  20           weight_alts = - 1.0 * num_alts / delta; /* discount alternative paths */
  21         else if (delta >= 1.0)
  22           weight_alts = 1.0 * num_alts * delta; /* overcount alternative paths */
  23         else
  24           weight_alts = 1.0 * num_alts; /* count alternative paths normally */
  25
  26         fprintf (stderr,
  27                  "Paths: %u  Offset: %u  Delta: %5.2f  SCORE: %f\n",
  28                  num_alts,
  29                  off,
  30                  delta,
  31                  ((off + 1.0) / (weight_alts * weight_alts)));
  32       }
  33
  34
  35 }