version 1.12, 2003/12/08 07:18:14 |
version 1.14, 2003/12/11 15:12:35 |
Line 249 def roundret(V){ |
|
Line 249 def roundret(V){ |
|
VN=size(V)[0]$ |
VN=size(V)[0]$ |
|
|
RET0=V$ |
RET0=V$ |
for(I=2;I<1000;I++){ |
for(I=1;I<1000;I++){ |
RET1=I*RET0$ |
RET1=I*RET0$ |
for(J=0;J<VN;J++){ |
for(J=0;J<VN;J++){ |
X=drint(RET1[J])$ |
X=drint(RET1[J])$ |
Line 450 def leastsq(NormMat,ExpMat,Vars,FLAG){ |
|
Line 450 def leastsq(NormMat,ExpMat,Vars,FLAG){ |
|
if(nonposdegchk(Res)){ |
if(nonposdegchk(Res)){ |
TMP1=makeret(Res,Vars,1)$ |
TMP1=makeret(Res,Vars,1)$ |
if(TMP1[0]==0){ |
if(TMP1[0]==0){ |
TMP=roundret(TMP1[1]*1.0)$ |
TMP=roundret(TMP1[1])$ |
if(TMP!=[]) |
if(TMP!=[]) |
RET=append(RET,wsort(TMP1[1],Vars,TMP,FLAG))$ |
RET=append(RET,wsort(TMP1[1],Vars,TMP,FLAG))$ |
|
|
Line 566 def weightr(ExpMat,Vars,PolyListNum,OneMat,FLAG){ |
|
Line 566 def weightr(ExpMat,Vars,PolyListNum,OneMat,FLAG){ |
|
|
|
TMP1=makeret(Res,Vars,1)$ |
TMP1=makeret(Res,Vars,1)$ |
if(TMP1[0]==0){ |
if(TMP1[0]==0){ |
TMP=roundret(TMP1[1]*1.0)$ |
TMP=roundret(TMP1[1])$ |
if(TMP!=[]) |
if(TMP!=[]) |
RET=append(RET,wsort(TMP1[1],Vars, |
RET=append(RET,wsort(TMP1[1],Vars, |
TMP,FLAG))$ |
TMP,FLAG))$ |