=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/lib/weight,v retrieving revision 1.1 retrieving revision 1.3 diff -u -p -r1.1 -r1.3 --- OpenXM_contrib2/asir2000/lib/weight 2003/10/15 07:06:02 1.1 +++ OpenXM_contrib2/asir2000/lib/weight 2003/11/05 08:26:57 1.3 @@ -1,232 +1,560 @@ -#include -load("solve")$ - -def nonposdegchk(Res){ - - for(I=0;IB[I]) - return -1$ - } - - return 0$ -} - -def weight(PolyList,Vars){ - - dp_ord(2)$ - - PolyListNum=length(PolyList)$ - - ExpMat=[]$ - for(I=0;IED) - continue$ - for(K=ST;K<=ED;K++){ - NormMat[I][J]=NormMat[I][J]+ - extmat(ExpMat,OneMat,ExpMatColNum,ExtMatColNum,K,I)* - extmat(ExpMat,OneMat,ExpMatColNum,ExtMatColNum,K,J)$ - } - } - - for(I=0;IED) - continue$ - - for(K=ST;K<=ED;K++){ - NormMat[I][ExtMatColNum-1]=NormMat[I][ExtMatColNum-1]+ - extmat(ExpMat,OneMat,ExpMatColNum,ExtMatColNum,K,I)* - extmat(ExpMat,OneMat,ExpMatColNum,ExtMatColNum,K,ExtMatColNum-1)$ - } - } - - ExtVars=Vars$ - for(I=0;Inmono(B) ? 1:0))$ +} + +def junban2(A,B){ + + for(I=0;IB[I]) + return -1$ + } + + return 0$ +} + +def roundret(V){ + + VN=length(V)$ + RET0=newvect(VN,V)$ + + for(I=1;I<1000;I++){ + RET1=I*RET0$ + for(J=0;J