aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2015-05-08 18:29:43 +0200
committerThomas White <taw@physics.org>2015-05-13 13:48:33 +0200
commit1871cb71b20b29ea39082e470ffb73bee208e4f5 (patch)
treed8d946d1270d1c10f54eab06d63fab0969144ef8
parent22429efd8922251c78006e8d53155c89adfa8f1f (diff)
Fix guide_dev()
-rw-r--r--src/post-refinement.c5
1 files changed, 3 insertions, 2 deletions
diff --git a/src/post-refinement.c b/src/post-refinement.c
index a31ed221..57d3efd3 100644
--- a/src/post-refinement.c
+++ b/src/post-refinement.c
@@ -486,7 +486,7 @@ static double guide_dev(Crystal *cr, const RefList *full)
refl != NULL;
refl = next_refl(refl, iter) ) {
- double G, p;
+ double G, p, L;
signed int h, k, l;
Reflection *full_version;
double I_full, I_partial;
@@ -505,13 +505,14 @@ static double guide_dev(Crystal *cr, const RefList *full)
G = crystal_get_osf(cr);
p = get_partiality(refl);
+ L = get_lorentz(refl);
I_partial = get_intensity(refl);
I_full = get_intensity(full_version);
//STATUS("%3i %3i %3i %5.2f %5.2f %5.2f %5.2f %5.2f\n",
// h, k, l, G, p, I_partial, I_full,
// I_partial - p*G*I_full);
- dev += pow(I_partial - p*G*I_full, 2.0);
+ dev += pow(I_partial - p*I_full/(G*L), 2.0);
}