From 18e61c8671d72f56f6203fd49740e501cb53b784 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Tue, 26 Oct 2010 14:57:51 +0200 Subject: Return zero if reflection is out of range --- src/diffraction.c | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) (limited to 'src/diffraction.c') diff --git a/src/diffraction.c b/src/diffraction.c index 29e30431..93ac7f8e 100644 --- a/src/diffraction.c +++ b/src/diffraction.c @@ -229,7 +229,10 @@ static double molecule_factor(const double *intensities,const double *phases, h = (signed int)rint(hd); k = (signed int)rint(kd); l = (signed int)rint(ld); - r = lookup_intensity(intensities, h, k, l); + if ( abs(h) > INDMAX ) r = 0.0; + else if ( abs(k) > INDMAX ) r = 0.0; + else if ( abs(l) > INDMAX ) r = 0.0; + else r = lookup_intensity(intensities, h, k, l); break; case GRADIENT_INTERPOLATE : r = interpolate_intensity(intensities, hd, kd, ld); -- cgit v1.2.3