From 500b2962f771103cbaccb7c12573f6a54dca0a93 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Thu, 22 Feb 2018 16:06:18 +0100
Subject: Take modulus of profile radius

A bit nasty, but avoids lots of NANs
---
 src/post-refinement.c | 6 ------
 1 file changed, 6 deletions(-)

(limited to 'src')

diff --git a/src/post-refinement.c b/src/post-refinement.c
index 1f867d35..5c6faeaa 100644
--- a/src/post-refinement.c
+++ b/src/post-refinement.c
@@ -324,12 +324,6 @@ static double residual_f(const gsl_vector *v, void *pp)
 	crystal_set_image(cr, &im);
 	apply_parameters(v, pv->initial, pv->rv, cr);
 
-	if ( crystal_get_profile_radius(cr) <= 0.0 ) {
-		crystal_free(cr);
-		if ( pv->verbose ) STATUS("R < 0\n");
-		return GSL_NAN;
-	}
-
 	if ( im.lambda <= 0.0 ) {
 		crystal_free(cr);
 		if ( pv->verbose ) STATUS("lambda < 0\n");
-- 
cgit v1.2.3