comparison bayesianArraySrc/BayesianArrayStructure.cpp @ 8:572564b7cb85

added calculation posterior into both onset and pitch processes

 	setNewDistributionOffsets(0);
 	
 	int zeroIndex = posterior.getRealTermsAsIndex(0);
 	printf("ZERO INDEX %i\n", zeroIndex);
 	
-	posterior.addGaussianShapeFromRealTime(0, 60, 1);//one way to add at x msec
+	posterior.addGaussianShapeFromRealTime(0, 300, 1);//one way to add at x msec
 //	posterior.addGaussianShape(posterior.getRealTermsAsIndex(10), 50, 1);//alternative way
 	
-	//posterior.addToIndex(0, 1);
 	likelihood.addConstant(1);
 	
 	updateCounter = 0;
 	
 

 	printPostOffset();
 	
 	updateBestEstimate(timeDifference);
 	lastBestEstimateUpdateTime = newEventTime;//getTimeNow(timePlayed);
 
-	
+	//set TARGETS - commented tenmporarily
 	setNewDistributionOffsets(max(0., bestEstimate - (prior.scalar*prior.arraySize/2)));
+	
 	crossUpdateArrays(posterior, relativeSpeedPosterior, timeDifference);
 	
 	//i.e. using the same offset as prior
 	posterior.offset = prior.offset;// 
 	

 	
 }
 
 
 void BayesianArrayStructure::crossUpdateArrays(DynamicVector& position, DynamicVector& speed, double timeDifference){
+
 	//set the cutoff for offset of position first! XXX
-	
-//	printf("time difference %f, ", timeDifference);
 	
 	double timeDifferenceInPositionVectorUnits = timeDifference / prior.scalar;
 	
 	printf("CROSS UPDATE time diff %f ms is %f units; ", timeDifference, timeDifferenceInPositionVectorUnits);
 	prior.zero();//kill prior

 	else
 		translateByMaximumSpeed(timeDifferenceInPositionVectorUnits);	
 			
 
 	updateCounter++;
-	prior.renormalise();
+	prior.renormalise();//not strictly necessary??
 
 }
 
 void BayesianArrayStructure::complexCrossUpdate(const double& timeDifferenceInPositionVectorUnits){
 

 
 
 
 void BayesianArrayStructure::translateByMaximumSpeed(const double& timeDifferenceInPositionVectorUnits){
 
+	
 	int distanceMoved, newPriorIndex;
-	
-	double speedValue = relativeSpeedPosterior.getIndexInRealTerms(relativeSpeedPosterior.integratedEstimate);
+	double speedIndex = getSpeedEstimateIndex();
+	double speedValue = relativeSpeedPosterior.getIndexInRealTerms(speedIndex);
+	
+//	double speedValue = relativeSpeedPosterior.getIndexInRealTerms(relativeSpeedPosterior.integratedEstimate);
+	
 	//so for scalar 0.01, 50 -> speed value of 0.5
-	double speedContribution = relativeSpeedPosterior.array[relativeSpeedPosterior.integratedEstimate];
+	double speedContribution = relativeSpeedPosterior.array[(int)round(speedIndex)];
 		//so we have moved 
 		distanceMoved = round(timeDifferenceInPositionVectorUnits * speedValue);//round the value
-					//	printf("speed [%i] gives %f moved %i in %f units \n", i, speedValue, distanceMoved, timeDifferenceInPositionVectorUnits);
+				
+	//	printf("speed [%i] gives %f moved %i in %f units \n", i, speedValue, distanceMoved, timeDifferenceInPositionVectorUnits);
 			
 			for (int postIndex = 0;postIndex < posterior.arraySize;postIndex++){
 				//old posterior contributing to new prior
 				newPriorIndex = postIndex + posterior.offset - prior.offset + distanceMoved;
 				if (newPriorIndex >= 0 && newPriorIndex < prior.arraySize){