diff options
author | Thomas White <taw@physics.org> | 2013-03-13 10:43:42 +0100 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2013-04-17 17:33:48 +0200 |
commit | 2b115b938d8bfba13f4193e329728b09a750b79f (patch) | |
tree | 76a0fd434021d3c5b19eab23b544d88e1495e40a /src | |
parent | 1acd2a02602a98a6286553868ed770427087c021 (diff) |
Add Lorentz factor
Diffstat (limited to 'src')
-rw-r--r-- | src/partial_sim.c | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/src/partial_sim.c b/src/partial_sim.c index 7c93da1b..32e6666a 100644 --- a/src/partial_sim.c +++ b/src/partial_sim.c @@ -100,12 +100,13 @@ static void calculate_partials(Crystal *cr, { signed int h, k, l; Reflection *rfull; - double p, Ip, If; + double L, p, Ip, If; int bin; get_indices(refl, &h, &k, &l); get_asymm(sym, h, k, l, &h, &k, &l); p = get_partiality(refl); + L = get_lorentz(refl); pthread_mutex_lock(full_lock); rfull = find_refl(full, h, k, l); @@ -139,7 +140,7 @@ static void calculate_partials(Crystal *cr, } } - Ip = crystal_get_osf(cr) * p * If; + Ip = crystal_get_osf(cr) * L * p * If; res = resolution(crystal_get_cell(cr), h, k, l); bin = NBINS*2.0*res/max_q; |