aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2011-05-11 16:33:33 +0200
committerThomas White <taw@physics.org>2012-02-22 15:27:26 +0100
commit1f32046f8c429752fe6443e0391044697be4bb08 (patch)
tree2ae0e8833248e5815d089161e045083bbc7e4dbf
parent3a073fe7427632674dc7b2e2fe7f4ba1534510cd (diff)
Don't complain if post refinement equations can't be solved
-rw-r--r--src/post-refinement.c6
1 files changed, 2 insertions, 4 deletions
diff --git a/src/post-refinement.c b/src/post-refinement.c
index 19462edf..4ebaffd9 100644
--- a/src/post-refinement.c
+++ b/src/post-refinement.c
@@ -315,9 +315,7 @@ static double pr_iterate(struct image *image, const RefList *full,
shifts = gsl_vector_alloc(NUM_PARAMS);
max_shift = 0.0;
- if ( gsl_linalg_HH_solve(M, v, shifts) ) {
- ERROR("Couldn't solve normal equations!\n");
- } else {
+ if ( gsl_linalg_HH_solve(M, v, shifts) == 0 ) {
for ( param=0; param<NUM_PARAMS; param++ ) {
double shift = gsl_vector_get(shifts, param);
@@ -325,7 +323,7 @@ static double pr_iterate(struct image *image, const RefList *full,
if ( fabs(shift) > max_shift ) max_shift = fabs(shift);
}
- }
+ } /* else problem with the equations, so leave things as they were */
gsl_matrix_free(M);
gsl_vector_free(v);