trunk/filtrez/eigen_sort.f


! $Header: /home/cvsroot/LMDZ4/libf/filtrez/eigen_sort.F,v 1.1.1.1 2004/05/19
! 12:53:09 lmdzadmin Exp $

SUBROUTINE eigen_sort(d, v, n, np)
  INTEGER n, np
  REAL d(np), v(np, np)
  INTEGER i, j, k
  REAL p

  DO i = 1, n - 1
    k = i
    p = d(i)
    DO j = i + 1, n
      IF (d(j)>=p) THEN
        k = j
        p = d(j)
      END IF
    END DO

    IF (k/=i) THEN
      d(k) = d(i)
      d(i) = p
      DO j = 1, n
        p = v(j, i)
        v(j, i) = v(j, k)
        v(j, k) = p
      END DO
    END IF
  END DO

  RETURN
END SUBROUTINE eigen_sort
1	guez	3
2	guez	81	! $Header: /home/cvsroot/LMDZ4/libf/filtrez/eigen_sort.F,v 1.1.1.1 2004/05/19
3			! 12:53:09 lmdzadmin Exp $
4	guez	3
5	guez	81	SUBROUTINE eigen_sort(d, v, n, np)
6			INTEGER n, np
7			REAL d(np), v(np, np)
8			INTEGER i, j, k
9			REAL p
10
11			DO i = 1, n - 1
12			k = i
13			p = d(i)
14			DO j = i + 1, n
15			IF (d(j)>=p) THEN
16			k = j
17			p = d(j)
18			END IF
19			END DO
20
21			IF (k/=i) THEN
22			d(k) = d(i)
23			d(i) = p
24			DO j = 1, n
25			p = v(j, i)
26			v(j, i) = v(j, k)
27			v(j, k) = p
28			END DO
29			END IF
30			END DO
31
32			RETURN
33			END SUBROUTINE eigen_sort