imroved Hecke

Author: grahamknockillaree

HowToUpdateHap.txt

@@ -0,0 +1,39 @@
+Change Test:=testquick.g, date and version in PackageInfo.g
+
+Change ~/Hap/date
+
+Change ~/Hap/version
+
+Run ./updateAll.sh
+
+Fix Undocumented.xml
+
+Run clean in ~/Hap/doc and in ~/Hap/doc/tutorial
+
+git clone https://github.com/gap-packages/hap.git
+
+cp -r pkg/Hap1.30/* hap
+
+rm hap/HowToUpdateHap.txt    THIS SEEMS TO BE IMPORTANT
+
+diff -r pkg/Hap1.30/ hap/
+
+cd hap and then
+
+ git add .
+ git commit -m "message"
+ git push origin master  [with username grahamknockillaree and passwd MumWon]
+
+ git clone --branch gh-pages  https://github.com/gap-packages/hap.git
+ 
+Rename the gh-pages directory from whatever it is called to gh-pages
+
+Populate gh-pages/doc  gh-pages/tutorial  gh-pages/www with most recent files (xml, html files, etc.)
+
+Make sure you have the latest Release Tools and then from within hap type
+
+../ReleaseTools/release-gap-package --token ghp_whateverthetokenis
+
+For some reason it may be necessary to perform the last command as root.
+
+That's it.

lib/ArithmeticGroups/cplGTree.gi

@@ -199,7 +199,6 @@ local t1,t2,
     ##
     Homotopy:=function(n,w)
     local d,path,i,h,k,g,pk,r,t;
-
     if not n=0 then 
         return [];
     fi;

lib/ArithmeticGroups/freeZGRes.gi

@@ -916,7 +916,6 @@ local
     Homotopy:=function(k,w)
     local 
         f,g,q,s,r,t,v,e,h;
-
     ### h([])=[] 
         if w=[] then 
             return [];

lib/Congruence/bianchi.gi

@@ -1145,7 +1145,11 @@ function(OQ,LL)
 local K,L;
 
 L:=QNeighbourhoodOfUnimodularPairs(OQ,LL);
+if OQ!.bianchiInteger=-5 then
 L:=Filtered(L,x->IsQUnimodularPair(OQ,x));
+else
+L:=Filtered(L,x->IsQQUnimodularPair(OQ,x));
+fi;
 K:=UnimodularPairsReduced_NN(OQ,L);
 K:=QQNeighbourhoodOfUnimodularPairs(OQ,K[1]);
 K:=UnimodularPairsReduced_NN(OQ,K);
@@ -1168,7 +1172,7 @@ ABI:=AbsInt(BI);
 ###############################
 
 HAPRECORD:=[
-[ 1, 1 ],[ 2, 1 ],[ 3, 1 ],[ 5, 20 ],[ 6, 24 ],[ 7, 1 ],[ 10, 40 ],[ 11, 1 ],[ 13, 52 ],[ 14, 56 ],[ 15, 15 ],[ 17, 68 ],[ 19, 4 ],[ 21, 84 ],[ 22, 88 ],[ 23, 18 ],[ 26, 104 ],[ 29, 169 ],[ 30, 121 ],[ 31, 16 ],[ 33, 132 ],[ 34, 9 ],[ 35, 35 ],[ 37, 25 ],[ 38, 4 ],[ 39, 39 ],[ 41, 25 ],[ 42, 25 ],[ 43, 9 ],[ 46, 25 ],[ 47, 24 ],[ 51, 51 ],[ 53, 25 ],[ 55, 56 ],[ 57, 49 ],[ 59, 27 ],[ 67, 16 ],[ 71, 48 ],[ 79, 40 ],[ 83, 33 ],[ 87, 87 ],[ 91, 91 ],[ 95, 95 ],[ 103, 56 ],[ 107, 36 ],[ 111, 111 ],[ 115, 115 ],[ 119, 25 ],[ 123, 123 ],[ 127, 64 ],[ 131, 45 ],[ 139, 55 ],[ 143, 9 ],[ 151, 80 ],[ 155, 155 ],[ 163, 53 ],
+[ 1, 1 ],[ 2, 1 ],[ 3, 1 ],[ 5, 20 ],[ 6, 24 ],[ 7, 1 ],[ 10, 40 ],[ 11, 1 ],[ 13, 52 ],[ 14, 56 ],[ 15, 15 ],[ 17, 68 ],[ 19, 4 ],[ 21, 84 ],[ 22, 88 ],[ 23, 18 ],[ 26, 104 ],[ 29, 169 ],[ 30, 121 ],[ 31, 16 ],[ 33, 132 ],[ 34, 9 ],[ 35, 35 ],[ 37, 25 ],[ 38, 4 ],[ 39, 39 ],[ 41, 25 ],[ 42, 25 ],[ 43, 9 ],[ 46, 25 ],[ 47, 24 ],[ 51, 51 ],[ 53, 25 ],[ 55, 56 ],[ 57, 49 ],[ 59, 27 ],[ 67, 16 ],[ 71, 48 ],[ 79, 40 ],[ 83, 33 ],[ 87, 87 ],[ 91, 91 ],[ 95, 95 ],[ 103, 56 ],[ 107, 36 ],[ 111, 111 ],[ 115, 115 ],[ 119, 25 ],[ 123, 123 ],[ 127, 64 ],[ 131, 45 ],[ 139, 55 ],[ 143, 9 ],[ 151, 80 ],[ 155, 155 ],[ 163, 53 ],[167,384],
 [ 179, 155]
 ];;
                  #HAPRECORD contains some precomuted [d,r] for which Swan's

lib/Congruence/bianchiComplex.gi

@@ -6,7 +6,74 @@ end;
 #######################################################
 #######################################################
 
-#RT:=0;
+#######################################################
+#######################################################
+HAP_QuadToCyclotomicMat:=function(A)
+local B,R,r,x;
+B:=[];
+
+for r in A do
+R:=[];
+for x in r do
+Add(R,HAP_QuadToCyclotomic(x));
+od;
+Add(B,R);
+od;
+
+return B;
+end;
+#######################################################
+#######################################################
+
+#######################################################
+#######################################################
+HAP_QuadToCyclotomicGroup:=function(G)
+local gens;
+
+gens:=GeneratorsOfGroup(G);
+gens:=List(gens,HAP_QuadToCyclotomicMat);
+
+return Group(gens);
+end;
+#######################################################
+#######################################################
+
+#######################################################
+#######################################################
+QuadraticToCyclotomicCoefficients:=function(RR)
+local R,BI,n,k,d,Q,OQ,I,G,L,pos;
+
+for n in [1..Length(RR)] do
+for k in [1..RR!.dimension(n)] do
+RR!.boundary(n,k);
+od;
+od;
+
+R:=Objectify(HapResolution,
+                rec(
+                dimension:=RR!.dimension,
+                boundary:=RR!.boundary,
+                homotopy:=RR!.homotopy,
+                elts:=List(RR!.elts,HAP_QuadToCyclotomicMat),
+                group:=fail,
+                properties:= RR!.properties));
+
+d:=One(RR!.elts[1][1][1]);;
+d:=d!.bianchiInteger;
+Q:=QuadraticNumberField(d);;
+OQ:=RingOfIntegers(Q);;
+I:=QuadraticIdeal(OQ,1);;
+G:=HAP_CongruenceSubgroupGamma0(I);;
+G!.tree:=true;
+R!.group:=G;
+
+return R;
+
+end;
+#######################################################
+#######################################################
+
+
 #######################################################
 #######################################################
 HAP_BianchiAction:=function(A,PP)
@@ -28,7 +95,6 @@ b:=A[1][2]; b:=HAP_QuadToCyclotomic(b);
 c:=A[2][1]; c:=HAP_QuadToCyclotomic(c); 
 d:=A[2][2]; d:=HAP_QuadToCyclotomic(d);
 
-#RT:=RT-Runtime();  #MOST OF THE TIME
 cc:=ComplexConjugate(c*z+d);
 
 #zz:=(a*z+b)*ComplexConjugate(c*z+d) +a*ComplexConjugate(c)*t^2;
@@ -44,7 +110,6 @@ zz:=zz/D;
 #tt:=t/((c*z+d)*ComplexConjugate(c*z+d) +(c*t)*ComplexConjugate(c*t));
 tt:=t/D;
 
-#RT:=RT+Runtime();
 
 
 ans:=[];