From f61a2210c6afd5ff222e271c450fdc08cf2b4475 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Wed, 22 Feb 2017 15:08:24 +0100 Subject: Horrible restraint to keep algorithm from drifting too far --- libcrystfel/src/geometry.c | 2 ++ 1 file changed, 2 insertions(+) (limited to 'libcrystfel/src/geometry.c') diff --git a/libcrystfel/src/geometry.c b/libcrystfel/src/geometry.c index 43cf25f2..7f5e8635 100644 --- a/libcrystfel/src/geometry.c +++ b/libcrystfel/src/geometry.c @@ -523,6 +523,8 @@ static double do_integral(double q2, double zl, double R, const double N = 1.5; /* Pointiness of spectrum */ FILE *fh = NULL; + assert(R*R < q2); + /* Range over which P is different from zero */ k0 = (R*R - q2)/(2.0*(zl+R)); k1 = (R*R - q2)/(2.0*(zl-R)); -- cgit v1.2.3