Return to weight CVS log | Up to [local] / OpenXM_contrib2 / asir2000 / lib |
version 1.18, 2004/01/07 06:33:31 | version 1.19, 2004/01/07 06:53:11 | ||
---|---|---|---|
|
|
||
#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){ | ||
|
|
||
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$ |