Add left coset decomposition

1 files changed, 24 insertions, 0 deletions

diff --git a/tests/symmetry_check.c b/tests/symmetry_check.c
index 5a2c91d5..6934a3d9 100644
--- a/tests/symmetry_check.c
+++ b/tests/symmetry_check.c
@@ -65,6 +65,28 @@ static void check_pg_props(const char *pg, int answer, int centro, int *fail)
 }
 
 
+static void check_subgroup(const char *ssource, const char *starget,
+                           int *fail)
+{
+	SymOpList *source;
+	SymOpList *target;
+	SymOpList *twins;
+
+	source = get_pointgroup(ssource);
+	target = get_pointgroup(starget);
+	if ( (source == NULL) || (target == NULL) ) {
+		*fail = 1;
+		return;
+	}
+
+	twins = get_ambiguities(source, target);
+	describe_symmetry(twins);
+
+	free_symoplist(target);
+	free_symoplist(source);
+}
+
+
 int main(int argc, char *argv[])
 {
 	int fail = 0;
@@ -127,5 +149,7 @@ int main(int argc, char *argv[])
 	check_pg_props( "m-3m",  48,  1, &fail);
 	STATUS("\n");
 
+	check_subgroup("2/m", "m", &fail);
+
 	return fail;
 }