stfinvfourier.h

Go to the documentation of this file.

// include guard
#ifndef STFINV_STFINVFOURIER_H_VERSION

#define STFINV_STFINVFOURIER_H_VERSION \
  "STFINV_STFINVFOURIER_H   V1.4"

#include <stfinv/stfinvbase.h>
#include <aff/array.h>
#include <fourier/fftwaffar.h>

namespace stfinv {

  class STFFourierDomainEngine: public stfinv::STFBaseEngine {
    public:
      typedef stfinv::STFBaseEngine Tbase;
      typedef fourier::fft::DRFFTWAFFArrayEngine Tfftengine;
      typedef Tfftengine::TAseries TAseries;
      typedef Tfftengine::TAspectrum TAspectrum;
    protected:
      STFFourierDomainEngine(const stfinv::Tvectoroftriples& triples,
                             const stfinv::Waveform& stf,
                             const std::string& parameters);
      STFFourierDomainEngine(const stfinv::Tvectoroftriples& triples,
                             const stfinv::Waveform& stf,
                             const stfinv::Tvectorofpairs& pairs,
                             const std::string& parameters);
    public:
      virtual ~STFFourierDomainEngine() { }
      virtual void exec() { STFINV_baseillegal; }
      virtual void help(std::ostream& os=std::cout) const;
      static void classhelp(std::ostream& os=std::cout);
      virtual void usage(std::ostream& os=std::cout) const;
      static void classusage(std::ostream& os=std::cout);
      virtual const char* name() const;
    protected:
      void fftinput();
      void fftoutput();
      TAspectrum recordingspec() const;
      TAspectrum syntheticspec() const;
      TAspectrum stfspec() const;
      TAspectrum recordingcoeff(const unsigned int& i) const;
      TAspectrum syntheticcoeff(const unsigned int& i) const;
      TAspectrum::Tvalue& stfcoeff(const unsigned int& i) const;
      double frequency(const unsigned int& i) const;
      unsigned int nfreq() const;
    private:
      void initialize();
      void getinput();
      void putoutput();
      void convolve();
      void stfshift();
      void taperstf();
      TAseries stfseries() const;
      
    // member data
    // -----------
    private:
      Tfftengine Mfftengineinput;
      Tfftengine Mfftengineoutput;
      Tfftengine Mfftenginestf;
      double Mtshift;
      bool Mapplyshift;
      bool Mapplystftaper;
      double Mtt1, Mtt2, Mtt3, Mtt4;
  }; // class STFFourierDomainEngine

}

#endif // STFINV_STFINVFOURIER_H_VERSION (includeguard)

/* ----- END OF stfinvfourier.h ----- */