Blame view

fvn_test/test_matcon.f90 547 Bytes
27d3b84d6   daniau   git-svn-id: https...
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
  program test_matcon
  use fvn
  implicit none
  real(8),dimension(3,3) :: a
  real(8) :: 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