primer.c

   1 /*
   2  * Generates the repetitive part of the gnuplot input
   3  * for generating the initial triangle vertices from SGT's
   4  * model.
   5  */
   6
   7 #include "common.h"
   8 #include "mgraph.h"
   9
  10 int main(int argc, const char **argv) {
  11   static const int prec= DBL_DIG+2;
  12   int vi;
  13
  14   if (argc!=1) { fputs("need no args\n",stderr); exit(8); }
  15
  16   printf("%d %d %d %d %d\n%%-%d.%dg\n",
  17          X*Y, N, X, Y, D3,
  18          prec+5,prec);
  19
  20   FOR_VERTEX(vi) {
  21     int x= vi & XMASK; /* distance along strip */
  22     int y= vi >> YSHIFT; /* distance across strip */
  23     double u= y * 1.0 / (Y-1); /* SGT's u runs 0..1 across the strip */
  24
  25     /* SGT's v runs 0..pi along the strip, where the join is at 0==pi.
  26      * So that corresponds to 0..X (since 0==X in our scheme).
  27      * Vertices with odd y coordinate are halfway to the next x coordinate.
  28      */
  29     double v= (x*2 + (y&1)) * 1.0 / (X*2);
  30     v += 0.5;
  31     v *= M_PI;
  32
  33     printf("%-*.*g %-*.*g # %03x %2d %2d\n",
  34            prec+5,prec,u, prec+5,prec,v,
  35            vi, x, y);
  36   }
  37   flushoutput();
  38   return 0;
  39 }