aboutsummaryrefslogtreecommitdiff
path: root/src/ewald.c
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2009-11-13 12:22:42 +0100
committerThomas White <taw@physics.org>2009-11-13 12:22:42 +0100
commitebe633dab856ae1ea857a3c07c6f284c30899874 (patch)
tree30dd0903a4118e67f16d3ef6a7a1c107510533ad /src/ewald.c
parent4e1519eeadbfbfe65d0b000a3af610bec77bbfc3 (diff)
No need to store the "phactors"
Diffstat (limited to 'src/ewald.c')
-rw-r--r--src/ewald.c16
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") */
-
}
}