lxsd / fvn

Blame view

fvn_test/test_matev.f90 535 Bytes
edit raw normal view history
27d3b84d6   daniau   git-svn-id: https... 
1
 
  program test_matev
2919a9e2d   daniau   git-svn-id: https... 
2
 
  use fvn_linear
27d3b84d6   daniau   git-svn-id: https... 
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
 
  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 (*,'(3(e12.5))') a
  write (*,*)
  do i=1,3
  write(*,'("Eigenvalue ",I3," : (",e12.5,",",e12.5,")  ")') i,evala(i)
  write(*,*) "Associated Eigenvector :"
  do j=1,3
  write(*,'("(",e12.5,",",e12.5,")  ")') eveca(j,i)
  end do
  write(*,*)
  end do
  end program