From 1f32046f8c429752fe6443e0391044697be4bb08 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Wed, 11 May 2011 16:33:33 +0200
Subject: Don't complain if post refinement equations can't be solved

---
 src/post-refinement.c | 6 ++----
 1 file changed, 2 insertions(+), 4 deletions(-)

(limited to 'src/post-refinement.c')

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);
-- 
cgit v1.2.3