From 65305f9f5064f0537f365208b6b750e4e6d7a9cf Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Thu, 9 Jun 2011 18:57:49 +0200
Subject: Don't give NaN for shift

---
 src/hrs-scaling.c | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

(limited to 'src/hrs-scaling.c')

diff --git a/src/hrs-scaling.c b/src/hrs-scaling.c
index cbd54ee5..1d11e402 100644
--- a/src/hrs-scaling.c
+++ b/src/hrs-scaling.c
@@ -325,10 +325,18 @@ static double iterate_scale(struct image *images, int n,
 	} else {
 		shifts = gsl_vector_alloc(n);
 		for ( frame=0; frame<n; frame++ ) {
-			double num, den;
+
+			double num, den, sh;
+
 			num = gsl_vector_get(v, frame);
 			den = gsl_matrix_get(M, frame, frame);
-			gsl_vector_set(shifts, frame, num/den);
+			sh = num/den;
+
+			if ( isnan(sh) ) {
+				gsl_vector_set(shifts, frame, 0.0);
+			} else {
+				gsl_vector_set(shifts, frame, sh);
+			}
 		}
 	}
 
-- 
cgit v1.2.3