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 |
|
|
|