main()

int main	(	int	iargc,
		char *	argv[]
	)

Definition at line 64 of file autocorr.cc.

References AUTOCORR_VERSION, Options::debug, Options::overwrite, and Options::verbose.

{

  // define usage information
  char usage_text[]=
  {
    AUTOCORR_VERSION "\n"
    "usage: autocorr -v -o -D file [file [...]] outfile" "\n"
    "   or: autocorr --help|-h" "\n"
  };

  // define full help text
  char help_text[]=
  {
    "-v           be verbose" "\n"
    "-D           debug mode" "\n"
    "-o           overwrite existing output file" "\n"
    "file ...     input file(s)" "\n"
    "outfile      output file" "\n"
  };

  // define commandline options
  using namespace tfxx::cmdline;
  static Declare options[]= 
  {
    // 0: print help
    {"help",arg_no,"-"},
    // 1: verbose mode
    {"v",arg_no,"-"},
    // 2: overwrite mode
    {"o",arg_no,"-"},
    // 3: debug mode
    {"D",arg_no,"-"},
    {NULL}
  };

  // no arguments? print usage...
  if (iargc<2) 
  {
    cerr << usage_text << endl;
    exit(0);
  }

  // collect options from commandline
  Commandline cmdline(iargc, argv, options);

  // help requested? print full help text...
  if (cmdline.optset(0))
  {
    cerr << usage_text << endl;
    cerr << help_text << endl;
    exit(0);
  }

  /*
  // dummy operation: print option settings
  for (int iopt=0; iopt<2; iopt++)
  {
    cout << "option: '" << options[iopt].opt_string << "'" << endl;
    if (cmdline.optset(iopt)) {  cout << "  option was set"; }
    else { cout << "option was not set"; }
    cout << endl;
    cout << "  argument (string): '" << cmdline.string_arg(iopt) << "'" << endl;
    cout << "     argument (int): '" << cmdline.int_arg(iopt) << "'" << endl;
    cout << "    argument (long): '" << cmdline.long_arg(iopt) << "'" << endl;
    cout << "   argument (float): '" << cmdline.float_arg(iopt) << "'" << endl;
    cout << "  argument (double): '" << cmdline.double_arg(iopt) << "'" << endl;
    cout << "    argument (bool): '";
    if (cmdline.bool_arg(iopt))
    { cout << "true"; } else { cout << "false"; }
    cout << "'" << endl;
  }
  while (cmdline.extra()) { cout << cmdline.next() << endl; }

  // dummy operation: print rest of command line
  while (cmdline.extra()) { cout << cmdline.next() << endl; }
  */

  Options opt;
  opt.verbose=cmdline.optset(1);
  opt.overwrite=cmdline.optset(2);
  opt.debug=cmdline.optset(3);

  TFXX_assert(cmdline.extra(), "missing input file!");
  std::string outfile=cmdline.next();
  TFXX_assert(cmdline.extra(), "missing output file!");
  Tnames infiles;
  while (cmdline.extra())
  {
    infiles.push_back(outfile);
    outfile=cmdline.next();
  }

  if (opt.verbose) { cout << "open output file " << outfile << endl; }
  std::ofstream ofs(outfile.c_str());
  datrw::osffstream os(ofs);

  sff::SRCE srce=sff::srce_reference();
  os << srce;

  Tnames::const_iterator I=infiles.begin();
  while (I!=infiles.end())
  {
    if (opt.verbose) { cout << "open " << std::string(*I) << endl; }
    std::ifstream ifs(I->c_str());
    datrw::isffstream is(ifs, opt.debug);
    while (is.good())
    {
      if (opt.verbose) { cout << "read trace" << endl; }
      Tts data;
      Tts::Tseries samples;
      is >> samples >> data.header;
      data=samples;
      sff::INFO info;
      is >> info;
      if (opt.verbose) { cout << data.header.line() << endl; }
      Tts result;
      result=ts::correlate(data,data);
      result /= Tts::Tvalue(data.size());
      result.header=data.header;
      result.header.auxid="corr";
      result.header.date=srce.date
        +(result.first()*libtime::double2time(result.header.dt));
      if (is.hasinfo()) { os << info; }
      os << result;
    }
    ++I;
  }
}