From 8683aef3880d17c64a10598ff9c2edd1532280cb Mon Sep 17 00:00:00 2001
From: Thomas White <taw@physics.org>
Date: Sat, 23 Nov 2013 08:55:58 +0100
Subject: get_hkl: Fix symmetric expansion

---
 src/get_hkl.c | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

diff --git a/src/get_hkl.c b/src/get_hkl.c
index b4025251..029f8d21 100644
--- a/src/get_hkl.c
+++ b/src/get_hkl.c
@@ -284,7 +284,7 @@ static RefList *expand_reflections(RefList *in, const SymOpList *initial,
 		for ( j=0; j<n; j++ ) {
 
 			signed int he, ke, le;
-			Reflection *new;
+			Reflection *copy;
 			int have_phase;
 			double ph;
 
@@ -294,14 +294,16 @@ static RefList *expand_reflections(RefList *in, const SymOpList *initial,
 			/* Put it into the asymmetric unit for the target */
 			get_asymm(target, he, ke, le, &he, &ke, &le);
 
+			if ( find_refl(out, he, ke, le) != NULL ) continue;
+
 			/* Make sure the intensity is in the right place */
-			new = add_refl(out, he, ke, le);
-			copy_data(new, refl);
+			copy = add_refl(out, he, ke, le);
+			copy_data(copy, refl);
 
 			/* FIXME: Make phase negative if the reflection is
 			 * separated from the original via an inversion */
 			ph = get_phase(refl, &have_phase);
-			if ( have_phase ) set_phase(new, -ph);
+			if ( have_phase ) set_phase(copy, -ph);
 
 		}
 
-- 
cgit v1.2.3