aboutsummaryrefslogtreecommitdiff
path: root/src/compare_hkl.c
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2013-06-29 20:39:07 -0700
committerThomas White <taw@physics.org>2013-06-29 20:39:07 -0700
commit2791ba6e657dc928c9e99ac9937df230af2f6875 (patch)
treed7911152a0bbacea800c3ad97d9a1410f97974b7 /src/compare_hkl.c
parent491295f2d71efae499ee0b4b87275260b971e25e (diff)
compare_hkl: For 'CRDano', use known mean of zero
Diffstat (limited to 'src/compare_hkl.c')
-rw-r--r--src/compare_hkl.c16
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);
}