#############################################################################
##
#W  nocentre.gi                 GrpConst                         Bettina Eick
#W                                                         Hans Ulrich Besche
##

#############################################################################
##
#F PermRepDP( D ) . . . . . . . . . . . . . permutation rep of direct product
##
BindGlobal( "PermRepDP", function( D )
    local A, B, operA, operB, G, emb1, emb2;

    A := Image( Projection( D, 1 ) );
    B := Image( Projection( D, 2 ) );
    operA := IsomorphismPermGroup( A );
    operB := IsomorphismPermGroup( B );
    G := DirectProduct( Image( operA ), Image( operB ) );
    emb1 := Embedding( G, 1 );
    emb2 := Embedding( G, 2 );
    return [operA * emb1, operB * emb2];

end );

#############################################################################
##
#F PermOper( oper, elms ) . . . . . . . . perm rep of direct product elements
##
BindGlobal( "PermOper", function( oper, gens )
    local new, g, h1, h2, h;
    new := [];
    for g in gens do
        h1 := Image( oper[1], g[1] );
        h2 := Image( oper[2], g[2] );
        h := h1 * h2;
        Add( new, h );
    od;
    return Group( new, () );
end );


#############################################################################
##
#F CosetReps . . . . . . . . . . . . . double coset reps for induced aut grps
##
BindGlobal( "CosetReps", function( C, hom, NU, NL )
    local U, H, CU, gens, oper, g, imgs, aut, CL, reps;

    if Size( C ) = 1 then return [Identity(C)]; fi;
    U := Image( hom );
    H := Source( hom );

    # automorphisms of U induced by NU
    CU := [];
    gens := GeneratorsOfGroup( U );
    oper := GeneratorsOfGroup( NU );
    for g in oper do
        imgs := List( gens, x -> x^g );
        if imgs <> gens then
            aut := GroupHomomorphismByImagesNC(U, U, gens, imgs);
            SetIsBijective( aut, true );
            Add( CU, aut );
        fi;
    od;
    CU := SubgroupNC( C, CU );

    # automorphisms of U induced by NL
    CL := [];
    gens := List(GeneratorsOfGroup(H), x -> Image(hom,x));
    oper := GeneratorsOfGroup( NL );
    for g in oper do
        imgs := List( GeneratorsOfGroup( H ), x -> Image( g, x ) );
        imgs := List( imgs, x -> Image( hom, x ) );
        aut := GroupHomomorphismByImagesNC(U, U, gens, imgs);
        SetIsBijective( aut, true );
        Add( CL, aut );
    od;
    CL := SubgroupNC( C, CL );

    # correct: C, CL, CU / bug: C, CU, CL
    reps := DoubleCosets( C, CL, CU );
    #reps := DoubleCosets( C, CU, CL );
    return List( reps, Representative );
end );
 
#############################################################################
##
## N is a center-free perfect group and H is soluble.
## classify extensions of N by H up to isomorphism
##
BindGlobal( "ExtensionsByGroupNoCentre", function( N, H )
    local A, I, hom, O, clU, B, clL, f, D, gensN, oper, pairs, U, L, nat, 
          F, iso, res, NU, NL, C, reps, r, new, gens, G, pair, g, h;

    # the automorphism group of N
    Info( InfoGrpCon, 2, " compute Aut N ");
    A := AutomorphismGroup(N);
    I := InnerAutomorphismsAutomorphismGroup(A);
    hom := NaturalHomomorphismByNormalSubgroupNC( A, I );
    O := Image(hom);

    # possible projections in O
    Info( InfoGrpCon, 2, " compute subgroups of Out N  ");
    clU := ConjugacyClassesSubgroups( O );
    clU := List( clU, Representative );

    # the automorphism group of H
    Info( InfoGrpCon, 2, " compute Aut H ");
    B := AutomorphismGroup( H );
    IsomorphismPermGroup(B);

    # possible kernels in H
    Info( InfoGrpCon, 2, " compute possible centralizers in H ");
    clL := NormalSubgroups( H );
    clL := Filtered( clL, x -> Index(H, x) <= Size(O) );
    f := function( pt, aut ) return Image( aut, pt ); end;
    clL := Orbits( B, clL, f );
    clL := List( clL, x -> x[1] );

    # compute direct product
    Info( InfoGrpCon, 2, " compute direct product and projections ");
    D := DirectProduct( A, H );
    gensN := List( GeneratorsOfGroup( I ), x -> DirectProductElement( [x, Identity(H)] ) );

    # compute perm reps
    Info( InfoGrpCon, 2, " compute perm rep of D ");
    oper := PermRepDP( D );

    # compute pairs
    Info( InfoGrpCon, 2, " computing pairs ");
    pairs := [];
    for U in clU do
        for L in clL do
            nat := NaturalHomomorphismByNormalSubgroupNC( H, L );
            F := Image( nat );
            if Size(U) = Size( F ) then
                if IdGroup( U ) = IdGroup( F ) then
                    iso := IsomorphismGroups( F, U );
                    iso := nat * iso;
                    SetIsSurjective( iso, true );
                    SetKernelOfMultiplicativeGeneralMapping( iso, L );
                    Add( pairs, iso );
                fi;
            fi;
        od;
    od;

    # classify subdirect products
    Info( InfoGrpCon, 2, " start to loop over ", Length(pairs), " pairs ");
    res := [];
    for pair in pairs do
        Info( InfoGrpCon, 3, "  start factors of id ", IdGroup(U),"");
        U := Image( pair );
        L := Kernel( pair );

        # operating groups
        Info( InfoGrpCon, 3, "  computing normalizer ");
        NU := Normalizer( O, U );
        Info( InfoGrpCon, 3, "  computing stabilizer ");
        NL := Stabilizer( B, L, f );
        Info( InfoGrpCon, 3, "  computing automorphism group ");
        C := AutomorphismGroup( U );

        # coset reps
        Info( InfoGrpCon, 3, "  computing coset reps ");
        reps := CosetReps( C, pair, NU, NL );

        # for each double coset rep create a group
        Info( InfoGrpCon, 3, "  loop over ", Length(reps), " coset reps ");
        for r in reps do
            new := pair * r;

            Info( InfoGrpCon, 4, 
                  "   compute generators of subdirect product "); 
            gens := [];
            for g in GeneratorsOfGroup( H ) do
                h := Image( new, g );
                h := PreImagesRepresentative( hom, h );
                h := DirectProductElement( [h, g ] );
                Add( gens, h );
            od;
            Append( gens, gensN );

            Info( InfoGrpCon, 4, "   compute perm rep ");
            Add( res, PermOper( oper, gens ) );
        od;
    od;
    
    return res;
end );

#############################################################################
##
#F UpwardsExtensionsNoCentre( N, stepsize )
##
BindGlobal( "UpwardsExtensionsNoCentre", function( N, stepsize )
    local all, res, G, tmp;
    if stepsize = 1 then
        return [N];
    fi;
    all := AllSmallGroups( stepsize, IsSolvableGroup, true );
    res := [];
    for G in all do
        IsPGroup( G );
        Info( InfoGrpCon, 1, "start extending by group with id ",IdGroup(G) );
        tmp := ExtensionsByGroupNoCentre( N, G );
        Append( res, tmp );
    od;
    return res;
end );