I have the following problem. I want to define expression which contains elements of several Weyl Character Rings. The problem comes from physics. I have expression for theory with several symmetries, which can be of the same or different groups which involves characters of these groups representations. Let's take a simplest example of such expression
$x = \chi_{R_1}(U)\chi_{R_2}(V)$
where $U\in G_1 $, $V\in G_2 $ are matrices in groups $G_1$ and $G_2$ and $R_{1,2}$ are corresponding representations. Let's say I want to find $x^2$ (or any other power) expanded in terms of irreps. If $G_1 \neq G_2$ I can write the following code that will work
G1 = WeylCharacterRing(['A',2], base_ring=QQ, style="coroots")
G2 = WeylCharacterRing(['A',1], base_ring=G1, style="coroots")
fund1 = G1(G1.fundamental_weights()[1]);
fund2 = G2(G2.fundamental_weights()[1]);
x = fund1*fund2
x**2
which gives me what I want
(A2(0,1)+A2(2,0))*A1(0) + (A2(0,1)+A2(2,0))*A1(2)
Here for simplicity I took $SU(2)$ and $SU(3)$ groups and fundamental representations for both of them. Now if $G_1 = G_2$ (for example $A_2$ root system in both cases) it also works but result is not very readable now since I can not understand which representation belongs to one group or another. For example
G1 = WeylCharacterRing(['A',2], base_ring=QQ, style="coroots")
G2 = WeylCharacterRing(['A',2], base_ring=G1, style="coroots")
fund1 = G1(G1.fundamental_weights()[1]);
fund2 = G2(G2.fundamental_weights()[1]);
ad1 = G1.adjoint_representation()
ad2 = G2.adjoint_representation()
x = fund1*ad2
x**2
gives me back
(A2(0,1)+A2(2,0))*A2(0,0) + (A2(0,1)+A2(2,0))*A2(0,3) + (2*A2(0,1)+2*A2(2,0))*A2(1,1) + (A2(0,1)+A2(2,0))*A2(3,0) + (A2(0,1)+A2(2,0))*A2(2,2)
which is hard to read. Is there a simple way to keep track which of the expressions belongs to one group or another?
