diff options
author | Thomas White <taw@physics.org> | 2011-04-26 15:13:42 +0200 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2012-02-22 15:27:24 +0100 |
commit | 6815bab588d3e995bd1319e23614372aa1a07f0c (patch) | |
tree | 14b2e12ff2cda12f6ed4baed992cdee6448a44ea /src/partialator.c | |
parent | 074034e2dcc34a9a0bab7c979f0deea060a3be21 (diff) |
Progress with post refinement...
Diffstat (limited to 'src/partialator.c')
-rw-r--r-- | src/partialator.c | 41 |
1 files changed, 4 insertions, 37 deletions
diff --git a/src/partialator.c b/src/partialator.c index 54022e9e..9f931496 100644 --- a/src/partialator.c +++ b/src/partialator.c @@ -136,43 +136,12 @@ static void refine_all(struct image *images, int n_total_patterns, qargs.task_defaults = task_defaults; qargs.n = 0; qargs.n_done = 0; - qargs.n_total_patterns = n_total_patterns; - qargs.images = images; + /* FIXME: Not refining the first image, for now */ + qargs.n_total_patterns = n_total_patterns-1; + qargs.images = images+1; run_threads(nthreads, refine_image, get_image, done_image, - &qargs, n_total_patterns, 0, 0, 0); -} - - -/* Decide which reflections can be scaled */ -static void select_scalable_reflections(struct image *images, int n) -{ - int m; - int n_scalable = 0; - - for ( m=0; m<n; m++ ) { - - Reflection *refl; - RefListIterator *iter; - - for ( refl = first_refl(images[m].reflections, &iter); - refl != NULL; - refl = next_refl(refl, iter) ) { - - int scalable = 1; - double v; - - if ( get_partiality(refl) < 0.1 ) scalable = 0; - v = fabs(get_intensity(refl)); - if ( v < 0.1 ) scalable = 0; - - set_scalable(refl, scalable); - if ( scalable ) n_scalable++; - - } - - } - STATUS("%i reflections selected as scalable.\n", n_scalable); + &qargs, n_total_patterns-1, 0, 0, 0); } @@ -369,7 +338,6 @@ int main(int argc, char *argv[]) /* Make initial estimates */ STATUS("Performing initial scaling.\n"); - select_scalable_reflections(images, n_total_patterns); full = scale_intensities(images, n_usable_patterns, sym, obs, cref); /* Iterate */ @@ -401,7 +369,6 @@ int main(int argc, char *argv[]) /* Re-estimate all the full intensities */ reflist_free(full); - select_scalable_reflections(images, n_usable_patterns); full = scale_intensities(images, n_usable_patterns, sym, obs, cref); |