From: Alexander Schmidt <alex@treefish.org>
Date: Mon, 5 Jan 2015 16:52:16 +0000 (+0100)
Subject: seems to be working.
X-Git-Url: http://git.treefish.org/~alex/phys/u1casc.git/commitdiff_plain

seems to be working.
---

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; imeas<O815->comargs.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);