program test_matcon
use fvn_linear
implicit none
real(kind=dp_kind),dimension(3,3) :: a
real(kind=dp_kind) :: rcond
integer :: status,i
call init_random_seed()
call random_number(a)
a=a*100
call fvn_matcon(3,a,rcond,status)
write(*,*) "Reasonnably conditionned matrix"
do i=1,3
    write (*,'(3(e12.5))') a(i,:)
end do
write (*,*)
write (*,*) "Cond = ",rcond
write (*,*)
write (*,*)
a(1,1)=a(1,1)*1d9
write(*,*) "Badly conditionned matrix"
do i=1,3
    write (*,'(3(e12.5))') a(i,:)
end do
call fvn_matcon(3,a,rcond,status)
write (*,*)
write (*,*) "Cond = ",rcond

end program