From 8604f4d110e244c8ae6ad1bc2aaad5169b8af116 Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Fri, 14 Aug 2015 17:45:57 +0200
Subject: Improve random partiality model

---
 libcrystfel/src/geometry.c | 21 ++++++++++++++++++---
 1 file changed, 18 insertions(+), 3 deletions(-)

(limited to 'libcrystfel/src/geometry.c')

diff --git a/libcrystfel/src/geometry.c b/libcrystfel/src/geometry.c
index cdd936e2..c9e30f4f 100644
--- a/libcrystfel/src/geometry.c
+++ b/libcrystfel/src/geometry.c
@@ -173,22 +173,37 @@ static double random_partiality(signed int h, signed int k, signed int l,
 	int i;
 
 	gsl_rng_set(rng, serial);
+	seed = gsl_rng_get(rng);
+	gsl_rng_set(rng, seed);
+
 	for ( i=0; i<abs(h)+1; i++ ) {
 		seed = gsl_rng_get(rng);
 	}
-	if ( h < 0 ) seed = gsl_rng_get(rng);
+	gsl_rng_set(rng, seed);
+	if ( h >= 0 ) {
+		seed = gsl_rng_get(rng);
+	}
+	seed = gsl_rng_get(rng);
 	gsl_rng_set(rng, seed);
 
 	for ( i=0; i<abs(k)+1; i++ ) {
 		seed = gsl_rng_get(rng);
 	}
-	if ( k < 0 ) seed = gsl_rng_get(rng);
+	gsl_rng_set(rng, seed);
+	if ( k >= 0 ) {
+		seed = gsl_rng_get(rng);
+	}
+	seed = gsl_rng_get(rng);
 	gsl_rng_set(rng, seed);
 
 	for ( i=0; i<abs(l)+1; i++ ) {
 		seed = gsl_rng_get(rng);
 	}
-	if ( l < 0 ) seed = gsl_rng_get(rng);
+	gsl_rng_set(rng, seed);
+	if ( l >= 0 ) {
+		seed = gsl_rng_get(rng);
+	}
+	seed = gsl_rng_get(rng);
 	gsl_rng_set(rng, seed);
 
 	p = gsl_rng_uniform(rng);
-- 
cgit v1.2.3