diff options
| -rw-r--r-- | src/symmetry.c | 13 |
1 files changed, 12 insertions, 1 deletions
diff --git a/src/symmetry.c b/src/symmetry.c index e304de48..4bc8d115 100644 --- a/src/symmetry.c +++ b/src/symmetry.c @@ -315,7 +315,8 @@ static ReflItemList *coset_decomp(signed int hs, signed int ks, signed int ls, ReflItemList *get_twins(ReflItemList *items, const char *holo, const char *mero) { int i; - ReflItemList *ops = new_items();; + ReflItemList *ops = new_items(); + int expected, actual; /* Run the coset decomposition for every reflection in the "pattern", * and see which gives the highest number of possibilities. This @@ -339,5 +340,15 @@ ReflItemList *get_twins(ReflItemList *items, const char *holo, const char *mero) } + /* Idiot check */ + actual = num_items(ops); + expected = num_general_equivs(holo) / num_general_equivs(mero); + if ( actual != expected ) { + ERROR("I couldn't find the number of twin laws that I expected." + " This is an internal error, and shouldn't happen. " + "Sorry.\n"); + abort(); + } + return ops; } |