Mercurial > octave-nkf

view libcruft/misc/r1mach.f @ 7789:82be108cc558

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
First attempt at single precision tyeps * * * corrections to qrupdate single precision routines * * * prefer demotion to single over promotion to double * * * Add single precision support to log2 function * * * Trivial PROJECT file update * * * Cache optimized hermitian/transpose methods * * * Add tests for tranpose/hermitian and ChangeLog entry for new transpose code
author David Bateman <dbateman@free.fr>
date Sun, 27 Apr 2008 22:34:17 +0200
parents
children 32b15d5c3147
line wrap: on
line source

      real function r1mach (i)
      integer i
      logical init
      real rmach(5)
      save init, rmach
      data init /.false./
      if (.not. init) then
        call smachar (rmach(1), rmach(2), rmach(3), rmach(4), rmach(5))
        init = .true.
      endif
      if (i .lt. 1  .or.  i .gt. 5) goto 999
      r1mach = rmach(i)
      return
  999 write(*,1999) i
 1999 format(' s1mach - i out of bounds', i10)
      call xstopx (' ')
      r1mach = 0
      end