[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.18 and 1.23

version 1.18, 2004/01/07 06:33:31 version 1.23, 2004/01/08 06:48:32
Line 4  load("gr")$
Line 4  load("gr")$
 #define EPS 1E-6  #define EPS 1E-6
 #define TINY 1E-20  #define TINY 1E-20
 #define MAX_ITER 100  #define MAX_ITER 100
   #define ROUND_THRESHOLD 0.4
   
 def rotate(A,I,J,K,L,C,S){  def rotate(A,I,J,K,L,C,S){
   
Line 352  def roundret(V){
Line 353  def roundret(V){
                 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])$
                         if(dabs(X-RET1[J])<0.2)                          if(dabs(X-RET1[J])<ROUND_THRESHOLD)
                                 RET1[J]=X$                                  RET1[J]=X$
                         else                          else
                                 break$                                  break$
Line 564  def leastsq(NormMat,ExpMat,Vars,FLAG,ID){
Line 565  def leastsq(NormMat,ExpMat,Vars,FLAG,ID){
                         return RET$                          return RET$
                 }                  }
                 else{                  else{
                         RET=append(RET,[[ID,Vars,vtol(TMP1[1]*1.0)]])$                          RET=append(RET,[[ID,Vars,vtol(TMP1[1])]])$
                         return RET$                          return RET$
                 }                  }
         }          }
Line 689  def weight(PolyList,Vars,FLAG){
Line 690  def weight(PolyList,Vars,FLAG){
   
         RET=append(RET,TMP[1])$          RET=append(RET,TMP[1])$
   
         RET=append(RET,leastsq(TMP[0],ExpMat,Vars,FLAG,3))$          TMP0=leastsq(TMP[0],ExpMat,Vars,FLAG,3)$
   
           RET=append(RET,TMP0)$
   
         ExpMat=qsort(ExpMat,junban)$          ExpMat=qsort(ExpMat,junban)$
   
         ExpMat2=[]$          ExpMat2=[]$
Line 701  def weight(PolyList,Vars,FLAG){
Line 704  def weight(PolyList,Vars,FLAG){
         if(size(ExpMat)[0]!=length(ExpMat2)){          if(size(ExpMat)[0]!=length(ExpMat2)){
                 ExpMat=newvect(length(ExpMat2),ExpMat2)$                  ExpMat=newvect(length(ExpMat2),ExpMat2)$
                 RET=append(RET,leastsq(0,ExpMat,Vars,FLAG,5))$                  RET=append(RET,leastsq(0,ExpMat,Vars,FLAG,5))$
           }
           else{
                   TMP0=map(ltov,TMP0)$
   
                   for(I=0;I<length(TMP0);I++)
                           TMP0[I][0]+=2$
   
                   TMP0=map(vtol,TMP0)$
   
                   RET=append(RET,TMP0)$
         }          }
   
         return RET$          return RET$
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.23
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>