? \p60
   realprecision = 67 significant digits (60 digits displayed)
? v(rr)=vector((2*n)+1,l,G[rr,l])
? t(ss)=vector(2*n,l,K[ss,l])
? f(ii)=sum(j=1,n,F[ii,j]*log(a[j]))+sum(j=1,n,F[ii,n+j]*log(1-a[j]))+F[ii,(2*n)+1]*Pi*I
? g(ii,jj)=(F[ii,jj]/a[jj])-(F[ii,n+jj]/(1-a[jj]))
delta(ii,jj)=if(abs(F[ii,jj])+abs(F[ii,jj+n]),(Eps/(2*n))*(1/((abs(F[ii,jj])/(abs(a[jj]))^2)+(abs(F[ii,n+jj])/(abs(1-a[jj]))^2))),((4*normb)/Eta)+1)
? vecdelta(jj)=vector(n,i,delta(i,jj))
? mindelta(jj)=vecmin(vecdelta(jj))
? vecmindelta(jj)=[mindelta(jj),imag(a[jj]),(norml2(a[jj]))/2,(norml2(1-a[jj]))/2]
? mmindelta(jj)=vecmin(vecmindelta(jj))
? c(i,j)=(abs(F[i,j])/(abs(a[j])-2*normhh)^2)+(abs(F[i,j+n])/(abs(1-a[j])-2*normhh)^2)
template(n,k,a,FF,GG,ttet,mmfld,compare=5)=local(F,G,H,K,r);F=FF;G=GG;H=matrix(k,2*n,i,j,F[i,j]);K=matrix(n,2*n,i,j,G[i,j]);r=1;while(n-matrank(H)&&(n+1-r),if((matrank(concat(F,v(r)))-matrank(F))&&(matrank(concat(H,t(r)))-matrank(H)),F=concat(F,v(r));H=concat(H,t(r)),r=r++));eval(F);eval(H);b=vector(n,i,f(i));normb=sqrt(norml2(b));A=matrix(n,n,i,j,g(i,j));matdet(A);M=conj(A)*mattranspose(A);wapprox=polroots(charpoly(M,x));w=real(wapprox);Eta=sqrt(vecmin(w))/2;Eps=(sqrt(vecmin(w))-Eta)/n;VECDELTA=vector(n,j,mmindelta(j));DELTA=vecmin(VECDELTA);(Eta*DELTA)/4;if(normb<((Eta*DELTA)/4),compare=2);B=matrix(n,1,j,i,b[j]);hhh=-(A)^(-1)*(B);hh=vector(n,j,hhh[j,1]);normhh=sqrt(norml2(hh));2*(normhh)<DELTA;atilde=a+hh;for(j=1,n,if(normhh<imag(atilde[j]),,error("failure at atilde[",j,"]")));for(j=1,n,if(normhh<(1/2)*abs(a[j]),,error("failure at atilde[",j,"]")));for(j=1,n,if(normhh<(1/2)*abs(1-a[j]),,error("failure at atilde[",j,"]")));Mlips=sqrt(sum(j=1,n,sum(i=1,n,c(i,j)^2)));normAinv=1/(2*Eta);1/(2*(normAinv)^2*Mlips);if(normb<=1/(2*(normAinv)^2*Mlips),compare=2);sqrt(norml2(A^(-1)));1/(2*norml2(A^(-1))*Mlips);if(normb<=1/(2*norml2(A^(-1))*Mlips),compare=3);if(compare<4,print("manifold ",ttet,mmfld," is complete hyperbolic"),print("manifold ",ttet,mmfld," can not be proven to be complete hyperbolic"))
? trap()
? counter=1
%1 = 1
? cancel=0
%2 = 0
? nummflds=4
%3 = 4
for(counter=1,2*nummflds,ID=counter;va=Vec(extern("sed -n "counter"p perloutfix"));counter=counter+1;vn=matsize(va)[2];soltype=2;for(imm=1,vn,if(imag(va[imm])<0,soltype=3));for(imm=1,vn,if(imag(va[imm])<10^(-28),soltype=4));FG=extern("sed -n "counter"p perloutfix");numalleq=matsize(FG)[1];vk=numalleq-vn;vF=matrix(vk,2*(vn)+1,ii,jj,FG[ii,jj]);vG=matrix(vn,2*(vn)+1,ii,jj,FG[vk+ii,jj]);if(counter==2,tet="5 ";mfld=168);if(counter==4,tet="6 ";mfld=297);if(counter==6,tet="7 ";mfld=1431);if(counter==8,tet="7 ";mfld=1927);if(soltype<3,template(vn,vk,va,vF,vG,tet,mfld),print("manifold ",tet,mfld," has param with neg or <10^(-28) imag component as approx sol"));)
? [0.7733011742417981161776267305+1.467711508710224270201778287*I,1.102784715200295155851014308+0.6654569511528134767061906116*I,1.102784715200295155851014308+0.6654569511528134767061906116*I,-0.6647417704792485198366937892+0.4011272787787053967477938379*I,0.7733011742417981161776267305+1.467711508710224270201778287*I]~
? [-1,1,-1,0,1,2,0,2,0,0,2;0,0,-1,1,1,0,-1,1,0,0,-1;0,0,0,-2,0,-1,1,1,2,-1,2;0,0,0,0,0,-1,-1,-1,0,-1,-2;1,-1,0,1,0,0,1,-1,0,0,-1;-1,1,1,0,-1,2,0,0,-2,2,2]
manifold 5 168 is complete hyperbolic
? [0.3478103847799310287081835500+1.028852254136692403433083008*I,-0.2948772561507291921394008735+0.8722716254613285720020490241*I,-0.5897545123014583842788017471+1.744543250922657144004098048*I,1.173905192389965514354091775+0.5144261270683462017165415041*I,1.285375281861320441695168847+0.3131612573507276628620679683*I,1.347810384779931028708183550+1.028852254136692403433083008*I]~
? [1,0,0,-1,-1,1,0,0,0,0,0,-1,-1;0,-1,0,-1,0,1,-1,1,-1,1,-1,0,0;1,1,0,-1,0,0,0,0,0,1,-1,0,-1;0,0,-1,1,1,0,0,0,1,0,1,1,2;0,1,1,0,-1,0,0,0,0,-1,1,0,-1;-1,0,0,0,0,-1,1,-1,1,0,0,0,1;0,-1,0,1,0,0,0,0,-1,-1,0,-1,-1]
manifold 6 297 is complete hyperbolic
? [-0.009920107820026848545925882970+1.301649705484044125585953369*I,1.024183670835232400150303549+0.5002629675593989803076184995*I,0.7808191399759772588889597495+1.195715513907902381291948852*I,1.719841218123934286808126798+0.2283954397470835105856209192*I,1.357831020696311948482045809+0.8361115400313675455103251122*I,-0.6290666437234602844422870577+1.186720257102513708028009943*I,0.7035353903287042017265504249+0.6906134794995176471576303239*I]~
? [0,0,-1,-1,0,0,-1,-1,-1,0,0,0,1,0,0;0,0,-2,0,-1,-1,0,-1,-1,2,0,1,1,-1,2;0,-1,0,1,2,0,-1,0,1,-2,0,0,0,1,0;0,0,0,0,0,0,-1,0,-1,0,-1,0,-1,1,-1;2,-1,0,0,0,1,0,0,1,0,0,0,0,0,-1;0,0,0,0,-1,1,1,0,0,0,-1,1,0,0,-1;0,0,1,-1,0,0,0,0,0,0,1,-2,0,-1,-1;-2,2,1,0,0,-1,1,1,0,0,1,0,0,0,2]
manifold 7 1431 is complete hyperbolic
? [0.5623154588478815762567920339+0.7628885866479953546056459654*I,-1.257038198390025394072240446+0.3810313780849329458294774414*I,0.4820538020406238057914206539+0.3218065600430791556256727294*I,1.152615340835564275199349334+0.6046773521493973046435975971*I,0.6990074247490748023324129366+0.8492676345240877909938075136*I,0.2921537766135374190026237150+1.401552862757634524005222736*I,0.6110832937929372542198651992+1.778043289911934249768273000*I]~
? [0,0,1,1,0,-1,-1,0,-1,0,0,0,2,2,2;0,1,0,-1,0,1,0,0,0,0,1,-1,0,0,-1;0,-1,0,0,-1,-2,1,-1,0,0,1,1,2,0,3;0,-1,0,0,0,0,0,0,1,-2,0,-1,-1,-2,-1;0,1,1,0,-1,0,0,0,0,0,-1,1,0,0,-1;0,0,-2,0,1,0,1,0,-1,2,0,0,0,0,0;2,0,1,1,1,1,-2,0,0,0,0,0,0,2,0;-2,0,0,0,0,0,0,1,0,0,-1,0,-1,0,0]
manifold 7 1927 is complete hyperbolic
? quit
Goodbye!