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.
cannam@0:     This file copyright 2005-2006 Christian Landone.
cannam@0:     All rights reserved.
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: }