PMNS_Sterile.cxx
//......................................................................
///
/// GSL Copy-constructor.
///
/// Allocate new memory for GSL objects when copying.
///
/// @param other - GSL object to copy
///
GSL Copy-constructor.

Allocate new memory for GSL objects when copying.
GSL Copy - constructor . Allocate new memory for GSL objects when copying .
GSL_EinSys::GSL_EinSys(const GSL_EinSys &other) :
fNumNus(other.fNumNus), 
fEvalGSL(0), fEvecGSL(0), H_GSL(0), W_GSL(0)
{
  fEvalGSL = gsl_vector_alloc(fNumNus);
  fEvecGSL = gsl_matrix_complex_alloc(fNumNus, fNumNus);
  H_GSL = gsl_matrix_complex_alloc(fNumNus, fNumNus);
  W_GSL = gsl_eigen_hermv_alloc(fNumNus);
  gsl_vector_memcpy(fEvalGSL, other.fEvalGSL);
  gsl_matrix_complex_memcpy(fEvecGSL, other.fEvecGSL);
  gsl_matrix_complex_memcpy(H_GSL, other.H_GSL);
}
GSL_EinSys :: GSL_EinSys ( const GSL_EinSys & other ) : fNumNus ( other . fNumNus ) , fEvalGSL ( 0 ) , fEvecGSL ( 0 ) , H_GSL ( 0 ) , W_GSL ( 0 ) { fEvalGSL = gsl_vector_alloc ( fNumNus ) ; fEvecGSL = gsl_matrix_complex_alloc ( fNumNus , fNumNus ) ; H_GSL = gsl_matrix_complex_alloc ( fNumNus , fNumNus ) ; W_GSL = gsl_eigen_hermv_alloc ( fNumNus ) ; gsl_vector_memcpy ( fEvalGSL , other . fEvalGSL ) ; gsl_matrix_complex_memcpy ( fEvecGSL , other . fEvecGSL ) ; gsl_matrix_complex_memcpy ( H_GSL , other . H_GSL ) ; }
......................................................................
GSL Copy-constructor.
...................................................................... GSL Copy - constructor .
// Allocate memory for the GSL objects