program test_matev
use fvn_linear
implicit none
complex(8),dimension(3,3) :: a
real(8),dimension(3,3) :: ra,ia
complex(8),dimension(3) :: evala
complex(8),dimension(3,3) :: eveca
integer :: status,i,j

call init_random_seed()
call random_number(ra)
call random_number(ia)
a=ra+fvn_i*ia
a=a*100
call fvn_matev(3,a,evala,eveca,status)

write(*,*) "The matrix :"
write (*,'(3("(",e12.5,",",e12.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

! 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