convolve.h

Go to the documentation of this file.

// include guard
#ifndef TF_CONVOLVE_H_VERSION

#define TF_CONVOLVE_H_VERSION \
  "TF_CONVOLVE_H   V1.0"

#include<aff/series.h>

namespace ts {

  template<class T>
    aff::Series<T> convolve(const aff::ConstSeries<T>& a, 
                            const aff::ConstSeries<T>& b)
    {
      aff::Series<T> retval(a.first()+b.first(),b.last()+a.last());
      for (long int k=retval.first(); k<=retval.last(); ++k)
      {
        retval(k)=T(0);
        long int lmin=a.first() > (k-b.last()) ? a.first() : (k-b.last());
        long int lmax=a.last() < (k-b.first()) ? a.last() : (k-b.first());
        for(long int l=lmin; l<=lmax; ++l)
        {
          retval(k)+=a(l)*b(k-l);
        }
      }
      return(retval);
    } // convolve

} // namespace ts

#endif // TF_CONVOLVE_H_VERSION (includeguard)

/* ----- END OF convolve.h ----- */