aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2011-06-22 17:29:55 +0200
committerThomas White <taw@physics.org>2012-02-22 15:27:30 +0100
commit25372dd121fcead4e52f71de422c74558b1900ad (patch)
tree336520c84e1d5a94fa3d312b9e97888dc705682b
parentc8464aea05863bceaf17aeb0640f5bc5d1004f13 (diff)
Get rid of solve_householder()
-rw-r--r--src/post-refinement.c24
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 ) {