From bd5d2ec01be4aef4731498424fe777e3cdb797ef Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Mon, 26 Feb 2018 17:21:50 +0100
Subject: Add --force-bandwidth and --force-radius

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

(limited to 'src/post-refinement.c')

diff --git a/src/post-refinement.c b/src/post-refinement.c
index 60810349..acafb8aa 100644
--- a/src/post-refinement.c
+++ b/src/post-refinement.c
@@ -330,6 +330,13 @@ static double residual_f(const gsl_vector *v, void *pp)
 		return GSL_NAN;
 	}
 
+	/* Can happen with grid scans and certain --force-radius values */
+	if ( fabs(crystal_get_profile_radius(cr)) < 0.0000001e9 ) {
+		crystal_free(cr);
+		if ( pv->verbose ) STATUS("radius very small\n");
+		return GSL_NAN;
+	}
+
 	if ( im.lambda <= 0.0 ) {
 		crystal_free(cr);
 		if ( pv->verbose ) STATUS("lambda < 0\n");
-- 
cgit v1.2.3