passagemath-sympow 10.6.39__cp314-cp314t-macosx_13_0_x86_64.whl

sage_wheels/share/sympow/datafiles/m02ES.txt

@@ -0,0 +1,101 @@
+PS 0
+1.0000000000000000000
+0
+0
+0
+PS 1
+0
+0
+0
+0
+PS 2
+-0.5000000000000000000
+0
+0
+0
+PS 3
+0.3333333333333333148
+1.8503717077085941313e-17
+0
+0
+PS 4
+-0.12500000000000000000
+0
+0
+0
+PS 5
+0.03333333333333333287
+4.625929269271485521e-19
+0
+0
+PS 6
+-0.006944444444444444059
+-3.854941057726237774e-19
+0
+0
+PS 7
+0.0011904761904761905836
+-1.0738764375094520475e-19
+0
+0
+PS 8
+-0.00017361111111111111774
+6.625679942966971324e-21
+0
+0
+PS 9
+2.2045855379188714009e-5
+-1.4867146192371776648e-21
+0
+0
+PS 10
+-2.480158730158730157e-6
+-2.1511947866775881608e-24
+0
+0
+PS 11
+2.505210838544171758e-7
+1.1981638892337766303e-23
+0
+0
+PS 12
+-2.2964432686654910205e-8
+1.3280795650245859229e-24
+0
+0
+PS 13
+1.927085260418593704e-9
+4.762595849589671153e-26
+0
+0
+PS 14
+-1.4911969277048642990e-10
+8.623893279380427754e-27
+0
+0
+PS 15
+1.0706029224547742766e-11
+3.004905946192624507e-28
+0
+0
+PS 16
+-7.169215998581077817e-13
+-1.2909695470519215890e-29
+0
+0
+PS 17
+4.498331606952832957e-14
+2.6414148369378292159e-30
+0
+0
+PS 18
+-2.6552651846596585180e-15
+1.9385575368755759074e-32
+0
+0
+PS 19
+1.4797143443923793400e-16
+9.040530305091900585e-34
+0
+0
+END

sage_wheels/share/sympow/datafiles/param_data

@@ -0,0 +1,13 @@
+M02H,285.3116706110000109,-7,11
+M02L,259.3742460099999789,-7,11
+m01E,75.13148009015777973,-8,10
+m02E,82.64462809917355912,-8,10
+A01O,35.04938994813925035,-5,6
+A012,31.863081771035680844,-5,6
+A013,28.966437973668799800,-4,4
+A014,26.333125430607999818,-4,4
+A015,21.762913579014877996,-4,4
+A016,19.784466890013526096,-4,4
+A017,16.350799082655804284,-4,4
+A018,14.864362802414369469,-3,2
+A019,12.284597357367246673,-3,2

sage_wheels/share/sympow/standard1.gp

@@ -0,0 +1,24 @@
+allocatemem(2^28); \\ need 2^28 for 64-bit machines; else 2^27
+L; X;
+SP=mx+2;
+default(seriesprecision,SP+1);
+H=vector(SP);
+{for(i=1,SP,H[i]=vector(N);
+ for(j=1,N,if(i==1,if(j==1,H[1][1]=1,H[i][j]=H[i][j-1]+1.0/j),
+           if(j==1,H[i][1]=H[i-1][1],H[i][j]=H[i][j-1]+H[i-1][j]/j*1.0))));}
+Hf(n,k)=if(k==0,0,H[n][k])
+ZETA(n)=if(n==1,Euler,zeta(n))
+J(k,v)=if(k<0,(v-k-1)*J(k+1,v),\
+                1.0*(-1)^k/k!/v*\
+                exp(sum(n=1,SP,(-1)^n*(ZETA(n)/n)*v^n))*\
+                (1+sum(n=1,SP,Hf(n,k)*v^n)))
+sinv(k,v)=if(k==0,1/v,-1/k-sum(l=1,SP,v^l/k^(l+1)))
+two1ms(k,v)=2^(1+k)*sum(l=0,SP,log(1/2)^l/l!*v^l)
+
+{DERIV(M)=
+ sum(i=0,poldegree(M,L),
+        (deriv(polcoeff(M,i,L),X)+(i+1)*polcoeff(M,i+1,L)/X)*L^i)}
+        
+{ev(T,v,C)=U=0;
+ for(i=0,mx,U+=(truncate(polcoeff(T,i,L)+O(X^C))*L^i));
+ subst(subst(U,L,log(v)),X,v);}

sage_wheels/share/sympow/standard2.gp

@@ -0,0 +1,24 @@
+
+P=sum(k=0,N,sum(l=1,mx+1,-(-1)^l/(l-1)!*polcoeff(F(k),-l,X)*L^(l-1))*X^k);
+K=40; Ds=vector(K); Ds[1]=DERIV(P); for(i=2,K,Ds[i]=DERIV(Ds[i-1]));
+
+ieee(x)=(round(x<<53)*1.0)>>53
+{IEEE(x)=if(x==0,return(0));if(length(x)<2,return(0));
+         y=ceil(log(abs(x))/log(2));ieee(x/2^y)*2^y;}
+{QD(x)=local(A);A=[IEEE(x),0,0,0];A[2]=IEEE(x-A[1]);A[3]=IEEE(x-A[1]-A[2]);
+ A[4]=IEEE(x-A[1]-A[2]-A[3]); A=precision(A,18);
+ print(A[1]); print(A[2]); print(A[3]); print(A[4]);}
+
+ 
+{doit(x,C)=print("AT ",x); QD(ev(P,x,C));
+ for(d=1,35,print("DERIV ",d); QD(ev(Ds[d]/d!,x,C)));}
+ch(x,C)=ev(P,x,C)-ev(P,x,C+20)
+
+setup(a,b)=for(i=0,31,doit(2^a+2^a*i/32,b));
+print("About to find TOO_BIG"); \\ ends up in /dev/null
+MM=100; while(log(10^(-70)+abs(ev(P,MM,N)))>-148,MM*=1.1);
+write1(PARAMDATAFILE,STR,",",precision(IEEE(MM),18));
+print("Now working backwards..."); \\ ends up in /dev/null
+l=-10; while(log(10^(-70)+abs(ch(2^l,20)))<-148,l+=1); l-=1; s=l;
+write1(PARAMDATAFILE,",",l-5); n=20;
+print("Starting to write mesh files"); \\ ends up in /dev/null

sage_wheels/share/sympow/standard3.gp

@@ -0,0 +1,7 @@
+while(2^l<MM,while(log(10^(-70)+abs(ch(2*2^l,n)))>-148,n+=5);setup(l,n);l+=1);
+write(PARAMDATAFILE,",",l-s);
+coeffs(j)=for(i=0,19,print("PS ",i);QD(polcoeff(polcoeff(P,j,L),i,X)));
+coeffO(j)=\
+for(i=0,9,print("PSL ",2*i+1);QD(polcoeff(polcoeff(P,j,L),2*i+1,X)));
+coeffE(j)=\
+for(i=0,9,print("PSL ",2*i);QD(polcoeff(polcoeff(P,j,L),2*i,X)));