@@ -1133,10 +1133,16 @@ N:=5;
11331133BI:= OQ!. bianchiInteger;
11341134ABI:= AbsInt(BI);
11351135# ##############################
1136- HAPRECORD:= [ [ 1 , 1 ] , [ 2 , 1 ] , [ 3 , 1 ] , [ 5 , 20 ] , [ 6 , 24 ] , [ 7 , 1 ] , [ 10 , 40 ] , [ 11 , 1 ] , [ 13 , 52 ] , [ 14 , 56 ] , [ 15 , 16 ] , [ 17 , 68 ] , [ 19 , 4 ] , [ 21 , 84 ] , [ 22 , 88 ] , [ 23 , 12 ] , [ 26 , 104 ] , [ 29 , 169 ] , [ 30 , 121 ] , [ 31 , 16 ] , [ 33 ,132 ] , [ 34 , 0 ] , [ 35 , 35 ] , [ 37 , 0 ] , [ 38 , 0 ] , [ 39 , 39 ] , [ 41 , 0 ] , [ 42 , 0 ] , [ 43 , 9 ] , [ 46 , 0 ] , [ 47 , 24 ] , [ 51 , 51 ] , [ 53 , 0 ] , [ 55 , 56 ] , [ 57 , 0 ] , [ 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 , 0 ] , [ 123 , 123 ] , [ 127 , 64 ] , [ 131 , 45 ] , [ 139 , 55 ] , [ 143 , 0 ] , [ 151 , 80 ] , [ 163 , 53 ] ] ;
11371136
1138- # HAPRECORD:=[]; #Uncomment this line to calculate
1139- # HAPRECORD entries from scratch.
1137+ HAPRECORD:= [
1138+ [ 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 , 12 ] ,[ 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 ] ,
1139+ [ 179 , 155 ]
1140+ ] ;;
1141+ # HAPRECORD contains some precomuted [d,r] for which Swan's
1142+ # criterion is satisfied for the given norm r.
1143+ # HAPRECORD:=[]; #Uncomment this line to calculate HAPRECORD entries from
1144+ # scratch. The value [179,155] was only partialy checked (up
1145+ # to norm r=5000).
11401146
11411147if Length(arg)= 2 then pos:= arg[ 2 ] ;
11421148else
@@ -1158,9 +1164,9 @@ while bool do
11581164N:= N+ 1 ;
11591165# if NRMS[N]>ABI+30 then NRMS[N]:=0; break; fi; #REMEBER TO DELETE THIS
11601166#
1161- # if Length(arg)>1 then
1167+ if Length(arg)> 1 then
11621168Print(" Adding hemispheres of squared radius " ,1 / NRMS[ N] ," \n " );
1163- # fi;
1169+ fi ;
11641170
11651171L:= UnimodularPairs(OQ,NRMS[ N] ,true ,L);
11661172# L:=QQNeighbourhoodOfUnimodularPairs(OQ,L);
@@ -1183,15 +1189,9 @@ L:=K[1]; K:=K[2];
11831189if pos= infinity then
11841190A:= BianchiPolyhedron(OQ!. bianchiInteger,L);
11851191A:= SwanBianchiCriterion(A);
1186- if Length(A)= 0 then bool:= false ; return [ L,K] ; fi ;
1187- # if A[1]=false then Print("YES\n"); return CoverOfUnimodularPairs(OQ,N+1,true);
1188- # else
1189- # A:=List(A,x->HAPNorm(OQ,x[1]));
1190- # A:=Minimum(A);
1191- # return CoverOfUnimodularPairs(OQ,A);
1192- # fi;
1193- N:= N+ 1 ;
1194- return CoverOfUnimodularPairs(OQ,N,true );
1192+ if not Length(A)= 0 then
1193+ return CoverOfUnimodularPairs(OQ,N+ 1 ,true );
1194+ fi ;
11951195fi ;
11961196# #################################################
11971197
0 commit comments