lxsd / fvn

Blame view

fvn_test/test_matcon.f90 576 Bytes
edit raw normal view history
27d3b84d6   daniau   git-svn-id: https... 
1
 
  program test_matcon
2919a9e2d   daniau   git-svn-id: https... 
2
 
  use fvn_linear
27d3b84d6   daniau   git-svn-id: https... 
3
 
  implicit none
f6bacaf83   cwaterkeyn   ChW 11/09: ANSI c... 
4
5
 
  real(kind=dp_kind),dimension(3,3) :: a
  real(kind=dp_kind) :: rcond
27d3b84d6   daniau   git-svn-id: https... 
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
 
  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