yarrg/rsmain.c

   1 /**/
   2
   3 #include "rscommon.h"
   4
   5 #include <ctype.h>
   6
   7 int o_quiet= 0;
   8 double max_mass=-1, max_volu=-1, max_capi=-1;
   9 double distance_loss_factor_per_league;
  10 int max_dist= -1;
  11
  12 FILE *debug_file;
  13 FILE *output;
  14
  15 #define tabdebugf printf
  16
  17
  18 #define CTR(x)    int ctr_##x;
  19 #define CTRA(x,n) int ctr_##x[n];
  20   COUNTER_LIST
  21 #undef CTR
  22 #undef CTRA
  23
  24 static Bucket ****results[GRANUS];
  25   /* results[GRANUS][start_isle_ix][finalisle][midisle]-> */
  26
  27 static pid_t debugoutpid;
  28
  29 int main(int argc, const char **argv) {
  30   const char *arg;
  31   int i, ap;
  32   int granui;
  33
  34 #ifndef debug_flags
  35   debug_flags= ~( dbg_sql2 );
  36 #endif
  37
  38   for (;;) {
  39     arg= *++argv;
  40     if (arg[0] != '-') break;
  41 #ifndef debug_flags
  42     if (!strcmp(arg,"-DN")) {
  43       debug_flags= 0;
  44     } else if (!strcmp(arg,"-D1")) {
  45       debug_flags= ~(dbg_sql2|dbg_lp|dbg_value2);
  46     } else
  47 #endif
  48     {
  49       abort();
  50     }
  51   }
  52
  53   if (debug_flags) {
  54     /* glpk insists on writing stuff to stdout, and it does buffering,
  55      * so we route all our debug through it this too */
  56     int realstdout;
  57     sysassert( (realstdout= dup(1)) > 2 );
  58     sysassert( output= fdopen(realstdout,"w") );
  59
  60     int pfd[2];
  61     sysassert(! pipe(pfd) );
  62     sysassert( (debugoutpid= fork()) >=0 );
  63     if (!debugoutpid) {
  64       sysassert( dup2(pfd[0],0)==0 );
  65       sysassert( dup2(2,1)==1 );
  66       sysassert(! close(pfd[0]) );
  67       sysassert(! close(pfd[1]) );
  68       sysassert(! execlp("cat","cat",(char*)0) );
  69     }
  70     sysassert( dup2(pfd[1],1)==1 );
  71     sysassert(! close(pfd[0]) );
  72     sysassert(! close(pfd[1]) );
  73
  74     debug_file= stdout;
  75   } else {
  76     output= stdout;
  77     debug_file= stderr;
  78   }
  79
  80   const char *database= *argv++;
  81
  82   sysassert( !setvbuf(debug,0,_IOLBF,0) );
  83
  84   max_mass= atof(*argv++);
  85   max_volu= atof(*argv++);
  86   max_capi= atof(*argv++);
  87   double loss_per_league= atof(*argv++);
  88   distance_loss_factor_per_league= 1.0 - loss_per_league;
  89
  90   setup_sql(database);
  91   setup_value();
  92   setup_search();
  93
  94   for (i=0; i<narches; i++)
  95     fprintf(output,"arch %d %s\n",i,archnames[i]);
  96   fprintf(output,"setup complete, starting search\n");
  97
  98   arg= *argv++;
  99   if (!strcmp(arg,"specific")) {
 100     int ia[argc], ni=0;
 101     while ((arg= *argv++))
 102       ia[ni++]= atoi(arg);
 103
 104     double val= value_route(ni, ia, 0);
 105     fprintf(output, "route value is %g\n", val);
 106   } else if (!strcmp(arg,"search")) {
 107     for (granui=0; granui<GRANUS; granui++)
 108       MCALLOC(results[granui], argc);
 109
 110     max_dist= atoi(*argv++);
 111
 112     for (ap=0; ap<AP; ap++) {
 113       int nhs= atoi(*argv++);
 114       for (granui=0; granui<GRANUS; granui++) {
 115         nhighscores[granui][ap]= nhs;
 116         MCALLOC(highscores[granui][ap], nhs);
 117       }
 118     }
 119     const char *final_isle_spec= *argv++;
 120
 121     int resultsix= 0;
 122     while ((arg= argv[resultsix])) {
 123       int init_isle= atoi(arg);
 124
 125       int final_isle;
 126       if (!strcmp(final_isle_spec,"circ")) final_isle= init_isle;
 127       else if (!strcmp(final_isle_spec,"any")) final_isle= -1;
 128       else final_isle= atoi(final_isle_spec);
 129       assert(final_isle);
 130
 131       Bucket ****buckets_base_io[GRANUS];
 132       for (granui=0; granui<GRANUS; granui++)
 133         buckets_base_io[granui]= &results[granui][resultsix];
 134
 135       search(init_isle, final_isle, buckets_base_io);
 136       resultsix++;
 137     }
 138
 139     int mid, fin;
 140     for (granui=0; granui<granus; granui++) {
 141       fprintf(output,"\n");
 142       for (i=0; i<resultsix; i++) {
 143         tabdebugf("========== start #%d granui%d %s ==========\n",
 144                   i, granui, argv[i]);
 145         Bucket ***buckets_resultsix= results[granui][i];
 146         if (!buckets_resultsix) continue;
 147         tabdebugf("    ");
 148         for (mid=0; mid<granusz_mid[granui]; mid++) {
 149           tabdebugf("|   m%-3d   ",mid);
 150         }
 151         tabdebugf("\n");
 152         for (fin=0; fin<granusz_fin[granui]; fin++) {
 153           Bucket **buckets_fin= buckets_resultsix[fin];
 154           if (!buckets_fin) continue;
 155           tabdebugf("f%-3d",fin);
 156           for (mid=0; mid<granusz_mid[granui]; mid++) {
 157             Bucket *result= buckets_fin[mid];
 158             if (!result) {
 159               tabdebugf("|          ");
 160             } else {
 161               tabdebugf("|%5d",(int)(result->prs[A].value[A]));
 162               tabdebugf(" ");
 163               tabdebugf("%4d",(int)(result->prs[P].value[P]));
 164             }
 165           }
 166           tabdebugf("\n");
 167         }
 168       } /* i */
 169
 170       for (ap=0; ap<AP; ap++) {
 171         int pos;
 172         fprintf(output,"============== granui%d ap=%d ==============\n",
 173                 granui, ap);
 174         for (pos=nhighscores[granui][ap]-1; pos>=0; pos--) {
 175           HighScoreEntry *hs= &highscores[granui][ap][pos];
 176           Bucket *bucket= hs->bucket;
 177           if (!bucket) continue;
 178           OnePotentialResult *pr= &bucket->prs[ap];
 179           const int *const ports= pr->ports;
 180           int nports;
 181           for (nports=0; nports<MAX_ROUTELEN && ports[nports]>=0; nports++);
 182           int finisle= ports[nports-1];
 183           int finarch= isle2arch(finisle);
 184           int midisle= ports[nports/2];
 185           int midarch= route2midarch(ports,nports);
 186           fprintf(output,
 187                   " @%2d %c#%2d | start%3d mid%d:%3d f%d:%3d | %5d %5d %4d |",
 188                   pos, "ap"[ap], nhighscores[granui][ap] - 1 - pos,
 189                   ports[0], midarch,midisle, finarch,finisle,
 190                   (int)hs->value, (int)pr->value[A], (int)pr->value[P]);
 191           for (i=0; i<nports; i++) fprintf(output," %d",ports[i]);
 192           fprintf(output,"\n");
 193         } /* pos */
 194       } /* ap */
 195     } /* granui */
 196     fprintf(output,"\n");
 197
 198   } else {
 199     abort();
 200   }
 201
 202 #define CTR(x) fprintf(output,"  %-30s %10d\n",#x,ctr_##x);
 203 #define CTRA(x,n) for (i=0;i<n;i++) \
 204   fprintf(output,"  %-27s[%d] %10d\n",#x,i,ctr_##x[i]);
 205   COUNTER_LIST
 206 #undef CTR
 207
 208   if (debug_flags) {
 209     sysassert(! fclose(debug) );
 210     waitpid_check_exitstatus(debugoutpid,"debug cat",1);
 211   }
 212   sysassert(! fclose(output) );
 213
 214   return 0;
 215 }