robert@464: /*
robert@464:  * spear_parser.cpp v1.2
robert@464:  *
robert@464:  *  Created on: May 6, 2014
robert@464:  *      Author: Victor Zappi
robert@464:  */
robert@464: 
robert@464: #include "spear_parser.h"
robert@464: 
robert@464: using namespace std;
robert@464: 
robert@464: //#define DO_CHECKS
robert@464: 
robert@464: //------------------------------------------------------------------------------------------------
robert@464: // partials
robert@464: //------------------------------------------------------------------------------------------------
robert@464: 
robert@464: Partials::Partials()
robert@464: {
robert@464: 	partialFrequencies	= NULL;
robert@464: //	partialAmplitudes	= NULL;
robert@464: //	partialNumFrames	= NULL;
robert@464: //	partialStartSample	= NULL;
robert@464: //	partialEndSample	= NULL;
robert@464: //	partialCurrentFrame	= NULL;
robert@464: //	partialFreqDelta	= NULL;
robert@464: //	partialAmpDelta	= NULL;
robert@464: 
robert@464: 
robert@464: 	activePartialNum	= NULL;
robert@464: //	activePartials		= NULL;
robert@464: 
robert@464: 	currentSample = -1;
robert@464: }
robert@464: 
robert@464: Partials::~Partials()
robert@464: {
robert@464: 	if(partialFrequencies != NULL)			// check on one is enough
robert@464: 	{
robert@464: 		if(partialFrequencies[0] != NULL)	// check on one is enough
robert@464: 		{
robert@464: 			for(unsigned int i=0; i<parNum; i++)
robert@464: 			{
robert@464: 				delete[] partialFrequencies[i];
robert@464: 				delete[] partialAmplitudes[i];
robert@464: 				delete[] partialFreqDelta[i];
robert@464: 				delete[] partialAmpDelta[i];
robert@464: 
robert@464: 			}
robert@464: 		}
robert@464: 
robert@464: 		delete[] partialFrequencies;
robert@464: 		delete[] partialAmplitudes;
robert@464: 		delete[] partialNumFrames;
robert@464: 		delete[] partialFreqDelta;
robert@464: 		delete[] partialAmpDelta;
robert@464: 		delete[] partialFreqMean;
robert@464: 	}
robert@464: 
robert@464: 	if(activePartialNum != NULL)
robert@464: 	{
robert@464: 		for(unsigned int i=0; i<hopNum+1; i++)
robert@464: 			delete[] activePartials[i];
robert@464: 
robert@464: 		delete[] activePartialNum;
robert@464: 		delete[] activePartials ;
robert@464: 	}
robert@464: }
robert@464: 
robert@464: void Partials::init(int parN, int hopS, bool isDBX)
robert@464: {
robert@464: 	if(!isDBX)
robert@464: 	{
robert@464: 		parNum	= parN;
robert@464: 		hopSize	= hopS;
robert@464: 
robert@464: 		partialFrequencies	= new float *[parNum];
robert@464: 		partialAmplitudes	= new float *[parNum];
robert@464: 		partialNumFrames	= new unsigned int[parNum];
robert@464: 		partialStartFrame	= new unsigned int[parNum];
robert@464: 		partialStartSample	= new unsigned int[parNum];
robert@464: 		partialEndSample	= new unsigned int[parNum];
robert@464: 		partialFreqDelta	= new float *[parNum];
robert@464: 		partialAmpDelta		= new float *[parNum];
robert@464: 		partialFreqMean		= new float[parNum];
robert@464: 
robert@464: 
robert@464: 
robert@464: 		// init in one shot
robert@464: 		fill(partialFreqMean, partialFreqMean+parNum, 0);			// mean is zero
robert@464: 
robert@464: 		partialFrequencies[0] 	= NULL;								// for free check
robert@464: 	}
robert@464: 	else
robert@464: 	{
robert@464: 		parNum	= parN;
robert@464: 		hopSize	= hopS;
robert@464: 
robert@464: 		partialFrequencies	= new float *[parNum];
robert@464: 		partialAmplitudes	= new float *[parNum];
robert@464: 		partialNumFrames	= new unsigned int[parNum];
robert@464: 		partialStartFrame	= new unsigned int[parNum];
robert@464: 		partialFreqDelta	= new float *[parNum];
robert@464: 		partialAmpDelta		= new float *[parNum];
robert@464: 		partialFreqMean		= new float[parNum];
robert@464: 
robert@464: 		partialFrequencies[0] 	= NULL;								// for free check
robert@464: 	}
robert@464: }
robert@464: 
robert@464: 
robert@464: 
robert@464: void Partials::update(int parIndex, int frameNum)
robert@464: {
robert@464: 	partialFrequencies[parIndex] = new float[frameNum];
robert@464: 	partialAmplitudes[parIndex]	 = new float[frameNum];
robert@464: 	partialFreqDelta[parIndex]	 = new float[frameNum];
robert@464: 	partialAmpDelta[parIndex]	 = new float[frameNum];
robert@464: 
robert@464: 	fill(partialFreqDelta[parIndex], partialFreqDelta[parIndex]+frameNum, 99999.0);	// in the end, only the last one will have 99999
robert@464: 	fill(partialAmpDelta[parIndex], partialAmpDelta[parIndex]+frameNum, 99999.0);	// in the end, only the last one will have 99999
robert@464: }
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: //------------------------------------------------------------------------------------------------
robert@464: // spear parser
robert@464: //------------------------------------------------------------------------------------------------
robert@464: Spear_parser::Spear_parser()
robert@464: {
robert@464: 	// some default values
robert@464: 	hopSize			= -1;
robert@464: 	fileSampleRate	= -1;
robert@464: }
robert@464: 
robert@464: Spear_parser::~Spear_parser()
robert@464: {
robert@464: }
robert@464: 
robert@464: void Spear_parser::calculateHopSize(char *filename)
robert@464: {
robert@464: 	int index 		= 0;
robert@464: 	bool prevWas_ 	= false;
robert@464: 	bool found_h	= false;
robert@464: 	int n 			= 0;
robert@464: 
robert@464: 	hopSize 		= 0;
robert@464: 
robert@464: 	do
robert@464: 	{
robert@464: 		// check if '_'
robert@464: 		if(filename[index] == '_')
robert@464: 			prevWas_ = true;
robert@464: 		else if( (filename[index] == 'h') && prevWas_) // if it is not, but it is 'h' and previous was '_', found "_h"!
robert@464: 		{
robert@464: 			found_h = true;
robert@464: 			while(filename[index] != '\0')
robert@464: 			{
robert@464: 				index++;
robert@464: 				if( (filename[index] == '.') || (filename[index] == '_'))
robert@464: 					break;
robert@464: 				else // i am not checking if char are digits...!
robert@464: 				{
robert@464: 					n = filename[index];
robert@464: 					hopSize =  hopSize*10+(n-48);
robert@464: 				}
robert@464: 			}
robert@464: 		}
robert@464: 		else	// else, nothing
robert@464: 			prevWas_ = false;
robert@464: 		index++;
robert@464: 	}
robert@464: 	while( (filename[index] != '\0') && !found_h );
robert@464: 
robert@464: 	if( !found_h || (hopSize<1) )
robert@464: 		hopSize = 551;	// default val
robert@464: 
robert@464: }
robert@464: 
robert@464: 
robert@464: bool Spear_parser::parser(char *filename, int hopsize, int samplerate)
robert@464: {
robert@464: 	string name = string(filename);
robert@464: 	int len		= name.length();
robert@464: 	// invoke correct parser according to the type of file...just checking the extension, crude but functional
robert@464: 	if( (name[len-4]=='.') && (name[len-3]=='d') && (name[len-2]=='b') && (name[len-1]=='x') )
robert@464: 		return DBXparser(filename, samplerate);				// .dbox
robert@464: 	else
robert@464: 		return TXTparser(filename, hopSize, samplerate);	// .txt, or whatever
robert@464: }
robert@464: 
robert@464: 
robert@464: bool Spear_parser::DBXparser(char *filename, int samplerate)
robert@464: {
robert@464: 	fileSampleRate 	= samplerate;
robert@464: 
robert@464: 	// working vars
robert@464: 	int parNum		= 0;	// total num of partials
robert@464: 	int hopNum		= 0;	// total num of hops
robert@464: 
robert@464: 	//----------------------------------------------------------------------------------------
robert@464: 	// open a file
robert@464: 	ifstream fin;
robert@464: 	fin.open(filename, ios::in | ios::binary);
robert@464: 	if (!fin.good())
robert@464: 	{
robert@464: 		cout << "Parser Error: file not found" << endl;	// exit if file not found
robert@464: 		return false;
robert@464: 	}
robert@464: 
robert@464: 	gettimeofday(&start, NULL);
robert@464: 	//----------------------------------------------------------------------------------------
robert@464: 	// general data
robert@464: 
robert@464: 	// look for partial count
robert@464: 	fin.read((char *) &parNum, sizeof(int));
robert@464: 	partials.parNum 		= parNum;
robert@464: 
robert@464: 	// look for hop count
robert@464: 	fin.read((char *) &hopNum, sizeof(int));
robert@464: 	partials.setHopNum(hopNum);
robert@464: 
robert@464: 	// look for hop size
robert@464: 	fin.read((char *) &hopSize, sizeof(int));
robert@464: 	partials.hopSize 		= hopSize;		// it's handy for both classes to know it
robert@464: 
robert@464: 	// init partials data structure
robert@464: 	partials.init(parNum, hopSize, true);
robert@464: 
robert@464: 	// look for max active par num
robert@464: 	fin.read((char *) &(partials.maxActiveParNum), sizeof(int));
robert@464: 
robert@464: 
robert@464: 
robert@464: 	// partial data
robert@464: 
robert@464: 	// start frame of each partial
robert@464: 	fin.read((char *) partials.partialStartFrame, sizeof(int)*parNum);
robert@464: 
robert@464: 	// num of frames of each partial
robert@464: 	fin.read((char *) partials.partialNumFrames, sizeof(int)*parNum);
robert@464: 
robert@464: 	// frequency mean of each partial
robert@464: 	fin.read((char *) partials.partialFreqMean, sizeof(int)*parNum);
robert@464: 
robert@464: 	for(int par=0; par<parNum; par++)
robert@464: 	{
robert@464: 		int frameNum = partials.partialNumFrames[par];
robert@464: 		partials.update(par, frameNum);
robert@464: 		fin.read((char *)partials.partialAmplitudes[par], sizeof(float)*frameNum);		// amplitude of each partial in each frame
robert@464: 		fin.read((char *)partials.partialFrequencies[par], sizeof(float)*frameNum);		// frequency of each partial in each frame
robert@464: 		fin.read((char *)partials.partialAmpDelta[par], sizeof(float)*frameNum);			// amplitude delta of each partial in each frame
robert@464: 		fin.read((char *)partials.partialFreqDelta[par], sizeof(float)*frameNum);			// frequency delta of each partial in each frame
robert@464: 	}
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 	// frame data
robert@464: 
robert@464: 	// number of active partial per each frame
robert@464: 	fin.read((char *) partials.activePartialNum, sizeof(short)*(hopNum+1));
robert@464: 	// init array
robert@464: 	for(int frame=0; frame<hopNum+1; frame++)
robert@464: 	{
robert@464: 		partials.activePartials[frame] = new unsigned int[partials.activePartialNum[frame]];
robert@464: 		fin.read((char *)partials.activePartials[frame], sizeof(int)*partials.activePartialNum[frame]);			// active partials per each frame
robert@464: 	}
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: 	gettimeofday(&stop, NULL);
robert@464: 	parserT = ( (stop.tv_sec*1000000+stop.tv_usec) - (start.tv_sec*1000000+start.tv_usec) );
robert@464: 
robert@464: 
robert@464: 	printf("\n-----------------------\n");
robert@464: 	printf("\nFile: %s\n", filename);
robert@464: 	printf("\n-----------------------\n");
robert@464: 	printf("Profiler\n");
robert@464: 	printf("-----------------------\n");
robert@464: 	printf("File parser:\t\t\t%lu usec\n", parserT);
robert@464: 	printf("\n\nTotal:\t\t%lu usec\n", parserT);
robert@464: 	printf("-----------------------\n");
robert@464: 
robert@464: 	fin.close();
robert@464: 
robert@464: 	return true;
robert@464: }
robert@464: 
robert@464: 
robert@464: 
robert@464: 
robert@464: bool Spear_parser::TXTparser(char *filename, int hopsize, int samplerate)
robert@464: {
robert@464: 	hopSize 		= hopsize;
robert@464: 	fileSampleRate 	= samplerate;
robert@464: 	if(hopsize<0)
robert@464: 	{
robert@464: 		gettimeofday(&start, NULL);
robert@464: 		calculateHopSize(filename);
robert@464: 		gettimeofday(&stop, NULL);
robert@464: 		hopSizeT = ( (stop.tv_sec*1000000+stop.tv_usec) - (start.tv_sec*1000000+start.tv_usec) );
robert@464: 	}
robert@464: 	else
robert@464: 		hopSizeT = 0;
robert@464: 
robert@464: 	calculateDeltaTime();
robert@464: 
robert@464: 	// working vars
robert@464: 	char * token;			// where to save single figures from file
robert@464: 	string s		= "";	// where to save lines from file
robert@464: 	int parNum		= 0;	// total num of partials
robert@464: 	int parIndex	= -1;	// index of current partial
robert@464: 	int frameNum	= 0;	// total num of frames
robert@464: 	int frameIndex	= -1;	// index of current frame
robert@464: 	int startSample	= -1;	// sample value for first frame of partials
robert@464: 	int endSample	= -1;	// sample value for last frame of partials
robert@464: 	int maxSample	= 0;	// to calculate total number of hops in file
robert@464: 	int missSampCnt = 0;	// number of mising samples
robert@464: 	double freq		= 0;	// to calculate frequency delta
robert@464: 	double prevFreq	= 0;	// to calculate frequency delta
robert@464: 	double amp		= 0;	// to calculate amplitude delta
robert@464: 	double prevAmp	= 0;	// to calculate amplitude delta
robert@464: 
robert@464: 
robert@464: 	//----------------------------------------------------------------------------------------
robert@464: 	// open a file
robert@464: 	ifstream fin;
robert@464: 	fin.open(filename);
robert@464: 	if (!fin.good())
robert@464: 	{
robert@464: 		cout << "Parser Error: file not found" << endl;	// exit if file not found
robert@464: 		return false;
robert@464: 	}
robert@464: 
robert@464: 	gettimeofday(&start, NULL);
robert@464: 	//----------------------------------------------------------------------------------------
robert@464: 	// init partials data structure
robert@464: 	getline(fin, s);
robert@464: 	getline(fin, s);
robert@464: 	getline(fin, s);	// third line is the first we are interested into
robert@464: 
robert@464: 	// look for partial count
robert@464: 	token = strtok((char *)s.c_str(), " ");
robert@464: 	// check if first token is there
robert@464: 	if(token)
robert@464: 	{
robert@464: 		token = strtok(0, " ");
robert@464: 		// check if second token is there
robert@464: 		if(token)
robert@464: 			parNum = atoi(token);
robert@464: 		#ifdef DO_CHECKS
robert@464: 		else
robert@464: 		{
robert@464: 			cout << "Parser Error: partial count not found, bad file format" << endl;	// exit if value not found
robert@464: 			return false;
robert@464: 		}
robert@464: 		#endif
robert@464: 	}
robert@464: 	#ifdef DO_CHECKS
robert@464: 	else
robert@464: 	{
robert@464: 		cout << "Parser Error: partial count not found, bad file format" << endl;		// exit if value not found
robert@464: 		return false;
robert@464: 	}
robert@464: 	#endif
robert@464: 	// from now on we take for granted that format is correct
robert@464: 
robert@464: 	// init partials data structure
robert@464: 	partials.init(parNum, hopSize);
robert@464: 
robert@464: 	//----------------------------------------------------------------------------------------
robert@464: 	// fill in partials data structure
robert@464: 	getline(fin, s);		// get rid of intro line "partials-data"
robert@464: 	getline(fin, s);		// first important line
robert@464: 
robert@464: 	while (!fin.eof())
robert@464: 	{
robert@464: 		//-------------------------------------
robert@464: 		// partial specific info
robert@464: 		token		= strtok((char *)s.c_str(), " ");
robert@464: 		parIndex	= atoi(token);						// partial index
robert@464: 
robert@464: 		token		= strtok(0, " ");					// num of frames, not used, cos we will do linear interpolation for missing frames
robert@464: //		frameNum	= atoi(token);
robert@464: //		partials.partialNumFrames[parIndex]	= frameNum;
robert@464: 
robert@464: 		token		= strtok(0, " ");					// time of first frame, still char *
robert@464: 		startSample = fromTimeToSamples(atof(token));	// convert time to samples
robert@464: 		partials.partialStartSample[parIndex]	= startSample;
robert@464: 
robert@464: 		token		= strtok(0, " ");					// time of last frame, still char *
robert@464: 		endSample	= fromTimeToSamples(atof(token)); 	// convert time to samples
robert@464: 		partials.partialEndSample[parIndex]		= endSample;
robert@464: 
robert@464: 		frameNum	= ((endSample-startSample)/hopSize) + 1;	// num of frames, including missing consecutive ones [+1 one cos we count frames, not hops]
robert@464: 		partials.partialNumFrames[parIndex]		= frameNum;
robert@464: 
robert@464: 
robert@464: 		// check if this is the highest sample value so far
robert@464: 		if(endSample > maxSample)
robert@464: 			maxSample = endSample;
robert@464: 
robert@464: 		// update data structure
robert@464: 		partials.update(parIndex, frameNum);
robert@464: 
robert@464: 
robert@464: 		//-------------------------------------
robert@464: 		// frames
robert@464: 		getline(fin, s);
robert@464: 		token		= strtok((char *)s.c_str(), " ");	// frame time
robert@464: 		frameIndex	= -1;
robert@464: 
robert@464: 		// unroll first iteration, so that in the following loop we save the check on the last frame to calculate increments
robert@464: 		if(token)						// all frames data are on one line, in groups of 3 entries
robert@464: 		{
robert@464: 			frameIndex++;
robert@464: 
robert@464: 			endSample	= fromTimeToSamples(atof(token));
robert@464: 
robert@464: 			token		= strtok(0, " ");	// frame frequency
robert@464: 			prevFreq	= atof(token);
robert@464: 			partials.partialFrequencies[parIndex][frameIndex]	= (float)prevFreq;
robert@464: 			partials.partialFreqMean[parIndex] 					+= prevFreq;		// for frequency mean
robert@464: 
robert@464: 			token	 	= strtok(0, " ");	// frame amplitude
robert@464: 			prevAmp  	= atof(token);
robert@464: 			partials.partialAmplitudes[parIndex][frameIndex]	= (float)prevAmp;
robert@464: 
robert@464: 			token 		= strtok(0, " ");	// next frame frequency, to be checked
robert@464: 		}
robert@464: 
robert@464: 		// here the loop starts
robert@464: 		while(token)						// all frames data are on one line, in groups of 3 entries
robert@464: 		{
robert@464: 			frameIndex++;
robert@464: 			missSampCnt 	= 0;
robert@464: 
robert@464: 			startSample		= fromTimeToSamples(atof(token));
robert@464: 
robert@464: 			token			= strtok(0, " ");	// frame frequency
robert@464: 			freq			= atof(token);
robert@464: 
robert@464: 			token			= strtok(0, " ");	// frame amplitude
robert@464: 			amp				= atof(token);
robert@464: 			// now we know all about the current frame, but we want to know if some frames are missing between this and the last one
robert@464: 
robert@464: 			// while current frame sample is farther than one hopsize...
robert@464: 			while(startSample > endSample+hopSize)
robert@464: 			{
robert@464: 				missSampCnt++;				// ...one sample is missing
robert@464: 				endSample += hopSize;		// move to next hop
robert@464: 			}
robert@464: 
robert@464: 			// if frames are missing do interpolation and update indices
robert@464: 			if(missSampCnt>0)
robert@464: 				startSample = interpolateSamples(parIndex, &frameIndex, missSampCnt, endSample+hopSize, freq, amp, &prevFreq, &prevAmp);
robert@464: 
robert@464: 			partials.partialFrequencies[parIndex][frameIndex]	= (float)freq;
robert@464: 			partials.partialFreqMean[parIndex] 					+= freq;			// for frequency mean
robert@464: 			partials.setFreqDelta(parIndex, frameIndex-1, (freq-prevFreq)/hopSize);	// freq delta between prev and current frame
robert@464: 			prevFreq 	= freq;
robert@464: 
robert@464: 			partials.partialAmplitudes[parIndex][frameIndex]	= (float)amp;
robert@464: 			partials.setAmpDelta(parIndex, frameIndex-1, (amp-prevAmp)/hopSize);	// amp delta between prev and current frame
robert@464: 			prevAmp		= amp;
robert@464: 
robert@464: 			endSample	= startSample;
robert@464: 			token		= strtok(0, " ");	// next frame frequency, to be checked
robert@464: 		}
robert@464: 		#ifdef DO_CHECKS
robert@464: 		if(frameIndex != (frameNum-1))
robert@464: 		{
robert@464: 			cout << "Parser Error: frame count mismatch on partial " << parIndex << ", bad file format"  << endl;	// exit if mismatch
robert@464: 			cout << "frameIndex: " << frameIndex << endl;
robert@464: 			cout << "frameNum: " << frameNum << endl;
robert@464: 			return false;
robert@464: 		}
robert@464: 		#endif
robert@464: 
robert@464: 		partials.partialFreqMean[parIndex] /= partials.partialNumFrames[parIndex];									// frequency mean
robert@464: 
robert@464: 		getline(fin, s);					// next partial line, to check
robert@464: 	}
robert@464: 	#ifdef DO_CHECKS
robert@464: 	if(parIndex != (parNum-1))
robert@464: 	{
robert@464: 		cout << "Parser Error: partial count mismatch, bad file format"  << endl;									// exit if mismatch
robert@464: 		return false;
robert@464: 	}
robert@464: 	#endif
robert@464: 
robert@464: 	partials.setHopNum(maxSample/hopSize);
robert@464: 
robert@464: 	gettimeofday(&stop, NULL);
robert@464: 	parserT = ( (stop.tv_sec*1000000+stop.tv_usec) - (start.tv_sec*1000000+start.tv_usec) );
robert@464: 
robert@464: 	gettimeofday(&start, NULL);
robert@464: 	staticCalculations();
robert@464: 	gettimeofday(&stop, NULL);
robert@464: 	staticT = ( (stop.tv_sec*1000000+stop.tv_usec) - (start.tv_sec*1000000+start.tv_usec) );
robert@464: 
robert@464: 	fin.close();
robert@464: 
robert@464: 
robert@464: 	printf("\n-----------------------\n");
robert@464: 	printf("\nFile: %s\n", filename);
robert@464: 	printf("\n-----------------------\n");
robert@464: 	printf("Profiler\n");
robert@464: 	printf("-----------------------\n");
robert@464: 	printf("Hop size parser:\t\t%lu usec\n", hopSizeT);
robert@464: 	printf("File parser:\t\t\t%lu usec\n", parserT);
robert@464: 	printf("Static calculations:\t\t%lu usec\n", staticT);
robert@464: 	printf("\n\nTotal:\t\t%lu usec\n", hopSizeT+parserT+staticT);
robert@464: 	printf("-----------------------\n");
robert@464: 
robert@464: 	return true;
robert@464: }
robert@464: 
robert@464: 
robert@464: int Spear_parser::interpolateSamples(int parIndex, int *frameIndex, int missCnt, int nextSample, double nextFreq, double nextAmp, double *prevFreq, double *prevAmp)
robert@464: {
robert@464: 	int frame			= *frameIndex;						// current frame index
robert@464: 	int sample			= nextSample - (hopSize*(missCnt)); // move from next real frame sample to first missing frame sample
robert@464: 	double freq			= *prevFreq;						// freq of the prev real frame
robert@464: 	double freqStep		= (nextFreq-*prevFreq)/(missCnt+1);	// fixed freq step between hops, for missing frames [linear interpolation]
robert@464: 	double deltaFreq	= freqStep/hopSize;					// fixed hop freq step in samples
robert@464: 	double amp			= *prevAmp;							// same for amp...
robert@464: 	double ampStep		= (nextAmp-*prevAmp)/(missCnt+1);
robert@464: 	double deltaAmp		= ampStep/hopSize;
robert@464: 
robert@464: 	// for each missing frame
robert@464: 	for(int i=0; i<missCnt; i++)
robert@464: 	{
robert@464: 		// calculate values for current missing frame
robert@464: 		freq	+= freqStep;
robert@464: 		amp		+= ampStep;
robert@464: 		// save values
robert@464: 		partials.partialFrequencies[parIndex][frame]	= freq;
robert@464: 		partials.partialAmplitudes[parIndex][frame]		= amp;
robert@464: 		partials.partialFreqMean[parIndex]				+= freq;	// for frequency mean
robert@464: 		// set deltas of previous frame [real or missing]
robert@464: 		partials.setFreqDelta(parIndex, frame-1, deltaFreq);
robert@464: 		partials.setAmpDelta(parIndex, frame-1, deltaAmp);
robert@464: 		// move to next frame [missing or real]
robert@464: 		sample += hopSize;
robert@464: 		frame++;
robert@464: 	}
robert@464: 
robert@464: 	// update global values
robert@464: 	*frameIndex	= frame;
robert@464: 	*prevFreq	= freq;
robert@464: 	*prevAmp	= amp;
robert@464: 
robert@464: 	return sample;	// return the frame sample of the next real frame
robert@464: }
robert@464: 
robert@464: 
robert@464: 
robert@464: // for each frame, statically calculate:
robert@464: // - which partial is active [and the total num of active partials]
robert@464: // - at which local frame each partial is
robert@464: void Spear_parser::staticCalculations()
robert@464: {
robert@464: 	partials.maxActiveParNum = 0;				// init to find maximum
robert@464: 
robert@464: 	unsigned short *indices	= new unsigned short[partials.parNum];	// temp array to store up to the maximum num of active partial indices
robert@464: 	unsigned int activeCnt	= 0;						// counts the num of active partials in each frame
robert@464: 
robert@464: 	unsigned int frameSample = 0;						// current frame in samples
robert@464: 
robert@464: 	char *partialStarted = new char [partials.parNum];	// index of the last local frame found per each partial
robert@464: 	fill(partialStarted, partialStarted+partials.parNum, 0);
robert@464: 
robert@464: 	for(unsigned int i=0; i<partials.hopNum+1; i++)		// for each frame [not hops, this explains the +1]
robert@464: 	{
robert@464: 		//partials.localPartialFrames[i] = new int[partials.parNum];	// init all local frames to -1
robert@464: 		//fill(partials.localPartialFrames[i], partials.localPartialFrames[i]+partials.parNum, -1);
robert@464: 
robert@464: 		frameSample = i*hopSize;	// current frame, expressed in samples
robert@464: 		activeCnt	  = 0;			// reset a each frame
robert@464: 
robert@464: 		for(unsigned int j=0; j<partials.parNum; j++)	// for each partial
robert@464: 		{
robert@464: 			// check if inside active time region [expressed in samples]
robert@464: 			if( (frameSample>=partials.partialStartSample[j]) && (frameSample<partials.partialEndSample[j]) )	// frame sample not equal to end sample, this filters out last frames and partials with one frame only
robert@464: 			{
robert@464: 				// activity
robert@464: 				indices[activeCnt] = j;	// save active index
robert@464: 				activeCnt++;			// increase counter
robert@464: 
robert@464: 				// partial local frames
robert@464: 				if(partialStarted[j]==0)	// this partial has just started, so current local frame is first frame
robert@464: 				{
robert@464: 					partialStarted[j] 		 		= 1;
robert@464: 					partials.partialStartFrame[j]	= i;	// here is the number of the first frame
robert@464: 				}
robert@464: 			}
robert@464: 		}
robert@464: 
robert@464: 		// activity
robert@464: 		partials.activePartialNum[i] = activeCnt;							// save number of active partials for this frame
robert@464: 		partials.activePartials[i]	 = new unsigned int[activeCnt];					// set correct size to save all indices
robert@464: 
robert@464: 		// look for maximum number of active partials at the same time
robert@464: 		if(activeCnt > partials.maxActiveParNum)
robert@464: 			partials.maxActiveParNum = activeCnt;
robert@464: 
robert@464: 		// copy indices
robert@464: 		for(unsigned int k=0; k<activeCnt; k++)
robert@464: 			partials.activePartials[i][k] = indices[k];
robert@464: 	}
robert@464: 
robert@464: 	delete[] indices;
robert@464: 	delete[] partialStarted;
robert@464: 
robert@464: 	delete[] partials.partialStartSample;
robert@464: 	delete[] partials.partialEndSample;
robert@464: }
robert@464: 
robert@464: 
robert@464: