From 1f83bf629093e45c2c16d91b9947158127265a86 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Fri, 17 Apr 2015 11:11:01 +0200
Subject: Add a simple rejection criterion to the prediction refinement

---
 src/process_image.c | 7 +++++++
 1 file changed, 7 insertions(+)

(limited to 'src')

diff --git a/src/process_image.c b/src/process_image.c
index c00976e1..c159465f 100644
--- a/src/process_image.c
+++ b/src/process_image.c
@@ -71,6 +71,13 @@ static void try_refine_autoR(struct image *image, Crystal *cr)
 	refine_radius(cr, image);
 	new_R = crystal_get_profile_radius(cr);
 
+	if ( (new_R > 0.02e9) || (new_R > 1.2*old_R) ) {
+		STATUS("Rejecting refinement solution with "
+		       "R %.3f -> %.3f nm^-1\n", old_R, new_R);
+		crystal_set_user_flag(cr, 1);
+		return;
+	}
+
 	snprintf(notes, 1024, "predict_refine/R old = %.5f new = %.5f nm^-1",
 	                      old_R/1e9, new_R/1e9);
 	crystal_add_notes(cr, notes);
-- 
cgit v1.2.3