diff options
author | Thomas White <taw@physics.org> | 2020-08-01 15:27:52 +0200 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2020-08-01 15:27:52 +0200 |
commit | 4942b91e76e84a296b79fa1f967a59d1cbfd8315 (patch) | |
tree | 376b4d16389199da2dbec6d3c9ba3c14ba7058b7 | |
parent | 189da15810deabd739d7c11c6e95fea55739fe60 (diff) |
Not used, and doesn't compile with recent GSL
-rw-r--r-- | src/refine-lmder.c | 5 |
1 files changed, 0 insertions, 5 deletions
diff --git a/src/refine-lmder.c b/src/refine-lmder.c index ba99eaf..5b9e935 100644 --- a/src/refine-lmder.c +++ b/src/refine-lmder.c @@ -321,7 +321,6 @@ void refine_lmder(AtomicModel *model, ReflectionList *reflections, RefinementSpe unsigned int j; RefinementPair pair; int iter_status, conv_status; - gsl_matrix *covar; unsigned int n_params, n_atoms, idx; double scale; ReflectionList *calc; @@ -372,7 +371,6 @@ void refine_lmder(AtomicModel *model, ReflectionList *reflections, RefinementSpe pair.model = model; pair.reflections = reflections; pair.spec = spec; f.params = &pair; - covar = gsl_matrix_alloc(f.p, f.p); gsl_multifit_fdfsolver_set(s, &f, coordinates); printf("initial: scale=%f, |f(x)|=%g\n", gsl_vector_get(s->x, 0), gsl_blas_dnrm2(s->f)); @@ -429,9 +427,6 @@ void refine_lmder(AtomicModel *model, ReflectionList *reflections, RefinementSpe printf("%i iterations performed\n", iter); if ( iter == MAX_REFINEMENT_ITERATIONS ) printf("Reached maximum allowed number of iterations"); - gsl_multifit_covar(s->J, 0.0, covar); - gsl_matrix_free(covar); - gsl_multifit_fdfsolver_free(s); gsl_vector_free(coordinates); |