| version 1.13, 2003/12/08 07:30:43 |
version 1.15, 2004/01/06 08:38:05 |
| 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 559 def weightr(ExpMat,Vars,PolyListNum,OneMat,FLAG){ |
|
| Line 559 def weightr(ExpMat,Vars,PolyListNum,OneMat,FLAG){ |
|
| if(member(Rea[I],Vars)) |
if(member(Rea[I],Vars)) |
| TMP=cons(Rea[I],TMP)$ |
TMP=cons(Rea[I],TMP)$ |
| |
|
| TMP=cons(ExtVars[F],TMP)$ |
if(member(ExtVars[F],Vars)) |
| |
TMP=cons(ExtVars[F],TMP)$ |
| |
|
| Res=getgcd(Res,TMP)$ |
Res=getgcd(Res,TMP)$ |
| |
|
| 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, |
RET=append(RET,wsort(TMP1[1],Vars, |
| TMP,FLAG))$ |
TMP,FLAG))$ |
![[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