lisousi.h

Go to the documentation of this file.

// include guard
#ifndef TF_LISOUSI_H_VERSION

#define TF_LISOUSI_H_VERSION \
  "TF_LISOUSI_H   V1.0"

#include <iostream>
#include <fstream>
#include <string>
#include <sstream>
#include <gsl/gsl_sf_bessel.h>
#include <tfxx/commandline.h>
#include <tfxx/xcmdline.h>
#include <tfxx/error.h>
#include <tfxx/stringfunc.h>
#include <tfxx/rangestring.h>
#include <tfxx/rangelist.h>
#include <tfxx/misc.h>
#include <tsxx/anyfilter.h>
#include <tsxx/convolve.h>
#include <datrwxx/readany.h>
#include <datrwxx/writeany.h>
#include <sffxx.h>
#include <sffostream.h>
#include <aff/dump.h>
#include <aff/seriesoperators.h>
#include <aff/subarray.h>
#include <tsxx/ovtaper.h>
#include <fourier/fftwaff.h>

using std::cout;
using std::cerr;
using std::endl;

typedef ts::filter::Ttimeseries Ttimeseries;
typedef Ttimeseries::Tseries Tseries;
typedef fourier::fft::DRFFTWAFF TFourier;

// the global Fourier processor will be instantiated in the main programmin
// unit
extern TFourier Fourier;

// global workspace and filters
extern Tseries filter, taper, workseries;

/*----------------------------------------------------------------------*/
// constants
  
// imaginary unit
const TFourier::Tcoeff IME=TFourier::Tcoeff(0.,1.);
// my Fourier transformation constant factor
const TFourier::Tcoeff CFTfac=exp(-IME*M_PI_4);

/*----------------------------------------------------------------------*/
// type of Fourier domain solution
enum Efdtype {
  Ffdfarfield,
  Ffdexplosion,
  Ffdzforce,
  Ffdwizforce,
  Ffdlamb
};

/*----------------------------------------------------------------------*/

struct IntegParam {
  int nsteps;
  double edge;
  double taper;
}; // struct IntegParam

/*----------------------------------------------------------------------*/

struct Options {
  bool verbose, debug;
  bool overwrite, taperfirst, sqrttaper, limitlength, radial;
  bool fredomain, fdfilter, tdfilter, nointeg, transition, tapsloset;
  bool spatialdistance;
  double tfac, tshift, tlim, integshift, tapdel;
  double tapslo, transition1, transition2;
  int npad;
  std::string inputformat, outputformat;
  // single-velocity switches:
  Efdtype fdtype;
  // propagation parameters:
  double velocity, vpvsratio, pquality, squality;
  // trapezoid integration:
  IntegParam ipa;
}; // struct Options

struct Parameters {
  Parameters():
    offsetfactor(0.), offset(0.), t0(0.), T(0.), dt(0.) { }
  double offsetfactor;
  double offset;
  double t0;
  double T;
  double dt;
  int nsamples;
}; // struct Parameters

#endif // TF_LISOUSI_H_VERSION (includeguard)

/* ----- END OF lisousi.h ----- */