From e9254176f8a59e74ff7413c5112a42d226ad1e04 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Tue, 10 May 2011 16:07:22 +0200
Subject: Add some sanity checks

---
 src/partialator.c | 33 +++++++++++++++++++++++++++++++++
 1 file changed, 33 insertions(+)

(limited to 'src/partialator.c')

diff --git a/src/partialator.c b/src/partialator.c
index a410c507..74ba95ac 100644
--- a/src/partialator.c
+++ b/src/partialator.c
@@ -347,6 +347,39 @@ int main(int argc, char *argv[])
 	full = scale_intensities(images, n_usable_patterns, sym,
 	                         scalable, cref);
 
+	for ( i=0; i<num_items(scalable); i++ ) {
+		Reflection *f;
+		struct refl_item *it = get_item(scalable, i);
+		f = find_refl(full, it->h, it->k, it->l);
+		if ( f == NULL ) {
+			ERROR("%3i %3i %3i was designated scalable, but no"
+			      " full intensity was recorded.\n",
+			      it->h, it->k, it->l);
+		}
+	}
+
+	for ( i=0; i<n_usable_patterns; 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;
+
+			if ( !get_scalable(refl) ) continue;
+			get_indices(refl, &h, &k, &l);
+
+			if ( find_item(scalable, h, k, l) == 0 ) {
+				ERROR("%3i %3i %3i in image %i is scalable"
+				      " but is not in the list of scalable"
+				      " reflections.\n", h, k, l, i);
+			}
+		}
+	}
+
 	/* Iterate */
 	for ( i=0; i<n_iter; i++ ) {
 
-- 
cgit v1.2.3