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