diff options
| author | Thomas White <taw@physics.org> | 2013-06-29 20:39:07 -0700 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2013-06-29 20:39:07 -0700 |
| commit | 2791ba6e657dc928c9e99ac9937df230af2f6875 (patch) | |
| tree | d7911152a0bbacea800c3ad97d9a1410f97974b7 | |
| parent | 491295f2d71efae499ee0b4b87275260b971e25e (diff) | |
compare_hkl: For 'CRDano', use known mean of zero
| -rw-r--r-- | src/compare_hkl.c | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/src/compare_hkl.c b/src/compare_hkl.c index 47b52419..d3185def 100644 --- a/src/compare_hkl.c +++ b/src/compare_hkl.c @@ -294,10 +294,10 @@ static double fom_overall(struct fom_context *fctx) overall_n++; } } - variance_signal = gsl_stats_variance(overall_along_diagonal, 1, - overall_n); - variance_error = gsl_stats_variance(overall_perpend_diagonal, 1, - overall_n); + variance_signal = gsl_stats_variance_m(overall_along_diagonal, + 1, overall_n, 0.0); + variance_error = gsl_stats_variance_m(overall_perpend_diagonal, + 1, overall_n, 0.0); cc = sqrt(variance_signal / variance_error ); free(overall_along_diagonal); @@ -372,10 +372,10 @@ static double fom_shell(struct fom_context *fctx, int i) perpend_diagonal[j] =( fctx->vec1[i][j] - fctx->vec2[i][j] ) / sqrt(2.0); } - variance_signal = gsl_stats_variance(along_diagonal, 1, - fctx->n[i]); - variance_error = gsl_stats_variance(perpend_diagonal, 1, - fctx->n[i]); + variance_signal = gsl_stats_variance_m(along_diagonal, 1, + fctx->n[i], 0.0); + variance_error = gsl_stats_variance_m(perpend_diagonal, 1, + fctx->n[i], 0.0); return sqrt(variance_signal / variance_error); } |