Get rid of solve_householder()

1 files changed, 0 insertions, 24 deletions

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 ) {