version 1.34, 2004/05/05 13:15:01 |
version 1.36, 2004/05/06 02:55:39 |
Line 780 def leastsq(NormMat,ExpMat,Vars,FLAG,ID){ |
|
Line 780 def leastsq(NormMat,ExpMat,Vars,FLAG,ID){ |
|
RET=append(RET,[[ID,Vars,vtol(TMP1[1])]])$ |
RET=append(RET,[[ID,Vars,vtol(TMP1[1])]])$ |
} |
} |
|
|
return RET$ |
return [NormMat0,RET]$ |
} |
} |
|
|
def weight(PolyList,Vars,FLAG){ |
def weight(PolyList,Vars,FLAG){ |
Line 827 def weight(PolyList,Vars,FLAG){ |
|
Line 827 def weight(PolyList,Vars,FLAG){ |
|
if(TMP[1]!=[]) |
if(TMP[1]!=[]) |
RET=append(RET,TMP[1])$ |
RET=append(RET,TMP[1])$ |
else{ |
else{ |
RET=append(RET,leastsq(0,ExpMat,Vars,FLAG,7))$ |
TMP=leastsq(0,ExpMat,Vars,FLAG,7)$ |
|
RET=append(RET,TMP[1])$ |
} |
} |
|
|
ExpMat=qsort(ExpMat,junban)$ |
ExpMat=qsort(ExpMat,junban)$ |
Line 843 def weight(PolyList,Vars,FLAG){ |
|
Line 844 def weight(PolyList,Vars,FLAG){ |
|
if(TMP[1]!=[]) |
if(TMP[1]!=[]) |
RET=append(RET,TMP[1])$ |
RET=append(RET,TMP[1])$ |
else{ |
else{ |
RET=append(RET,leastsq(0,ExpMat,Vars,FLAG,9))$ |
TMP=leastsq(0,ExpMat,Vars,FLAG,9)$ |
|
RET=append(RET,TMP[1])$ |
} |
} |
|
|
} |
} |
else{ |
else{ |
TMP0=map(ltov,TMP0)$ |
TMP=map(ltov,TMP[1])$ |
|
|
for(I=0;I<length(TMP0);I++) |
for(I=0;I<length(TMP);I++){ |
if(TMP0[I][0]==3) |
if(TMP[I][0]==3) |
TMP0[I][0]=5$ |
TMP[I][0]=5$ |
else if(TMP0[I][0]==4) |
else if(TMP[I][0]==4) |
TMP0[I][0]=6$ |
TMP[I][0]=6$ |
|
else if(TMP[I][0]==7) |
|
TMP[I][0]=9$ |
|
else if(TMP[I][0]==8) |
|
TMP[I][0]=10$ |
|
} |
|
|
TMP0=map(vtol,TMP0)$ |
TMP=map(vtol,TMP)$ |
|
|
RET=append(RET,TMP0)$ |
RET=append(RET,TMP)$ |
} |
} |
|
|
return RET$ |
return RET$ |