Rev | Line | |
---|
[8] | 1 | program check |
---|
| 2 | |
---|
| 3 | implicit none |
---|
| 4 | |
---|
| 5 | real(kind=8), parameter :: tol = 1.e-9 |
---|
| 6 | real(kind=8) :: x, y |
---|
| 7 | integer :: ios1, ios2 |
---|
| 8 | integer :: irc |
---|
| 9 | integer :: nl = 0 |
---|
| 10 | |
---|
| 11 | open(10,file='temp1',form='formatted',status='old') |
---|
| 12 | open(11,file='temp2',form='formatted',status='old') |
---|
| 13 | |
---|
| 14 | irc = 0 |
---|
| 15 | |
---|
| 16 | do |
---|
| 17 | |
---|
| 18 | nl = nl + 1 |
---|
| 19 | |
---|
| 20 | read(10,*,iostat=ios1) x |
---|
| 21 | read(11,*,iostat=ios2) y |
---|
| 22 | |
---|
| 23 | if (ios1 /= 0 .and. ios2 /= 0) then |
---|
| 24 | exit |
---|
| 25 | else if (ios1 /= 0 .xor. ios2 /= 0) then |
---|
| 26 | irc = 1 |
---|
| 27 | exit |
---|
| 28 | ! Mizture de test relatif et de test absolu |
---|
| 29 | else if ( abs(x-y) <= max(tol*(abs(x)+abs(y)),tol) ) then |
---|
| 30 | continue |
---|
| 31 | else |
---|
| 32 | irc = 1 |
---|
| 33 | exit |
---|
| 34 | end if |
---|
| 35 | |
---|
| 36 | end do |
---|
| 37 | |
---|
| 38 | if (irc == 0) then |
---|
| 39 | print *,'Test OK' |
---|
| 40 | else |
---|
| 41 | print *,'Problemes ligne', nl |
---|
| 42 | stop 1 |
---|
| 43 | end if |
---|
| 44 | |
---|
| 45 | end program check |
---|
Note: See
TracBrowser
for help on using the repository browser.