Various improvements

1 files changed, 7 insertions, 5 deletions

diff --git a/src/statistics.c b/src/statistics.c
index 9219b069..50bfd513 100644
--- a/src/statistics.c
+++ b/src/statistics.c
@@ -29,7 +29,8 @@ static double stat_scale_intensity(double *obs, double *calc, unsigned int *c,
 	double bot = 0.0;
 	int i;
 
-	for ( i=0; i<size; i++ ) {
+	/* Start from i=1 -> skip central beam */
+	for ( i=1; i<size; i++ ) {
 
 		if ( c[i] > 0 ) {
 			double obsi;
@@ -54,16 +55,17 @@ double stat_r2(double *obs, double *calc, unsigned int *c, int size,
 	scale = stat_scale_intensity(obs, calc, c, size);
 	*scalep = scale;
 
-	for ( i=0; i<size; i++ ) {
+	/* Start from i=1 -> skip central beam */
+	for ( i=1; i<size; i++ ) {
 
 		if ( c[i] > 0 ) {
 			double obsi;
 			obsi = obs[i] / (double)c[i];
-			top += fabs(obsi/scale - calc[i]);
-			bot += obsi/scale;
+			top += pow(fabs(obsi - scale*calc[i]), 2.0);
+			bot += pow(obsi, 2.0);
 		}
 
 	} /* else reflection not measured so don't worry about it */
 
-	return top/bot;
+	return sqrt(top/bot);
 }