comparison src/DynamicVector.cpp @ 9:75dcd1308658

looked at tempo process, likelihood function changed and improved. Now drawing using constrained vector function. Good working version.

 	arraySize = 0;
 	maximumValue = 0;
 	MAPestimate = 0;
 	offset = 0;
 	scalar = 1;
+	integratedEstimate = length/2;
 	
 	gaussianLookupMean = (double) GAUSSIAN_LOOKUP_LENGTH/2;
 	gaussianLookupStdDev = (double)(GAUSSIAN_LOOKUP_LENGTH/16);
 	double factor = 1.0;//(1.0 / (gaussianLookupStdDev*sqrt(2*PI)) );//1.0;//-1.0/(2*PI*sqrt(gaussianLookupStdDev));
 	for (int i = 0;i < GAUSSIAN_LOOKUP_LENGTH;i++){

 	for (int i = 0; i < len;i++){
 		array.push_back(0);
 	}
 	length = len;
 	arraySize = array.size();
+	integratedEstimate = length/2;
 }
 
 
 double DynamicVector::getMaximum(){
 	int i;

 	return max;
 }
 
 double DynamicVector::getIntegratedEstimate(){
 	//returns the index of the integrated average - where the probability distribution is centred
-	double estimate = 0;
+	integratedEstimate = 0;
 	double integratedTotal = 0;
 	for (int i = 0;i < length;i++){
-		estimate += array[i]*i;
+		integratedEstimate += array[i]*i;
 		integratedTotal += array[i];
 	}
 	if (integratedTotal > 0){
-		estimate /= integratedTotal;
-	}
-	return estimate;
-}
+		integratedEstimate /= integratedTotal;
+	}
+	return integratedEstimate;
+}
+
+void DynamicVector::updateIntegratedEstimate(){
+	//returns the index of the integrated average - where the probability distribution is centred
+	integratedEstimate = 0;
+	double integratedTotal = 0;
+	for (int i = 0;i < length;i++){
+		integratedEstimate += array[i]*i;
+		integratedTotal += array[i];
+	}
+	if (integratedTotal > 0){
+		integratedEstimate /= integratedTotal;
+	}
+
+}
+
 
 void DynamicVector::zero(){
 	for (int i = 0;i < array.size();i++)
 		array[i] = 0;
 }

 }
 
 void DynamicVector::addGaussianShape(double mean, double StdDev, double factor){
 	
 	int i;
+	double std_dev_factor = (2*StdDev*StdDev);
 	factor *= (1/(StdDev*sqrt(2*PI)));
-	for (i=0;i<array.size();i++){
-		array[i] += factor*exp(-1*(i-mean)*(i-mean)/(2*StdDev*StdDev));
-	}
-	
-	//addGaussianShapeByLookupTable(mean, StdDev, factor);
+	int maxVal = min((int) array.size(), (int)(mean + 4.8*StdDev));
+	int minVal = max(0, (int)(mean - 4.8*StdDev));
+	
+	for (i=minVal;i < maxVal;i++){
+		array[i] += factor*exp(-1*(i-mean)*(i-mean)/(std_dev_factor));
+	}
+	
+//	addGaussianShapeByLookupTable(mean, StdDev, factor);
 }
 
 void DynamicVector::addGaussianShapeByLookupTable(double& mean, double& StdDev, double& factor){
 	int i;
 	int lookupIndex ;

 
 void DynamicVector::drawVector(const int& minIndex, const int& maxIndex){
 
 		
 		double stepSize = ofGetWidth() / (double)(maxIndex - minIndex);
-		double screenHeight = ofGetHeight();
+		double screenHeight = (double) ofGetHeight();
 		double maxVal = getMaximum();
+	
 		int startInt = max(1,minIndex+1);
 		int endInt = min(maxIndex, (int)array.size());
 		double heightConstant = screenHeight / maxVal;
 		int lastHeightPixel = heightConstant * (1 - array[startInt-1]);
 		int newHeightPixel;