fed2_var_clr.c

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

    else if ( strcmp (argv[i], "-gamma2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c2->gamma));
       }
    }
    else if ( strcmp (argv[i], "-metric2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c2->metric));
       }
    }
    else if ( strcmp (argv[i], "-o5smax2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c2->macarg.stepmax));
       }
    }
    else if ( strcmp (argv[i], "-ftol2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c2->ftol));
       }
    }
    else if ( strcmp (argv[i], "-nl2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c2->NL));
       }
    }
    else if ( strcmp (argv[i], "-opt2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c2->opt));
       }
    }
    else if ( strcmp (argv[i], "-itmax2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c2->itmax));
       }
    }
    else if ( strcmp (argv[i], "-betap2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c2->betap));
       }
    }
    else if ( strcmp (argv[i], "-sdx2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c2->sdx));
       }
    }
    else if ( strcmp (argv[i], "-seed2") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%ld", &(c2->seed));
       }
    }
    else if ( strcmp (argv[i], "-b2") == 0 ) {
      if ( i + 3 >= argc ) { ERROR2;      }
      else {
	  cs *= sscanf(argv[++i], "%d", &(c2->betastyle)); 
	  cs *= sscanf(argv[++i], "%lf", &(c2->beta0)); 
	  cs *= sscanf(argv[++i], "%lf", &(c2->beta1)); 
      }
    }
    else if ( strcmp (argv[i], "-gamma3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c3->gamma));
       }
    }
    else if ( strcmp (argv[i], "-metric3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c3->metric));
       }
    }
    else if ( strcmp (argv[i], "-o5smax3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c3->macarg.stepmax));
       }
    }
    else if ( strcmp (argv[i], "-ftol3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c3->ftol));
       }
    }
    else if ( strcmp (argv[i], "-nl3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c3->NL));
       }
    }
    else if ( strcmp (argv[i], "-opt3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c3->opt));
       }
    }
    else if ( strcmp (argv[i], "-itmax3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%d", &(c3->itmax));
       }
    }
    else if ( strcmp (argv[i], "-betap3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c3->betap));
       }
    }
    else if ( strcmp (argv[i], "-sdx3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%lf", &(c3->sdx));
       }
    }
    else if ( strcmp (argv[i], "-seed3") == 0 ) {
      if ( i + 1 == argc ) { ERROR1; }
      else {
        cs *= sscanf(argv[++i], "%ld", &(c3->seed));
       }
    }
    else if ( strcmp (argv[i], "-b3") == 0 ) {
      if ( i + 3 >= argc ) { ERROR2;      }
      else {
	  cs *= sscanf(argv[++i], "%d", &(c3->betastyle)); 
	  cs *= sscanf(argv[++i], "%lf", &(c3->beta0)); 
	  cs *= sscanf(argv[++i], "%lf", &(c3->beta1)); 
      }
    }