diff options
| author | Thomas White <taw@bitwiz.org.uk> | 2010-12-17 18:03:14 -0800 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:10 +0100 |
| commit | a764d0a9512a7bd5c0e17af70b3c55e5fe7f9f71 (patch) | |
| tree | c8515c03eb8ffa00e15bc3cfb379c0dd63994d12 /src | |
| parent | ec152b28b6405af9ac0bb651fe645873ad5bfa3f (diff) | |
Scaling matrix is symmetric
Diffstat (limited to 'src')
| -rw-r--r-- | src/hrs-scaling.c | 2 |
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); } |