#############################################################################
##  TestDihedral
#############################################################################
##  Description:
##
##  For an integer `N`,
##  this tests all normal subgroups of the dihedral group D(2n)
##
##  The normal subgroups of D(2n) are classified as follows:
##  - 1. D(2n)
##  - 2. C(k) for k dividing n
##  - 3. if n is even, two times D(n)
#############################################################################

TestDihedral := function(N, UseLIS...)
    local G, L, Current, n, k;

    if Length(UseLIS) > 1 then
        Error("Unknown number of arguments!");
    elif Length(UseLIS) = 1 then
        UseLIS := UseLIS[1];
    else
        UseLIS := false;
    fi;

    if not IsEvenInt(N) then
        Error("Dihedral group has even size!");
    fi;

    n := N / 2;
    G := DihedralGroup(N);
    L := List(LowIndexNormalSubgroupsSearch(G, N, rec(UseLIS := UseLIS)));
    Current := 2;
    if IsEvenInt(n) then
        for k in [1 .. 2] do
            if Order(L[Current]!.Grp) <> n then
                Error("LINS did not find D(", n, ")!");
            fi;
            Current := Current + 1;
        od;
    fi;

    for k in Reversed(DivisorsInt(n)) do
        if Order(L[Current]!.Grp) <> k then
            Error("LINS did not find C(", k, ")!");
        fi;
        Current := Current + 1;
    od;

    if Current <= Length(L) then
        Error("LINS found too many subgroups!");
    fi;

    return true;
end;