diff options
| -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); |