From 0d9370a585d27aedd2cdf0d3b082be1cfca8e5e8 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Thu, 1 Aug 2013 16:07:42 +0200 Subject: Count filtered eigenvalues --- src/post-refinement.c | 3 +++ 1 file changed, 3 insertions(+) diff --git a/src/post-refinement.c b/src/post-refinement.c index 3716d728..c444ae31 100644 --- a/src/post-refinement.c +++ b/src/post-refinement.c @@ -321,6 +321,7 @@ static void check_eigen(gsl_vector *e_val) const int n = e_val->size; const double max_condition = 1e6; const int verbose = 0; + int n_filt = 0; if ( verbose ) STATUS("Eigenvalues:\n"); vmin = +INFINITY; @@ -336,6 +337,7 @@ static void check_eigen(gsl_vector *e_val) double val = gsl_vector_get(e_val, i); if ( val < vmax/max_condition ) { gsl_vector_set(e_val, i, 0.0); + n_filt++; } } @@ -350,6 +352,7 @@ static void check_eigen(gsl_vector *e_val) if ( verbose ) { STATUS("Condition number: %e / %e = %5.2f\n", vmax, vmin, vmax/vmin); + STATUS("%i out of %i eigenvalues filtered.\n", n_filt, n); } } -- cgit v1.2.3