Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/filtrez/jacobi.F90
- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/filtrez/jacobi.F90
r1907 r5082 38 38 ENDDO 39 39 ENDDO 40 IF(SM .EQ.0.)RETURN41 IF(I .LT.4)THEN40 IF(SM==0.)RETURN 41 IF(I<4)THEN 42 42 TRESH=0.2*SM/N**2 43 43 ELSE … … 47 47 DO IQ=IP+1,N 48 48 G=100.*ABS(A(IP,IQ)) 49 IF((I .GT.4).AND.(ABS(D(IP))+G.EQ.ABS(D(IP))) &50 .AND.(ABS(D(IQ))+G .EQ.ABS(D(IQ))))THEN49 IF((I>4).AND.(ABS(D(IP))+G==ABS(D(IP))) & 50 .AND.(ABS(D(IQ))+G==ABS(D(IQ))))THEN 51 51 A(IP,IQ)=0. 52 ELSE IF(ABS(A(IP,IQ)) .GT.TRESH)THEN52 ELSE IF(ABS(A(IP,IQ))>TRESH)THEN 53 53 H=D(IQ)-D(IP) 54 IF(ABS(H)+G .EQ.ABS(H))THEN54 IF(ABS(H)+G==ABS(H))THEN 55 55 T=A(IP,IQ)/H 56 56 ELSE 57 57 THETA=0.5*H/A(IP,IQ) 58 58 T=1./(ABS(THETA)+SQRT(1.+THETA**2)) 59 IF(THETA .LT.0.)T=-T59 IF(THETA<0.)T=-T 60 60 ENDIF 61 61 C=1./SQRT(1+T**2)
Note: See TracChangeset
for help on using the changeset viewer.