Add more wordrep tests

tst/testinstall/wordrep.tst

@@ -1,97 +1,208 @@
+#
+gap> START_TEST("wordrep.tst");
+
+#
 gap> f := FreeGroup(IsSyllableWordsFamily,4);;
 gap> fam := FamilyObj(f.1);
 <Family: "FreeGroupElementsFamily">
 
 #
-gap> w8:=ObjByExtRep(fam,[1,1]); Is8BitsAssocWord(w8);
-f1
+gap> w8:=ObjByExtRep(fam,[1,1,2,-10,1,4]); Is8BitsAssocWord(w8);
+f1*f2^-10*f1^4
 true
 gap> ExtRepOfObj(w8);
-[ 1, 1 ]
-gap> w16:=ObjByExtRep(fam,[1,2^10]); Is16BitsAssocWord(w16);
-f1^1024
+[ 1, 1, 2, -10, 1, 4 ]
+gap> w16:=ObjByExtRep(fam,[1,2^10,2,-10,1,4]); Is16BitsAssocWord(w16);
+f1^1024*f2^-10*f1^4
 true
 gap> ExtRepOfObj(w16);
-[ 1, 1024 ]
-gap> w32:=ObjByExtRep(fam,[1,2^20]); Is32BitsAssocWord(w32);
-f1^1048576
+[ 1, 1024, 2, -10, 1, 4 ]
+gap> w32:=ObjByExtRep(fam,[1,2^20,2,-10,1,4]); Is32BitsAssocWord(w32);
+f1^1048576*f2^-10*f1^4
 true
 gap> ExtRepOfObj(w32);
-[ 1, 1048576 ]
-gap> winf:=ObjByExtRep(fam,[1,2^40]);; IsInfBitsAssocWord(winf);
+[ 1, 1048576, 2, -10, 1, 4 ]
+gap> winf:=ObjByExtRep(fam,[1,2^40,2,-10,1,4]);; IsInfBitsAssocWord(winf);
 true
 gap> ExtRepOfObj(winf);
-[ 1, 1099511627776 ]
+[ 1, 1099511627776, 2, -10, 1, 4 ]
+
+#
+# ExponentSums
+#
+
+#
+gap> ExponentSums(w8);
+[ 5, -10, 0, 0 ]
+gap> ExponentSums(w8, 3, 4);
+[ 0, 0 ]
+gap> ExponentSums(w8, 4, 1);
+[ ]
+gap> ExponentSums(w8, 0, 1);
+Error, <start> must be a positive integer
+gap> ExponentSums(w8, 1, 0);
+Error, <end> must be a positive integer
+
+#
+gap> ExponentSums(w16);
+[ 1028, -10, 0, 0 ]
+gap> ExponentSums(w16, 3, 4);
+[ 0, 0 ]
+gap> ExponentSums(w16, 4, 1);
+[ ]
+gap> ExponentSums(w16, 0, 1);
+Error, <start> must be a positive integer
+gap> ExponentSums(w16, 1, 0);
+Error, <end> must be a positive integer
+
+#
+gap> ExponentSums(w32);
+[ 1048580, -10, 0, 0 ]
+gap> ExponentSums(w32, 3, 4);
+[ 0, 0 ]
+gap> ExponentSums(w32, 4, 1);
+[ ]
+gap> ExponentSums(w32, 0, 1);
+Error, <start> must be a positive integer
+gap> ExponentSums(w32, 1, 0);
+Error, <end> must be a positive integer
+
+#
+gap> ExponentSums(winf);
+[ 1099511627780, -10, 0, 0 ]
+gap> ExponentSums(winf, 3, 4);
+[ 0, 0 ]
+gap> ExponentSums(winf, 4, 1);
+[ ]
+gap> ExponentSums(winf, 0, 1);
+Error, <from> must be a positive integer
+gap> ExponentSums(winf, 1, 0);
+Error, <to> must be a positive integer
+
+#
+# syllables
+#
+
+#
+gap> words := [w8, w16, w32, winf];;
+gap> ForAll(words, w -> GeneratorSyllable(w, 2) = 2);
+true
+gap> ForAll(words, w -> ExponentSyllable(w, 2) = -10);
+true
+
+#
+gap> 8Bits_ExponentSyllable(w8, 0);
+Error, <i> must be an integer between 1 and 3
+gap> 16Bits_ExponentSyllable(w16, 0);
+Error, <i> must be an integer between 1 and 3
+gap> 32Bits_ExponentSyllable(w32, 0);
+Error, <i> must be an integer between 1 and 3
+
+#
+gap> 8Bits_GeneratorSyllable(w8, 0);
+Error, <i> must be an integer between 1 and 3
+gap> 16Bits_GeneratorSyllable(w16, 0);
+Error, <i> must be an integer between 1 and 3
+gap> 32Bits_GeneratorSyllable(w32, 0);
+Error, <i> must be an integer between 1 and 3
 
 #
 # test powering (esp. 8Bits_Power), with various shapes
 #
 
 # special case: w * gi^n * w^-1
-gap> w:=ObjByExtRep(fam,[1,1, 2,1, 1,-1]); w^0; w^1; w^-1; w^3; w^100;
+gap> u8:=ObjByExtRep(fam,[1,1, 2,1, 1,-1]); u8^0; u8^1; u8^-1; u8^3; u8^-3; u8^100;
 f1*f2*f1^-1
 <identity ...>
 f1*f2*f1^-1
 f1*f2^-1*f1^-1
 f1*f2^3*f1^-1
+f1*f2^-3*f1^-1
 f1*f2^100*f1^-1
-gap> w:=ObjByExtRep(fam,[1,1, 2,2^10, 1,-1]); w^0; w^1; w^-1; w^3; w^100;
+gap> u16:=ObjByExtRep(fam,[1,1, 2,2^10, 1,-1]); u16^0; u16^1; u16^-1; u16^3; u16^-3; u16^100;
 f1*f2^1024*f1^-1
 <identity ...>
 f1*f2^1024*f1^-1
 f1*f2^-1024*f1^-1
 f1*f2^3072*f1^-1
+f1*f2^-3072*f1^-1
 f1*f2^102400*f1^-1
-gap> w:=ObjByExtRep(fam,[1,1, 2,2^10, 1,-1]); w^0; w^1; w^-1; w^3; w^100;
-f1*f2^1024*f1^-1
+gap> u32:=ObjByExtRep(fam,[1,1, 2,2^20, 1,-1]); u32^0; u32^1; u32^-1; u32^3; u32^-3; u32^100;;
+f1*f2^1048576*f1^-1
 <identity ...>
-f1*f2^1024*f1^-1
-f1*f2^-1024*f1^-1
-f1*f2^3072*f1^-1
-f1*f2^102400*f1^-1
-gap> w:=ObjByExtRep(fam,[1,1, 2,2^40, 1,-1]);; w^0; w^1;; w^-1;; w^3;; w^100;;
+f1*f2^1048576*f1^-1
+f1*f2^-1048576*f1^-1
+f1*f2^3145728*f1^-1
+f1*f2^-3145728*f1^-1
+gap> uinf:=ObjByExtRep(fam,[1,1, 2,2^40, 1,-1]);; uinf^0; uinf^1;; uinf^-1;; uinf^3;; uinf^-3;; uinf^100;;
 <identity ...>
 
 # special case: w * gj^x * t * gj^y * w^-1, x != -y
-gap> w:=ObjByExtRep(fam,[1,1, 2,1, 3,1, 4,-1, 2,2, 1,-1]); w^0; w^1; w^-1; w^3; w^100;
+gap> v8:=ObjByExtRep(fam,[1,1, 2,1, 3,1, 4,-1, 2,2, 1,-1]); v8^0; v8^1; v8^-1; v8^3; v8^-3; v8^100;
 f1*f2*f3*f4^-1*f2^2*f1^-1
 <identity ...>
 f1*f2*f3*f4^-1*f2^2*f1^-1
 f1*f2^-2*f4*f3^-1*f2^-1*f1^-1
 f1*(f2*f3*f4^-1*f2^2)^3*f1^-1
+f1*(f2^-2*f4*f3^-1*f2^-1)^3*f1^-1
 f1*(f2*f3*f4^-1*f2^2)^100*f1^-1
-gap> w:=ObjByExtRep(fam,[1,1, 2,1, 3,2^10, 4,-1, 2,2, 1,-1]); w^0; w^1; w^-1; w^3; w^100;
+gap> v16:=ObjByExtRep(fam,[1,1, 2,1, 3,2^10, 4,-1, 2,2, 1,-1]); v16^0; v16^1; v16^-1; v16^3; v16^-3; v16^100;;
 f1*f2*f3^1024*f4^-1*f2^2*f1^-1
 <identity ...>
 f1*f2*f3^1024*f4^-1*f2^2*f1^-1
 f1*f2^-2*f4*f3^-1024*f2^-1*f1^-1
 f1*(f2*f3^1024*f4^-1*f2^2)^3*f1^-1
-f1*(f2*f3^1024*f4^-1*f2^2)^100*f1^-1
-gap> w:=ObjByExtRep(fam,[1,1, 2,1, 3,2^20, 4,-1, 2,2, 1,-1]); w^0; w^1; w^-1; w^3;
+f1*(f2^-2*f4*f3^-1024*f2^-1)^3*f1^-1
+gap> v32:=ObjByExtRep(fam,[1,1, 2,1, 3,2^20, 4,-1, 2,2, 1,-1]); v32^0; v32^1; v32^-1; v32^3; v32^-3; v32^100;;
 f1*f2*f3^1048576*f4^-1*f2^2*f1^-1
 <identity ...>
 f1*f2*f3^1048576*f4^-1*f2^2*f1^-1
 f1*f2^-2*f4*f3^-1048576*f2^-1*f1^-1
 f1*(f2*f3^1048576*f4^-1*f2^2)^3*f1^-1
+f1*(f2^-2*f4*f3^-1048576*f2^-1)^3*f1^-1
 
 # general case: w * t * w^-1
-gap> w:=ObjByExtRep(fam,[1,1, 2,1, 3,1, 4,-1]); w^0; w^1; w^-1; w^3; w^100;
+gap> x8:=ObjByExtRep(fam,[1,1, 2,1, 3,1, 4,-1]); x8^0; x8^1; x8^-1; x8^3; x8^-3; x8^100;
 f1*f2*f3*f4^-1
 <identity ...>
 f1*f2*f3*f4^-1
 f4*f3^-1*f2^-1*f1^-1
 (f1*f2*f3*f4^-1)^3
+(f4*f3^-1*f2^-1*f1^-1)^3
 (f1*f2*f3*f4^-1)^100
-gap> w:=ObjByExtRep(fam,[1,1, 2,1, 3,2^10, 4,-1]); w^0; w^1; w^-1; w^3; w^100;
+gap> x16:=ObjByExtRep(fam,[1,1, 2,1, 3,2^10, 4,-1]); x16^0; x16^1; x16^-1; x16^3; x16^-3; x16^100;
 f1*f2*f3^1024*f4^-1
 <identity ...>
 f1*f2*f3^1024*f4^-1
 f4*f3^-1024*f2^-1*f1^-1
 (f1*f2*f3^1024*f4^-1)^3
+(f4*f3^-1024*f2^-1*f1^-1)^3
 (f1*f2*f3^1024*f4^-1)^100
-gap> w:=ObjByExtRep(fam,[1,1, 2,1, 3,2^20, 4,-1]); w^0; w^1; w^-1; w^3;
+gap> x32:=ObjByExtRep(fam,[1,1, 2,1, 3,2^20, 4,-1]); x32^0; x32^1; x32^-1; x32^3; x32^-3;
 f1*f2*f3^1048576*f4^-1
 <identity ...>
 f1*f2*f3^1048576*f4^-1
 f4*f3^-1048576*f2^-1*f1^-1
 (f1*f2*f3^1048576*f4^-1)^3
+(f4*f3^-1048576*f2^-1*f1^-1)^3
+
+#
+# quotients
+#
+gap> words8 := [u8,v8,w8,x8];; ForAll(words8, Is8BitsAssocWord);
+true
+gap> words16 := [u16,v16,w16,x16];; ForAll(words16, Is16BitsAssocWord);
+true
+gap> words32 := [u32,v32,w32,x32];; ForAll(words32, Is32BitsAssocWord);
+true
+
+#
+gap> SetX(words8, words8, {a,b} -> (a/b) = 8Bits_Quotient(a,b));
+[ true ]
+gap> SetX(words16, words16, {a,b} -> (a/b) = 16Bits_Quotient(a,b));
+[ true ]
+gap> SetX(words32, words32, {a,b} -> (a/b) = 32Bits_Quotient(a,b));
+[ true ]
+
+#
+gap> STOP_TEST("wordrep.tst", 1);