n:=4;
load "affineHecke.txt";
Bas := rec<RFBasis|>;
BasC := rec<RFBasis|>;
BigBas := rec<RFBasis|>;
IndBasTC := rec<RFBasis|>;
IndBasC := rec<RFBasis|>;
IndType := rec<IndElementType|>;
FromTrivialC := rec<RFBasis|>;
FromTrivial := rec<RFBasis|>;
InitInd([We],[],~IndType);

// Bas and BasC encode trivial modules for the Hecke algebra H of type A_3
InitBasis(0,0,~Bas,["C"],IndType: trivial := true);
ReorderBasis(~Bas);
BasC := Bas;

// FromTrivial encodes the module that is Ind_{H_J}^H BasC, where J = { s_2 }.
time FromTrivial := Induce(BasC,0,0,{2},We);
print "is the time for Induce from BasC";
ReorderBasis(~FromTrivial);

// Replaces the basis encoding FromTrivial with the canonical basis
time FromTrivialC := ConvertBasisType("C",We,FromTrivial,BasC);
print "is the time for computing canonical bases FromTrivial";

Mycells:=<>;
// Computes the cells of Res_{H_J} FromTrivialC, where J = {1, 2} and some other data about these cells.
time Mycells := PrintCells(FromTrivialC,We:FullLabel:=false,PCJ:=[1..n-2]);

Mycells[1];
[<Mycells[2][i], Mycells[3][i]>:i in [1..#(Mycells[2])]];

Mycells[2];
/* Should output the following. These are labels for the cells of Res_{H_J} FromTrivialC, where J = {1, 2}.
[
    [
        [ 2, 3, 4 ]
    ],
    [
        [ 2, 3 ],
        [ 4 ]
    ],
    [
        [ 1, 2, 4 ]
    ],
    [
        [ 1, 2 ],
        [ 4 ]
    ],
    [
        [ 1, 4 ],
        [ 2 ]
    ],
    [
        [ 1 ],
        [ 2 ],
        [ 4 ]
    ],
    [
        [ 1, 2, 3 ]
    ],
    [
        [ 1, 3 ],
        [ 2 ]
    ]
]
*/