Changeset 5246 for LMDZ6/trunk/libf/filtrez/eigen_sort.f90
- Timestamp:
- Oct 21, 2024, 2:58:45 PM (23 hours ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/filtrez/eigen_sort.f90
r5245 r5246 2 2 ! $Header$ 3 3 ! 4 5 INTEGERn,np6 REALd(np),v(np,np)7 INTEGERi,j,k8 REALp4 SUBROUTINE eigen_sort(d,v,n,np) 5 INTEGER :: n,np 6 REAL :: d(np),v(np,np) 7 INTEGER :: i,j,k 8 REAL :: p 9 9 10 DO i=1,n-1 11 k=i 12 p=d(i) 13 DO j=i+1,n 14 IF(d(j).ge.p) THEN 15 k=j 16 p=d(j) 17 ENDIF 18 ENDDO 19 20 IF(k.ne.i) THEN 21 d(k)=d(i) 22 d(i)=p 23 DO j=1,n 24 p=v(j,i) 25 v(j,i)=v(j,k) 26 v(j,k)=p 27 ENDDO 28 ENDIF 29 ENDDO 10 DO i=1,n-1 11 k=i 12 p=d(i) 13 DO j=i+1,n 14 IF(d(j).ge.p) THEN 15 k=j 16 p=d(j) 17 ENDIF 18 ENDDO 30 19 31 RETURN 32 END 20 IF(k.ne.i) THEN 21 d(k)=d(i) 22 d(i)=p 23 DO j=1,n 24 p=v(j,i) 25 v(j,i)=v(j,k) 26 v(j,k)=p 27 ENDDO 28 ENDIF 29 ENDDO 30 31 RETURN 32 END SUBROUTINE eigen_sort
Note: See TracChangeset
for help on using the changeset viewer.