diff options
| author | Thomas White <taw@physics.org> | 2013-03-12 12:26:14 +0100 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2013-04-17 17:33:48 +0200 |
| commit | 1acd2a02602a98a6286553868ed770427087c021 (patch) | |
| tree | af6e75ebc7b8e162370596c8b95bb1df68d85519 /src/post-refinement.c | |
| parent | d3b353ad71a4d53e108d8e1dfeaff09ebaf14e17 (diff) | |
Fix radius gradient
Diffstat (limited to 'src/post-refinement.c')
| -rw-r--r-- | src/post-refinement.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/post-refinement.c b/src/post-refinement.c index 32a93f38..9e4649a2 100644 --- a/src/post-refinement.c +++ b/src/post-refinement.c @@ -81,7 +81,7 @@ static double partiality_rgradient(double r, double profile_radius) dpdq = 6.0*(q-pow(q, 2.0)); /* dq/drad */ - dqdrad = 0.5 * (1.0 - r * pow(profile_radius, -2.0)); + dqdrad = -0.5 * r * pow(profile_radius, -2.0); return dpdq * dqdrad; } @@ -162,7 +162,7 @@ double gradient(Crystal *cr, int k, Reflection *refl, PartialityModel pmodel) case REF_R : gr = partiality_rgradient(rlow, r); - gr += partiality_rgradient(rhigh, r); + gr -= partiality_rgradient(rhigh, r); return gr; /* Cell parameters and orientation */ |