From a11ebc01ac2c3dbddb1974a171f5755cf3016922 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Thu, 21 Oct 2010 18:26:33 -0700 Subject: compare_hkl: Display numbers of reflections and measurements --- src/compare_hkl.c | 4 ++++ 1 file changed, 4 insertions(+) (limited to 'src/compare_hkl.c') diff --git a/src/compare_hkl.c b/src/compare_hkl.c index 84759cd7..1100a05e 100644 --- a/src/compare_hkl.c +++ b/src/compare_hkl.c @@ -70,6 +70,7 @@ static void plot_shells(const double *ref1, const double *ref2, FILE *fh; double snr_total = 0; int nmeas = 0; + int nmeastot = 0; if ( cell == NULL ) { ERROR("Need the unit cell to plot resolution shells.\n"); @@ -232,9 +233,12 @@ static void plot_shells(const double *ref1, const double *ref2, snr_total += (lookup_intensity(ref1, h, k, l) / lookup_intensity(sigma, h, k, l)); nmeas++; + nmeastot += lookup_count(char_counts, h, k, l); } STATUS("overall = %f\n", snr_total/(double)nmeas); + STATUS("%i measurements in total.\n", nmeastot); + STATUS("%i reflections in total.\n", nmeas); den = 0.0; ctot = 0; -- cgit v1.2.3