comparison jnmr/midiEventHolder.cpp @ 46:43edc8abe2a7

Fixed bug in complex update due to mismatch between ms and vector units

  */
 
 
 //Main file to look at here is newNoteEvent() - this calls everything else to update the Bayesian array
 
+//Relative speed is RECORDED / PLAYED
+//so for timedifference t
+//this is t/scalar in terms of position units of played time
+//and by relative speed r, r*t/scalar
+
+
 #include "midiEventHolder.h"
 
 #include <iostream>
 #include <fstream>
 #include <assert.h>
 
 midiEventHolder::midiEventHolder(){
 	
 
+	double scalarForPositionVectors = 20;
 	
 //	recordedNoteOnIndex = 0;
 	alignmentPosition = 0;
 	
 	useTempoPrior = false;//puts sine wave round tempo

 
 
 	
 	speedPriorValue = 1.0;
 	
-	
-	bayesStruct.resetSize(matchWindowWidth);
-	bayesStruct.setPositionDistributionScalar(1);
+	int tmpArraySize = (int)(matchWindowWidth/scalarForPositionVectors);
+	bayesStruct.resetSize(tmpArraySize);
+	bayesStruct.setPositionDistributionScalar(scalarForPositionVectors);
+	printf("ARRAY SIZE IS %i\n", tmpArraySize);
 	
 	bayesStruct.resetSpeedSize(200);
 	bayesStruct.setRelativeSpeedScalar(0.01);
 	bayesStruct.relativeSpeedPrior.getMaximum();
 	//bayesStruct.simpleExample();

 //	relativeSpeedForSmooth = 1.0;
 //	storedSmoothPlayPosition = smoothPlayPosition;
 //	lastSmoothUpdateTime = getTimeNow(0);
 	
 	printf("reset speed prior is %f\n", speedPriorValue);
-	bayesStruct.resetSpeedToOne();
+//	bayesStruct.resetSpeedToOne();
 	bayesStruct.setSpeedPrior(speedPriorValue);
 	setMatchedNotesBackToFalse();
 
-	periodCounter = 0;
-	for (int i = 0;i < periodValues.size();i++){
+//	periodCounter = 0;
+//	for (int i = 0;i < periodValues.size();i++){
 	//	printf("period at %f is %f\n", periodValues[i][2], periodValues[i][1]);
-	}
+//	}
 /*	if (periodValues.size() > 0){
 		updatePeriodValue(0);// periodValues[0][2];
 		printf("Resetting period to %f , size is %i\n", period,  (int)periodValues.size());
 	}
 */

 	return (period * ticks / (double) pulsesPerQuarternote);
 }
 
 void midiEventHolder::newNoteOnEvent(int pitch, int velocity, double timePlayed){
 //	tempoSpeedString = "";
-	
+
 	//MOVE INTO BAYESSTRUCT?? XXX
 	//bayesStruct.copyPriorToPosterior();
 	//why was this here??
 	bayesStruct.prior.copyFromDynamicVector(bayesStruct.posterior);//try the otehr way
 	//bayesStruct.copyPriorToPosterior();

 //	timeString += " :  Predicted MAP time" + ofToString(newMAPestimateTime,0);
 
 	//then we recalculate the window start based on MAP being central
 	//then we do the matches on these and the likelihood on these.
 	
-	bayesStruct.setNewDistributionOffsets(max(0., bayesStruct.bestEstimate - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2)));
-//	bayesStruct.prior.offset = max(0.,newMAPestimateTime - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2));
+	bayesStruct.setNewDistributionOffsets(max(bayesStruct.posterior.offset, bayesStruct.bestEstimate - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2)));
+	//bayesStruct.setNewDistributionOffsets(max(0.,newMAPestimateTime - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2));
 	
 	timeString += " \n :  new offset " + ofToString(bayesStruct.prior.offset , 0);
 	timeString += " \n best estimate "+ofToString(bayesStruct.bestEstimate, 1);
 	timeString += " error "+ofToString(minimumMatchError, 0);
 	timeString += " map "+ofToString(bayesStruct.relativeSpeedPosterior.integratedEstimate, 1);

 	
 	timeString += " new OFF "+ofToString(bayesStruct.bestEstimate - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2), 1);
 	timeString += " notearrayindex "+ofToString(noteArrayIndex, 0);
 	//when this is off teh screen there is a problem somehow XXX
 	
-	bayesStruct.posterior.offset = max(0., bayesStruct.bestEstimate - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2));//	bayesStruct.prior.offset = max(0., bayesStruct.bestEstimate - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2));
+	bayesStruct.posterior.offset = bayesStruct.prior.offset;
+	//max(0., bayesStruct.bestEstimate - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2));//	bayesStruct.prior.offset = max(0., bayesStruct.bestEstimate - (bayesStruct.prior.scalar*bayesStruct.prior.arraySize/2));
 
 	//trying to switch to prior
 	
 	
 	bayesStruct.lastEventTime = timePlayed;//bayesStruct.lastEventTime = ofGetElapsedTimeMillis();

 	if (recordedEventTimes.size() > 0){
 	updateTempo();
 	//calcuateNewInterNoteIntervals();
 	}
 	
+
 	//storedSmoothPlayPosition = smoothPlayPosition;
 	 
 }
 
 void midiEventHolder::updateTempo(){

 	double quantity = likelihoodToNoiseRatio / numberOfMatches;
 	
 	for (int i = 0;i < numberOfMatches && matchesFound[i] >= 0 && matchesFound[i] < recordedEventTimes.size();i++){
 	//	printf("match times %i of %i::%f adding likelihood to %f\n", i, numberOfMatches, recordedEventTimes[matchesFound[i]], recordedEventTimes[matchesFound[i]] - bayesStruct.likelihood.offset);
 		//this is the vent time since start of file
-		if (recordedEventTimes[matchesFound[i]] - bayesStruct.likelihood.offset < bayesStruct.likelihood.arraySize){
+		if (recordedEventTimes[matchesFound[i]] - bayesStruct.likelihood.offset < bayesStruct.likelihood.arraySize*bayesStruct.likelihood.scalar){
 	//		double confidenceMeasure = 0;
 	//		if (totalConfidence > 0)
 	//			confidenceMeasure =	bayesStruct.posterior.getValueAtMillis(recordedEventTimes[matchesFound[i]])/totalConfidence;
 
 			//bayesStruct.likelihood.addGaussianShape(recordedEventTimes[matchesFound[i]] - bayesStruct.likelihood.offset, likelihoodWidth,  quantity);//* confidenceMeasure

 		
 		
 		if (recordedTimeDifference > minimumTimeIntervalForTempoUpdate
 			&& speedRatio < maximumMatchSpeed && speedRatio > minimumMatchSpeed){
 			
-			/*		 printf("(%i)", previousIndex);
+		/*	 printf("(%i)", previousIndex);
 			 printf("[%i] :: ", recordedPreviousIndex);
 			 //	printf(" conf %f & %f ", currentMatchConfidence, previousMatchConfidence);
 			 printf(" rec{%f} vs play(%f) ", recordedTimeDifference, playedTimeDifference);
 			 printf("update on speed ratio %f\n", speedRatio);
-			 */	 
+		*/	 	 
 			//	matchString += " speed: "+ofToString(speedRatio, 3);
 			//	commented for debug
 			
 			
 			double priorWeighting = 1;

 	}
 	
 	
 	if (needToUpdate)
 		bayesStruct.updateTempoDistribution();
+	
 	//printf("current speed is %f\n", bayesStruct.relativeSpeedPosterior.getIndexInRealTerms(bayesStruct.relativeSpeedPosterior.MAPestimate));
 }
 
 
 void midiEventHolder::calcuateNewInterNoteIntervals(){

 		
 	smoothPlayPosition = newPosition;
 		
 }
 
-
+/*
 void midiEventHolder::updatePeriodValue(const double& millis){
 
 	double tmp = period;
-/*	
-	while (periodCounter >= 0 && periodCounter < periodValues.size()-1 && periodValues[periodCounter][2] < millis){
-		periodCounter++;
-	}
-	while (periodCounter > 0 && periodValues[periodCounter][2] > millis){
-		periodCounter--;
-	}
- */
+	
+//	while (periodCounter >= 0 && periodCounter < periodValues.size()-1 && periodValues[periodCounter][2] < millis){
+//		periodCounter++;
+//	}
+//	while (periodCounter > 0 && periodValues[periodCounter][2] > millis){
+//		periodCounter--;
+//	}
+ //
 	//period = periodValues[periodCounter][1];
 	
 	if (period != tmp){
 		printf("new period at %f of %f\n", millis, period);
 	}
 }
+*/
 
 void  midiEventHolder::updateNoteCounter(){
 	while (totalNoteCounterIndex < bestMatchIndex){
 		int tmpPitch = recordedNoteOnMatrix[totalNoteCounterIndex][1];
 		recordedTotalNoteCounterByPitch[tmpPitch] += 1;

 	}
 }
 
 
 void midiEventHolder::drawMidiFile(){
+	//ofBackground(80,80,80);
 	
 	//draws midi file on scrolling screen
 	int size = recordedNoteOnMatrix.size();
 	if (size > 0){
 		

 		ofDrawBitmapString(indexString, 20, 40);
 	 */
 	}	
 	
 	//ofDrawBitmapString(ofToString(timeOffsetForScreen, 1), 20,20);
-	
-	//ofDrawBitmapString(timeString, 20, 60);
+	ofSetColor(255,255,255);
+	ofDrawBitmapString(timeString, 20, 60);
 	
 //last played piutch
 	ofSetColor(0,200,0,50);
 	int yLocation = (*screenHeight) - ((lastPlayedPitch - noteMinimum )*(*screenHeight)/ (float)(noteMaximum - noteMinimum));						 
 	ofRect(0,yLocation, 100,  noteHeight);

 				   printf("\n");
 	}			   
 	
 }
 
-int midiEventHolder::getLocationFromTicks(double tickPosition){
-	return 0;
+//int midiEventHolder::getLocationFromTicks(double tickPosition){
+//	return 0;
 //	return (int)((float)(tickPosition - numberOfScreensIn*ticksPerScreen)*(*screenWidth)/(float)ticksPerScreen);
 //not used
-}
+//}
 
 int midiEventHolder::getLocationFromMillis(double millisPosition){
 	//(getEventTimeTicks(windowStartTime+matchWindowWidth) - numberOfScreensIn*ticksPerScreen)*(*screenWidth) / (double)ticksPerScreen
 	return (millisPosition - timeOffsetForScreen)*(*screenWidth)/getEventTimeMillis(ticksPerScreen);
 }
 
-
+/*
 void midiEventHolder::exampleCrossUpdate(){
 	
 	bayesStruct.crossUpdateArrays(bayesStruct.posterior, bayesStruct.relativeSpeedPosterior, 200);
 	
 }
-
+*/
 
 void midiEventHolder::setStartPlayingTimes(){
 	startPlayingTime = getTimeNow(0);//ofGetElapsedTimeMillis();
 	//startTime = startPlayingTime;
 	printf("starting playing at time %f\n", startPlayingTime);