Mercurial > octave-nkf
view libcruft/misc/d1mach.f @ 12119:e320928eeb3a release-3-2-x release-3-2-4
version 3.2.4
author | Jaroslav Hajek <highegg@gmail.com> |
---|---|
date | Fri, 22 Jan 2010 12:43:12 +0100 |
parents | 587d268cf64e |
children |
line wrap: on
line source
double precision function d1mach (i) integer i logical init double precision dmach(5) double precision dlamch external dlamch save init, dmach data init /.false./ if (.not. init) then dmach(1) = dlamch ('u') dmach(2) = dlamch ('o') dmach(3) = dlamch ('e') dmach(4) = dlamch ('p') dmach(5) = log10 (dlamch ('b')) init = .true. endif if (i .lt. 1 .or. i .gt. 5) goto 999 d1mach = dmach(i) return 999 write (*, 1999) i 1999 format (' d1mach - i out of bounds', i10) call xstopx (' ') d1mach = 0 end