From ebe633dab856ae1ea857a3c07c6f284c30899874 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Fri, 13 Nov 2009 12:22:42 +0100 Subject: No need to store the "phactors" --- src/ewald.c | 16 ---------------- 1 file changed, 16 deletions(-) (limited to 'src/ewald.c') 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; xwidth; x++ ) { for ( y=0; yheight; 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") */ - } } -- cgit v1.2.3