KM3::TOALIB Namespace Reference

Functions
void	vectorMultiply (complex_type *_vectorSrcA, complex_type *_vectorSrcB, complex_type *_vectorDest, size_t _length)

Function Documentation

vectorMultiply()

void KM3::TOALIB::vectorMultiply	(	complex_type *	_vectorSrcA,
		complex_type *	_vectorSrcB,
		complex_type *	_vectorDest,
		size_t	_length )

Definition at line 13 of file km3_toa_utils.cc.

{
 
        // TODO: there is space for optimization
        // for instance: how do we manage denormals?
        // do we use SSE/2?
        // so: search for a ready-to-use lib and use it
 
        /*
         * given two complex numbers: a = a_real + j*a_imag, b = b_real + j*b_imag
         * the product c = a*b can be calculated as:
         * c = (d1-d2) + j*(d1+d3)
         * with:
         * d1 = a_real*(b_real+b_imag)
         * d2 = b_imag*(a_real+b_real)
         * d3 = b_real*(a_imag-a_real)
         *
         * see: http://www.embedded.com/design/real-time-and-performance/4007256/Digital-Signal-Processing-Tricks--Fast-multiplication-of-complex-numbers
         */
 
        real_type a_real, b_real, a_imag, b_imag;
        real_type k1, k2, k3;
        for(unsigned int k = 0; k < _length; k++) {
                a_real = _vectorSrcA[k][0];
                a_imag = _vectorSrcA[k][1];
                b_real = _vectorSrcB[k][0];
                b_imag = _vectorSrcB[k][1];
 
                k1 = a_real*(b_real+b_imag);
                k2 = b_imag*(a_real+a_imag);
                k3 = b_real*(a_imag-a_real);
 
                _vectorDest[k][0] = k1-k2;
                _vectorDest[k][1] = k1+k3;
        }
 
}