diff options
| author | Thomas White <taw@physics.org> | 2015-05-08 18:29:43 +0200 |
|---|---|---|
| committer | Thomas White <taw@physics.org> | 2015-05-13 13:48:33 +0200 |
| commit | 1871cb71b20b29ea39082e470ffb73bee208e4f5 (patch) | |
| tree | d8d946d1270d1c10f54eab06d63fab0969144ef8 | |
| parent | 22429efd8922251c78006e8d53155c89adfa8f1f (diff) | |
Fix guide_dev()
| -rw-r--r-- | src/post-refinement.c | 5 |
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); } |