aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorThomas White <taw@bitwiz.org.uk>2010-12-17 18:03:14 -0800
committerThomas White <taw@physics.org>2012-02-22 15:27:10 +0100
commita764d0a9512a7bd5c0e17af70b3c55e5fe7f9f71 (patch)
treec8515c03eb8ffa00e15bc3cfb379c0dd63994d12
parentec152b28b6405af9ac0bb651fe645873ad5bfa3f (diff)
Scaling matrix is symmetric
-rw-r--r--src/hrs-scaling.c2
1 files changed, 2 insertions, 0 deletions
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);
}