aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2010-10-25 18:51:22 +0200
committerThomas White <taw@physics.org>2012-02-22 15:27:03 +0100
commita57e91ffc5b9aed57b2257d6b5e9b5cb59eb77ad (patch)
tree0bce34b84d14d678399b1cc95baa51e6cb95c9bb
parent89fc199308230dd77757f7b02c26974e46712ee9 (diff)
Fix wavelength
-rw-r--r--src/diffraction-gpu.c4
-rw-r--r--src/diffraction.c4
2 files changed, 4 insertions, 4 deletions
diff --git a/src/diffraction-gpu.c b/src/diffraction-gpu.c
index 825f6372..7f05471f 100644
--- a/src/diffraction-gpu.c
+++ b/src/diffraction-gpu.c
@@ -136,8 +136,8 @@ void get_diffraction_gpu(struct gpu_context *gctx, struct image *image,
cell.s[6] = cx; cell.s[7] = cy; cell.s[8] = cz;
/* Calculate wavelength */
- klow = 1.0/(image->lambda + image->beam->bandwidth/2.0);
- khigh = 1.0/(image->lambda - image->beam->bandwidth/2.0);
+ klow = 1.0/(image->lambda*(1.0 + image->beam->bandwidth/2.0));
+ khigh = 1.0/(image->lambda*(1.0 - image->beam->bandwidth/2.0));
bwstep = (khigh-klow) / BWSAMPLING;
/* Orientation */
diff --git a/src/diffraction.c b/src/diffraction.c
index 20490c99..29e30431 100644
--- a/src/diffraction.c
+++ b/src/diffraction.c
@@ -311,8 +311,8 @@ void get_diffraction(struct image *image, int na, int nb, int nc,
/* Needed later for Lorentz calculation */
image->twotheta = malloc(image->width * image->height * sizeof(double));
- klow = 1.0/(image->lambda + image->beam->bandwidth/2.0);
- khigh = 1.0/(image->lambda - image->beam->bandwidth/2.0);
+ klow = 1.0/(image->lambda*(1.0 + image->beam->bandwidth/2.0));
+ khigh = 1.0/(image->lambda*(1.0 - image->beam->bandwidth/2.0));
bwstep = (khigh-klow) / BWSAMPLING;
for ( xs=0; xs<image->width*SAMPLING; xs++ ) {