tapers.h

Go to the documentation of this file.

// include guard
#ifndef TF_TAPERS_H_VERSION

#define TF_TAPERS_H_VERSION \
  "TF_TAPERS_H   V1.3"

namespace ts {

  namespace tapers {

    class Taper {
      public:
        virtual ~Taper() { }
        template<class C>
          void apply(C c) const;
        template<class C>
          void psdapply(C c) const;
        virtual double psdnorm() const =0;
      private:
        virtual void init(const int& f, const int& l) const =0;
        virtual double value(const int& i) const =0;
    }; // class Taper

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

    class Cosine: public Taper {
      public:
        Cosine(const double& f=0.1);
        ~Cosine() { }
      private:
        void init(const int& f, const int& l) const;
        double value(const int& i) const;
        double psdnorm() const;
        mutable int Mf,Ml; //<! index of first and last sample
        mutable int Msf; //<! sample index fraction
        mutable double Mfac; //<! sine function argument scaling
        double Mfrac; //<! taper fraction
    }; // class Hanning

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

    class Hanning: public Taper {
      public:
        ~Hanning() { }
      private:
        void init(const int& f, const int& l) const;
        double value(const int& i) const;
        double psdnorm() const;
        mutable int Mf; //<! first index
        mutable double Mfac; //<! sine function argument scaling
    }; // class Hanning

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

    class FourPoint: public Taper {
      public:
        FourPoint(const double& t1=0.,
                  const double& t2=0.1,
                  const double& t3=0.9,
                  const double& t4=1.);
        ~FourPoint() { }
      private:
        void init(const int& f, const int& l) const;
        double value(const int& i) const;
        double psdnorm() const;
        double Mt1, Mt2, Mt3, Mt4;
        mutable double Mti1, Mti2, Mti3, Mti4;
        mutable double Mfac1, Mfac2;
    }; // class FourPoint

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

    template<class C>
      void Taper::apply(C c) const
      {
        this->init(c.f(), c.l());
        for (int i=c.f(); i<=c.l(); ++i)
        { c(i) *= this->value(i); }
      } // template<class C> void Taper::apply(C c)

    template<class C>
      void Taper::psdapply(C c) const
      {
        this->init(c.f(), c.l());
        for (int i=c.f(); i<=c.l(); ++i)
        { c(i) *= this->value(i)/this->psdnorm(); }
      } // template<class C> void Taper::psdapply(C c)
    
  } // namespace tapers

} // namespace ts

#endif // TF_TAPERS_H_VERSION (includeguard)

/* ----- END OF tapers.h ----- */