comparison dsp/tempotracking/TempoTrack.cpp @ 46:4f1870dbab2c

* Add lots of debug output

 
 #include "maths/MathAliases.h"
 #include "maths/MathUtilities.h"
 
 #include <iostream>
+
+#include <cassert>
 
 
 #define RAY43VAL
 
 //////////////////////////////////////////////////////////////////////

     else
     {
 	numelem = tsig;
     }
 
+    std::cerr << "tempoMM: m_winLength = " << m_winLength << ", m_lagLength = " << m_lagLength << ", numelem = " << numelem << std::endl;
+
     for(i=1;i<m_lagLength-1;i++)
     {
 	//first and last output values are left intentionally as zero
 	for (a=1;a<=numelem;a++)
 	{

 	m_RCFConditioning->process( m_tempoScratch, m_smoothRCF);
 
 	if (tsig != 0) // i.e. in context dependent state
 	{	
 //     NOW FIND MAX INDEX OF ACFOUT
-    	for( i = 0; i < m_lagLength; i++)
-    	{
-			if( m_tempoScratch[ i ] > maxValRCF)
-			{
-	    		maxValRCF = m_tempoScratch[ i ];
-	    		maxIndexRCF = i;
-			}
- 	   }
+            for( i = 0; i < m_lagLength; i++)
+            {
+                if( m_tempoScratch[ i ] > maxValRCF)
+                {
+                    maxValRCF = m_tempoScratch[ i ];
+                    maxIndexRCF = i;
+                }
+            }
 	}
 	else // using rayleigh weighting
 	{
 		vector <vector<double> > rcfMat;
 	

 			beatPeriod = p_i;
 		}
 		
 		// now write the output
 		maxIndexRCF = static_cast<int>(beatPeriod);
-
 	}
 
 
     double locked = 5168.f / maxIndexRCF;
     if (locked >= 30 && locked <= 180) {

     }
 
     if( tsig == 0 )
 	tsig = 4;
 
+
+std::cerr << "tempoMM: maxIndexRCF = " << maxIndexRCF << std::endl;
 	
     if( tsig == 4 )
     {
 	pdPeaks = new double[ 4 ];
 	for( i = 0; i < 4; i++ ){ pdPeaks[ i ] = 0.0;}

 
 void TempoTrack::createPhaseExtractor(double *Filter, unsigned int winLength, double period, unsigned int fsp, unsigned int lastBeat)
 {	
     int p = (int)MathUtilities::round( period );
     int predictedOffset = 0;
+
+    std::cerr << "TempoTrack::createPhaseExtractor: period = " << period << ", p = " << p << std::endl;
+
+    assert(p < 10000);
 
     double* phaseScratch = new double[ p*2 ];
 
 	
     if( lastBeat != 0 )

 	{
 	    temp = phaseScratch[ i ];
 	    phaseScratch[ i ] = (temp - PhaseMin)/PhaseMax;
 	}
 
+        std::cerr << "predictedOffset = " << predictedOffset << std::endl;
+
 	unsigned int index = 0;
 	for(int i = p - ( predictedOffset - 1); i < p + ( p - predictedOffset) + 1; i++)
 	{
+            std::cerr << "assigning to filter index " << index << " (size = " << p*2 << ")" << std::endl;
 	    Filter[ index++ ] = phaseScratch[ i ];
 	}
     }
     else
     {

 			
 	    periodG[ TTLoopIndex ] = tempoMM( m_frameACF, GW, tsig ); 
 
 	    period = periodG[ TTLoopIndex ];
 
+            std::cerr << "TempoTrack::process(2): constFlag == " << constFlag << ", TTLoopIndex = " << TTLoopIndex << ", period from periodG = " << period << std::endl;
+
 	    createPhaseExtractor( PW, m_winLength, period, FSP, 0 ); 
 
 	    constFlag = 0;
 
 	}
 	else
 	{
 	    if( GW[ 0 ] != 0 )
 	    {
 		period = periodG[ TTLoopIndex ];
+
+                std::cerr << "TempoTrack::process(2): GW[0] == " << GW[0] << ", TTLoopIndex = " << TTLoopIndex << ", period from periodG = " << period << std::endl;
+
+                if (period > 10000) {
+                    std::cerr << "WARNING!  Highly implausible period value!" << std::endl;
+                    std::cerr << "periodG contains (of " << TTFrames << " frames): " << std::endl;
+                    for (int i = 0; i < TTLoopIndex + 3 && i < TTFrames; ++i) {
+                        std::cerr << i << " -> " << periodG[i] << std::endl;
+                    }
+                    std::cerr << "periodP contains (of " << TTFrames << " frames): " << std::endl;
+                    for (int i = 0; i < TTLoopIndex + 3 && i < TTFrames; ++i) {
+                        std::cerr << i << " -> " << periodP[i] << std::endl;
+                    }
+                }
+
 		createPhaseExtractor( PW, m_winLength, period, FSP, lastBeat ); 
 
 	    }
 	    else
 	    {
 		period = periodP[ TTLoopIndex ];
+
+                std::cerr << "TempoTrack::process(2): GW[0] == " << GW[0] << ", TTLoopIndex = " << TTLoopIndex << ", period from periodP = " << period << std::endl;
+
 		createPhaseExtractor( PW, m_winLength, period, FSP, 0 ); 
 	    }
 	}
 
 	alignment[ TTLoopIndex ] = phaseMM( m_rawDFFrame, PW, m_winLength, period );