wnintegration.h

Go to the documentation of this file.

// include guard
#ifndef TF_WNINTEGRATION_H_VERSION

#define TF_WNINTEGRATION_H_VERSION \
  "TF_WNINTEGRATION_H   V1.0"

#include "lisousi.h"

struct Model {
  double Vp;
  double Vs;
  double Qp;
  double Qs;
}; // class Model

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

class ExpCoefficients {
  protected:
    ExpCoefficients(const Model& m): 
      Mm(m) 
  { 
    // square of complex velocities according to recipe I
    // Forbriger and Friederich (2000) eq. (7)
    // for the calculations to be carried out, we will almost exclusively
    // require the squares of the complex frequencies
    Malphaq=Mm.Vp*Mm.Vp*(1.+IME/Mm.Qp);
    Mbetaq=Mm.Vs*Mm.Vs*(1.+IME/Mm.Qs);
  }
  public:
    virtual ~ExpCoefficients();
    double maxp() const;
    virtual TFourier::Tcoeff coeff(const double& p) const=0;
  protected:
    Model Mm;
    TFourier::Tcoeff Malphaq;
    TFourier::Tcoeff Mbetaq;
}; // class ExpCoefficients

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

class FSECZ: public ExpCoefficients {
  public:
    FSECZ(const Model& m): ExpCoefficients(m) { }
    virtual ~FSECZ();
    virtual TFourier::Tcoeff coeff(const double& p) const;
}; // class FSECZ

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

class HSEC: public ExpCoefficients {
  public:
    HSEC(const Model& m): ExpCoefficients(m) { }
    virtual ~HSEC();
    virtual TFourier::Tcoeff coeff(const double& p) const=0;
  protected:
    TFourier::Tcoeff rayleigh(const double& p) const;
    mutable TFourier::Tcoeff Ma, Mb, Maq, Mbq, Mpq;
}; // class HSEC

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

class HSECZ: public HSEC {
  public:
    HSECZ(const Model& m): HSEC(m) { }
    virtual ~HSECZ();
    virtual TFourier::Tcoeff coeff(const double& p) const;
}; // class HSECZ

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

class HSECR: public HSEC {
  public:
    HSECR(const Model& m): HSEC(m) { }
    virtual ~HSECR();
    virtual TFourier::Tcoeff coeff(const double& p) const;
}; // class HSECR

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

class Exco: public TFourier::Tspectrum {
  public:
    Exco(): Mdp(0.) { }
    Exco(const ExpCoefficients& ec, const IntegParam& par);
    double dp() const { return(Mdp); }
  private:
    double Mdp;
}; // class Exco

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

enum Ebasis {
  Fbessel0,
  Fbessel1,
  Fcos,
  Fsin 
}; // enum Ebasis

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

TFourier::Tcoeff wnintegration(const Exco& ec,
                               const double& f,
                               const double& offset,
                               const Ebasis& fb);

#endif // TF_WNINTEGRATION_H_VERSION (includeguard)

/* ----- END OF wnintegration.h ----- */