From ee2c11a10b39210c75a9c9002a83ec0dda95418e Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Wed, 22 Jul 2015 17:23:30 +0200
Subject: Add weightings to PR

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

diff --git a/src/post-refinement.c b/src/post-refinement.c
index a10ee781..517ed6a5 100644
--- a/src/post-refinement.c
+++ b/src/post-refinement.c
@@ -553,7 +553,7 @@ static double pr_iterate(Crystal *cr, const RefList *full,
 		s = resolution(crystal_get_cell(cr), ha, ka, la);
 
 		/* Calculate the weight for this reflection */
-		w = 1.0;
+		w = (s/1e9)*(s/1e9) / (esd*esd);
 
 		/* Calculate all gradients for this reflection */
 		for ( k=0; k<num_params; k++ ) {
@@ -725,7 +725,7 @@ static double residual(Crystal *cr, const RefList *full, int verbose, int free,
 
 		fx = exp(G)*p*exp(-B*s*s)*I_full/L;
 		dc = I_partial - fx;
-		w = 1.0;
+		w = (s/1e9)*(s/1e9)/(esd*esd);
 		dev += w*dc*dc;
 		n_used++;
 
-- 
cgit v1.2.3