innerproduct.h

Go to the documentation of this file.

// include guard
#ifndef TF_INNERPRODUCT_H_VERSION

#define TF_INNERPRODUCT_H_VERSION \
  "TF_INNERPRODUCT_H   V1.0   "

#include<aff/series.h>

namespace ts {

  template<class T>
    T innerproduct(const aff::ConstSeries<T>& a, 
                   const aff::ConstSeries<T>& b)
    {
      T retval=0;
      long int first=(a.first() > b.first() ? a.first() : b.first());
      long int last=(a.last() < b.last() ? a.last() : b.last());
      if (last >= first) 
      {
        for (int i=first; i<=last; i++)
        {
          retval += a(i)*b(i);
        }
      }
      return(retval);
    } // innerproduct

  template<class T>
    T rms(const aff::ConstSeries<T>& a) 
    { return(sqrt(innerproduct(a,a)/a.size())); }

} // namespace ts

#endif // TF_INNERPRODUCT_H_VERSION (includeguard)

/* ----- END OF innerproduct.h ----- */