--- ../../dev2/gift/FeatureExtraction/gabor.c	2006-09-20 16:30:58.000000000 -0500
+++ FeatureExtraction/gabor.c	2006-09-26 11:20:57.000000000 -0500
@@ -125,11 +125,22 @@
 	for (x = 0; x < width; x++) {
 	for (y = 0; y < height; y++) {
 		target_conv=&conv[((width*height)-1)-(x*height+y+(kernal_size[filter_scale]/2))];
+		if (((y>=kernal_size[filter_scale]/2)) && ((y+kernal_size[filter_scale]/2)<height))
+		  {
+		    /* first do the matrix multiply, then do the summation, so the matrix multiply can be vectored. */
+		    for (k = 0; k < kernal_size[filter_scale]; k++)
+		      temparray[k] = target_kernal[k]*target_conv[k];
+		    for (k = 0; k < kernal_size[filter_scale]; k++)
+		      output[y*width+x] += temparray[k];
+		  }
+		else
+		  {
 		for (k = 0; k < kernal_size[filter_scale]; k++) {
 			if (((y+(kernal_size[filter_scale]/2))>=k) && (y+(kernal_size[filter_scale]/2)<height+k))
 				output[y*width + x] +=
 					target_kernal[k]*target_conv[k];
 		}
+		  }
 	}
 	}
 
@@ -164,11 +175,22 @@
 	for (x = 0; x < width; x++) {
 	for (y = 0; y < height; y++) {
 		target_conv=&conv[((width*height)-1)-(x*height+y+(kernal_size[filter_scale]/2))];
+		if (((y>=kernal_size[filter_scale]/2)) && ((y+kernal_size[filter_scale]/2)<height))
+		  {
+		    /* first do the matrix multiply, then do the summation, so the matrix multiply can be vectored. */
+		    for (k = 0; k < kernal_size[filter_scale]; k++)
+		      temparray[k] = target_kernal[k]*target_conv[k];
+		    for (k = 0; k < kernal_size[filter_scale]; k++)
+		      output[y*width+x] -= temparray[k];
+		  }
+		else
+		  {
 		for (k = 0; k < kernal_size[filter_scale]; k++) {
 			if (((y+(kernal_size[filter_scale]/2))>=k) && (y+(kernal_size[filter_scale]/2)<height+k))
 				output[y*width + x] -=
 					target_kernal[k]*target_conv[k];
 		}
+		}
 	}
 	}