aboutsummaryrefslogtreecommitdiff
path: root/src/partialator.c
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2011-04-26 15:13:42 +0200
committerThomas White <taw@physics.org>2012-02-22 15:27:24 +0100
commit6815bab588d3e995bd1319e23614372aa1a07f0c (patch)
tree14b2e12ff2cda12f6ed4baed992cdee6448a44ea /src/partialator.c
parent074034e2dcc34a9a0bab7c979f0deea060a3be21 (diff)
Progress with post refinement...
Diffstat (limited to 'src/partialator.c')
-rw-r--r--src/partialator.c41
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);