fe_var6_clr.c

快速傅立叶变换程序代码,学信号的同学,可要注意了

    else if ( strcmp (argv[i], "-gamma") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->gamma));
       }
    }
    else if ( strcmp (argv[i], "-metric") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->metric));
       }
    }
    else if ( strcmp (argv[i], "-log") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->writelog));
       }
    }
    else if ( strcmp (argv[i], "-logfile") == 0 ) {
       if ( i + 1 == argc ) { ERROR1; }
       else {
         strcpy(c->logfile, argv[++i]);
       }
    }
    else if ( strcmp (argv[i], "-o5smax") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c->macarg.stepmax));
       }
    }
    else if ( strcmp (argv[i], "-second") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->second));
       }
    }
    else if ( strcmp (argv[i], "-third") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->third));
       }
    }
    else if ( strcmp (argv[i], "-sdx") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c->sdx));
       }
    }
    else if ( strcmp (argv[i], "-fv") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->verbose));
       }
    }
    else if ( strcmp (argv[i], "-sv") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->seqverbose));
       }
    }
    else if ( strcmp (argv[i], "-fcg") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->CG));
       }
    }
    else if ( strcmp (argv[i], "-oappend") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->outappend));
       }
    }
    else if ( strcmp (argv[i], "-o") == 0 ) {
       if ( i + 1 == argc ) { ERROR1; }
       else {
         strcpy(c->outfile, argv[++i]);
         if(!(c->printout))c->printout=1;
       }
    }
    else if ( strcmp (argv[i], "-to") == 0 ) {
       if ( i + 1 == argc ) { ERROR1; }
       else {
         strcpy(c->totoutfile, argv[++i]);
         if(!(c->printtot))c->printtot=1;
       }
    }
    else if ( strcmp (argv[i], "-ho") == 0 ) {
       if ( i + 1 == argc ) { ERROR1; }
       else {
         strcpy(c->histofile, argv[++i]);
         if(!(c->printhisto))c->printhisto=1;
       }
    }
    else if ( strcmp (argv[i], "-fseed") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%ld", &(c->seed));
       }
    }
    else if ( strcmp (argv[i], "-e") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c->epsilon));
       }
    }
    else if ( strcmp (argv[i], "-fdemo") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->DEMO));
       }
    }
    else if ( strcmp (argv[i], "-norder") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->NORDER));
       }
    }
    else if ( strcmp (argv[i], "-betap") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c->betap));
       }
    }
    else if ( strcmp (argv[i], "-nl") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->NL));
       }
    }
    else if ( strcmp (argv[i], "-rich") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->rich));
       }
    }
    else if ( strcmp (argv[i], "-opt") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->opt));
       }
    }
    else if ( strcmp (argv[i], "-itmax") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->itmax));
       }
    }
    else if ( strcmp (argv[i], "-ftol") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c->ftol));
       }
    }
    else if ( strcmp (argv[i], "-flintol") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c->flintol));
       }
    }
    else if ( strcmp (argv[i], "-seqp") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c->seq_period));
       }
    }