From a0c10794ee1276a61a221a398f96df89658149dd Mon Sep 17 00:00:00 2001 From: Alexander Schmidt Date: Mon, 5 Jan 2015 17:52:16 +0100 Subject: [PATCH] seems to be working. --- u1casc-ordinary/obs_diagcorr.hpp | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) diff --git a/u1casc-ordinary/obs_diagcorr.hpp b/u1casc-ordinary/obs_diagcorr.hpp index a8d88b5..d52a6cd 100644 --- a/u1casc-ordinary/obs_diagcorr.hpp +++ b/u1casc-ordinary/obs_diagcorr.hpp @@ -103,7 +103,6 @@ void obs_diagcorr::cdiag (gsl_vector *v, gsl_matrix_complex *evec, gsl_matrix_co gsl_eigen_hermv_sort(v, evec, GSL_EIGEN_SORT_VAL_ASC); gsl_eigen_hermv_free(wspace); - gsl_matrix_complex_free(evec); } void obs_diagcorr::_finish() { @@ -146,7 +145,7 @@ void obs_diagcorr::_finish() { ); } cdiag(jackres, evecres, tmpmatrix); - + for (int imeas=0; imeascomargs.nmeas; imeas++) { gsl_matrix_complex_memcpy (tmpmatrix, totalval); gsl_matrix_complex_sub (tmpmatrix, measurements[imeas][itsep]); @@ -164,7 +163,7 @@ void obs_diagcorr::_finish() { discopart ) ); - } + } cdiag(tmpvec, tmpmatrix2, tmpmatrix); gsl_vector_sub(tmpvec, jackres); -- 2.39.5