gcc/testsuite/gfortran.dg/read_float_2.f03

   1 ! { dg-do run }
   2 ! Contributed by Dominique Dhumieres <dominiq@lps.ens.fr>
   3
   4 character(15) :: str="+ .339  567+2"
   5 real, parameter :: should_be = .339567e2
   6 real, parameter :: eps = 10 * epsilon (should_be)
   7 real :: x, y
   8
   9 read(str,'(BN,F15.6)') x
  10 print *, x
  11 read(str,'(G15.7)') y
  12 print *, y
  13
  14 if (abs (x - should_be) > eps .or. abs (y - should_be) > eps) then
  15   STOP 1
  16 end if
  17
  18 end