Mercurial > geeqie.yaz

changeset 462:6a2934cd0883

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
histogram_read(): speed up calculations by 20%.
author zas_
date Mon, 21 Apr 2008 13:39:17 +0000
parents 65475285957d
children 544934a1ff89
files src/histogram.c
diffstat 1 files changed, 28 insertions(+), 13 deletions(-) [+]
line wrap: on
line diff
--- a/src/histogram.c	Mon Apr 21 11:44:03 2008 +0000
+++ b/src/histogram.c	Mon Apr 21 13:39:17 2008 +0000
@@ -104,7 +104,7 @@
 
 gulong histogram_read(Histogram *histogram, GdkPixbuf *imgpixbuf)
 {
-	gint w, h, i, j, srs, has_alpha;
+	gint w, h, i, j, srs, has_alpha, step;
 	guchar *s_pix;
 
 	if (!histogram) return 0;
@@ -117,25 +117,40 @@
 
 	memset(histogram->histmap, 0, sizeof(histogram->histmap));
 
-	for (i = 0; i < h; i++)
+	/* code duplication is here to speed up the calculation */
+	step = 3 + !!(has_alpha);
+	if (histogram->histogram_chan == HCHAN_MAX)
 		{
-		guchar *sp = s_pix + (i * srs); /* 8bit */
-		for (j = 0; j < w; j++)
+		for (i = 0; i < h; i++)
 			{
-			histogram->histmap[sp[0] + 0 * HISTOGRAM_SIZE]++;
-			histogram->histmap[sp[1] + 1 * HISTOGRAM_SIZE]++;
-			histogram->histmap[sp[2] + 2 * HISTOGRAM_SIZE]++;
-			if (histogram->histogram_chan == HCHAN_MAX)
+			guchar *sp = s_pix + (i * srs); /* 8bit */
+			for (j = 0; j < w; j++)
 				{
 				guchar t = sp[0];
 				if (sp[1]>t) t = sp[1];
 				if (sp[2]>t) t = sp[2];
-  				histogram->histmap[t + 3 * HISTOGRAM_SIZE]++;
+
+				histogram->histmap[sp[0] + 0 * HISTOGRAM_SIZE]++;
+				histogram->histmap[sp[1] + 1 * HISTOGRAM_SIZE]++;
+				histogram->histmap[sp[2] + 2 * HISTOGRAM_SIZE]++;
+				histogram->histmap[t + 3 * HISTOGRAM_SIZE]++;
+				sp += step;
 				}
-			else
-  				histogram->histmap[3 * HISTOGRAM_SIZE + (sp[0]+sp[1]+sp[2])/3]++;
-			sp += 3;
-			if (has_alpha) sp++;
+			}
+		}
+	else
+		{
+		for (i = 0; i < h; i++)
+			{
+			guchar *sp = s_pix + (i * srs); /* 8bit */
+			for (j = 0; j < w; j++)
+				{
+				histogram->histmap[sp[0] + 0 * HISTOGRAM_SIZE]++;
+				histogram->histmap[sp[1] + 1 * HISTOGRAM_SIZE]++;
+				histogram->histmap[sp[2] + 2 * HISTOGRAM_SIZE]++;
+				histogram->histmap[3 * HISTOGRAM_SIZE + (sp[0]+sp[1]+sp[2])/3]++;
+				sp += step;
+				}
 			}
 		}