aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorThomas White <taw@bitwiz.org.uk>2010-10-21 18:26:33 -0700
committerThomas White <taw@physics.org>2012-02-22 15:27:03 +0100
commita11ebc01ac2c3dbddb1974a171f5755cf3016922 (patch)
treeee9581cb10922ed100fddbe894198df1ea8f4d48
parentf7d612eb193e718490276b05b55af72cc3e70aac (diff)
compare_hkl: Display numbers of reflections and measurements
-rw-r--r--src/compare_hkl.c4
1 files changed, 4 insertions, 0 deletions
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 <snr> = %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;