diff options
author | Thomas White <taw@physics.org> | 2009-11-13 12:22:42 +0100 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2009-11-13 12:22:42 +0100 |
commit | ebe633dab856ae1ea857a3c07c6f284c30899874 (patch) | |
tree | 30dd0903a4118e67f16d3ef6a7a1c107510533ad /src/ewald.c | |
parent | 4e1519eeadbfbfe65d0b000a3af610bec77bbfc3 (diff) |
No need to store the "phactors"
Diffstat (limited to 'src/ewald.c')
-rw-r--r-- | src/ewald.c | 16 |
1 files changed, 0 insertions, 16 deletions
diff --git a/src/ewald.c b/src/ewald.c index e9f36596..d543d98a 100644 --- a/src/ewald.c +++ b/src/ewald.c @@ -19,10 +19,6 @@ #include "cell.h" -/* Pulse energy density in J/m^2 */ -#define PULSE_ENERGY_DENSITY (30.0e7) - - void get_ewald(struct image *image) { int x, y; @@ -31,18 +27,9 @@ void get_ewald(struct image *image) k = 1/image->lambda; - /* How many photons are scattered per electron? */ - i0fac = PULSE_ENERGY_DENSITY * pow(THOMSON_LENGTH, 2.0) - / image->xray_energy; - - printf("%e photons are scattered per electron\n", i0fac); - image->qvecs = malloc(image->width * image->height * sizeof(struct threevec)); - image->phactors = malloc(image->width * image->height - * sizeof(double)); - for ( x=0; x<image->width; x++ ) { for ( y=0; y<image->height; y++ ) { @@ -66,9 +53,6 @@ void get_ewald(struct image *image) image->qvecs[x + image->width*y].u = qx; image->qvecs[x + image->width*y].v = qy; image->qvecs[x + image->width*y].w = qz; - - /* Calculate photon factor ("phactor") */ - } } |