diff options
| author | Thomas White <taw@physics.org> | 2011-05-11 16:33:33 +0200 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:26 +0100 |
| commit | 1f32046f8c429752fe6443e0391044697be4bb08 (patch) | |
| tree | 2ae0e8833248e5815d089161e045083bbc7e4dbf /src/post-refinement.c | |
| parent | 3a073fe7427632674dc7b2e2fe7f4ba1534510cd (diff) | |
Don't complain if post refinement equations can't be solved
Diffstat (limited to 'src/post-refinement.c')
| -rw-r--r-- | src/post-refinement.c | 6 |
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); |