From efd6562f9156ddff2fe073c97b2ddbf25c45688e Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Mon, 4 Jul 2011 17:50:56 +0200
Subject: Separate "refinable" and "scalable" concepts

---
 src/partialator.c | 45 +++++++++++++++++++++++++++++++++++++++++----
 1 file changed, 41 insertions(+), 4 deletions(-)

(limited to 'src/partialator.c')

diff --git a/src/partialator.c b/src/partialator.c
index 167f736d..43c4bda0 100644
--- a/src/partialator.c
+++ b/src/partialator.c
@@ -191,6 +191,35 @@ static int select_scalable_reflections(RefList *list, ReflItemList *sc_l)
 }
 
 
+static void select_reflections_for_refinement(struct image *images, int n,
+                                              ReflItemList *scalable)
+{
+	int i;
+
+	for ( i=0; i<n; i++ ) {
+
+		Reflection *refl;
+		RefListIterator *iter;
+
+		for ( refl = first_refl(images[i].reflections, &iter);
+		      refl != NULL;
+		      refl = next_refl(refl, iter) )
+		{
+			signed int h, k, l;
+			int sc;
+
+			get_indices(refl, &h, &k, &l);
+			sc = get_scalable(refl);
+
+			if ( sc && find_item(scalable, h, k, l) ) {
+				set_refinable(refl, 1);
+			}
+		}
+
+	}
+}
+
+
 int main(int argc, char *argv[])
 {
 	int c;
@@ -397,8 +426,8 @@ int main(int argc, char *argv[])
 		reflist_free(cur->reflections);
 		cur->reflections = as;
 
-		update_partialities(cur, sym,
-		                    &n_expected, &n_found, &n_notfound);
+		predict_corresponding_reflections(cur, sym, &n_expected,
+		                                 &n_found, &n_notfound);
 
 		nobs += select_scalable_reflections(cur->reflections, scalable);
 
@@ -419,6 +448,8 @@ int main(int argc, char *argv[])
 	full = scale_intensities(images, n_usable_patterns, sym,
 	                         scalable, cref, reference);
 
+	select_reflections_for_refinement(images, n_usable_patterns, scalable);
+
 	for ( i=0; i<num_items(scalable); i++ ) {
 		Reflection *f;
 		struct refl_item *it = get_item(scalable, i);
@@ -459,6 +490,7 @@ int main(int argc, char *argv[])
 		FILE *fhg;
 		FILE *fhp;
 		char filename[1024];
+		int j;
 
 		STATUS("Post refinement cycle %i of %i\n", i+1, n_iter);
 
@@ -484,9 +516,14 @@ int main(int argc, char *argv[])
 
 		nobs = 0;
 		clear_items(scalable);
-		for ( i=0; i<n_usable_patterns; i++ ) {
+		for ( j=0; j<n_usable_patterns; j++ ) {
+
+			struct image *cur = &images[j];
+
+			predict_corresponding_reflections(cur, sym, &n_expected,
+		                                          &n_found,
+		                                          &n_notfound);
 
-			struct image *cur = &images[i];
 			nobs += select_scalable_reflections(cur->reflections,
 			                                    scalable);
 
-- 
cgit v1.2.3