| Line |
Branch |
Exec |
Source |
| 1 |
|
|
! |
| 2 |
|
|
! $Header$ |
| 3 |
|
|
! |
| 4 |
|
✗ |
SUBROUTINE eigen( e,d) |
| 5 |
|
|
IMPLICIT NONE |
| 6 |
|
|
!----------------------------------------------------------------------- |
| 7 |
|
|
! INCLUDE 'dimensions.h' |
| 8 |
|
|
! |
| 9 |
|
|
! dimensions.h contient les dimensions du modele |
| 10 |
|
|
! ndm est tel que iim=2**ndm |
| 11 |
|
|
!----------------------------------------------------------------------- |
| 12 |
|
|
|
| 13 |
|
|
INTEGER iim,jjm,llm,ndm |
| 14 |
|
|
|
| 15 |
|
|
PARAMETER (iim= 32,jjm=32,llm=39,ndm=1) |
| 16 |
|
|
|
| 17 |
|
|
!----------------------------------------------------------------------- |
| 18 |
|
|
real :: e( iim,iim ), d( iim ) |
| 19 |
|
|
real :: asm( iim ) |
| 20 |
|
|
integer :: im,i,j |
| 21 |
|
|
im=iim |
| 22 |
|
|
c |
| 23 |
|
✗ |
DO 48 i = 1,im |
| 24 |
|
✗ |
asm( i ) = d( im-i+1 ) |
| 25 |
|
✗ |
48 CONTINUE |
| 26 |
|
✗ |
DO 49 i = 1,iim |
| 27 |
|
✗ |
d( i ) = asm( i ) |
| 28 |
|
✗ |
49 CONTINUE |
| 29 |
|
|
c |
| 30 |
|
|
c PRINT 70,d |
| 31 |
|
|
70 FORMAT(5x,'Valeurs propres',/,8(1x,8f10.4,/),/) |
| 32 |
|
✗ |
print * |
| 33 |
|
|
c |
| 34 |
|
✗ |
DO 51 i = 1,im |
| 35 |
|
✗ |
DO 52 j = 1,im |
| 36 |
|
✗ |
asm( j ) = e( i , im-j+1 ) |
| 37 |
|
✗ |
52 CONTINUE |
| 38 |
|
✗ |
DO 50 j = 1,im |
| 39 |
|
✗ |
e( i,j ) = asm( j ) |
| 40 |
|
✗ |
50 CONTINUE |
| 41 |
|
✗ |
51 CONTINUE |
| 42 |
|
|
|
| 43 |
|
✗ |
RETURN |
| 44 |
|
|
END |
| 45 |
|
|
|