[BACK]Return to weight CVS log [TXT][DIR] Up to [local] / OpenXM_contrib2 / asir2000 / lib

Diff for /OpenXM_contrib2/asir2000/lib/weight between version 1.12 and 1.14

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))$

Legend:
Removed from v.1.12  
changed lines
  Added in v.1.14

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>