aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2020-08-01 15:27:52 +0200
committerThomas White <taw@physics.org>2020-08-01 15:27:52 +0200
commit4942b91e76e84a296b79fa1f967a59d1cbfd8315 (patch)
tree376b4d16389199da2dbec6d3c9ba3c14ba7058b7 /src
parent189da15810deabd739d7c11c6e95fea55739fe60 (diff)
Remove covariance stuffHEADmaster
Not used, and doesn't compile with recent GSL
Diffstat (limited to 'src')
-rw-r--r--src/refine-lmder.c5
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);