aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--tests/integration_check.c79
1 files changed, 45 insertions, 34 deletions
diff --git a/tests/integration_check.c b/tests/integration_check.c
index 51725634..9d77b562 100644
--- a/tests/integration_check.c
+++ b/tests/integration_check.c
@@ -54,6 +54,8 @@ int main(int argc, char *argv[])
const int ir_inn = 2;
const int ir_mid = 4;
const int ir_out = 6;
+ double int_sum = 0.0;
+ int i;
fh = fopen("/dev/urandom", "r");
fread(&seed, sizeof(seed), 1, fh);
@@ -101,43 +103,52 @@ int main(int argc, char *argv[])
image.n_crystals = 0;
image.crystals = NULL;
- for ( fs=0; fs<w; fs++ ) {
- for ( ss=0; ss<h; ss++ ) {
- image.dp[0][fs+w*ss] = 10.0*poisson_noise(40);
- if ( (fs-64)*(fs-64) + (ss-64)*(ss-64) > 2 ) continue;
- //image.dp[0][fs+w*ss] += 10.0*poisson_noise(10);
- }
- }
+ for ( i=0; i<1000; i++ ) {
+
+ for ( fs=0; fs<w; fs++ ) {
+ for ( ss=0; ss<h; ss++ ) {
+ image.dp[0][fs+w*ss] = 10.0*poisson_noise(40);
+ if ( (fs-64)*(fs-64) + (ss-64)*(ss-64) > 2 ) continue;
+ //image.dp[0][fs+w*ss] += 10.0*poisson_noise(10);
+ }
+ }
+
+ list = reflist_new();
+ refl = add_refl(list, 0, 0, 0);
+ set_detector_pos(refl, 0.0, 64, 64);
+ cell = cell_new();
+ cell_set_lattice_type(cell, L_CUBIC);
+ cell_set_centering(cell, 'P');
+ cell_set_parameters(cell, 800.0e-10, 800.0e-10, 800.0e-10,
+ deg2rad(90.0), deg2rad(90.0), deg2rad(90.0));
+ cell = cell_rotate(cell, random_quaternion());
+
+ ic.halfw = ir_out;
+ ic.image = &image;
+ ic.k = 1.0/image.lambda;
+ ic.n_saturated = 0;
+ ic.n_implausible = 0;
+ ic.cell = cell;
+ ic.ir_inn = ir_inn;
+ ic.ir_mid = ir_mid;
+ ic.ir_out = ir_out;
+ ic.limit = 0.0;
+ ic.meth = INTEGRATION_RINGS;
+ if ( init_intcontext(&ic) ) {
+ ERROR("Failed to initialise integration.\n");
+ return 1;
+ }
+ setup_ring_masks(&ic, ir_inn, ir_mid, ir_out);
+
+ integrate_rings_once(refl, &image, &ic, cell);
+
+ cell_free(cell);
+
+ int_sum += get_intensity(refl);
- list = reflist_new();
- refl = add_refl(list, 0, 0, 0);
- set_detector_pos(refl, 0.0, 64, 64);
- cell = cell_new();
- cell_set_lattice_type(cell, L_CUBIC);
- cell_set_centering(cell, 'P');
- cell_set_parameters(cell, 800.0e-10, 800.0e-10, 800.0e-10,
- deg2rad(90.0), deg2rad(90.0), deg2rad(90.0));
- cell = cell_rotate(cell, random_quaternion());
- cell_print(cell);
-
- ic.halfw = ir_out;
- ic.image = &image;
- ic.k = 1.0/image.lambda;
- ic.n_saturated = 0;
- ic.n_implausible = 0;
- ic.cell = cell;
- ic.ir_inn = ir_inn;
- ic.ir_mid = ir_mid;
- ic.ir_out = ir_out;
- ic.limit = 0.0;
- ic.meth = INTEGRATION_RINGS;
- if ( init_intcontext(&ic) ) {
- ERROR("Failed to initialise integration.\n");
- return 1;
}
- setup_ring_masks(&ic, ir_inn, ir_mid, ir_out);
- integrate_rings_once(refl, &image, &ic, cell);
+ STATUS("mean value = %f\n", int_sum/1000.0);
free(image.beam);
free(image.det->panels);