-
+ programme de test pour matev séparé en un test réel et un test complexe
test_matev_r.f90
1.27 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
program test_matev
use fvn_linear
implicit none
real(8),dimension(3,3) :: a
complex(8),dimension(3) :: evala
complex(8),dimension(3,3) :: eveca
integer :: status,i,j
call init_random_seed()
call random_number(a)
a=a*100
call fvn_matev(3,a,evala,eveca,status)
write(*,*) "The matrix :"
write (*,'(3e12.5)') a
write (*,*)
do i=1,3
write(*,'("Eigenvalue ",I3," : (",e12.5,",",e12.5,") ")') i,evala(i)
write(*,'("Modulus : ",e12.5)') abs(evala(i))
write(*,*) "Associated Eigenvector :"
do j=1,3
write(*,'("(",e12.5,",",e12.5,") ")') eveca(j,i)
end do
write(*,*)
end do
write(*,*) "Computation using matevx"
call fvn_d_matevx(3,a,evala,eveca,status)
do i=1,3
write(*,'("Eigenvalue ",I3," : (",e12.5,",",e12.5,") ")') i,evala(i)
write(*,'("Modulus : ",e12.5)') abs(evala(i))
write(*,*) "Associated Eigenvector :"
do j=1,3
write(*,'("(",e12.5,",",e12.5,") ")') eveca(j,i)
end do
write(*,*)
end do
! tri
write(*,*) "With sort option"
call fvn_matev(3,a,evala,eveca,status,.true.)
do i=1,3
write(*,'("Eigenvalue ",I3," : (",e12.5,",",e12.5,") ")') i,evala(i)
write(*,'("Modulus : ",e12.5)') abs(evala(i))
write(*,*) "Associated Eigenvector :"
do j=1,3
write(*,'("(",e12.5,",",e12.5,") ")') eveca(j,i)
end do
write(*,*)
end do
end program