| version 1.30, 2004/01/14 09:29:39 |
version 1.31, 2004/01/15 07:20:31 |
| Line 199 def fixedpoint(A,FLAG){ |
|
| Line 199 def fixedpoint(A,FLAG){ |
|
| return fixpointmain(F,Vars)$ |
return fixpointmain(F,Vars)$ |
| } |
} |
| |
|
| def nonzerovec(A){ |
|
| |
|
| for(I=0;I<size(A)[0];I++) |
|
| if(A[I]!=0) |
|
| return 1$ |
|
| |
|
| return 0$ |
|
| } |
|
| |
|
| def junban(A,B){ |
def junban(A,B){ |
| return (A<B ? 1:(A>B ? -1:0))$ |
return (A<B ? 1:(A>B ? -1:0))$ |
| } |
} |
| |
|
| def worder(A,B){ |
|
| return (A[0]<B[0] ? 1:(A[0]>B[0] ? -1:0))$ |
|
| } |
|
| |
|
| def bsort(A){ |
def bsort(A){ |
| |
|
| K=size(A)[0]-1$ |
K=size(A)[0]-1$ |
| Line 263 def nonposdegchk(Res){ |
|
| Line 250 def nonposdegchk(Res){ |
|
| |
|
| def getgcd(A,B){ |
def getgcd(A,B){ |
| |
|
| VarsNumA=length(A)$ |
Anum=length(A)$ |
| VarsNumB=length(B)$ |
|
| |
|
| C=newvect(VarsNumB,B)$ |
TMP=[]$ |
| |
for(I=0;I<length(B);I++){ |
| |
|
| for(I=0;I<VarsNumA;I++){ |
for(J=0;J<Anum;J++) |
| |
if(B[I]==A[J][0]) |
| |
break; |
| |
|
| for(J=0;J<VarsNumB;J++) |
if(J==Anum) |
| if(B[J]==A[I][0]) |
TMP=cons([B[I],B[I]],TMP)$ |
| break$ |
|
| |
|
| if(J<VarsNumB) |
|
| C[J]=A[I][1]$ |
|
| } |
} |
| |
A=append(A,TMP)$ |
| |
|
| D=0$ |
Anum=length(A)$ |
| for(I=0;I<VarsNumB;I++) |
A=map(ltov,A)$ |
| D=gcd(D,C[I])$ |
|
| |
|
| |
for(D=0,I=0;I<Anum;I++) |
| |
D=gcd(D,A[I][1])$ |
| |
|
| if(D!=0){ |
if(D!=0){ |
| C=C/D$ |
for(I=0;I<Anum;I++) |
| C=map(red,C)$ |
A[I][1]=red(A[I][1]/D)$ |
| } |
} |
| |
|
| for(L=1,D=0,I=0;I<VarsNumB;I++){ |
for(L=1,D=0,I=0;I<Anum;I++){ |
| if(type(TMP=dn(C[I]))==1) |
if(type(TMP=dn(A[I][1]))==1) |
| L=ilcm(L,TMP)$ |
L=ilcm(L,TMP)$ |
| |
|
| if(type(TMP=nm(C[I]))==1) |
if(type(TMP=nm(A[I][1]))==1) |
| D=igcd(D,TMP)$ |
D=igcd(D,TMP)$ |
| } |
} |
| |
|
| C=C*L$ |
for(I=0;I<Anum;I++) |
| if(D!=0) |
A[I][1]=A[I][1]*L$ |
| C=C/D$ |
|
| |
|
| RET=[]$ |
if(D!=0){ |
| for(I=0;I<VarsNumB;I++) |
|
| RET=cons([B[I],C[I]],RET)$ |
|
| |
|
| return RET$ |
for(I=0;I<Anum;I++) |
| |
A[I][1]=A[I][1]/D$ |
| |
} |
| |
|
| |
return map(vtol,A)$ |
| } |
} |
| |
|
| def makeret(Res,Vars,FLAG){ |
def makeret(Res,Vars,FLAG){ |
| Line 351 def makeret(Res,Vars,FLAG){ |
|
| Line 340 def makeret(Res,Vars,FLAG){ |
|
| RET[J]=ResVec[I]$ |
RET[J]=ResVec[I]$ |
| } |
} |
| |
|
| |
|
| for(J=0;J<length(Vars);J++) |
|
| RET=map(subst,RET,Vars[J], |
|
| strtov(rtostr(Vars[J])+"_deg"))$ |
|
| |
|
| for(I=0;I<VarsNum;I++) |
for(I=0;I<VarsNum;I++) |
| if(type(RET[I])!=1) |
if(type(RET[I])!=1) |
| return [1,RET]$ |
return [1,RET]$ |
| Line 477 def value2(Vars,Ans){ |
|
| Line 461 def value2(Vars,Ans){ |
|
| Res[I][0]=Vars[I]$ |
Res[I][0]=Vars[I]$ |
| Res[I][1]=Ans[I]$ |
Res[I][1]=Ans[I]$ |
| } |
} |
| |
Res=map(vtol,Res)$ |
| |
Res=vtol(Res)$ |
| |
|
| Res=getgcd(Res,Vars)$ |
Res=getgcd(Res,Vars)$ |
| |
|
| Line 511 def qcheck(PolyList,Vars,FLAG){ |
|
| Line 497 def qcheck(PolyList,Vars,FLAG){ |
|
| |
|
| VarsList=[]$ |
VarsList=[]$ |
| for(I=0;I<VarsNum;I++) |
for(I=0;I<VarsNum;I++) |
| if(member(Vars[CHAGORD[I]],Rea)) |
if(member(TMP0=Vars[CHAGORD[I]],Rea)) |
| VarsList=cons(Vars[CHAGORD[I]],VarsList)$ |
VarsList=cons(Vars[CHAGORD[I]],VarsList)$ |
| |
|
| Res=solve(reverse(SolveList),reverse(VarsList))$ |
Res=solve(reverse(SolveList),reverse(VarsList))$ |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to weight CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / lib