diff options
Diffstat (limited to 'src/rejection.c')
| -rw-r--r-- | src/rejection.c | 9 |
1 files changed, 8 insertions, 1 deletions
diff --git a/src/rejection.c b/src/rejection.c index fbe631ab..83738657 100644 --- a/src/rejection.c +++ b/src/rejection.c @@ -139,6 +139,8 @@ static double calculate_cchalf(RefList *template, RefList *full, c = get_contributions(refl); assert(c != NULL); + STATUS("reflection %4i %4i %4i\n", h, k, l); + /* Mean of contributions */ refl_sum = 0.0; for ( j=0; j<c->n_contrib; j++ ) { @@ -152,14 +154,15 @@ static double calculate_cchalf(RefList *template, RefList *full, /* FIXME: Apply corrections */ Ii = get_intensity(c->contribs[j]); + STATUS("contrib %f (crystal %p)\n", Ii, c->contrib_crystals[j]); refl_sum += Ii; nc++; } if ( nc < 2 ) continue; - refl_mean = refl_sum / nc; + STATUS(" refl_mean = %f\n", refl_mean); /* Variance of contributions */ refl_sumsq = 0.0; @@ -180,6 +183,7 @@ static double calculate_cchalf(RefList *template, RefList *full, } refl_var = refl_sumsq/(nc-1.0); + STATUS(" refl_var = %f\n", refl_var); all_sum_var += refl_var; n++; total_contribs += nc; @@ -194,7 +198,10 @@ static double calculate_cchalf(RefList *template, RefList *full, } hsig2E = all_sum_var / n; + STATUS("mean of reflection variances = %f = 0.5*sig2E\n", hsig2E); sig2Y = S / (wSum - 1.0); + STATUS("variance of reflection means = %f = sig2Y\n", sig2Y); + STATUS("%lli total contribs\n", hsig2E, sig2Y, total_contribs); if ( pnref != NULL ) { *pnref = n; |