home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #26 / NN_1992_26.iso / spool / comp / graphics / gnuplot / 364 < prev    next >
Encoding:
Text File  |  1992-11-12  |  3.8 KB  |  123 lines
  1. Newsgroups: comp.graphics.gnuplot
  2. Path: sparky!uunet!ukma!wupost!bcm!convex!seas.smu.edu!mustafa
  3. From: mustafa@seas.smu.edu (Mustafa Kocaturk)
  4. Subject: Histograms in Gnuplot
  5. Message-ID: <1992Nov13.014817.3123@seas.smu.edu>
  6. Summary: C filter for producing a histogram in Gnuplot
  7. Keywords: histogram insertion sort filter
  8. Sender: news@seas.smu.edu (USENET News System)
  9. Nntp-Posting-Host: turbo_f.seas.smu.edu
  10. Organization: SMU School Of Engineering and Applied Science
  11. References: <9211121849.AA13763@aulne.lifl.fr>
  12. Date: Fri, 13 Nov 1992 01:48:17 GMT
  13. Lines: 108
  14.  
  15. In article <9211121849.AA13763@aulne.lifl.fr> hemery@lifl.fr writes:
  16. >I am looking for tool like gnuplot, for drawing histogram graphics.
  17. >
  18. >Do you know where I can find this kind of tool ?
  19.  
  20. The short program below is a filter that calculates a histogram
  21. from a sequence of numbers and prints the output in such a format
  22. that it can be plotted in Gnuplot by the command sequence
  23.  
  24.     !histogram < datain > tmp
  25.     plot "tmp" with impulses
  26.  
  27. where
  28.     `histogram' is the program's name
  29.     `datain' is a file containing input data samples
  30.     `tmp' is a file to hold histogram's output
  31.  
  32. The number of points generated is 100 by default. A different
  33. number can be given as the first command line argument.
  34. Supplying a non-integer argument will cause a help message to be
  35. printed to the standard error.
  36.  
  37. Another usage example is
  38.  
  39.      !myprog | histogram 50 > tmp
  40.      plot "tmp" with impulses
  41.  
  42. The maximum number of input data points that can be handled
  43. was chosen arbitrarily to be 50000.
  44.  
  45. Yours respectfully,
  46.  
  47. Mustafa Kocaturk     mustafa@seas.smu.edu
  48.  
  49. /* ------- cut here and save the lines below as histogram.c --------- */
  50. /*
  51. This source code is placed in the public domain with no guarantee
  52. of fitness for a particular purpose or of portability, but
  53. only with the hope that it may be useful.
  54.  
  55. I wrote `histogram.c' as a part of my effort to contribute
  56. to a paper by two of my colleagues at SMU. The paper
  57. discussed various programming environments such as Gnuplot,
  58. awk, spice, and xfig.
  59.  
  60. Mustafa Kocaturk    mustafa@seas.smu.edu  EE Dept., SEAS, SMU
  61. */
  62. #include <stdio.h>
  63. #include <math.h>
  64. #define HELP \
  65. "This program calculates the histogram of a sequence of numbers\n\
  66. read from the standard input until an end of file is encountered.\n\
  67. The results are printed on the standard output as lines of the format\n\
  68. <subinterval midpoint> <number of sample points in subinterval>\n\
  69. \n\
  70. Usage: %s [number of subintervals (default==100)]\n" /* Message string */
  71.  
  72. #define CALSIZ    50000 /* Number of doubles to be allocated for sort */
  73.  
  74. main(argc,argv)
  75.  int argc;
  76.  char **argv;
  77.  {
  78.  double x,t,z,*a,*b,*f,*fe,*fl;
  79.  unsigned i=0,k=100;
  80.  int j=1;
  81.  char c;
  82.  if(argc>1)j=sscanf(argv[1],"%u",&k);
  83.  if(j!=1){fprintf(stderr,HELP,argv[0]);return 0;}
  84.  if(NULL==(a=(double *)calloc(CALSIZ,sizeof x)))
  85.   {
  86.   fprintf(stderr,"%s: Unable to allocate memory\n",(argv[0]));
  87.   return 0; /* let other work continue */
  88.   }
  89.  fl=a+CALSIZ-1;
  90.  fe=a;
  91.  b=a+1;
  92.  for(i=0;j>0&&fe<fl;)
  93.   {
  94.   if((j=scanf("%lf%*[^\n]%*c",&x))>0) /* If number was recognized */
  95.    {
  96.    f=fe++;
  97.    while(x<*f&&f>a) { f[1]=*f; f--; }
  98.    f[1]=x;
  99.    i++;
  100.    }
  101.   else
  102.    j=scanf("%*[^\n]%*c");
  103.    /* Read data until end-of-line,
  104.    if end-of-stdin (i. e. ctrl-D) then j<0 (i. e. break the for loop) */
  105.   }
  106.  x=*b; /* This was Insertion sort (increasing) */
  107.  t=*fe;
  108.  printf("# %d points in [%lg, %lg]; Frequencies in %d subintervals:\n",
  109.   i,x,t,k);
  110.  t=(t-x)/k;
  111.  x+=t;
  112.  for(f=b,i=0;i<k;i++,x+=t)
  113.   {
  114.   for(j=0;x>*f&&f<=fe;j++,f++);/* Count points in the subinterval */
  115.   printf(" %l g %d\n",x-t*.5,j); /* Print subinterval center and
  116.             number of points in the subinterval to stdout */
  117.   }
  118.  }
  119. /* ------- cut here and save the lines above as histogram.c --------- */
  120. -- 
  121. Mustafa KOCATURK    mustafa@seas.smu.edu  EE Dept., Room 305A, Caruth Bldg.
  122. Home: 214-706-5954  Office: 214-768-1475  SMU Box 753190, Dallas, TX 75275
  123.