home *** CD-ROM | disk | FTP | other *** search
/ Programmer 7500 / MAX_PROGRAMMERS.iso / INFO / C / FFT142.ZIP / RI2M.C < prev    next >
Encoding:
C/C++ Source or Header  |  1991-06-06  |  2.9 KB  |  134 lines
  1. /*    name...
  2.         ri2mag - convert frequencies with real & imaginary parts to frequencies
  3.                 with magnitudes
  4.     bugs...
  5.         should read in all data, then write, so disk head doesn't have
  6.         to move between input and output file.
  7. */
  8. #include <stdio.h>
  9. #include <math.h>
  10.  
  11. #define VERSION "1.1"
  12. #define MAXLABELS 50
  13. #define BUFSIZE 200
  14.  
  15. FILE *ifile=stdin, *ofile=stdout;
  16.  
  17. char buffer[128], oname[35];
  18. char buf[BUFSIZE];
  19.  
  20. int automatic_abscissas, /* nonzero if abscissas to be calculated */
  21. abscissa_arguments;
  22.  
  23. double
  24. abscissa,
  25. abscissa_step;
  26.  
  27.  
  28. main(argc, argv) int argc; char **argv;
  29. {    int i;
  30.  
  31.     double amp, phase, factor, q, pi, time, dt;
  32.     double this, big, sweep, scale;
  33.  
  34.     {double junk;
  35.     sscanf("3.4","%lf",&junk);        /* workaround for DESMET sscanf bug */
  36.     }
  37.  
  38.     argc = args(argc, argv);
  39.     read_data(argc, argv);
  40.  
  41.     if(ofile!=stdout)
  42.         {fclose(ofile);
  43.         }
  44. }
  45.  
  46.  
  47.  
  48. read_data(argc, argv) int argc; char **argv;
  49. {    int i, j, nums, ac;
  50.     double ff, xx, yy;
  51.     char *s, *t, *stop, **av;
  52.     char *strchr();
  53.  
  54.     argc--; argv++;
  55.     if(strchr(argv[0], '?')) help();
  56.  
  57. /*            open input file         */
  58.     if(argc)
  59.         {ifile=fopen(argv[0], "r");
  60.         if(ifile==0) {printf("file %s not found\n", argv[0]); exit(1);}
  61.         argc--; argv++;
  62.         }
  63. /*            open output file         */
  64.     if(argc)
  65.         {strcpy(oname, argv[0]);
  66.         unlink(oname);
  67.         if((ofile=fopen(oname, "w"))==0)
  68.             {fprintf(stderr, "can\'t open output file");
  69.             exit(1);
  70.             }        
  71.         }
  72.  
  73.     fprint_cmd(stdout, "; ri2m%s\n");
  74.  
  75.     while(fgets(buf, BUFSIZE, ifile))
  76.         {t=buf;
  77.         while(*t && isspace(*t)) t++;
  78.         if(*t == 0) continue;        /* skip blank lines */
  79.         buf[strlen(buf)-1]=0;        /* zero out the line feed */
  80.         if(buf[0]==';')                /* copy comment */
  81.             {fprintf(ofile, "%s\n", buf); continue;
  82.             }
  83.         if(t = strchr(buf, ';')) *t=0;    /* zap same-line comment */
  84.         if(automatic_abscissas)
  85.             {ff=abscissa;
  86.             abscissa+=abscissa_step;
  87.             sscanf(buf, "%lf %lf", &xx, &yy);
  88.             }
  89.         else
  90.             {sscanf(buf, "%lf %lf %lf", &ff, &xx, &yy);
  91.             }
  92.         s=buf;                      /* start looking for label */
  93.         nums=3;
  94.         if(automatic_abscissas) nums--;
  95.         while(nums--)                    /* skip the numbers */
  96.             {while(*s==' ')s++;
  97.             while(*s && (*s!=' '))s++;
  98.             }
  99.         while(*s==' ')s++;
  100.         fprintf(ofile, "%15.8g %15.8g %s\n", ff, sqrt(xx*xx + yy*yy), s);
  101.         }
  102. }
  103.  
  104. int streq(a,b) char *a,*b;
  105. {    while(*a)
  106.         {if(*a!=*b)return 0;
  107.         a++; b++;
  108.         }
  109.     return (*b==0);
  110. }
  111.  
  112. /* get_parameter - process one command line option
  113.         (return # parameters used) */
  114. get_parameter(argc, argv) int argc; char **argv;
  115. {    int i;
  116.     if(streq(*argv, "-a"))
  117.         {i=get_double(argc, argv, 2, &abscissa_step, &abscissa, &abscissa);
  118.         abscissa_arguments=i-1;
  119.         automatic_abscissas=1;
  120.         return i;
  121.         }
  122.     else gripe(argv);
  123. }
  124.  
  125. char *message[]={
  126. "ri2m   version ", VERSION,
  127. " - convert real & imag parts to magnitude\n",
  128. "usage: ri2m  infile  [outfile]  [options]\n",
  129. "options are:\n",
  130. "  -a  [step [start]] automatic abscissas \n",
  131. 0
  132. };
  133.  
  134.