cannam@0: /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
cannam@0: 
cannam@0: /*
cannam@0:     QM DSP Library
cannam@0: 
cannam@0:     Centre for Digital Music, Queen Mary, University of London.
Chris@84:     This file 2005-2006 Christian Landone.
Chris@84: 
Chris@84:     This program is free software; you can redistribute it and/or
Chris@84:     modify it under the terms of the GNU General Public License as
Chris@84:     published by the Free Software Foundation; either version 2 of the
Chris@84:     License, or (at your option) any later version.  See the file
Chris@84:     COPYING included with this distribution for more information.
cannam@0: */
cannam@0: 
cannam@0: #include "Correlation.h"
cannam@0: 
cannam@0: //////////////////////////////////////////////////////////////////////
cannam@0: // Construction/Destruction
cannam@0: //////////////////////////////////////////////////////////////////////
cannam@0: 
cannam@0: Correlation::Correlation()
cannam@0: {
cannam@0: 
cannam@0: }
cannam@0: 
cannam@0: Correlation::~Correlation()
cannam@0: {
cannam@0: 
cannam@0: }
cannam@0: 
cannam@0: void Correlation::doAutoUnBiased(double *src, double *dst, unsigned int length)
cannam@0: {
cannam@0:     double tmp = 0.0;
cannam@0:     double outVal = 0.0;
cannam@0: 
cannam@0:     unsigned int i,j;
cannam@0: 
cannam@0:     for( i = 0; i <  length; i++)
cannam@0:     {
cannam@0: 	for( j = i; j < length; j++)
cannam@0: 	{
cannam@0: 	    tmp += src[ j-i ] * src[ j ]; 
cannam@0: 	}
cannam@0: 
cannam@0: 
cannam@0: 	outVal = tmp / ( length - i );
cannam@0: 
cannam@0: 	if( outVal <= 0 )
cannam@0: 	    dst[ i ] = EPS;
cannam@0: 	else
cannam@0: 	    dst[ i ] = outVal;
cannam@0: 
cannam@0: 	tmp = 0.0;
cannam@0:     }
cannam@0: }