From 854f85f90d6cc8424bc0459372720ed951165211 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Fri, 17 Jul 2015 15:35:02 +0200
Subject: partialator: Detect (in a terminal fashion) NAN partialities

---
 libcrystfel/src/geometry.c | 11 +++++++++++
 1 file changed, 11 insertions(+)

(limited to 'libcrystfel/src/geometry.c')

diff --git a/libcrystfel/src/geometry.c b/libcrystfel/src/geometry.c
index 048dc9f7..3325923c 100644
--- a/libcrystfel/src/geometry.c
+++ b/libcrystfel/src/geometry.c
@@ -236,6 +236,17 @@ static Reflection *check_reflection(struct image *image, Crystal *cryst,
 	/* Calculate partiality */
 	part = partiality(pmodel, rlow, rhigh, pr);
 
+	if ( isnan(part) ) {
+		ERROR("Assigning NAN partiality!\n");
+		ERROR("%3i %3i %3i  rlow = %e, rhigh = %e\n",
+		      h, k, l, rlow, rhigh);
+		ERROR("div + m = %e, R = %e, bw = %e\n", del, pr, image->bw);
+		ERROR("D = %e\n", rlow - rhigh);
+		ERROR("klow = %e  khigh = %e\n", klow, khigh);
+		ERROR("distance(..) = %e\n", distance(cet, cez, tl, zl));
+		abort();
+	}
+
 	if ( updateme == NULL ) {
 		refl = reflection_new(h, k, l);
 	} else {
-- 
cgit v1.2.3