◆ innerproduct()

template<class T >

T ts::innerproduct	(	const aff::ConstSeries< T > &	a,
		const aff::ConstSeries< T > &	b
	)

Calculate inner product of two series.

If both series are of unequal index range, they are padded with zeros.

The function evaluates

$c_k = \sum\limits_{l=F}^{L} a_l b_l,$

where the index ranges of the input series are $a_l: l \in [F_a,L_a]$ and $b_l: l \in [F_b,L_b].$ and $F=\textrm{max}(F_a,F_b)$ as well as $L=\textrm{min}(L_a,L_b).$

Parameters

a
b

Returns: c

Definition at line 66 of file innerproduct.h.

References ts::seife::first().

Referenced by rms().

     {
       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

Here is the call graph for this function:

Here is the caller graph for this function: