![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to weight CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / lib|
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$ |