Go to the documentation of this file.00001 double precision function r1mach (i)
00002 integer i
00003 logical init
00004 double precision rmach(5)
00005 double precision slamch
00006 external slamch
00007 save init, rmach
00008 data init /.false./
00009 if (.not. init) then
00010 rmach(1) = slamch ('u')
00011 rmach(2) = slamch ('o')
00012 rmach(3) = slamch ('e')
00013 rmach(4) = slamch ('p')
00014 rmach(5) = log10 (slamch ('b'))
00015 init = .true.
00016 endif
00017 if (i .lt. 1 .or. i .gt. 5) goto 999
00018 r1mach = rmach(i)
00019 return
00020 999 write (*, 1999) i
00021 1999 format (' r1mach - i out of bounds', i10)
00022 call xstopx (' ')
00023 d1mach = 0
00024 end
1.7.1
