home *** CD-ROM | disk | FTP | other *** search
/ C/C++ Interactive Guide / c-cplusplus-interactive-guide.iso / c_ref / csource3 / 126_01 / facsplot.c < prev    next >
Text File  |  1985-03-11  |  24KB  |  854 lines
  1. /*********************************************************************\
  2. ** .---------------------------------------------------------------. **
  3. ** |                                                               | **
  4. ** |                                                               | **
  5. ** |         Copyright (c) 1981, 1982, 1983 by Eric Martz.         | **
  6. ** |                                                               | **
  7. ** |                                                               | **
  8. ** |       Permission is hereby granted to use this source         | **
  9. ** |       code only for non-profit purposes. Publication of       | **
  10. ** |       all or any part of this source code, as well as         | **
  11. ** |       use for business purposes is forbidden without          | **
  12. ** |       written permission of the author and copyright          | **
  13. ** |       holder:                                                 | **
  14. ** |                                                               | **
  15. ** |                          Eric Martz                           | **
  16. ** |                         POWER  TOOLS                          | **
  17. ** |                    48 Hunter's Hill Circle                    | **
  18. ** |                      Amherst MA 01002 USA                     | **
  19. ** |                                                               | **
  20. ** |                                                               | **
  21. ** `---------------------------------------------------------------' **
  22. \*********************************************************************/
  23.  
  24. /*------------------------------------------------------------------*/
  25. /* This source provides an example of how HIPLOT.C has been used
  26. in producing actual plots. The application concerns distributions
  27. of 256 integers, each representing a count in a channel, from an
  28. instrument called a flow cytofluorometer, or
  29. fluorescence-activated cell sorter. Since use of the facsplot
  30. program as it stands is hardware dependent, you will not be able
  31. to run this program (or even compile and link it in its present
  32. form). However, it provides some functions which may serve as
  33. useful models for anyone wishing to write a more general purpose
  34. plotting package. In particular, see
  35.  
  36.     axis() which draws an axis
  37.     axis_scale() which puts tick marks and numbers on it
  38.         tick()
  39.     vlabel() and hlabel() which write labels next to the vertical
  40.         and horizontal axes
  41.     
  42.     scalepoint()
  43. */
  44. /*------------------------------------------------------------------*/
  45. #include "bdscio.h"
  46.  
  47. #define NIOBUFS 1
  48. #include "fio1.h"
  49.  
  50. /* FACS */
  51. #define STDERR 1
  52. #define NDISTS 5 /* NUMBER OF DISTRIBUTION BUFFERS IN MEMORY */
  53. #define DISTLEN 255
  54. char Buf[MAXLINE], **Pp[32];
  55. char Memdist[NDISTS+1][(DISTLEN+1)*5]; /* BDSC LACKS 3-DIMENSIONAL ARRAYS */
  56. char Dfinmem[NDISTS+1][16]; /* NAMES OF DISTRIBUTION FILES IN MEMORY */
  57. char Tot[NDISTS+1][5], Av[NDISTS+1][5], Max[NDISTS+1][5],
  58.     Condist[NDISTS+1][15+1], Pos[NDISTS+1][5];
  59. int Location; /* INDEX LOCATION OF CURRENT DISTRIBUTION IN MEMORY */
  60. int Nexttogo;    /* NEXT MEMORY ARRAY TO BE REFILLED BY NEEDED DISTRIBUTION */
  61. int Verbose;
  62. char *Arr;    /* POINTER TO CURRENT DISTRIBUTION IN MEMORY */
  63.  
  64. char Fpone[5], Fpzero[5], Fpmone[5],
  65.     Fphundred[5], Fpthousand[5]; /* FLOAT */
  66. char F1[5], F2[5], F3[5], F4[5]; /* FLOAT */
  67. FILE *Fpin, *Fpout;
  68.  
  69. /* MOSC */
  70. int Begregn, Endregn;
  71.  
  72. */
  73. /*------------------------------------------------------------------*/
  74. /* FACS PLOT */
  75. #define LABEL 2
  76. #define SCALE 2
  77. #define REBUF 4000
  78. #define VHALF 755
  79. #define HHALF 1015
  80. char Hmin_f[5], Vmin_f[5], Hfact_f[5], Vfact_f[5];
  81. char Xmin_f[5], Ymin_f[5];
  82. char Buf_repeat[REBUF], *At_repeat;
  83. int Lsize, Nsize; /* SIZES OF CHARACTERS FOR LABEL, SCALE */
  84. int Fill_repeat;
  85. /*------------------------------------------------------------------*/
  86.  
  87. #define FACSPLOT 1
  88.  
  89. /*------------------------------------------------------------------*/
  90. main(argc, argv)
  91.     int argc;
  92.     char **argv;
  93.     {
  94.     int option;
  95. #include "fio2.h"
  96.     if (0) putchar('a');
  97.     if (0) getchar();
  98.     if (0) fopen("x",'r');
  99.     if (0) fclose(Fpin);
  100.     if (argc >1) Verbose = YES;
  101.     else Verbose = NO;
  102.     init_facs();
  103.  
  104.     option = '1';
  105.     FOREVER {
  106. /*        fprintf(STDERR,
  107.            "\n\t1. MOSC\n%s%s%s%s%s",
  108.             "\t2. Regions\n",
  109.             "\t3. Display statistics\n",
  110.             "\t4. List distribution\n",
  111.             "\t5. Normalize\n",
  112.             "\t6. Quit\n\n");
  113.         option = charq("Option","123456");
  114. */
  115.  
  116.         switch(option) {
  117.             case '1':
  118.                 fplot();
  119.                 break;
  120.             case '2':
  121.                 break;
  122.             case '3':
  123.                 break;
  124.             case '4':
  125.                 exit(0);
  126.         }
  127.         if (!ynq("Another plot")) exit(0);
  128.     }
  129. }
  130. /*------------------------------------------------------------------*/
  131. init_facs() {
  132.     int i;
  133.     for (i=1; i <=NDISTS; i++)
  134.         Dfinmem[i][0] = NULL;
  135.     itof(Fpone,1);
  136.     itof(Fpzero,0);
  137.     itof(Fpmone,-1);
  138.     itof(Fphundred,100);
  139.     itof(Fpthousand,1000);
  140.     Nexttogo = 1;
  141. }
  142. /*------------------------------------------------------------------*/
  143.  
  144. #include "fio3.h"
  145.  
  146. /*------------------------------------------------------------------*/
  147. /* END OF FACS.C */
  148. /*------------------------------------------------------------------*/
  149. /*------------------------------------------------------------------*/
  150. fplot() {
  151.     int i, qmax, max, h, v, more, check, nofile, dstyle, atdist,
  152.         left_margin, bottom_margin, distcnt, logplot;
  153.     char dnames[128];
  154.     int line_type[10];
  155.     int curve_thickness[10];
  156.     set4baud(1, 300);
  157.     while(!ynq("Is paper ready in plotter (1200 BAUD)"));
  158.     init_plot();
  159.  
  160.     /* SELECT PLOT SIZE */
  161.     Pp[1] = "\t-------------\n";
  162.     fprintf(STDERR,"\n%s\t|  1  |  2  |\n%s\t|  3  |  4  |\n%s\n",
  163.         Pp[1],Pp[1],Pp[1]);
  164.     checkint("\nQuarter page position 1-4 (Default is full page)? ",
  165.         &i, Buf, 0, 4, 0);
  166.     init_quarter();
  167.     switch(i) {
  168.         case 0:
  169.             init_full();
  170.             break;
  171.         case 1:
  172.             Vmin = VHALF;
  173.             break;
  174.         case 2:
  175.             Hmin = HHALF;
  176.             Vmin = VHALF;
  177.             break;
  178.         case 3:
  179.             break;
  180.         case 4:
  181.             Hmin = HHALF;
  182.             break;
  183.     }
  184.     fix_origin();
  185.     
  186.     /* DEFAULT STYLE? */
  187.     dstyle = ynq("\nDefault style");
  188.  
  189.     /* CHANGE LINE THICKNESS OR CHARACTER SIZE? */
  190.     if (!dstyle AND !ynq("Use default character sizes and line thickness")) {
  191.         checkint("Curve thickness (default %d)? ",&Plot_thickness,
  192.             Buf, 1, 10, 1);
  193.         checkint("Offset (default %d)? ",&Offset,
  194.             Buf, 1, 3, 2);
  195.         checkint("Character size:\n\tfor axis labels (default %d)? ",&Lsize,
  196.             Buf, 1, 5, Lsize);
  197.         checkint("\tfor scale numbers (default %d)? ",&Nsize,
  198.             Buf, 1, 5, Nsize);
  199.         if (Lsize EQ Nsize) {
  200.             checkint("Character thickness (default %d)? ",
  201.                 &Cthick[Lsize], Buf, 1, 10, Cthick[Lsize]);
  202.         }
  203.         else {
  204.             checkint(
  205.             "Character thickness:\n\tfor axis labels (default %d)? ",
  206.                 &Cthick[Lsize], Buf, 1, 10, Cthick[Lsize]);
  207.             checkint("\tfor scale numbers (default %d)? ",
  208.                 &Cthick[Nsize], Buf, 1, 10, Cthick[Nsize]);
  209.         }
  210.     }
  211.     
  212.     /* ALLOW SPACE FOR LABELS, SCALES, AND AXES */
  213.  
  214.     /* RAISE BOTTOM TO ALLOW FOR 1 LINE OF LABEL, 1 BLANK LINE,
  215.         1 LINE OF SCALE, AND ONE BETLIN BELOW AXIS */
  216.     bottom_margin = (2 * Cheight[Lsize]) + Cheight[Nsize] + Cbetlin[Nsize];
  217.     Vmin += bottom_margin;
  218.  
  219.     /* SHIFT TO RIGHT TO ALLOW FOR:
  220.         1 HEIGHT FOR LABEL,
  221.         1 BLANK HEIGHT BETWEEN LABEL AND SCALE,
  222.         4 WIDTHS FOR SCALE,
  223.         0.5 BLANK WIDTHS BETWEEN SCALE AND AXIS. */
  224.     left_margin = (2 * Cheight[Lsize]) +
  225.         (4 * Cwidth[Lsize]) + (Cwidth[Lsize]/2);
  226.     Hmin += left_margin;
  227.  
  228.     /* GET LIST OF DISTRIBUTIONS TO BE PLOTTED */
  229.     distcnt = 0;
  230.     while (!distcnt)
  231.         distcnt = getwrds(
  232. "List on one line distributions to be plotted, separated by commas.\n",
  233.             dnames, Pp);
  234.  
  235.     /* GET CURVE THICKNESS AND LINE TYPES */
  236.     for (i=1; i <= distcnt; i++) {
  237.         line_type[i] = 0;
  238.         curve_thickness[i] = Plot_thickness;
  239.     }
  240.     if (ynq("Vary line type or curve thickness")) {
  241.         fprintf(STDERR,
  242. "Line types:\n\t0  Solid\n\t2  Broken, fine\n\t4  Broken, medium\n%s",
  243.             "\t6  Broken, coarse\n");
  244.         for (i=1; i <= distcnt; i++) {
  245.             fprintf(STDERR,"\nCurve %s:\n",Pp[i]);
  246.             checkint("\tLine type (default %d)? ",
  247.                 &line_type[i], Buf, 0, 8, 0);
  248.             checkint("\tThickness (default %d)? ",
  249.                 &curve