From b9df4629055f939cfbc4d0d61a4290e7a6974998 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Thu, 23 Jul 2015 11:25:10 +0200
Subject: Fix PR gradients (again - reverts
 22dee8348b165d33262112ed1bbaf7f574b96c7e)

The gradients are correct, the test way wrong.
---
 src/post-refinement.c | 6 ++++++
 1 file changed, 6 insertions(+)

(limited to 'src/post-refinement.c')

diff --git a/src/post-refinement.c b/src/post-refinement.c
index 6630cfc9..16fc1791 100644
--- a/src/post-refinement.c
+++ b/src/post-refinement.c
@@ -449,6 +449,12 @@ static double pr_iterate(Crystal *cr, const RefList *full,
 		/* Calculate all gradients for this reflection */
 		for ( k=0; k<num_params; k++ ) {
 			gradients[k] = gradient(cr, rv[k], refl, pmodel);
+
+			/* Gradients affecting the "p" part of the residual
+			 * need to account for it being ln(p) */
+			if ( (rv[k] != GPARAM_OSF) && (rv[k] != GPARAM_BFAC) ) {
+				gradients[k] /= p;
+			}
 		}
 
 		for ( k=0; k<num_params; k++ ) {
-- 
cgit v1.2.3