From b7540780d46b92f9ce96de74eb4c2e05b2369983 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Wed, 13 May 2015 13:43:42 +0200
Subject: Squash NaN shifts in prediction refinement

---
 src/predict-refine.c | 3 +++
 1 file changed, 3 insertions(+)

(limited to 'src/predict-refine.c')

diff --git a/src/predict-refine.c b/src/predict-refine.c
index 97371172..2d845feb 100644
--- a/src/predict-refine.c
+++ b/src/predict-refine.c
@@ -618,6 +618,9 @@ static int iterate(struct reflpeak *rps, int n, UnitCell *cell,
 
 	for ( i=0; i<num_params; i++ ) {
 	//	STATUS("Shift %i = %e\n", i, gsl_vector_get(shifts, i));
+		if ( isnan(gsl_vector_get(shifts, i)) ) {
+			gsl_vector_set(shifts, i, 0.0);
+		}
 	}
 
 	/* Apply shifts */
-- 
cgit v1.2.3