From a764d0a9512a7bd5c0e17af70b3c55e5fe7f9f71 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@bitwiz.org.uk>
Date: Fri, 17 Dec 2010 18:03:14 -0800
Subject: Scaling matrix is symmetric

---
 src/hrs-scaling.c | 2 ++
 1 file changed, 2 insertions(+)

(limited to 'src')

diff --git a/src/hrs-scaling.c b/src/hrs-scaling.c
index 80894131..3f03b8d8 100644
--- a/src/hrs-scaling.c
+++ b/src/hrs-scaling.c
@@ -172,6 +172,7 @@ static double iterate_scale(struct image *images, int n,
 			int mcomp;
 			struct image *imagem = &images[m];
 
+			if ( m > l ) continue;  /* Matrix is symmetric */
 			if ( m == crossed ) continue;
 
 			for ( h=0; h<n_ref; h++ ) {
@@ -206,6 +207,7 @@ static double iterate_scale(struct image *images, int n,
 			mcomp = m;
 			if ( m > crossed ) mcomp--;
 			gsl_matrix_set(M, lcomp, mcomp, mc_tot);
+			gsl_matrix_set(M, mcomp, lcomp, mc_tot);
 
 		}
 
-- 
cgit v1.2.3