From 25372dd121fcead4e52f71de422c74558b1900ad Mon Sep 17 00:00:00 2001 From: Thomas White Date: Wed, 22 Jun 2011 17:29:55 +0200 Subject: Get rid of solve_householder() --- src/post-refinement.c | 24 ------------------------ 1 file changed, 24 deletions(-) (limited to 'src/post-refinement.c') diff --git a/src/post-refinement.c b/src/post-refinement.c index 1782e7da..8f7b9862 100644 --- a/src/post-refinement.c +++ b/src/post-refinement.c @@ -298,29 +298,6 @@ static gsl_vector *solve_svd(gsl_vector *v, gsl_matrix *M) } -static gsl_vector *solve_householder(gsl_vector *v, gsl_matrix *M) -{ - int n, err; - gsl_vector *shifts; - - n = v->size; - if ( v->size != M->size1 ) return NULL; - if ( v->size != M->size2 ) return NULL; - - shifts = gsl_vector_alloc(n); - - err = gsl_linalg_HH_solve(M, v, shifts); - - if ( err != 0 ) { - ERROR("Failed to solve equations: %s\n", gsl_strerror(err)); - gsl_vector_free(shifts); - return NULL; - } - - return shifts; -} - - /* Perform one cycle of post refinement on 'image' against 'full' */ static double pr_iterate(struct image *image, const RefList *full, const char *sym) @@ -415,7 +392,6 @@ static double pr_iterate(struct image *image, const RefList *full, } max_shift = 0.0; - //shifts = solve_householder(v, M); shifts = solve_svd(v, M); if ( shifts != NULL ) { -- cgit v1.2.3