Make generation work (hopefully) and add unit test

1 files changed, 41 insertions, 90 deletions

diff --git a/src/symmetry.c b/src/symmetry.c
index 3b995967..b0e430b5 100644
--- a/src/symmetry.c
+++ b/src/symmetry.c
@@ -235,71 +235,41 @@ static void combine_ops(signed int *h1, signed int *k1, signed int *l1,
 }
 
 
-static void expand_all_ops(signed int *hs, signed int *ks, signed int *ls,
-                           int skip, SymOpList *s, SymOpList *expanded)
+static void combine_and_add_symop(struct sym_op *opi, int oi,
+                                  struct sym_op *opj,
+                                  SymOpList *s)
 {
-	int i, n;
-
-	n = num_ops(s);
-
-	for ( i=0; i<n; i++ ) {
-
-		int oi;
-		struct sym_op *opi = &s->ops[i];
-
-		for ( oi=0; oi<opi->order; oi++ ) {
-
-			int oi_it;
-			signed int *h_ordered, *k_ordered, *l_ordered;
-
-			h_ordered = malloc(3*sizeof(signed int));
-			k_ordered = malloc(3*sizeof(signed int));
-			l_ordered = malloc(3*sizeof(signed int));
-			assert(h_ordered != NULL);
-			assert(k_ordered != NULL);
-			assert(l_ordered != NULL);
-
-			memcpy(h_ordered, hs, 3*sizeof(signed int));
-			memcpy(k_ordered, ks, 3*sizeof(signed int));
-			memcpy(l_ordered, ls, 3*sizeof(signed int));
+	int i;
+	signed int *h, *k, *l;
 
-			/* Apply "opi" this number of times */
-			for ( oi_it=0; oi_it<oi; oi_it++ ) {
+	h = malloc(3*sizeof(signed int));
+	k = malloc(3*sizeof(signed int));
+	l = malloc(3*sizeof(signed int));
+	assert(h != NULL);
+	assert(k != NULL);
+	assert(l != NULL);
 
-				signed int hfs[3], kfs[3], lfs[3];
+	memcpy(h, opj->h, 3*sizeof(signed int));
+	memcpy(k, opj->k, 3*sizeof(signed int));
+	memcpy(l, opj->l, 3*sizeof(signed int));
 
-				combine_ops(h_ordered, k_ordered, l_ordered,
-				            opi->h, opi->k, opi->l,
-				            hfs, kfs, lfs);
+	for ( i=0; i<oi; i++ ) {
 
-				if ( i != skip ) {
+		signed int hfs[3], kfs[3], lfs[3];
 
-					memcpy(h_ordered, hfs,
-					       3*sizeof(signed int));
-					memcpy(k_ordered, kfs,
-					       3*sizeof(signed int));
-					memcpy(l_ordered, lfs,
-					       3*sizeof(signed int));
+		combine_ops(h, k, l, opi->h, opi->k, opi->l, hfs, kfs, lfs);
 
-				}
-			}
+		memcpy(h, hfs, 3*sizeof(signed int));
+		memcpy(k, kfs, 3*sizeof(signed int));
+		memcpy(l, lfs, 3*sizeof(signed int));
 
-			STATUS("i=%i, oi=%i\n", i, oi);
-			STATUS("Creating %3i %3i %3i\n", h_ordered[0],
-			                                 h_ordered[1],
-			                                 h_ordered[2]);
-			STATUS("         %3i %3i %3i\n", k_ordered[0],
-			                                 k_ordered[1],
-			                                 k_ordered[2]);
-			STATUS("         %3i %3i %3i\n", l_ordered[0],
-			                                 l_ordered[1],
-			                                 l_ordered[2]);
-			STATUS("\n");
-			add_symop(expanded, h_ordered, k_ordered, l_ordered, 1);
+	}
 
-		}
+//	STATUS("Creating %3i %3i %3i\n", h[0], h[1], h[2]);
+//	STATUS("         %3i %3i %3i\n", k[0], k[1], k[2]);
+//	STATUS("         %3i %3i %3i\n", l[0], l[1], l[2]);
 
-	}
+	add_symop(s, h, k, l, 1);
 }
 
 
@@ -308,55 +278,37 @@ static void expand_all_ops(signed int *hs, signed int *ks, signed int *ls,
 static SymOpList *expand_ops(SymOpList *s)
 {
 	int n, i;
-	SymOpList *expanded;
+	SymOpList *e;
 
-	n = num_ops(s);
-	expanded = new_symoplist();
-	if ( expanded == NULL ) return NULL;
-	expanded->name = strdup(symmetry_name(s));
-	STATUS("%i ops to expand.\n", n);
+	e = new_symoplist();
+	if ( e == NULL ) return NULL;
+	e->name = strdup(symmetry_name(s));
+
+	add_symop(e, v(1,0,0,0), v(0,1,0,0), v(0,0,0,1), 1);  /* I */
 
+	n = num_ops(s);
 	for ( i=0; i<n; i++ ) {
 
-		int oi;
+		int j, nj;
 		struct sym_op *opi = &s->ops[i];
-		STATUS("Op %i, order=%i\n", i, opi->order);
 
-		for ( oi=0; oi<opi->order; oi++ ) {
+		/* Apply op 'i' to all the current ops in the list */
+		nj = num_ops(e);
+		for ( j=0; j<nj; j++ ) {
 
-			int oi_it;
-			signed int h_ordered[3], k_ordered[3], l_ordered[3];
-
-			h_ordered[0] = 1;  h_ordered[1] = 0;  h_ordered[2] = 0;
-			k_ordered[0] = 0;  k_ordered[1] = 1;  k_ordered[2] = 0;
-			l_ordered[0] = 0;  l_ordered[1] = 0;  l_ordered[2] = 1;
-
-			/* Apply "opi" this number of times */
-			for ( oi_it=0; oi_it<oi; oi_it++ ) {
-
-				signed int hfs[3], kfs[3], lfs[3];
-
-				combine_ops(h_ordered, k_ordered, l_ordered,
-				            opi->h, opi->k, opi->l,
-				            hfs, kfs, lfs);
-
-				memcpy(h_ordered, hfs, 3*sizeof(signed int));
-				memcpy(k_ordered, kfs, 3*sizeof(signed int));
-				memcpy(l_ordered, lfs, 3*sizeof(signed int));
+			int oi;
 
+			for ( oi=0; oi<opi->order-1; oi++ ) {
+				combine_and_add_symop(opi, oi+1, &e->ops[j], e);
 			}
 
-			STATUS("Upper: i=%i, oi=%i\n", i, oi);
-			expand_all_ops(h_ordered, k_ordered, l_ordered, i,
-			               s, expanded);
-
 		}
 
 	}
 
 	free_symoplist(s);
 
-	return expanded;
+	return e;
 }
 
 
@@ -365,7 +317,7 @@ static SymOpList *expand_ops(SymOpList *s)
 static SymOpList *make_1bar()
 {
 	SymOpList *new = new_symoplist();
-	add_symop(new, v(1,0,0,0), v(0,1,0,0), v(0,0,0,1), 2);  /* -I */
+	add_symop(new, v(-1,0,0,0), v(0,-1,0,0), v(0,0,0,-1), 2);  /* -I */
 	new->name = strdup("-1");
 	return expand_ops(new);
 }
@@ -374,7 +326,6 @@ static SymOpList *make_1bar()
 static SymOpList *make_1()
 {
 	SymOpList *new = new_symoplist();
-	add_symop(new, v(1,0,0,0), v(0,1,0,0), v(0,0,0,1), 1);  /* I */
 	new->name = strdup("1");
 	return expand_ops(new);
 }