Commit ee729abd authored by cortes_c's avatar cortes_c
Browse files

Merge branch 'matched-gauss-fixes' of gitlab.psi.ch:OPAL/src into matched-gauss-fixes

parents 509c3ec4 3073a6d9
......@@ -2534,9 +2534,12 @@ void Distribution::generateGaussZ(size_t numberOfParticles) {
*gmsg << " \\\\" << endl;
}
#endif
/*
//Sets the GSL error handler off, exception will be handled internaly with a renormalization method
gsl_set_error_handler_off();
*/
int errcode = gsl_linalg_cholesky_decomp(corMat);
/*
double rn = 1e-12;
while (errcode == GSL_EDOM) {
......@@ -2562,6 +2565,7 @@ void Distribution::generateGaussZ(size_t numberOfParticles) {
}
//Sets again the standard GSL error handler on
gsl_set_error_handler(NULL);
*/
//Just to be sure
if (errcode == GSL_EDOM) {
throw OpalException("Distribution::GenerateGaussZ",
......@@ -2573,7 +2577,6 @@ void Distribution::generateGaussZ(size_t numberOfParticles) {
gsl_matrix_set (corMat, i, j, 0.0);
}
}
#define DISTDBG2
#ifdef DISTDBG2
*gmsg << "* m after gsl_linalg_cholesky_decomp" << endl;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment