From c274adc00703bf725137dd846a1d0f32aaa8bb1f Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Thu, 20 Jun 2019 16:28:53 +0200
Subject: Modify residual to give clearer minimum, particularly as function of
 profile radius

---
 src/merge.c | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

(limited to 'src')

diff --git a/src/merge.c b/src/merge.c
index 33567ecd..e56c3378 100644
--- a/src/merge.c
+++ b/src/merge.c
@@ -402,8 +402,8 @@ double residual(Crystal *cr, const RefList *full, int free,
 			w = 0.0;
 		}
 
-		num += fabs(int1 - int2) * w;
-		den += fabs(int1 + int2) * w/2.0;
+		num += w*fabs(int1 - int2)/((int1+int2)/2.0);
+		den += w;
 
 		n_used++;
 
-- 
cgit v1.2.3