summarylogtreecommitdiffstats
diff options
context:
space:
mode:
authorPhil Ruffwind2015-05-14 08:09:00 +0000
committerPhil Ruffwind2015-06-03 19:49:13 -0400
commit2787fe3a075a521e31a3bdacc252c57c24029164 (patch)
tree1d0aed46a01902d0cee03b9f11e27e273c61441d
downloadaur-2787fe3a075a521e31a3bdacc252c57c24029164.tar.gz
Initial version from AUR3
Diffstat
-rw-r--r--.SRCINFO20
-rw-r--r--LICENSE12
-rw-r--r--Makefile297
-rw-r--r--PKGBUILD33
-rw-r--r--slatec.patch2188
5 files changed, 2550 insertions, 0 deletions
diff --git a/.SRCINFO b/.SRCINFO
new file mode 100644
index 000000000000..ae8feb711bba
--- /dev/null
+++ b/.SRCINFO
@@ -0,0 +1,20 @@
+pkgbase = slatec
+ pkgdesc = SLATEC Common Mathematical Library
+ pkgver = 4.1
+ pkgrel = 2
+ url = http://netlib.org/slatec
+ arch = i686
+ arch = x86_64
+ license = PD
+ makedepends = gcc-fortran
+ source = http://netlib.org/slatec/slatec_src.tgz
+ source = LICENSE
+ source = Makefile
+ source = slatec.patch
+ sha256sums = df009d9ef9c18aae06ce68711b1ae108d3533b4f174582c3cbea0915c4fdfe01
+ sha256sums = d83a2ee63ee4bbd56d59072a0778e06685fa001047f153a3f0293d30849c8046
+ sha256sums = 7d8fb187324596d4d882d895c31db0afddf2984a933ab6e623de576d5f36558e
+ sha256sums = 4fb39658c6da25618d20e6aec0680c48d4666ab0eb1322f6e7f93a34bb51ce4f
+
+pkgname = slatec
+
diff --git a/LICENSE b/LICENSE
new file mode 100644
index 000000000000..4dc38da53a18
--- /dev/null
+++ b/LICENSE
@@ -0,0 +1,12 @@
+Public Domain
+=============
+
+This material was prepared as an account of work sponsored by the
+United States Government. Neither the United States, nor the
+Department of Energy, nor the Department of Defense, nor any of
+their employees, nor any of their contractors, subcontractors, or
+their employees, makes any warranty, expressed or implied, or
+assumes any legal liability or responsibility for the accuracy,
+completeness, or usefulness of any information, apparatus, product,
+or process disclosed, or represents that its use would not infringe
+upon privately owned rights.
diff --git a/Makefile b/Makefile
new file mode 100644
index 000000000000..8f00f4ef6310
--- /dev/null
+++ b/Makefile
@@ -0,0 +1,297 @@
+PREFIX=/usr/local
+ARFLAGS=-cru
+FFLAGS=-O2 -fPIC
+
+major=4
+version=$(major).1
+
+objs=src/aaaaaa.o src/acosh.o src/ai.o src/aie.o src/albeta.o src/algams.o \
+src/ali.o src/alngam.o src/alnrel.o src/asinh.o src/asyik.o src/asyjy.o \
+src/atanh.o src/avint.o src/bakvec.o src/balanc.o src/balbak.o src/bandr.o \
+src/bandv.o src/bcrh.o src/bdiff.o src/besi.o src/besi0.o src/besi0e.o \
+src/besi1.o src/besi1e.o src/besj.o src/besj0.o src/besj1.o src/besk.o \
+src/besk0.o src/besk0e.o src/besk1.o src/besk1e.o src/beskes.o src/besknu.o \
+src/besks.o src/besy.o src/besy0.o src/besy1.o src/besynu.o src/beta.o \
+src/betai.o src/bfqad.o src/bi.o src/bie.o src/binom.o src/bint4.o \
+src/bintk.o src/bisect.o src/bkias.o src/bkisr.o src/bksol.o src/blktr1.o \
+src/blktri.o src/bndacc.o src/bndsol.o src/bnfac.o src/bnslv.o src/bqr.o \
+src/bsgq8.o src/bskin.o src/bspdoc.o src/bspdr.o src/bspev.o src/bsplvd.o \
+src/bsplvn.o src/bsppp.o src/bspvd.o src/bspvn.o src/bsqad.o src/bsrh.o \
+src/bvalu.o src/c0lgmc.o src/c1merg.o src/c9lgmc.o src/c9ln2r.o src/cacai.o \
+src/cacon.o src/cacos.o src/cacosh.o src/cairy.o src/carg.o src/casin.o \
+src/casinh.o src/casyi.o src/catan.o src/catan2.o src/catanh.o src/caxpy.o \
+src/cbabk2.o src/cbal.o src/cbesh.o src/cbesi.o src/cbesj.o src/cbesk.o \
+src/cbesy.o src/cbeta.o src/cbinu.o src/cbiry.o src/cbknu.o src/cblkt1.o \
+src/cblktr.o src/cbrt.o src/cbuni.o src/cbunk.o src/ccbrt.o src/cchdc.o \
+src/cchdd.o src/cchex.o src/cchud.o src/ccmpb.o src/ccopy.o src/ccosh.o \
+src/ccot.o src/cdcdot.o src/cdcor.o src/cdcst.o src/cdiv.o src/cdntl.o \
+src/cdntp.o src/cdotc.o src/cdotu.o src/cdpsc.o src/cdpst.o src/cdriv1.o \
+src/cdriv2.o src/cdriv3.o src/cdscl.o src/cdstp.o src/cdzro.o src/cexprl.o \
+src/cfftb.o src/cfftb1.o src/cfftf.o src/cfftf1.o src/cffti.o src/cffti1.o \
+src/cfod.o src/cg.o src/cgamma.o src/cgamr.o src/cgbco.o src/cgbdi.o \
+src/cgbfa.o src/cgbmv.o src/cgbsl.o src/cgeco.o src/cgedi.o src/cgeev.o \
+src/cgefa.o src/cgefs.o src/cgeir.o src/cgemm.o src/cgemv.o src/cgerc.o \
+src/cgeru.o src/cgesl.o src/cgtsl.o src/ch.o src/chbmv.o src/chemm.o \
+src/chemv.o src/cher.o src/cher2.o src/cher2k.o src/cherk.o src/chfcm.o \
+src/chfdv.o src/chfev.o src/chfie.o src/chico.o src/chidi.o src/chiev.o \
+src/chifa.o src/chisl.o src/chkder.o src/chkpr4.o src/chkprm.o src/chksn4.o \
+src/chksng.o src/chpco.o src/chpdi.o src/chpfa.o src/chpmv.o src/chpr.o \
+src/chpr2.o src/chpsl.o src/chu.o src/cinvit.o src/ckscl.o src/clbeta.o \
+src/clngam.o src/clnrel.o src/clog10.o src/cmgnbn.o src/cmlri.o src/cmpcsg.o \
+src/cmposd.o src/cmposn.o src/cmposp.o src/cmptr3.o src/cmptrx.o src/cnbco.o \
+src/cnbdi.o src/cnbfa.o src/cnbfs.o src/cnbir.o src/cnbsl.o src/combak.o \
+src/comhes.o src/comlr.o src/comlr2.o src/compb.o src/comqr.o src/comqr2.o \
+src/cortb.o src/corth.o src/cosdg.o src/cosgen.o src/cosqb.o src/cosqb1.o \
+src/cosqf.o src/cosqf1.o src/cosqi.o src/cost.o src/costi.o src/cot.o \
+src/cpadd.o src/cpbco.o src/cpbdi.o src/cpbfa.o src/cpbsl.o src/cpevl.o \
+src/cpevlr.o src/cpoco.o src/cpodi.o src/cpofa.o src/cpofs.o src/cpoir.o \
+src/cposl.o src/cppco.o src/cppdi.o src/cppfa.o src/cppsl.o src/cpqr79.o \
+src/cproc.o src/cprocp.o src/cprod.o src/cprodp.o src/cpsi.o src/cptsl.o \
+src/cpzero.o src/cqrdc.o src/cqrsl.o src/crati.o src/crotg.o src/cs1s2.o \
+src/cscal.o src/cseri.o src/csevl.o src/cshch.o src/csico.o src/csidi.o \
+src/csifa.o src/csinh.o src/csisl.o src/cspco.o src/cspdi.o src/cspfa.o \
+src/cspsl.o src/csroot.o src/csrot.o src/csscal.o src/csvdc.o src/cswap.o \
+src/csymm.o src/csyr2k.o src/csyrk.o src/ctan.o src/ctanh.o src/ctbmv.o \
+src/ctbsv.o src/ctpmv.o src/ctpsv.o src/ctrco.o src/ctrdi.o src/ctrmm.o \
+src/ctrmv.o src/ctrsl.o src/ctrsm.o src/ctrsv.o src/cuchk.o src/cunhj.o \
+src/cuni1.o src/cuni2.o src/cunik.o src/cunk1.o src/cunk2.o src/cuoik.o \
+src/cv.o src/cwrsk.o src/d1mach.o src/d1merg.o src/d1mpyq.o src/d1updt.o \
+src/d9aimp.o src/d9atn1.o src/d9b0mp.o src/d9b1mp.o src/d9chu.o src/d9gmic.o \
+src/d9gmit.o src/d9knus.o src/d9lgic.o src/d9lgit.o src/d9lgmc.o src/d9ln2r.o \
+src/d9pak.o src/d9upak.o src/dacosh.o src/dai.o src/daie.o src/dasinh.o \
+src/dasum.o src/dasyik.o src/dasyjy.o src/datanh.o src/davint.o src/daws.o \
+src/daxpy.o src/dbcg.o src/dbdiff.o src/dbesi.o src/dbesi0.o src/dbesi1.o \
+src/dbesj.o src/dbesj0.o src/dbesj1.o src/dbesk.o src/dbesk0.o src/dbesk1.o \
+src/dbesks.o src/dbesy.o src/dbesy0.o src/dbesy1.o src/dbeta.o src/dbetai.o \
+src/dbfqad.o src/dbhin.o src/dbi.o src/dbie.o src/dbinom.o src/dbint4.o \
+src/dbintk.o src/dbkias.o src/dbkisr.o src/dbndac.o src/dbndsl.o src/dbnfac.o \
+src/dbnslv.o src/dbocls.o src/dbols.o src/dbolsm.o src/dbsgq8.o src/dbsi0e.o \
+src/dbsi1e.o src/dbsk0e.o src/dbsk1e.o src/dbskes.o src/dbskin.o src/dbsknu.o \
+src/dbspdr.o src/dbspev.o src/dbsppp.o src/dbspvd.o src/dbspvn.o src/dbsqad.o \
+src/dbsynu.o src/dbvalu.o src/dcbrt.o src/dcdot.o src/dcfod.o src/dcg.o \
+src/dcgn.o src/dcgs.o src/dchdc.o src/dchdd.o src/dchex.o src/dchfcm.o \
+src/dchfdv.o src/dchfev.o src/dchfie.o src/dchkw.o src/dchu.o src/dchud.o \
+src/dckder.o src/dcopy.o src/dcopym.o src/dcosdg.o src/dcot.o src/dcov.o \
+src/dcpplt.o src/dcsevl.o src/dcv.o src/ddaini.o src/ddajac.o src/ddanrm.o \
+src/ddaslv.o src/ddassl.o src/ddastp.o src/ddatrp.o src/ddaws.o src/ddawts.o \
+src/ddcor.o src/ddcst.o src/ddeabm.o src/ddebdf.o src/dderkf.o src/ddes.o \
+src/ddntl.o src/ddntp.o src/ddoglg.o src/ddot.o src/ddpsc.o src/ddpst.o \
+src/ddriv1.o src/ddriv2.o src/ddriv3.o src/ddscl.o src/ddstp.o src/ddzro.o \
+src/de1.o src/deabm.o src/debdf.o src/defc.o src/defcmn.o src/defe4.o \
+src/defehl.o src/defer.o src/dei.o src/denorm.o src/derf.o src/derfc.o \
+src/derkf.o src/derkfs.o src/des.o src/dexint.o src/dexprl.o src/dfac.o \
+src/dfc.o src/dfcmn.o src/dfdjc1.o src/dfdjc3.o src/dfehl.o src/dfspvd.o \
+src/dfspvn.o src/dfulmt.o src/dfzero.o src/dgami.o src/dgamic.o src/dgamit.o \
+src/dgamlm.o src/dgamln.o src/dgamma.o src/dgamr.o src/dgamrn.o src/dgaus8.o \
+src/dgbco.o src/dgbdi.o src/dgbfa.o src/dgbmv.o src/dgbsl.o src/dgeco.o \
+src/dgedi.o src/dgefa.o src/dgefs.o src/dgemm.o src/dgemv.o src/dger.o \
+src/dgesl.o src/dglss.o src/dgmres.o src/dgtsl.o src/dh12.o src/dhels.o \
+src/dheqr.o src/dhfti.o src/dhkseq.o src/dhstrt.o src/dhvnrm.o src/dintp.o \
+src/dintrv.o src/dintyd.o src/dir.o src/djairy.o src/dlbeta.o src/dlgams.o \
+src/dli.o src/dllsia.o src/dllti2.o src/dlngam.o src/dlnrel.o src/dlpdoc.o \
+src/dlpdp.o src/dlsei.o src/dlsi.o src/dlsod.o src/dmout.o src/dmpar.o \
+src/dnbco.o src/dnbdi.o src/dnbfa.o src/dnbfs.o src/dnbsl.o src/dnls1.o \
+src/dnls1e.o src/dnrm2.o src/dnsq.o src/dnsqe.o src/dogleg.o src/domn.o \
+src/dorth.o src/dp1vlu.o src/dpbco.o src/dpbdi.o src/dpbfa.o src/dpbsl.o \
+src/dpchbs.o src/dpchce.o src/dpchci.o src/dpchcm.o src/dpchcs.o src/dpchdf.o \
+src/dpchfd.o src/dpchfe.o src/dpchia.o src/dpchic.o src/dpchid.o src/dpchim.o \
+src/dpchkt.o src/dpchng.o src/dpchsp.o src/dpchst.o src/dpchsw.o src/dpcoef.o \
+src/dpfqad.o src/dpigmr.o src/dpincw.o src/dpinit.o src/dpintm.o src/dpjac.o \
+src/dplint.o src/dplpce.o src/dplpdm.o src/dplpfe.o src/dplpfl.o src/dplpmn.o \
+src/dplpmu.o src/dplpup.o src/dpnnzr.o src/dpoch.o src/dpoch1.o src/dpoco.o \
+src/dpodi.o src/dpofa.o src/dpofs.o src/dpolcf.o src/dpolft.o src/dpolvl.o \
+src/dpopt.o src/dposl.o src/dppco.o src/dppdi.o src/dpperm.o src/dppfa.o \
+src/dppgq8.o src/dppqad.o src/dppsl.o src/dppval.o src/dprwpg.o src/dprwvr.o \
+src/dpsi.o src/dpsifn.o src/dpsixn.o src/dpsort.o src/dptsl.o src/dqag.o \
+src/dqage.o src/dqagi.o src/dqagie.o src/dqagp.o src/dqagpe.o src/dqags.o \
+src/dqagse.o src/dqawc.o src/dqawce.o src/dqawf.o src/dqawfe.o src/dqawo.o \
+src/dqawoe.o src/dqaws.o src/dqawse.o src/dqc25c.o src/dqc25f.o src/dqc25s.o \
+src/dqcheb.o src/dqdota.o src/dqdoti.o src/dqelg.o src/dqform.o src/dqk15.o \
+src/dqk15i.o src/dqk15w.o src/dqk21.o src/dqk31.o src/dqk41.o src/dqk51.o \
+src/dqk61.o src/dqmomo.o src/dqnc79.o src/dqng.o src/dqpsrt.o src/dqrdc.o \
+src/dqrfac.o src/dqrsl.o src/dqrslv.o src/dqwgtc.o src/dqwgtf.o src/dqwgts.o \
+src/drc.o src/drc3jj.o src/drc3jm.o src/drc6j.o src/drd.o src/dreadp.o \
+src/drf.o src/drj.o src/drkfs.o src/drlcal.o src/drot.o src/drotg.o \
+src/drotm.o src/drotmg.o src/drsco.o src/ds2lt.o src/ds2y.o src/dsbmv.o \
+src/dscal.o src/dsd2s.o src/dsdbcg.o src/dsdcg.o src/dsdcgn.o src/dsdcgs.o \
+src/dsdgmr.o src/dsdi.o src/dsdomn.o src/dsdot.o src/dsds.o src/dsdscl.o \
+src/dsgs.o src/dsiccg.o src/dsico.o src/dsics.o src/dsidi.o src/dsifa.o \
+src/dsilur.o src/dsilus.o src/dsindg.o src/dsisl.o src/dsjac.o src/dsli.o \
+src/dsli2.o src/dsllti.o src/dslubc.o src/dslucn.o src/dslucs.o src/dslugm.o \
+src/dslui.o src/dslui2.o src/dslui4.o src/dsluom.o src/dsluti.o src/dslvs.o \
+src/dsmmi2.o src/dsmmti.o src/dsmtv.o src/dsmv.o src/dsort.o src/dsos.o \
+src/dsoseq.o src/dsossl.o src/dspco.o src/dspdi.o src/dspenc.o src/dspfa.o \
+src/dsplp.o src/dspmv.o src/dspr.o src/dspr2.o src/dspsl.o src/dsteps.o \
+src/dstod.o src/dsvco.o src/dsvdc.o src/dswap.o src/dsymm.o src/dsymv.o \
+src/dsyr.o src/dsyr2.o src/dsyr2k.o src/dsyrk.o src/dtbmv.o src/dtbsv.o \
+src/dtin.o src/dtout.o src/dtpmv.o src/dtpsv.o src/dtrco.o src/dtrdi.o \
+src/dtrmm.o src/dtrmv.o src/dtrsl.o src/dtrsm.o src/dtrsv.o src/du11ls.o \
+src/du11us.o src/du12ls.o src/du12us.o src/dulsia.o src/dusrmt.o src/dvnrms.o \
+src/dvout.o src/dwnlit.o src/dwnlsm.o src/dwnlt1.o src/dwnlt2.o src/dwnlt3.o \
+src/dwnnls.o src/dwritp.o src/dwupdt.o src/dx.o src/dx4.o src/dxadd.o \
+src/dxadj.o src/dxc210.o src/dxcon.o src/dxlcal.o src/dxlegf.o src/dxnrmp.o \
+src/dxpmu.o src/dxpmup.o src/dxpnrm.o src/dxpqnu.o src/dxpsi.o src/dxqmu.o \
+src/dxqnu.o src/dxred.o src/dxset.o src/dy.o src/dy4.o src/dyairy.o src/e1.o \
+src/efc.o src/efcmn.o src/ei.o src/eisdoc.o src/elmbak.o src/elmhes.o \
+src/eltran.o src/enorm.o src/erf.o src/erfc.o src/exint.o src/exprel.o \
+src/ezfft1.o src/ezfftb.o src/ezfftf.o src/ezffti.o src/fac.o src/fc.o \
+src/fcmn.o src/fdjac1.o src/fdjac3.o src/fdump.o src/fftdoc.o src/figi.o \
+src/figi2.o src/fulmat.o src/fundoc.o src/fzero.o src/gami.o src/gamic.o \
+src/gamit.o src/gamlim.o src/gamln.o src/gamma.o src/gamr.o src/gamrn.o \
+src/gaus8.o src/genbun.o src/h12.o src/hfti.o src/hkseq.o src/hpperm.o \
+src/hpsort.o src/hqr.o src/hqr2.o src/hstart.o src/hstcrt.o src/hstcs1.o \
+src/hstcsp.o src/hstcyl.o src/hstplr.o src/hstssp.o src/htrib3.o src/htribk.o \
+src/htrid3.o src/htridi.o src/hvnrm.o src/hw3crt.o src/hwscrt.o src/hwscs1.o \
+src/hwscsp.o src/hwscyl.o src/hwsplr.o src/hwsss1.o src/hwsssp.o src/i1mach.o \
+src/i1merg.o src/icamax.o src/icopy.o src/idamax.o src/idloc.o src/imtql1.o \
+src/imtql2.o src/imtqlv.o src/indxa.o src/indxb.o src/indxc.o src/initds.o \
+src/inits.o src/intrv.o src/intyd.o src/invit.o src/inxca.o src/inxcb.o \
+src/inxcc.o src/iploc.o src/ipperm.o src/ipsort.o src/isamax.o src/isdbcg.o \
+src/isdcg.o src/isdcgn.o src/isdcgs.o src/isdgmr.o src/isdir.o src/isdomn.o \
+src/isort.o src/issbcg.o src/isscg.o src/isscgn.o src/isscgs.o src/issgmr.o \
+src/issir.o src/issomn.o src/iswap.o src/ivout.o src/j4save.o src/jairy.o \
+src/la05ad.o src/la05as.o src/la05bd.o src/la05bs.o src/la05cd.o src/la05cs.o \
+src/la05ed.o src/la05es.o src/llsia.o src/lmpar.o src/lpdp.o src/lsame.o \
+src/lsei.o src/lsi.o src/lsod.o src/mc20ad.o src/mc20as.o src/minfit.o \
+src/minso4.o src/minsol.o src/mpadd.o src/mpadd2.o src/mpadd3.o src/mpblas.o \
+src/mpcdm.o src/mpchk.o src/mpcmd.o src/mpdivi.o src/mperr.o src/mpmaxr.o \
+src/mpmlp.o src/mpmul.o src/mpmul2.o src/mpmuli.o src/mpnzr.o src/mpovfl.o \
+src/mpstr.o src/mpunfl.o src/numxer.o src/ortbak.o src/orthes.o src/ortho4.o \
+src/orthog.o src/ortran.o src/passb.o src/passb2.o src/passb3.o src/passb4.o \
+src/passb5.o src/passf.o src/passf2.o src/passf3.o src/passf4.o src/passf5.o \
+src/pchbs.o src/pchce.o src/pchci.o src/pchcm.o src/pchcs.o src/pchdf.o \
+src/pchdoc.o src/pchfd.o src/pchfe.o src/pchia.o src/pchic.o src/pchid.o \
+src/pchim.o src/pchkt.o src/pchngs.o src/pchsp.o src/pchst.o src/pchsw.o \
+src/pcoef.o src/pfqad.o src/pgsf.o src/pimach.o src/pinitm.o src/pjac.o \
+src/pnnzrs.o src/poch.o src/poch1.o src/pois3d.o src/poisd2.o src/poisn2.o \
+src/poisp2.o src/poistg.o src/polcof.o src/polfit.o src/polint.o src/polyvl.o \
+src/pos3d1.o src/postg2.o src/ppadd.o src/ppgq8.o src/ppgsf.o src/pppsf.o \
+src/ppqad.o src/ppsgf.o src/ppspf.o src/ppval.o src/proc.o src/procp.o \
+src/prod.o src/prodp.o src/prwpge.o src/prwvir.o src/psgf.o src/psi.o \
+src/psifn.o src/psixn.o src/pvalue.o src/pythag.o src/qag.o src/qage.o \
+src/qagi.o src/qagie.o src/qagp.o src/qagpe.o src/qags.o src/qagse.o \
+src/qawc.o src/qawce.o src/qawf.o src/qawfe.o src/qawo.o src/qawoe.o \
+src/qaws.o src/qawse.o src/qc25c.o src/qc25f.o src/qc25s.o src/qcheb.o \
+src/qelg.o src/qform.o src/qk15.o src/qk15i.o src/qk15w.o src/qk21.o \
+src/qk31.o src/qk41.o src/qk51.o src/qk61.o src/qmomo.o src/qnc79.o src/qng.o \
+src/qpdoc.o src/qpsrt.o src/qrfac.o src/qrsolv.o src/qs2i1d.o src/qs2i1r.o \
+src/qwgtc.o src/qwgtf.o src/qwgts.o src/qzhes.o src/qzit.o src/qzval.o \
+src/qzvec.o src/r1mach.o src/r1mpyq.o src/r1updt.o src/r9aimp.o src/r9atn1.o \
+src/r9chu.o src/r9gmic.o src/r9gmit.o src/r9knus.o src/r9lgic.o src/r9lgit.o \
+src/r9lgmc.o src/r9ln2r.o src/r9pak.o src/r9upak.o src/radb2.o src/radb3.o \
+src/radb4.o src/radb5.o src/radbg.o src/radf2.o src/radf3.o src/radf4.o \
+src/radf5.o src/radfg.o src/rand.o src/ratqr.o src/rc.o src/rc3jj.o \
+src/rc3jm.o src/rc6j.o src/rd.o src/rebak.o src/rebakb.o src/reduc.o \
+src/reduc2.o src/rf.o src/rfftb.o src/rfftb1.o src/rfftf.o src/rfftf1.o \
+src/rffti.o src/rffti1.o src/rg.o src/rgauss.o src/rgg.o src/rj.o \
+src/rpqr79.o src/rpzero.o src/rs.o src/rsb.o src/rsco.o src/rsg.o src/rsgab.o \
+src/rsgba.o src/rsp.o src/rst.o src/rt.o src/runif.o src/rwupdt.o \
+src/s1merg.o src/sasum.o src/saxpy.o src/sbcg.o src/sbhin.o src/sbocls.o \
+src/sbols.o src/sbolsm.o src/scasum.o src/scg.o src/scgn.o src/scgs.o \
+src/schdc.o src/schdd.o src/schex.o src/schkw.o src/schud.o src/sclosm.o \
+src/scnrm2.o src/scopy.o src/scopym.o src/scov.o src/scpplt.o src/sdaini.o \
+src/sdajac.o src/sdanrm.o src/sdaslv.o src/sdassl.o src/sdastp.o src/sdatrp.o \
+src/sdawts.o src/sdcor.o src/sdcst.o src/sdntl.o src/sdntp.o src/sdot.o \
+src/sdpsc.o src/sdpst.o src/sdriv1.o src/sdriv2.o src/sdriv3.o src/sdscl.o \
+src/sdsdot.o src/sdstp.o src/sdzro.o src/sepeli.o src/sepx4.o src/sgbco.o \
+src/sgbdi.o src/sgbfa.o src/sgbmv.o src/sgbsl.o src/sgeco.o src/sgedi.o \
+src/sgeev.o src/sgefa.o src/sgefs.o src/sgeir.o src/sgemm.o src/sgemv.o \
+src/sger.o src/sgesl.o src/sglss.o src/sgmres.o src/sgtsl.o src/shels.o \
+src/sheqr.o src/sindg.o src/sinqb.o src/sinqf.o src/sinqi.o src/sint.o \
+src/sinti.o src/sintrp.o src/sir.o src/sllti2.o src/slpdoc.o src/slvs.o \
+src/smout.o src/snbco.o src/snbdi.o src/snbfa.o src/snbfs.o src/snbir.o \
+src/snbsl.o src/snls1.o src/snls1e.o src/snrm2.o src/snsq.o src/snsqe.o \
+src/somn.o src/sopenm.o src/sorth.o src/sos.o src/soseqs.o src/sossol.o \
+src/spbco.o src/spbdi.o src/spbfa.o src/spbsl.o src/speli4.o src/spelip.o \
+src/spenc.o src/spigmr.o src/spincw.o src/spinit.o src/splp.o src/splpce.o \
+src/splpdm.o src/splpfe.o src/splpfl.o src/splpmn.o src/splpmu.o src/splpup.o \
+src/spoco.o src/spodi.o src/spofa.o src/spofs.o src/spoir.o src/spopt.o \
+src/sposl.o src/sppco.o src/sppdi.o src/spperm.o src/sppfa.o src/sppsl.o \
+src/spsort.o src/sptsl.o src/sqrdc.o src/sqrsl.o src/sreadp.o src/srlcal.o \
+src/srot.o src/srotg.o src/srotm.o src/srotmg.o src/ss2lt.o src/ss2y.o \
+src/ssbmv.o src/sscal.o src/ssd2s.o src/ssdbcg.o src/ssdcg.o src/ssdcgn.o \
+src/ssdcgs.o src/ssdgmr.o src/ssdi.o src/ssdomn.o src/ssds.o src/ssdscl.o \
+src/ssgs.o src/ssiccg.o src/ssico.o src/ssics.o src/ssidi.o src/ssiev.o \
+src/ssifa.o src/ssilur.o src/ssilus.o src/ssisl.o src/ssjac.o src/ssli.o \
+src/ssli2.o src/ssllti.o src/sslubc.o src/sslucn.o src/sslucs.o src/sslugm.o \
+src/sslui.o src/sslui2.o src/sslui4.o src/ssluom.o src/ssluti.o src/ssmmi2.o \
+src/ssmmti.o src/ssmtv.o src/ssmv.o src/ssort.o src/sspco.o src/sspdi.o \
+src/sspev.o src/sspfa.o src/sspmv.o src/sspr.o src/sspr2.o src/sspsl.o \
+src/ssvdc.o src/sswap.o src/ssymm.o src/ssymv.o src/ssyr.o src/ssyr2.o \
+src/ssyr2k.o src/ssyrk.o src/stbmv.o src/stbsv.o src/steps.o src/stin.o \
+src/stod.o src/stout.o src/stpmv.o src/stpsv.o src/strco.o src/strdi.o \
+src/strmm.o src/strmv.o src/strsl.o src/strsm.o src/strsv.o src/svco.o \
+src/svd.o src/svout.o src/swritp.o src/sxlcal.o src/tevlc.o src/tevls.o \
+src/tinvit.o src/tql1.o src/tql2.o src/tqlrat.o src/trbak1.o src/trbak3.o \
+src/tred1.o src/tred2.o src/tred3.o src/tri3.o src/tridib.o src/tridq.o \
+src/tris4.o src/trisp.o src/trix.o src/tsturm.o src/u11ls.o src/u11us.o \
+src/u12ls.o src/u12us.o src/ulsia.o src/usrmat.o src/vnwrms.o src/wnlit.o \
+src/wnlsm.o src/wnlt1.o src/wnlt2.o src/wnlt3.o src/wnnls.o src/xadd.o \
+src/xadj.o src/xc210.o src/xcon.o src/xerbla.o src/xerclr.o src/xercnt.o \
+src/xerdmp.o src/xerhlt.o src/xermax.o src/xermsg.o src/xerprn.o src/xersve.o \
+src/xgetf.o src/xgetua.o src/xgetun.o src/xlegf.o src/xnrmp.o src/xpmu.o \
+src/xpmup.o src/xpnrm.o src/xpqnu.o src/xpsi.o src/xqmu.o src/xqnu.o \
+src/xred.o src/xset.o src/xsetf.o src/xsetua.o src/xsetun.o src/yairy.o \
+src/zabs.o src/zacai.o src/zacon.o src/zairy.o src/zasyi.o src/zbesh.o \
+src/zbesi.o src/zbesj.o src/zbesk.o src/zbesy.o src/zbinu.o src/zbiry.o \
+src/zbknu.o src/zbuni.o src/zbunk.o src/zdiv.o src/zexp.o src/zkscl.o \
+src/zlog.o src/zmlri.o src/zmlt.o src/zrati.o src/zs1s2.o src/zseri.o \
+src/zshch.o src/zsqrt.o src/zuchk.o src/zunhj.o src/zuni1.o src/zuni2.o \
+src/zunik.o src/zunk1.o src/zunk2.o src/zuoik.o src/zwrsk.o
+
+objs_dbvp=src/dbksol.o src/dbvder.o src/dbvpor.o src/dbvsup.o src/dcoef.o \
+src/dcscal.o src/dexbvp.o src/dlssud.o src/dmacon.o src/dmgsbv.o src/dohtrl.o \
+src/dorthr.o src/dprvec.o src/dreort.o src/drkfab.o src/dstor1.o src/dstway.o \
+src/dsuds.o src/dvecs.o $(objs)
+
+objs_sbvp=src/bvder.o src/bvpor.o src/bvsup.o src/cscale.o src/exbvp.o \
+src/lssods.o src/lssuds.o src/macon.o src/mgsbv.o src/ohtrol.o src/ohtror.o \
+src/orthol.o src/orthor.o src/prvec.o src/reort.o src/rkfab.o src/scoef.o \
+src/sods.o src/stor1.o src/stway.o src/suds.o src/svecs.o $(objs)
+
+all: lib/libslatec.so lib/libslatec.a lib/libslatec-dbvp.a lib/libslatec-sbvp.a
+
+clean:
+ rm -fr lib src/*.o
+
+install:
+ install -d $(DESTDIR)$(PREFIX)/lib
+ install -m644 -t $(DESTDIR)$(PREFIX)/lib \
+ lib/libslatec.a lib/libslatec-dbvp.a lib/libslatec-sbvp.a
+ install -m755 -t $(DESTDIR)$(PREFIX)/lib lib/libslatec.so.$(version)
+ cp -P lib/libslatec.so lib/libslatec.so.$(major) \
+ $(DESTDIR)$(PREFIX)/lib
+
+uninstall:
+ rm -f \
+ $(DESTDIR)$(PREFIX)/lib/libslatec.a \
+ $(DESTDIR)$(PREFIX)/lib/libslatec-sbvp.a \
+ $(DESTDIR)$(PREFIX)/lib/libslatec-dbvp.a \
+ $(DESTDIR)$(PREFIX)/lib/libslatec.so \
+ $(DESTDIR)$(PREFIX)/lib/libslatec.so.$(major) \
+ $(DESTDIR)$(PREFIX)/lib/libslatec.so.$(version)
+
+.PHONY: all clean install uninstall
+
+lib/libslatec.a: $(objs)
+ mkdir -p lib
+ $(AR) $(ARFLAGS) $@ $(objs)
+
+lib/libslatec-dbvp.a: $(objs_dbvp)
+ mkdir -p lib
+ $(AR) $(ARFLAGS) $@ $(objs_dbvp)
+
+lib/libslatec-sbvp.a: $(objs_sbvp)
+ mkdir -p lib
+ $(AR) $(ARFLAGS) $@ $(objs_sbvp)
+
+lib/libslatec.so: lib/libslatec.so.$(major)
+ ln -fs libslatec.so.$(major) $@
+
+lib/libslatec.so.$(major): lib/libslatec.so.$(version)
+ ln -fs libslatec.so.$(version) $@
+
+lib/libslatec.so.$(version): $(objs)
+ mkdir -p lib
+ $(FC) $(LDFLAGS) -shared -Wl,-soname,libslatec.so.$(major) \
+ -o $@ $(objs)
diff --git a/PKGBUILD b/PKGBUILD
new file mode 100644
index 000000000000..2fdcf5b950c6
--- /dev/null
+++ b/PKGBUILD
@@ -0,0 +1,33 @@
+# Maintainer: Phil Ruffwind <rf@rufflewind.com>
+pkgname=slatec
+pkgver=4.1
+pkgrel=2
+pkgdesc="SLATEC Common Mathematical Library"
+arch=(i686 x86_64)
+url=http://netlib.org/slatec
+license=(PD)
+depends=()
+makedepends=(gcc-fortran)
+source=(http://netlib.org/slatec/slatec_src.tgz LICENSE Makefile slatec.patch)
+sha256sums=('df009d9ef9c18aae06ce68711b1ae108d3533b4f174582c3cbea0915c4fdfe01'
+ 'd83a2ee63ee4bbd56d59072a0778e06685fa001047f153a3f0293d30849c8046'
+ '7d8fb187324596d4d882d895c31db0afddf2984a933ab6e623de576d5f36558e'
+ '4fb39658c6da25618d20e6aec0680c48d4666ab0eb1322f6e7f93a34bb51ce4f')
+
+prepare() {
+ cd "$srcdir"
+ patch -p 0 <slatec.patch
+}
+
+build() {
+ cd "$srcdir"
+ make FC=gfortran all
+}
+
+package() {
+ cd "$srcdir"
+ make DESTDIR="$pkgdir" PREFIX=/usr install
+ install -Dm644 LICENSE "$pkgdir/usr/share/licenses/$pkgname/LICENSE"
+ install -m755 -d "$pkgdir/usr/share/slatec/src"
+ install -m755 -t "$pkgdir/usr/share/slatec/src" src/*.f
+}
diff --git a/slatec.patch b/slatec.patch
new file mode 100644
index 000000000000..d3de2553fd1e
--- /dev/null
+++ b/slatec.patch
@@ -0,0 +1,2188 @@
+--- src/d1mach.f
++++ src/d1mach.f
+@@ -67,436 +67,189 @@
+ C 920625 Added CONVEX -p8 and -pd8 compiler option constants.
+ C (BKS, WRB)
+ C 930201 Added DEC Alpha and SGI constants. (RWC and WRB)
++C 150501 Replaced with alternative implementation from BLAS.
+ C***END PROLOGUE D1MACH
+ C
++ INTEGER SMALL(2)
++ INTEGER LARGE(2)
++ INTEGER RIGHT(2)
++ INTEGER DIVER(2)
++ INTEGER LOG10(2)
++ INTEGER SC, CRAY1(38), J
++ COMMON /D9MACH/ CRAY1
++ SAVE SMALL, LARGE, RIGHT, DIVER, LOG10, SC
+- INTEGER SMALL(4)
+- INTEGER LARGE(4)
+- INTEGER RIGHT(4)
+- INTEGER DIVER(4)
+- INTEGER LOG10(4)
+-C
+ DOUBLE PRECISION DMACH(5)
+- SAVE DMACH
+-C
+ EQUIVALENCE (DMACH(1),SMALL(1))
+ EQUIVALENCE (DMACH(2),LARGE(1))
+ EQUIVALENCE (DMACH(3),RIGHT(1))
+ EQUIVALENCE (DMACH(4),DIVER(1))
+ EQUIVALENCE (DMACH(5),LOG10(1))
++C THIS VERSION ADAPTS AUTOMATICALLY TO MOST CURRENT MACHINES.
++C R1MACH CAN HANDLE AUTO-DOUBLE COMPILING, BUT THIS VERSION OF
++C D1MACH DOES NOT, BECAUSE WE DO NOT HAVE QUAD CONSTANTS FOR
++C MANY MACHINES YET.
++C TO COMPILE ON OLDER MACHINES, ADD A C IN COLUMN 1
++C ON THE NEXT LINE
++ DATA SC/0/
++C AND REMOVE THE C FROM COLUMN 1 IN ONE OF THE SECTIONS BELOW.
++C CONSTANTS FOR EVEN OLDER MACHINES CAN BE OBTAINED BY
++C mail netlib@research.bell-labs.com
++C send old1mach from blas
++C PLEASE SEND CORRECTIONS TO dmg OR ehg@bell-labs.com.
++C
++C MACHINE CONSTANTS FOR THE HONEYWELL DPS 8/70 SERIES.
++C DATA SMALL(1),SMALL(2) / O402400000000, O000000000000 /
++C DATA LARGE(1),LARGE(2) / O376777777777, O777777777777 /
++C DATA RIGHT(1),RIGHT(2) / O604400000000, O000000000000 /
++C DATA DIVER(1),DIVER(2) / O606400000000, O000000000000 /
++C DATA LOG10(1),LOG10(2) / O776464202324, O117571775714 /, SC/987/
++C
++C MACHINE CONSTANTS FOR PDP-11 FORTRANS SUPPORTING
++C 32-BIT INTEGERS.
++C DATA SMALL(1),SMALL(2) / 8388608, 0 /
++C DATA LARGE(1),LARGE(2) / 2147483647, -1 /
++C DATA RIGHT(1),RIGHT(2) / 612368384, 0 /
++C DATA DIVER(1),DIVER(2) / 620756992, 0 /
++C DATA LOG10(1),LOG10(2) / 1067065498, -2063872008 /, SC/987/
++C
++C MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES.
++C DATA SMALL(1),SMALL(2) / O000040000000, O000000000000 /
++C DATA LARGE(1),LARGE(2) / O377777777777, O777777777777 /
++C DATA RIGHT(1),RIGHT(2) / O170540000000, O000000000000 /
++C DATA DIVER(1),DIVER(2) / O170640000000, O000000000000 /
++C DATA LOG10(1),LOG10(2) / O177746420232, O411757177572 /, SC/987/
++C
++C ON FIRST CALL, IF NO DATA UNCOMMENTED, TEST MACHINE TYPES.
++ IF (SC .NE. 987) THEN
++ DMACH(1) = 1.D13
++ IF ( SMALL(1) .EQ. 1117925532
++ * .AND. SMALL(2) .EQ. -448790528) THEN
++* *** IEEE BIG ENDIAN ***
++ SMALL(1) = 1048576
++ SMALL(2) = 0
++ LARGE(1) = 2146435071
++ LARGE(2) = -1
++ RIGHT(1) = 1017118720
++ RIGHT(2) = 0
++ DIVER(1) = 1018167296
++ DIVER(2) = 0
++ LOG10(1) = 1070810131
++ LOG10(2) = 1352628735
++ ELSE IF ( SMALL(2) .EQ. 1117925532
++ * .AND. SMALL(1) .EQ. -448790528) THEN
++* *** IEEE LITTLE ENDIAN ***
++ SMALL(2) = 1048576
++ SMALL(1) = 0
++ LARGE(2) = 2146435071
++ LARGE(1) = -1
++ RIGHT(2) = 1017118720
++ RIGHT(1) = 0
++ DIVER(2) = 1018167296
++ DIVER(1) = 0
++ LOG10(2) = 1070810131
++ LOG10(1) = 1352628735
++ ELSE IF ( SMALL(1) .EQ. -2065213935
++ * .AND. SMALL(2) .EQ. 10752) THEN
++* *** VAX WITH D_FLOATING ***
++ SMALL(1) = 128
++ SMALL(2) = 0
++ LARGE(1) = -32769
++ LARGE(2) = -1
++ RIGHT(1) = 9344
++ RIGHT(2) = 0
++ DIVER(1) = 9472
++ DIVER(2) = 0
++ LOG10(1) = 546979738
++ LOG10(2) = -805796613
++ ELSE IF ( SMALL(1) .EQ. 1267827943
++ * .AND. SMALL(2) .EQ. 704643072) THEN
++* *** IBM MAINFRAME ***
++ SMALL(1) = 1048576
++ SMALL(2) = 0
++ LARGE(1) = 2147483647
++ LARGE(2) = -1
++ RIGHT(1) = 856686592
++ RIGHT(2) = 0
++ DIVER(1) = 873463808
++ DIVER(2) = 0
++ LOG10(1) = 1091781651
++ LOG10(2) = 1352628735
++ ELSE IF ( SMALL(1) .EQ. 1120022684
++ * .AND. SMALL(2) .EQ. -448790528) THEN
++* *** CONVEX C-1 ***
++ SMALL(1) = 1048576
++ SMALL(2) = 0
++ LARGE(1) = 2147483647
++ LARGE(2) = -1
++ RIGHT(1) = 1019215872
++ RIGHT(2) = 0
++ DIVER(1) = 1020264448
++ DIVER(2) = 0
++ LOG10(1) = 1072907283
++ LOG10(2) = 1352628735
++ ELSE IF ( SMALL(1) .EQ. 815547074
++ * .AND. SMALL(2) .EQ. 58688) THEN
++* *** VAX G-FLOATING ***
++ SMALL(1) = 16
++ SMALL(2) = 0
++ LARGE(1) = -32769
++ LARGE(2) = -1
++ RIGHT(1) = 15552
++ RIGHT(2) = 0
++ DIVER(1) = 15568
++ DIVER(2) = 0
++ LOG10(1) = 1142112243
++ LOG10(2) = 2046775455
++ ELSE
++ DMACH(2) = 1.D27 + 1
++ DMACH(3) = 1.D27
++ LARGE(2) = LARGE(2) - RIGHT(2)
++ IF (LARGE(2) .EQ. 64 .AND. SMALL(2) .EQ. 0) THEN
++ CRAY1(1) = 67291416
++ DO 10 J = 1, 20
++ CRAY1(J+1) = CRAY1(J) + CRAY1(J)
++ 10 CONTINUE
++ CRAY1(22) = CRAY1(21) + 321322
++ DO 20 J = 22, 37
++ CRAY1(J+1) = CRAY1(J) + CRAY1(J)
++ 20 CONTINUE
++ IF (CRAY1(38) .EQ. SMALL(1)) THEN
++* *** CRAY ***
++ CALL I1MCRY(SMALL(1), J, 8285, 8388608, 0)
++ SMALL(2) = 0
++ CALL I1MCRY(LARGE(1), J, 24574, 16777215, 16777215)
++ CALL I1MCRY(LARGE(2), J, 0, 16777215, 16777214)
++ CALL I1MCRY(RIGHT(1), J, 16291, 8388608, 0)
++ RIGHT(2) = 0
++ CALL I1MCRY(DIVER(1), J, 16292, 8388608, 0)
++ DIVER(2) = 0
++ CALL I1MCRY(LOG10(1), J, 16383, 10100890, 8715215)
++ CALL I1MCRY(LOG10(2), J, 0, 16226447, 9001388)
++ ELSE
++ WRITE(*,9000)
++ STOP 779
++ END IF
++ ELSE
++ WRITE(*,9000)
++ STOP 779
++ END IF
++ END IF
++ SC = 987
++ END IF
++* SANITY CHECK
++ IF (DMACH(4) .GE. 1.0D0) STOP 778
++ IF (I .LT. 1 .OR. I .GT. 5) THEN
++ WRITE(*,*) 'D1MACH(I): I =',I,' is out of bounds.'
++ STOP
++ END IF
+-C
+-C MACHINE CONSTANTS FOR THE AMIGA
+-C ABSOFT FORTRAN COMPILER USING THE 68020/68881 COMPILER OPTION
+-C
+-C DATA SMALL(1), SMALL(2) / Z'00100000', Z'00000000' /
+-C DATA LARGE(1), LARGE(2) / Z'7FEFFFFF', Z'FFFFFFFF' /
+-C DATA RIGHT(1), RIGHT(2) / Z'3CA00000', Z'00000000' /
+-C DATA DIVER(1), DIVER(2) / Z'3CB00000', Z'00000000' /
+-C DATA LOG10(1), LOG10(2) / Z'3FD34413', Z'509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE AMIGA
+-C ABSOFT FORTRAN COMPILER USING SOFTWARE FLOATING POINT
+-C
+-C DATA SMALL(1), SMALL(2) / Z'00100000', Z'00000000' /
+-C DATA LARGE(1), LARGE(2) / Z'7FDFFFFF', Z'FFFFFFFF' /
+-C DATA RIGHT(1), RIGHT(2) / Z'3CA00000', Z'00000000' /
+-C DATA DIVER(1), DIVER(2) / Z'3CB00000', Z'00000000' /
+-C DATA LOG10(1), LOG10(2) / Z'3FD34413', Z'509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE APOLLO
+-C
+-C DATA SMALL(1), SMALL(2) / 16#00100000, 16#00000000 /
+-C DATA LARGE(1), LARGE(2) / 16#7FFFFFFF, 16#FFFFFFFF /
+-C DATA RIGHT(1), RIGHT(2) / 16#3CA00000, 16#00000000 /
+-C DATA DIVER(1), DIVER(2) / 16#3CB00000, 16#00000000 /
+-C DATA LOG10(1), LOG10(2) / 16#3FD34413, 16#509F79FF /
+-C
+-C MACHINE CONSTANTS FOR THE BURROUGHS 1700 SYSTEM
+-C
+-C DATA SMALL(1) / ZC00800000 /
+-C DATA SMALL(2) / Z000000000 /
+-C DATA LARGE(1) / ZDFFFFFFFF /
+-C DATA LARGE(2) / ZFFFFFFFFF /
+-C DATA RIGHT(1) / ZCC5800000 /
+-C DATA RIGHT(2) / Z000000000 /
+-C DATA DIVER(1) / ZCC6800000 /
+-C DATA DIVER(2) / Z000000000 /
+-C DATA LOG10(1) / ZD00E730E7 /
+-C DATA LOG10(2) / ZC77800DC0 /
+-C
+-C MACHINE CONSTANTS FOR THE BURROUGHS 5700 SYSTEM
+-C
+-C DATA SMALL(1) / O1771000000000000 /
+-C DATA SMALL(2) / O0000000000000000 /
+-C DATA LARGE(1) / O0777777777777777 /
+-C DATA LARGE(2) / O0007777777777777 /
+-C DATA RIGHT(1) / O1461000000000000 /
+-C DATA RIGHT(2) / O0000000000000000 /
+-C DATA DIVER(1) / O1451000000000000 /
+-C DATA DIVER(2) / O0000000000000000 /
+-C DATA LOG10(1) / O1157163034761674 /
+-C DATA LOG10(2) / O0006677466732724 /
+-C
+-C MACHINE CONSTANTS FOR THE BURROUGHS 6700/7700 SYSTEMS
+-C
+-C DATA SMALL(1) / O1771000000000000 /
+-C DATA SMALL(2) / O7770000000000000 /
+-C DATA LARGE(1) / O0777777777777777 /
+-C DATA LARGE(2) / O7777777777777777 /
+-C DATA RIGHT(1) / O1461000000000000 /
+-C DATA RIGHT(2) / O0000000000000000 /
+-C DATA DIVER(1) / O1451000000000000 /
+-C DATA DIVER(2) / O0000000000000000 /
+-C DATA LOG10(1) / O1157163034761674 /
+-C DATA LOG10(2) / O0006677466732724 /
+-C
+-C MACHINE CONSTANTS FOR THE CDC 170/180 SERIES USING NOS/VE
+-C
+-C DATA SMALL(1) / Z"3001800000000000" /
+-C DATA SMALL(2) / Z"3001000000000000" /
+-C DATA LARGE(1) / Z"4FFEFFFFFFFFFFFE" /
+-C DATA LARGE(2) / Z"4FFE000000000000" /
+-C DATA RIGHT(1) / Z"3FD2800000000000" /
+-C DATA RIGHT(2) / Z"3FD2000000000000" /
+-C DATA DIVER(1) / Z"3FD3800000000000" /
+-C DATA DIVER(2) / Z"3FD3000000000000" /
+-C DATA LOG10(1) / Z"3FFF9A209A84FBCF" /
+-C DATA LOG10(2) / Z"3FFFF7988F8959AC" /
+-C
+-C MACHINE CONSTANTS FOR THE CDC 6000/7000 SERIES
+-C
+-C DATA SMALL(1) / 00564000000000000000B /
+-C DATA SMALL(2) / 00000000000000000000B /
+-C DATA LARGE(1) / 37757777777777777777B /
+-C DATA LARGE(2) / 37157777777777777777B /
+-C DATA RIGHT(1) / 15624000000000000000B /
+-C DATA RIGHT(2) / 00000000000000000000B /
+-C DATA DIVER(1) / 15634000000000000000B /
+-C DATA DIVER(2) / 00000000000000000000B /
+-C DATA LOG10(1) / 17164642023241175717B /
+-C DATA LOG10(2) / 16367571421742254654B /
+-C
+-C MACHINE CONSTANTS FOR THE CELERITY C1260
+-C
+-C DATA SMALL(1), SMALL(2) / Z'00100000', Z'00000000' /
+-C DATA LARGE(1), LARGE(2) / Z'7FEFFFFF', Z'FFFFFFFF' /
+-C DATA RIGHT(1), RIGHT(2) / Z'3CA00000', Z'00000000' /
+-C DATA DIVER(1), DIVER(2) / Z'3CB00000', Z'00000000' /
+-C DATA LOG10(1), LOG10(2) / Z'3FD34413', Z'509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -fn OR -pd8 COMPILER OPTION
+-C
+-C DATA DMACH(1) / Z'0010000000000000' /
+-C DATA DMACH(2) / Z'7FFFFFFFFFFFFFFF' /
+-C DATA DMACH(3) / Z'3CC0000000000000' /
+-C DATA DMACH(4) / Z'3CD0000000000000' /
+-C DATA DMACH(5) / Z'3FF34413509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -fi COMPILER OPTION
+-C
+-C DATA DMACH(1) / Z'0010000000000000' /
+-C DATA DMACH(2) / Z'7FEFFFFFFFFFFFFF' /
+-C DATA DMACH(3) / Z'3CA0000000000000' /
+-C DATA DMACH(4) / Z'3CB0000000000000' /
+-C DATA DMACH(5) / Z'3FD34413509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -p8 COMPILER OPTION
+-C
+-C DATA DMACH(1) / Z'00010000000000000000000000000000' /
+-C DATA DMACH(2) / Z'7FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF' /
+-C DATA DMACH(3) / Z'3F900000000000000000000000000000' /
+-C DATA DMACH(4) / Z'3F910000000000000000000000000000' /
+-C DATA DMACH(5) / Z'3FFF34413509F79FEF311F12B35816F9' /
+-C
+-C MACHINE CONSTANTS FOR THE CRAY
+-C
+-C DATA SMALL(1) / 201354000000000000000B /
+-C DATA SMALL(2) / 000000000000000000000B /
+-C DATA LARGE(1) / 577767777777777777777B /
+-C DATA LARGE(2) / 000007777777777777774B /
+-C DATA RIGHT(1) / 376434000000000000000B /
+-C DATA RIGHT(2) / 000000000000000000000B /
+-C DATA DIVER(1) / 376444000000000000000B /
+-C DATA DIVER(2) / 000000000000000000000B /
+-C DATA LOG10(1) / 377774642023241175717B /
+-C DATA LOG10(2) / 000007571421742254654B /
+-C
+-C MACHINE CONSTANTS FOR THE DATA GENERAL ECLIPSE S/200
+-C NOTE - IT MAY BE APPROPRIATE TO INCLUDE THE FOLLOWING CARD -
+-C STATIC DMACH(5)
+-C
+-C DATA SMALL / 20K, 3*0 /
+-C DATA LARGE / 77777K, 3*177777K /
+-C DATA RIGHT / 31420K, 3*0 /
+-C DATA DIVER / 32020K, 3*0 /
+-C DATA LOG10 / 40423K, 42023K, 50237K, 74776K /
+-C
+-C MACHINE CONSTANTS FOR THE DEC ALPHA
+-C USING G_FLOAT
+-C
+-C DATA DMACH(1) / '0000000000000010'X /
+-C DATA DMACH(2) / 'FFFFFFFFFFFF7FFF'X /
+-C DATA DMACH(3) / '0000000000003CC0'X /
+-C DATA DMACH(4) / '0000000000003CD0'X /
+-C DATA DMACH(5) / '79FF509F44133FF3'X /
+-C
+-C MACHINE CONSTANTS FOR THE DEC ALPHA
+-C USING IEEE_FORMAT
+-C
+-C DATA DMACH(1) / '0010000000000000'X /
+-C DATA DMACH(2) / '7FEFFFFFFFFFFFFF'X /
+-C DATA DMACH(3) / '3CA0000000000000'X /
+-C DATA DMACH(4) / '3CB0000000000000'X /
+-C DATA DMACH(5) / '3FD34413509F79FF'X /
+-C
+-C MACHINE CONSTANTS FOR THE DEC RISC
+-C
+-C DATA SMALL(1), SMALL(2) / Z'00000000', Z'00100000'/
+-C DATA LARGE(1), LARGE(2) / Z'FFFFFFFF', Z'7FEFFFFF'/
+-C DATA RIGHT(1), RIGHT(2) / Z'00000000', Z'3CA00000'/
+-C DATA DIVER(1), DIVER(2) / Z'00000000', Z'3CB00000'/
+-C DATA LOG10(1), LOG10(2) / Z'509F79FF', Z'3FD34413'/
+-C
+-C MACHINE CONSTANTS FOR THE DEC VAX
+-C USING D_FLOATING
+-C (EXPRESSED IN INTEGER AND HEXADECIMAL)
+-C THE HEX FORMAT BELOW MAY NOT BE SUITABLE FOR UNIX SYSTEMS
+-C THE INTEGER FORMAT SHOULD BE OK FOR UNIX SYSTEMS
+-C
+-C DATA SMALL(1), SMALL(2) / 128, 0 /
+-C DATA LARGE(1), LARGE(2) / -32769, -1 /
+-C DATA RIGHT(1), RIGHT(2) / 9344, 0 /
+-C DATA DIVER(1), DIVER(2) / 9472, 0 /
+-C DATA LOG10(1), LOG10(2) / 546979738, -805796613 /
+-C
+-C DATA SMALL(1), SMALL(2) / Z00000080, Z00000000 /
+-C DATA LARGE(1), LARGE(2) / ZFFFF7FFF, ZFFFFFFFF /
+-C DATA RIGHT(1), RIGHT(2) / Z00002480, Z00000000 /
+-C DATA DIVER(1), DIVER(2) / Z00002500, Z00000000 /
+-C DATA LOG10(1), LOG10(2) / Z209A3F9A, ZCFF884FB /
+-C
+-C MACHINE CONSTANTS FOR THE DEC VAX
+-C USING G_FLOATING
+-C (EXPRESSED IN INTEGER AND HEXADECIMAL)
+-C THE HEX FORMAT BELOW MAY NOT BE SUITABLE FOR UNIX SYSTEMS
+-C THE INTEGER FORMAT SHOULD BE OK FOR UNIX SYSTEMS
+-C
+-C DATA SMALL(1), SMALL(2) / 16, 0 /
+-C DATA LARGE(1), LARGE(2) / -32769, -1 /
+-C DATA RIGHT(1), RIGHT(2) / 15552, 0 /
+-C DATA DIVER(1), DIVER(2) / 15568, 0 /
+-C DATA LOG10(1), LOG10(2) / 1142112243, 2046775455 /
+-C
+-C DATA SMALL(1), SMALL(2) / Z00000010, Z00000000 /
+-C DATA LARGE(1), LARGE(2) / ZFFFF7FFF, ZFFFFFFFF /
+-C DATA RIGHT(1), RIGHT(2) / Z00003CC0, Z00000000 /
+-C DATA DIVER(1), DIVER(2) / Z00003CD0, Z00000000 /
+-C DATA LOG10(1), LOG10(2) / Z44133FF3, Z79FF509F /
+-C
+-C MACHINE CONSTANTS FOR THE ELXSI 6400
+-C (ASSUMING REAL*8 IS THE DEFAULT DOUBLE PRECISION)
+-C
+-C DATA SMALL(1), SMALL(2) / '00100000'X,'00000000'X /
+-C DATA LARGE(1), LARGE(2) / '7FEFFFFF'X,'FFFFFFFF'X /
+-C DATA RIGHT(1), RIGHT(2) / '3CB00000'X,'00000000'X /
+-C DATA DIVER(1), DIVER(2) / '3CC00000'X,'00000000'X /
+-C DATA LOG10(1), LOG10(2) / '3FD34413'X,'509F79FF'X /
+-C
+-C MACHINE CONSTANTS FOR THE HARRIS 220
+-C
+-C DATA SMALL(1), SMALL(2) / '20000000, '00000201 /
+-C DATA LARGE(1), LARGE(2) / '37777777, '37777577 /
+-C DATA RIGHT(1), RIGHT(2) / '20000000, '00000333 /
+-C DATA DIVER(1), DIVER(2) / '20000000, '00000334 /
+-C DATA LOG10(1), LOG10(2) / '23210115, '10237777 /
+-C
+-C MACHINE CONSTANTS FOR THE HONEYWELL 600/6000 SERIES
+-C
+-C DATA SMALL(1), SMALL(2) / O402400000000, O000000000000 /
+-C DATA LARGE(1), LARGE(2) / O376777777777, O777777777777 /
+-C DATA RIGHT(1), RIGHT(2) / O604400000000, O000000000000 /
+-C DATA DIVER(1), DIVER(2) / O606400000000, O000000000000 /
+-C DATA LOG10(1), LOG10(2) / O776464202324, O117571775714 /
+-C
+-C MACHINE CONSTANTS FOR THE HP 730
+-C
+-C DATA DMACH(1) / Z'0010000000000000' /
+-C DATA DMACH(2) / Z'7FEFFFFFFFFFFFFF' /
+-C DATA DMACH(3) / Z'3CA0000000000000' /
+-C DATA DMACH(4) / Z'3CB0000000000000' /
+-C DATA DMACH(5) / Z'3FD34413509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE HP 2100
+-C THREE WORD DOUBLE PRECISION OPTION WITH FTN4
+-C
+-C DATA SMALL(1), SMALL(2), SMALL(3) / 40000B, 0, 1 /
+-C DATA LARGE(1), LARGE(2), LARGE(3) / 77777B, 177777B, 177776B /
+-C DATA RIGHT(1), RIGHT(2), RIGHT(3) / 40000B, 0, 265B /
+-C DATA DIVER(1), DIVER(2), DIVER(3) / 40000B, 0, 276B /
+-C DATA LOG10(1), LOG10(2), LOG10(3) / 46420B, 46502B, 77777B /
+-C
+-C MACHINE CONSTANTS FOR THE HP 2100
+-C FOUR WORD DOUBLE PRECISION OPTION WITH FTN4
+-C
+-C DATA SMALL(1), SMALL(2) / 40000B, 0 /
+-C DATA SMALL(3), SMALL(4) / 0, 1 /
+-C DATA LARGE(1), LARGE(2) / 77777B, 177777B /
+-C DATA LARGE(3), LARGE(4) / 177777B, 177776B /
+-C DATA RIGHT(1), RIGHT(2) / 40000B, 0 /
+-C DATA RIGHT(3), RIGHT(4) / 0, 225B /
+-C DATA DIVER(1), DIVER(2) / 40000B, 0 /
+-C DATA DIVER(3), DIVER(4) / 0, 227B /
+-C DATA LOG10(1), LOG10(2) / 46420B, 46502B /
+-C DATA LOG10(3), LOG10(4) / 76747B, 176377B /
+-C
+-C MACHINE CONSTANTS FOR THE HP 9000
+-C
+-C DATA SMALL(1), SMALL(2) / 00040000000B, 00000000000B /
+-C DATA LARGE(1), LARGE(2) / 17737777777B, 37777777777B /
+-C DATA RIGHT(1), RIGHT(2) / 07454000000B, 00000000000B /
+-C DATA DIVER(1), DIVER(2) / 07460000000B, 00000000000B /
+-C DATA LOG10(1), LOG10(2) / 07764642023B, 12047674777B /
+-C
+-C MACHINE CONSTANTS FOR THE IBM 360/370 SERIES,
+-C THE XEROX SIGMA 5/7/9, THE SEL SYSTEMS 85/86, AND
+-C THE PERKIN ELMER (INTERDATA) 7/32.
+-C
+-C DATA SMALL(1), SMALL(2) / Z00100000, Z00000000 /
+-C DATA LARGE(1), LARGE(2) / Z7FFFFFFF, ZFFFFFFFF /
+-C DATA RIGHT(1), RIGHT(2) / Z33100000, Z00000000 /
+-C DATA DIVER(1), DIVER(2) / Z34100000, Z00000000 /
+-C DATA LOG10(1), LOG10(2) / Z41134413, Z509F79FF /
+-C
+-C MACHINE CONSTANTS FOR THE IBM PC
+-C ASSUMES THAT ALL ARITHMETIC IS DONE IN DOUBLE PRECISION
+-C ON 8088, I.E., NOT IN 80 BIT FORM FOR THE 8087.
+-C
+-C DATA SMALL(1) / 2.23D-308 /
+-C DATA LARGE(1) / 1.79D+308 /
+-C DATA RIGHT(1) / 1.11D-16 /
+-C DATA DIVER(1) / 2.22D-16 /
+-C DATA LOG10(1) / 0.301029995663981195D0 /
+-C
+-C MACHINE CONSTANTS FOR THE IBM RS 6000
+-C
+-C DATA DMACH(1) / Z'0010000000000000' /
+-C DATA DMACH(2) / Z'7FEFFFFFFFFFFFFF' /
+-C DATA DMACH(3) / Z'3CA0000000000000' /
+-C DATA DMACH(4) / Z'3CB0000000000000' /
+-C DATA DMACH(5) / Z'3FD34413509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE INTEL i860
+-C
+-C DATA DMACH(1) / Z'0010000000000000' /
+-C DATA DMACH(2) / Z'7FEFFFFFFFFFFFFF' /
+-C DATA DMACH(3) / Z'3CA0000000000000' /
+-C DATA DMACH(4) / Z'3CB0000000000000' /
+-C DATA DMACH(5) / Z'3FD34413509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE PDP-10 (KA PROCESSOR)
+-C
+-C DATA SMALL(1), SMALL(2) / "033400000000, "000000000000 /
+-C DATA LARGE(1), LARGE(2) / "377777777777, "344777777777 /
+-C DATA RIGHT(1), RIGHT(2) / "113400000000, "000000000000 /
+-C DATA DIVER(1), DIVER(2) / "114400000000, "000000000000 /
+-C DATA LOG10(1), LOG10(2) / "177464202324, "144117571776 /
+-C
+-C MACHINE CONSTANTS FOR THE PDP-10 (KI PROCESSOR)
+-C
+-C DATA SMALL(1), SMALL(2) / "000400000000, "000000000000 /
+-C DATA LARGE(1), LARGE(2) / "377777777777, "377777777777 /
+-C DATA RIGHT(1), RIGHT(2) / "103400000000, "000000000000 /
+-C DATA DIVER(1), DIVER(2) / "104400000000, "000000000000 /
+-C DATA LOG10(1), LOG10(2) / "177464202324, "476747767461 /
+-C
+-C MACHINE CONSTANTS FOR PDP-11 FORTRAN SUPPORTING
+-C 32-BIT INTEGERS (EXPRESSED IN INTEGER AND OCTAL).
+-C
+-C DATA SMALL(1), SMALL(2) / 8388608, 0 /
+-C DATA LARGE(1), LARGE(2) / 2147483647, -1 /
+-C DATA RIGHT(1), RIGHT(2) / 612368384, 0 /
+-C DATA DIVER(1), DIVER(2) / 620756992, 0 /
+-C DATA LOG10(1), LOG10(2) / 1067065498, -2063872008 /
+-C
+-C DATA SMALL(1), SMALL(2) / O00040000000, O00000000000 /
+-C DATA LARGE(1), LARGE(2) / O17777777777, O37777777777 /
+-C DATA RIGHT(1), RIGHT(2) / O04440000000, O00000000000 /
+-C DATA DIVER(1), DIVER(2) / O04500000000, O00000000000 /
+-C DATA LOG10(1), LOG10(2) / O07746420232, O20476747770 /
+-C
+-C MACHINE CONSTANTS FOR PDP-11 FORTRAN SUPPORTING
+-C 16-BIT INTEGERS (EXPRESSED IN INTEGER AND OCTAL).
+-C
+-C DATA SMALL(1), SMALL(2) / 128, 0 /
+-C DATA SMALL(3), SMALL(4) / 0, 0 /
+-C DATA LARGE(1), LARGE(2) / 32767, -1 /
+-C DATA LARGE(3), LARGE(4) / -1, -1 /
+-C DATA RIGHT(1), RIGHT(2) / 9344, 0 /
+-C DATA RIGHT(3), RIGHT(4) / 0, 0 /
+-C DATA DIVER(1), DIVER(2) / 9472, 0 /
+-C DATA DIVER(3), DIVER(4) / 0, 0 /
+-C DATA LOG10(1), LOG10(2) / 16282, 8346 /
+-C DATA LOG10(3), LOG10(4) / -31493, -12296 /
+-C
+-C DATA SMALL(1), SMALL(2) / O000200, O000000 /
+-C DATA SMALL(3), SMALL(4) / O000000, O000000 /
+-C DATA LARGE(1), LARGE(2) / O077777, O177777 /
+-C DATA LARGE(3), LARGE(4) / O177777, O177777 /
+-C DATA RIGHT(1), RIGHT(2) / O022200, O000000 /
+-C DATA RIGHT(3), RIGHT(4) / O000000, O000000 /
+-C DATA DIVER(1), DIVER(2) / O022400, O000000 /
+-C DATA DIVER(3), DIVER(4) / O000000, O000000 /
+-C DATA LOG10(1), LOG10(2) / O037632, O020232 /
+-C DATA LOG10(3), LOG10(4) / O102373, O147770 /
+-C
+-C MACHINE CONSTANTS FOR THE SILICON GRAPHICS
+-C
+-C DATA SMALL(1), SMALL(2) / Z'00100000', Z'00000000' /
+-C DATA LARGE(1), LARGE(2) / Z'7FEFFFFF', Z'FFFFFFFF' /
+-C DATA RIGHT(1), RIGHT(2) / Z'3CA00000', Z'00000000' /
+-C DATA DIVER(1), DIVER(2) / Z'3CB00000', Z'00000000' /
+-C DATA LOG10(1), LOG10(2) / Z'3FD34413', Z'509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE SUN
+-C
+-C DATA DMACH(1) / Z'0010000000000000' /
+-C DATA DMACH(2) / Z'7FEFFFFFFFFFFFFF' /
+-C DATA DMACH(3) / Z'3CA0000000000000' /
+-C DATA DMACH(4) / Z'3CB0000000000000' /
+-C DATA DMACH(5) / Z'3FD34413509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE SUN
+-C USING THE -r8 COMPILER OPTION
+-C
+-C DATA DMACH(1) / Z'00010000000000000000000000000000' /
+-C DATA DMACH(2) / Z'7FFEFFFFFFFFFFFFFFFFFFFFFFFFFFFF' /
+-C DATA DMACH(3) / Z'3F8E0000000000000000000000000000' /
+-C DATA DMACH(4) / Z'3F8F0000000000000000000000000000' /
+-C DATA DMACH(5) / Z'3FFD34413509F79FEF311F12B35816F9' /
+-C
+-C MACHINE CONSTANTS FOR THE SUN 386i
+-C
+-C DATA SMALL(1), SMALL(2) / Z'FFFFFFFD', Z'000FFFFF' /
+-C DATA LARGE(1), LARGE(2) / Z'FFFFFFB0', Z'7FEFFFFF' /
+-C DATA RIGHT(1), RIGHT(2) / Z'000000B0', Z'3CA00000' /
+-C DATA DIVER(1), DIVER(2) / Z'FFFFFFCB', Z'3CAFFFFF'
+-C DATA LOG10(1), LOG10(2) / Z'509F79E9', Z'3FD34413' /
+-C
+-C MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES FTN COMPILER
+-C
+-C DATA SMALL(1), SMALL(2) / O000040000000, O000000000000 /
+-C DATA LARGE(1), LARGE(2) / O377777777777, O777777777777 /
+-C DATA RIGHT(1), RIGHT(2) / O170540000000, O000000000000 /
+-C DATA DIVER(1), DIVER(2) / O170640000000, O000000000000 /
+-C DATA LOG10(1), LOG10(2) / O177746420232, O411757177572 /
+-C
+-C***FIRST EXECUTABLE STATEMENT D1MACH
+- IF (I .LT. 1 .OR. I .GT. 5) CALL XERMSG ('SLATEC', 'D1MACH',
+- + 'I OUT OF BOUNDS', 1, 2)
+-C
+ D1MACH = DMACH(I)
+ RETURN
++ 9000 FORMAT(/' Adjust D1MACH by uncommenting data statements'/
++ *' appropriate for your machine.')
++ END
++ SUBROUTINE I1MCRY(A, A1, B, C, D)
++**** SPECIAL COMPUTATION FOR OLD CRAY MACHINES ****
++ INTEGER A, A1, B, C, D
++ A1 = 16777216*B + C
++ A = 16777216*A1 + D
+-C
+ END
+--- src/i1mach.f
++++ src/i1mach.f
+@@ -88,801 +88,242 @@
+ C 930201 Added DEC Alpha and SGI constants. (RWC and WRB)
+ C 930618 Corrected I1MACH(5) for Convex -p8 and -pd8 compiler
+ C options. (DWL, RWC and WRB).
++C 150501 Replaced with alternative implementation from BLAS.
+ C***END PROLOGUE I1MACH
+ C
++ INTEGER IMACH(16), OUTPUT, SC, SMALL(2)
++ SAVE IMACH, SC
++ REAL RMACH
++ EQUIVALENCE (IMACH(4),OUTPUT), (RMACH,SMALL(1))
++ INTEGER I3, J, K, T3E(3)
++ DATA T3E(1) / 9777664 /
++ DATA T3E(2) / 5323660 /
++ DATA T3E(3) / 46980 /
++C THIS VERSION ADAPTS AUTOMATICALLY TO MOST CURRENT MACHINES,
++C INCLUDING AUTO-DOUBLE COMPILERS.
++C TO COMPILE ON OLDER MACHINES, ADD A C IN COLUMN 1
++C ON THE NEXT LINE
++ DATA SC/0/
++C AND REMOVE THE C FROM COLUMN 1 IN ONE OF THE SECTIONS BELOW.
++C CONSTANTS FOR EVEN OLDER MACHINES CAN BE OBTAINED BY
++C mail netlib@research.bell-labs.com
++C send old1mach from blas
++C PLEASE SEND CORRECTIONS TO dmg OR ehg@bell-labs.com.
++C
++C MACHINE CONSTANTS FOR THE HONEYWELL DPS 8/70 SERIES.
++C
++C DATA IMACH( 1) / 5 /
++C DATA IMACH( 2) / 6 /
++C DATA IMACH( 3) / 43 /
++C DATA IMACH( 4) / 6 /
++C DATA IMACH( 5) / 36 /
++C DATA IMACH( 6) / 4 /
++C DATA IMACH( 7) / 2 /
++C DATA IMACH( 8) / 35 /
++C DATA IMACH( 9) / O377777777777 /
++C DATA IMACH(10) / 2 /
++C DATA IMACH(11) / 27 /
++C DATA IMACH(12) / -127 /
++C DATA IMACH(13) / 127 /
++C DATA IMACH(14) / 63 /
++C DATA IMACH(15) / -127 /
++C DATA IMACH(16) / 127 /, SC/987/
++C
++C MACHINE CONSTANTS FOR PDP-11 FORTRANS SUPPORTING
+- INTEGER IMACH(16),OUTPUT
+- SAVE IMACH
+- EQUIVALENCE (IMACH(4),OUTPUT)
+-C
+-C MACHINE CONSTANTS FOR THE AMIGA
+-C ABSOFT COMPILER
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 5 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -126 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1022 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE APOLLO
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 129 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1025 /
+-C
+-C MACHINE CONSTANTS FOR THE BURROUGHS 1700 SYSTEM
+-C
+-C DATA IMACH( 1) / 7 /
+-C DATA IMACH( 2) / 2 /
+-C DATA IMACH( 3) / 2 /
+-C DATA IMACH( 4) / 2 /
+-C DATA IMACH( 5) / 36 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 33 /
+-C DATA IMACH( 9) / Z1FFFFFFFF /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -256 /
+-C DATA IMACH(13) / 255 /
+-C DATA IMACH(14) / 60 /
+-C DATA IMACH(15) / -256 /
+-C DATA IMACH(16) / 255 /
+-C
+-C MACHINE CONSTANTS FOR THE BURROUGHS 5700 SYSTEM
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 48 /
+-C DATA IMACH( 6) / 6 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 39 /
+-C DATA IMACH( 9) / O0007777777777777 /
+-C DATA IMACH(10) / 8 /
+-C DATA IMACH(11) / 13 /
+-C DATA IMACH(12) / -50 /
+-C DATA IMACH(13) / 76 /
+-C DATA IMACH(14) / 26 /
+-C DATA IMACH(15) / -50 /
+-C DATA IMACH(16) / 76 /
+-C
+-C MACHINE CONSTANTS FOR THE BURROUGHS 6700/7700 SYSTEMS
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 48 /
+-C DATA IMACH( 6) / 6 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 39 /
+-C DATA IMACH( 9) / O0007777777777777 /
+-C DATA IMACH(10) / 8 /
+-C DATA IMACH(11) / 13 /
+-C DATA IMACH(12) / -50 /
+-C DATA IMACH(13) / 76 /
+-C DATA IMACH(14) / 26 /
+-C DATA IMACH(15) / -32754 /
+-C DATA IMACH(16) / 32780 /
+-C
+-C MACHINE CONSTANTS FOR THE CDC 170/180 SERIES USING NOS/VE
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 64 /
+-C DATA IMACH( 6) / 8 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 63 /
+-C DATA IMACH( 9) / 9223372036854775807 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 47 /
+-C DATA IMACH(12) / -4095 /
+-C DATA IMACH(13) / 4094 /
+-C DATA IMACH(14) / 94 /
+-C DATA IMACH(15) / -4095 /
+-C DATA IMACH(16) / 4094 /
+-C
+-C MACHINE CONSTANTS FOR THE CDC 6000/7000 SERIES
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6LOUTPUT/
+-C DATA IMACH( 5) / 60 /
+-C DATA IMACH( 6) / 10 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 48 /
+-C DATA IMACH( 9) / 00007777777777777777B /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 47 /
+-C DATA IMACH(12) / -929 /
+-C DATA IMACH(13) / 1070 /
+-C DATA IMACH(14) / 94 /
+-C DATA IMACH(15) / -929 /
+-C DATA IMACH(16) / 1069 /
+-C
+-C MACHINE CONSTANTS FOR THE CELERITY C1260
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 0 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / Z'7FFFFFFF' /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -126 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1022 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -fn COMPILER OPTION
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1023 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -fi COMPILER OPTION
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 128 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1024 /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -p8 COMPILER OPTION
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 64 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 63 /
+-C DATA IMACH( 9) / 9223372036854775807 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 53 /
+-C DATA IMACH(12) / -1023 /
+-C DATA IMACH(13) / 1023 /
+-C DATA IMACH(14) / 113 /
+-C DATA IMACH(15) / -16383 /
+-C DATA IMACH(16) / 16383 /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -pd8 COMPILER OPTION
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 64 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 63 /
+-C DATA IMACH( 9) / 9223372036854775807 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 53 /
+-C DATA IMACH(12) / -1023 /
+-C DATA IMACH(13) / 1023 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1023 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE CRAY
+-C USING THE 46 BIT INTEGER COMPILER OPTION
+-C
+-C DATA IMACH( 1) / 100 /
+-C DATA IMACH( 2) / 101 /
+-C DATA IMACH( 3) / 102 /
+-C DATA IMACH( 4) / 101 /
+-C DATA IMACH( 5) / 64 /
+-C DATA IMACH( 6) / 8 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 46 /
+-C DATA IMACH( 9) / 1777777777777777B /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 47 /
+-C DATA IMACH(12) / -8189 /
+-C DATA IMACH(13) / 8190 /
+-C DATA IMACH(14) / 94 /
+-C DATA IMACH(15) / -8099 /
+-C DATA IMACH(16) / 8190 /
+-C
+-C MACHINE CONSTANTS FOR THE CRAY
+-C USING THE 64 BIT INTEGER COMPILER OPTION
+-C
+-C DATA IMACH( 1) / 100 /
+-C DATA IMACH( 2) / 101 /
+-C DATA IMACH( 3) / 102 /
+-C DATA IMACH( 4) / 101 /
+-C DATA IMACH( 5) / 64 /
+-C DATA IMACH( 6) / 8 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 63 /
+-C DATA IMACH( 9) / 777777777777777777777B /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 47 /
+-C DATA IMACH(12) / -8189 /
+-C DATA IMACH(13) / 8190 /
+-C DATA IMACH(14) / 94 /
+-C DATA IMACH(15) / -8099 /
+-C DATA IMACH(16) / 8190 /
+-C
+-C MACHINE CONSTANTS FOR THE DATA GENERAL ECLIPSE S/200
+-C
+-C DATA IMACH( 1) / 11 /
+-C DATA IMACH( 2) / 12 /
+-C DATA IMACH( 3) / 8 /
+-C DATA IMACH( 4) / 10 /
+-C DATA IMACH( 5) / 16 /
+-C DATA IMACH( 6) / 2 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 15 /
+-C DATA IMACH( 9) / 32767 /
+-C DATA IMACH(10) / 16 /
+-C DATA IMACH(11) / 6 /
+-C DATA IMACH(12) / -64 /
+-C DATA IMACH(13) / 63 /
+-C DATA IMACH(14) / 14 /
+-C DATA IMACH(15) / -64 /
+-C DATA IMACH(16) / 63 /
+-C
+-C MACHINE CONSTANTS FOR THE DEC ALPHA
+-C USING G_FLOAT
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 5 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1023 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE DEC ALPHA
+-C USING IEEE_FLOAT
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 128 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1024 /
+-C
+-C MACHINE CONSTANTS FOR THE DEC RISC
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 128 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1024 /
+-C
+-C MACHINE CONSTANTS FOR THE DEC VAX
+-C USING D_FLOATING
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 5 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 56 /
+-C DATA IMACH(15) / -127 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR THE DEC VAX
+-C USING G_FLOATING
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 5 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1023 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE ELXSI 6400
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 32 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -126 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1022 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE HARRIS 220
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 0 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 24 /
+-C DATA IMACH( 6) / 3 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 23 /
+-C DATA IMACH( 9) / 8388607 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 23 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 38 /
+-C DATA IMACH(15) / -127 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR THE HONEYWELL 600/6000 SERIES
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 43 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 36 /
+-C DATA IMACH( 6) / 6 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 35 /
+-C DATA IMACH( 9) / O377777777777 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 27 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 63 /
+-C DATA IMACH(15) / -127 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR THE HP 730
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 128 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1024 /
+-C
+-C MACHINE CONSTANTS FOR THE HP 2100
+-C 3 WORD DOUBLE PRECISION OPTION WITH FTN4
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 4 /
+-C DATA IMACH( 4) / 1 /
+-C DATA IMACH( 5) / 16 /
+-C DATA IMACH( 6) / 2 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 15 /
+-C DATA IMACH( 9) / 32767 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 23 /
+-C DATA IMACH(12) / -128 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 39 /
+-C DATA IMACH(15) / -128 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR THE HP 2100
+-C 4 WORD DOUBLE PRECISION OPTION WITH FTN4
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 4 /
+-C DATA IMACH( 4) / 1 /
+-C DATA IMACH( 5) / 16 /
+-C DATA IMACH( 6) / 2 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 15 /
+-C DATA IMACH( 9) / 32767 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 23 /
+-C DATA IMACH(12) / -128 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 55 /
+-C DATA IMACH(15) / -128 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR THE HP 9000
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 7 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 32 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -126 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1015 /
+-C DATA IMACH(16) / 1017 /
+-C
+-C MACHINE CONSTANTS FOR THE IBM 360/370 SERIES,
+-C THE XEROX SIGMA 5/7/9, THE SEL SYSTEMS 85/86, AND
+-C THE PERKIN ELMER (INTERDATA) 7/32.
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 7 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / Z7FFFFFFF /
+-C DATA IMACH(10) / 16 /
+-C DATA IMACH(11) / 6 /
+-C DATA IMACH(12) / -64 /
+-C DATA IMACH(13) / 63 /
+-C DATA IMACH(14) / 14 /
+-C DATA IMACH(15) / -64 /
+-C DATA IMACH(16) / 63 /
+-C
+-C MACHINE CONSTANTS FOR THE IBM PC
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 0 /
+-C DATA IMACH( 4) / 0 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE IBM RS 6000
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 0 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 128 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1024 /
+-C
+-C MACHINE CONSTANTS FOR THE INTEL i860
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 128 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1024 /
+-C
+-C MACHINE CONSTANTS FOR THE PDP-10 (KA PROCESSOR)
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 5 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 36 /
+-C DATA IMACH( 6) / 5 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 35 /
+-C DATA IMACH( 9) / "377777777777 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 27 /
+-C DATA IMACH(12) / -128 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 54 /
+-C DATA IMACH(15) / -101 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR THE PDP-10 (KI PROCESSOR)
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 5 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 36 /
+-C DATA IMACH( 6) / 5 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 35 /
+-C DATA IMACH( 9) / "377777777777 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 27 /
+-C DATA IMACH(12) / -128 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 62 /
+-C DATA IMACH(15) / -128 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR PDP-11 FORTRAN SUPPORTING
+ C 32-BIT INTEGER ARITHMETIC.
+ C
++C DATA IMACH( 1) / 5 /
++C DATA IMACH( 2) / 6 /
++C DATA IMACH( 3) / 7 /
++C DATA IMACH( 4) / 6 /
++C DATA IMACH( 5) / 32 /
++C DATA IMACH( 6) / 4 /
++C DATA IMACH( 7) / 2 /
++C DATA IMACH( 8) / 31 /
++C DATA IMACH( 9) / 2147483647 /
++C DATA IMACH(10) / 2 /
++C DATA IMACH(11) / 24 /
++C DATA IMACH(12) / -127 /
++C DATA IMACH(13) / 127 /
++C DATA IMACH(14) / 56 /
++C DATA IMACH(15) / -127 /
++C DATA IMACH(16) / 127 /, SC/987/
++C
++C MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES.
++C
++C NOTE THAT THE PUNCH UNIT, I1MACH(3), HAS BEEN SET TO 7
++C WHICH IS APPROPRIATE FOR THE UNIVAC-FOR SYSTEM.
++C IF YOU HAVE THE UNIVAC-FTN SYSTEM, SET IT TO 1.
++C
++C DATA IMACH( 1) / 5 /
++C DATA IMACH( 2) / 6 /
++C DATA IMACH( 3) / 7 /
++C DATA IMACH( 4) / 6 /
++C DATA IMACH( 5) / 36 /
++C DATA IMACH( 6) / 6 /
++C DATA IMACH( 7) / 2 /
++C DATA IMACH( 8) / 35 /
++C DATA IMACH( 9) / O377777777777 /
++C DATA IMACH(10) / 2 /
++C DATA IMACH(11) / 27 /
++C DATA IMACH(12) / -128 /
++C DATA IMACH(13) / 127 /
++C DATA IMACH(14) / 60 /
++C DATA IMACH(15) /-1024 /
++C DATA IMACH(16) / 1023 /, SC/987/
++C
++ IF (SC .NE. 987) THEN
++* *** CHECK FOR AUTODOUBLE ***
++ SMALL(2) = 0
++ RMACH = 1E13
++ IF (SMALL(2) .NE. 0) THEN
++* *** AUTODOUBLED ***
++ IF ( (SMALL(1) .EQ. 1117925532
++ * .AND. SMALL(2) .EQ. -448790528)
++ * .OR. (SMALL(2) .EQ. 1117925532
++ * .AND. SMALL(1) .EQ. -448790528)) THEN
++* *** IEEE ***
++ IMACH(10) = 2
++ IMACH(14) = 53
++ IMACH(15) = -1021
++ IMACH(16) = 1024
++ ELSE IF ( SMALL(1) .EQ. -2065213935
++ * .AND. SMALL(2) .EQ. 10752) THEN
++* *** VAX WITH D_FLOATING ***
++ IMACH(10) = 2
++ IMACH(14) = 56
++ IMACH(15) = -127
++ IMACH(16) = 127
++ ELSE IF ( SMALL(1) .EQ. 1267827943
++ * .AND. SMALL(2) .EQ. 704643072) THEN
++* *** IBM MAINFRAME ***
++ IMACH(10) = 16
++ IMACH(14) = 14
++ IMACH(15) = -64
++ IMACH(16) = 63
++ ELSE
++ WRITE(*,9010)
++ STOP 777
++ END IF
++ IMACH(11) = IMACH(14)
++ IMACH(12) = IMACH(15)
++ IMACH(13) = IMACH(16)
++ ELSE
++ RMACH = 1234567.
++ IF (SMALL(1) .EQ. 1234613304) THEN
++* *** IEEE ***
++ IMACH(10) = 2
++ IMACH(11) = 24
++ IMACH(12) = -125
++ IMACH(13) = 128
++ IMACH(14) = 53
++ IMACH(15) = -1021
++ IMACH(16) = 1024
++ SC = 987
++ ELSE IF (SMALL(1) .EQ. -1271379306) THEN
++* *** VAX ***
++ IMACH(10) = 2
++ IMACH(11) = 24
++ IMACH(12) = -127
++ IMACH(13) = 127
++ IMACH(14) = 56
++ IMACH(15) = -127
++ IMACH(16) = 127
++ SC = 987
++ ELSE IF (SMALL(1) .EQ. 1175639687) THEN
++* *** IBM MAINFRAME ***
++ IMACH(10) = 16
++ IMACH(11) = 6
++ IMACH(12) = -64
++ IMACH(13) = 63
++ IMACH(14) = 14
++ IMACH(15) = -64
++ IMACH(16) = 63
++ SC = 987
++ ELSE IF (SMALL(1) .EQ. 1251390520) THEN
++* *** CONVEX C-1 ***
++ IMACH(10) = 2
++ IMACH(11) = 24
++ IMACH(12) = -128
++ IMACH(13) = 127
++ IMACH(14) = 53
++ IMACH(15) = -1024
++ IMACH(16) = 1023
++ ELSE
++ DO 10 I3 = 1, 3
++ J = SMALL(1) / 10000000
++ K = SMALL(1) - 10000000*J
++ IF (K .NE. T3E(I3)) GO TO 20
++ SMALL(1) = J
++ 10 CONTINUE
++* *** CRAY T3E ***
++ IMACH( 1) = 5
++ IMACH( 2) = 6
++ IMACH( 3) = 0
++ IMACH( 4) = 0
++ IMACH( 5) = 64
++ IMACH( 6) = 8
++ IMACH( 7) = 2
++ IMACH( 8) = 63
++ CALL I1MCR1(IMACH(9), K, 32767, 16777215, 16777215)
++ IMACH(10) = 2
++ IMACH(11) = 53
++ IMACH(12) = -1021
++ IMACH(13) = 1024
++ IMACH(14) = 53
++ IMACH(15) = -1021
++ IMACH(16) = 1024
++ GO TO 35
++ 20 CALL I1MCR1(J, K, 16405, 9876536, 0)
++ IF (SMALL(1) .NE. J) THEN
++ WRITE(*,9020)
++ STOP 777
++ END IF
++* *** CRAY 1, XMP, 2, AND 3 ***
++ IMACH(1) = 5
++ IMACH(2) = 6
++ IMACH(3) = 102
++ IMACH(4) = 6
++ IMACH(5) = 46
++ IMACH(6) = 8
++ IMACH(7) = 2
++ IMACH(8) = 45
++ CALL I1MCR1(IMACH(9), K, 0, 4194303, 16777215)
++ IMACH(10) = 2
++ IMACH(11) = 47
++ IMACH(12) = -8188
++ IMACH(13) = 8189
++ IMACH(14) = 94
++ IMACH(15) = -8141
++ IMACH(16) = 8189
++ GO TO 35
++ END IF
++ END IF
++ IMACH( 1) = 5
++ IMACH( 2) = 6
++ IMACH( 3) = 7
++ IMACH( 4) = 6
++ IMACH( 5) = 32
++ IMACH( 6) = 4
++ IMACH( 7) = 2
++ IMACH( 8) = 31
++ IMACH( 9) = 2147483647
++ 35 SC = 987
++ END IF
++ 9010 FORMAT(/' Adjust autodoubled I1MACH by uncommenting data'/
++ * ' statements appropriate for your machine and setting'/
++ * ' IMACH(I) = IMACH(I+3) for I = 11, 12, and 13.')
++ 9020 FORMAT(/' Adjust I1MACH by uncommenting data statements'/
++ * ' appropriate for your machine.')
++ IF (I .LT. 1 .OR. I .GT. 16) GO TO 40
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 5 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 56 /
+-C DATA IMACH(15) / -127 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR PDP-11 FORTRAN SUPPORTING
+-C 16-BIT INTEGER ARITHMETIC.
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 5 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 16 /
+-C DATA IMACH( 6) / 2 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 15 /
+-C DATA IMACH( 9) / 32767 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 56 /
+-C DATA IMACH(15) / -127 /
+-C DATA IMACH(16) / 127 /
+-C
+-C MACHINE CONSTANTS FOR THE SILICON GRAPHICS
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 128 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1024 /
+-C
+-C MACHINE CONSTANTS FOR THE SUN
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -125 /
+-C DATA IMACH(13) / 128 /
+-C DATA IMACH(14) / 53 /
+-C DATA IMACH(15) / -1021 /
+-C DATA IMACH(16) / 1024 /
+-C
+-C MACHINE CONSTANTS FOR THE SUN
+-C USING THE -r8 COMPILER OPTION
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 6 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 32 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 31 /
+-C DATA IMACH( 9) / 2147483647 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 53 /
+-C DATA IMACH(12) / -1021 /
+-C DATA IMACH(13) / 1024 /
+-C DATA IMACH(14) / 113 /
+-C DATA IMACH(15) / -16381 /
+-C DATA IMACH(16) / 16384 /
+-C
+-C MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES FTN COMPILER
+-C
+-C DATA IMACH( 1) / 5 /
+-C DATA IMACH( 2) / 6 /
+-C DATA IMACH( 3) / 1 /
+-C DATA IMACH( 4) / 6 /
+-C DATA IMACH( 5) / 36 /
+-C DATA IMACH( 6) / 4 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 35 /
+-C DATA IMACH( 9) / O377777777777 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 27 /
+-C DATA IMACH(12) / -128 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 60 /
+-C DATA IMACH(15) / -1024 /
+-C DATA IMACH(16) / 1023 /
+-C
+-C MACHINE CONSTANTS FOR THE Z80 MICROPROCESSOR
+-C
+-C DATA IMACH( 1) / 1 /
+-C DATA IMACH( 2) / 1 /
+-C DATA IMACH( 3) / 0 /
+-C DATA IMACH( 4) / 1 /
+-C DATA IMACH( 5) / 16 /
+-C DATA IMACH( 6) / 2 /
+-C DATA IMACH( 7) / 2 /
+-C DATA IMACH( 8) / 15 /
+-C DATA IMACH( 9) / 32767 /
+-C DATA IMACH(10) / 2 /
+-C DATA IMACH(11) / 24 /
+-C DATA IMACH(12) / -127 /
+-C DATA IMACH(13) / 127 /
+-C DATA IMACH(14) / 56 /
+-C DATA IMACH(15) / -127 /
+-C DATA IMACH(16) / 127 /
+-C
+-C***FIRST EXECUTABLE STATEMENT I1MACH
+- IF (I .LT. 1 .OR. I .GT. 16) GO TO 10
+-C
+ I1MACH = IMACH(I)
+ RETURN
++ 40 WRITE(*,*) 'I1MACH(I): I =',I,' is out of bounds.'
+-C
+- 10 CONTINUE
+- WRITE (UNIT = OUTPUT, FMT = 9000)
+- 9000 FORMAT ('1ERROR 1 IN I1MACH - I OUT OF BOUNDS')
+-C
+-C CALL FDUMP
+-C
+ STOP
+ END
++ SUBROUTINE I1MCR1(A, A1, B, C, D)
++**** SPECIAL COMPUTATION FOR OLD CRAY MACHINES ****
++ INTEGER A, A1, B, C, D
++ A1 = 16777216*B + C
++ A = 16777216*A1 + D
++ END
+--- src/r1mach.f
++++ src/r1mach.f
+@@ -66,6 +66,7 @@
+ C 920625 Added CONVEX -p8 and -pd8 compiler option constants.
+ C (BKS, WRB)
+ C 930201 Added DEC Alpha and SGI constants. (RWC and WRB)
++C 150501 Replaced with alternative implementation from BLAS.
+ C***END PROLOGUE R1MACH
+ C
+ INTEGER SMALL(2)
+@@ -73,347 +74,194 @@
+ INTEGER RIGHT(2)
+ INTEGER DIVER(2)
+ INTEGER LOG10(2)
++C needs to be (2) for AUTODOUBLE, HARRIS SLASH 6, ...
++ INTEGER SC
++ SAVE SMALL, LARGE, RIGHT, DIVER, LOG10, SC
+-C
+ REAL RMACH(5)
+- SAVE RMACH
+-C
+ EQUIVALENCE (RMACH(1),SMALL(1))
+ EQUIVALENCE (RMACH(2),LARGE(1))
+ EQUIVALENCE (RMACH(3),RIGHT(1))
+ EQUIVALENCE (RMACH(4),DIVER(1))
+ EQUIVALENCE (RMACH(5),LOG10(1))
++ INTEGER J, K, L, T3E(3)
++ DATA T3E(1) / 9777664 /
++ DATA T3E(2) / 5323660 /
++ DATA T3E(3) / 46980 /
++C THIS VERSION ADAPTS AUTOMATICALLY TO MOST CURRENT MACHINES,
++C INCLUDING AUTO-DOUBLE COMPILERS.
++C TO COMPILE ON OLDER MACHINES, ADD A C IN COLUMN 1
++C ON THE NEXT LINE
++ DATA SC/0/
++C AND REMOVE THE C FROM COLUMN 1 IN ONE OF THE SECTIONS BELOW.
++C CONSTANTS FOR EVEN OLDER MACHINES CAN BE OBTAINED BY
++C mail netlib@research.bell-labs.com
++C send old1mach from blas
++C PLEASE SEND CORRECTIONS TO dmg OR ehg@bell-labs.com.
++C
++C MACHINE CONSTANTS FOR THE HONEYWELL DPS 8/70 SERIES.
++C DATA RMACH(1) / O402400000000 /
++C DATA RMACH(2) / O376777777777 /
++C DATA RMACH(3) / O714400000000 /
++C DATA RMACH(4) / O716400000000 /
++C DATA RMACH(5) / O776464202324 /, SC/987/
++C
++C MACHINE CONSTANTS FOR PDP-11 FORTRANS SUPPORTING
+-C
+-C MACHINE CONSTANTS FOR THE AMIGA
+-C ABSOFT FORTRAN COMPILER USING THE 68020/68881 COMPILER OPTION
+-C
+-C DATA SMALL(1) / Z'00800000' /
+-C DATA LARGE(1) / Z'7F7FFFFF' /
+-C DATA RIGHT(1) / Z'33800000' /
+-C DATA DIVER(1) / Z'34000000' /
+-C DATA LOG10(1) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE AMIGA
+-C ABSOFT FORTRAN COMPILER USING SOFTWARE FLOATING POINT
+-C
+-C DATA SMALL(1) / Z'00800000' /
+-C DATA LARGE(1) / Z'7EFFFFFF' /
+-C DATA RIGHT(1) / Z'33800000' /
+-C DATA DIVER(1) / Z'34000000' /
+-C DATA LOG10(1) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE APOLLO
+-C
+-C DATA SMALL(1) / 16#00800000 /
+-C DATA LARGE(1) / 16#7FFFFFFF /
+-C DATA RIGHT(1) / 16#33800000 /
+-C DATA DIVER(1) / 16#34000000 /
+-C DATA LOG10(1) / 16#3E9A209B /
+-C
+-C MACHINE CONSTANTS FOR THE BURROUGHS 1700 SYSTEM
+-C
+-C DATA RMACH(1) / Z400800000 /
+-C DATA RMACH(2) / Z5FFFFFFFF /
+-C DATA RMACH(3) / Z4E9800000 /
+-C DATA RMACH(4) / Z4EA800000 /
+-C DATA RMACH(5) / Z500E730E8 /
+-C
+-C MACHINE CONSTANTS FOR THE BURROUGHS 5700/6700/7700 SYSTEMS
+-C
+-C DATA RMACH(1) / O1771000000000000 /
+-C DATA RMACH(2) / O0777777777777777 /
+-C DATA RMACH(3) / O1311000000000000 /
+-C DATA RMACH(4) / O1301000000000000 /
+-C DATA RMACH(5) / O1157163034761675 /
+-C
+-C MACHINE CONSTANTS FOR THE CDC 170/180 SERIES USING NOS/VE
+-C
+-C DATA RMACH(1) / Z"3001800000000000" /
+-C DATA RMACH(2) / Z"4FFEFFFFFFFFFFFE" /
+-C DATA RMACH(3) / Z"3FD2800000000000" /
+-C DATA RMACH(4) / Z"3FD3800000000000" /
+-C DATA RMACH(5) / Z"3FFF9A209A84FBCF" /
+-C
+-C MACHINE CONSTANTS FOR THE CDC 6000/7000 SERIES
+-C
+-C DATA RMACH(1) / 00564000000000000000B /
+-C DATA RMACH(2) / 37767777777777777776B /
+-C DATA RMACH(3) / 16414000000000000000B /
+-C DATA RMACH(4) / 16424000000000000000B /
+-C DATA RMACH(5) / 17164642023241175720B /
+-C
+-C MACHINE CONSTANTS FOR THE CELERITY C1260
+-C
+-C DATA SMALL(1) / Z'00800000' /
+-C DATA LARGE(1) / Z'7F7FFFFF' /
+-C DATA RIGHT(1) / Z'33800000' /
+-C DATA DIVER(1) / Z'34000000' /
+-C DATA LOG10(1) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -fn COMPILER OPTION
+-C
+-C DATA RMACH(1) / Z'00800000' /
+-C DATA RMACH(2) / Z'7FFFFFFF' /
+-C DATA RMACH(3) / Z'34800000' /
+-C DATA RMACH(4) / Z'35000000' /
+-C DATA RMACH(5) / Z'3F9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -fi COMPILER OPTION
+-C
+-C DATA RMACH(1) / Z'00800000' /
+-C DATA RMACH(2) / Z'7F7FFFFF' /
+-C DATA RMACH(3) / Z'33800000' /
+-C DATA RMACH(4) / Z'34000000' /
+-C DATA RMACH(5) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE CONVEX
+-C USING THE -p8 OR -pd8 COMPILER OPTION
+-C
+-C DATA RMACH(1) / Z'0010000000000000' /
+-C DATA RMACH(2) / Z'7FFFFFFFFFFFFFFF' /
+-C DATA RMACH(3) / Z'3CC0000000000000' /
+-C DATA RMACH(4) / Z'3CD0000000000000' /
+-C DATA RMACH(5) / Z'3FF34413509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE CRAY
+-C
+-C DATA RMACH(1) / 200034000000000000000B /
+-C DATA RMACH(2) / 577767777777777777776B /
+-C DATA RMACH(3) / 377224000000000000000B /
+-C DATA RMACH(4) / 377234000000000000000B /
+-C DATA RMACH(5) / 377774642023241175720B /
+-C
+-C MACHINE CONSTANTS FOR THE DATA GENERAL ECLIPSE S/200
+-C NOTE - IT MAY BE APPROPRIATE TO INCLUDE THE FOLLOWING CARD -
+-C STATIC RMACH(5)
+-C
+-C DATA SMALL / 20K, 0 /
+-C DATA LARGE / 77777K, 177777K /
+-C DATA RIGHT / 35420K, 0 /
+-C DATA DIVER / 36020K, 0 /
+-C DATA LOG10 / 40423K, 42023K /
+-C
+-C MACHINE CONSTANTS FOR THE DEC ALPHA
+-C USING G_FLOAT
+-C
+-C DATA RMACH(1) / '00000080'X /
+-C DATA RMACH(2) / 'FFFF7FFF'X /
+-C DATA RMACH(3) / '00003480'X /
+-C DATA RMACH(4) / '00003500'X /
+-C DATA RMACH(5) / '209B3F9A'X /
+-C
+-C MACHINE CONSTANTS FOR THE DEC ALPHA
+-C USING IEEE_FLOAT
+-C
+-C DATA RMACH(1) / '00800000'X /
+-C DATA RMACH(2) / '7F7FFFFF'X /
+-C DATA RMACH(3) / '33800000'X /
+-C DATA RMACH(4) / '34000000'X /
+-C DATA RMACH(5) / '3E9A209B'X /
+-C
+-C MACHINE CONSTANTS FOR THE DEC RISC
+-C
+-C DATA RMACH(1) / Z'00800000' /
+-C DATA RMACH(2) / Z'7F7FFFFF' /
+-C DATA RMACH(3) / Z'33800000' /
+-C DATA RMACH(4) / Z'34000000' /
+-C DATA RMACH(5) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE DEC VAX
+-C (EXPRESSED IN INTEGER AND HEXADECIMAL)
+-C THE HEX FORMAT BELOW MAY NOT BE SUITABLE FOR UNIX SYSTEMS
+-C THE INTEGER FORMAT SHOULD BE OK FOR UNIX SYSTEMS
+-C
+-C DATA SMALL(1) / 128 /
+-C DATA LARGE(1) / -32769 /
+-C DATA RIGHT(1) / 13440 /
+-C DATA DIVER(1) / 13568 /
+-C DATA LOG10(1) / 547045274 /
+-C
+-C DATA SMALL(1) / Z00000080 /
+-C DATA LARGE(1) / ZFFFF7FFF /
+-C DATA RIGHT(1) / Z00003480 /
+-C DATA DIVER(1) / Z00003500 /
+-C DATA LOG10(1) / Z209B3F9A /
+-C
+-C MACHINE CONSTANTS FOR THE ELXSI 6400
+-C (ASSUMING REAL*4 IS THE DEFAULT REAL)
+-C
+-C DATA SMALL(1) / '00800000'X /
+-C DATA LARGE(1) / '7F7FFFFF'X /
+-C DATA RIGHT(1) / '33800000'X /
+-C DATA DIVER(1) / '34000000'X /
+-C DATA LOG10(1) / '3E9A209B'X /
+-C
+-C MACHINE CONSTANTS FOR THE HARRIS 220
+-C
+-C DATA SMALL(1), SMALL(2) / '20000000, '00000201 /
+-C DATA LARGE(1), LARGE(2) / '37777777, '00000177 /
+-C DATA RIGHT(1), RIGHT(2) / '20000000, '00000352 /
+-C DATA DIVER(1), DIVER(2) / '20000000, '00000353 /
+-C DATA LOG10(1), LOG10(2) / '23210115, '00000377 /
+-C
+-C MACHINE CONSTANTS FOR THE HONEYWELL 600/6000 SERIES
+-C
+-C DATA RMACH(1) / O402400000000 /
+-C DATA RMACH(2) / O376777777777 /
+-C DATA RMACH(3) / O714400000000 /
+-C DATA RMACH(4) / O716400000000 /
+-C DATA RMACH(5) / O776464202324 /
+-C
+-C MACHINE CONSTANTS FOR THE HP 730
+-C
+-C DATA RMACH(1) / Z'00800000' /
+-C DATA RMACH(2) / Z'7F7FFFFF' /
+-C DATA RMACH(3) / Z'33800000' /
+-C DATA RMACH(4) / Z'34000000' /
+-C DATA RMACH(5) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE HP 2100
+-C 3 WORD DOUBLE PRECISION WITH FTN4
+-C
+-C DATA SMALL(1), SMALL(2) / 40000B, 1 /
+-C DATA LARGE(1), LARGE(2) / 77777B, 177776B /
+-C DATA RIGHT(1), RIGHT(2) / 40000B, 325B /
+-C DATA DIVER(1), DIVER(2) / 40000B, 327B /
+-C DATA LOG10(1), LOG10(2) / 46420B, 46777B /
+-C
+-C MACHINE CONSTANTS FOR THE HP 2100
+-C 4 WORD DOUBLE PRECISION WITH FTN4
+-C
+-C DATA SMALL(1), SMALL(2) / 40000B, 1 /
+-C DATA LARGE(1), LARGE(2) / 77777B, 177776B /
+-C DATA RIGHT(1), RIGHT(2) / 40000B, 325B /
+-C DATA DIVER(1), DIVER(2) / 40000B, 327B /
+-C DATA LOG10(1), LOG10(2) / 46420B, 46777B /
+-C
+-C MACHINE CONSTANTS FOR THE HP 9000
+-C
+-C DATA SMALL(1) / 00004000000B /
+-C DATA LARGE(1) / 17677777777B /
+-C DATA RIGHT(1) / 06340000000B /
+-C DATA DIVER(1) / 06400000000B /
+-C DATA LOG10(1) / 07646420233B /
+-C
+-C MACHINE CONSTANTS FOR THE IBM 360/370 SERIES,
+-C THE XEROX SIGMA 5/7/9, THE SEL SYSTEMS 85/86 AND
+-C THE PERKIN ELMER (INTERDATA) 7/32.
+-C
+-C DATA RMACH(1) / Z00100000 /
+-C DATA RMACH(2) / Z7FFFFFFF /
+-C DATA RMACH(3) / Z3B100000 /
+-C DATA RMACH(4) / Z3C100000 /
+-C DATA RMACH(5) / Z41134413 /
+-C
+-C MACHINE CONSTANTS FOR THE IBM PC
+-C
+-C DATA SMALL(1) / 1.18E-38 /
+-C DATA LARGE(1) / 3.40E+38 /
+-C DATA RIGHT(1) / 0.595E-07 /
+-C DATA DIVER(1) / 1.19E-07 /
+-C DATA LOG10(1) / 0.30102999566 /
+-C
+-C MACHINE CONSTANTS FOR THE IBM RS 6000
+-C
+-C DATA RMACH(1) / Z'00800000' /
+-C DATA RMACH(2) / Z'7F7FFFFF' /
+-C DATA RMACH(3) / Z'33800000' /
+-C DATA RMACH(4) / Z'34000000' /
+-C DATA RMACH(5) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE INTEL i860
+-C
+-C DATA RMACH(1) / Z'00800000' /
+-C DATA RMACH(2) / Z'7F7FFFFF' /
+-C DATA RMACH(3) / Z'33800000' /
+-C DATA RMACH(4) / Z'34000000' /
+-C DATA RMACH(5) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE PDP-10 (KA OR KI PROCESSOR)
+-C
+-C DATA RMACH(1) / "000400000000 /
+-C DATA RMACH(2) / "377777777777 /
+-C DATA RMACH(3) / "146400000000 /
+-C DATA RMACH(4) / "147400000000 /
+-C DATA RMACH(5) / "177464202324 /
+-C
+-C MACHINE CONSTANTS FOR PDP-11 FORTRAN SUPPORTING
+ C 32-BIT INTEGERS (EXPRESSED IN INTEGER AND OCTAL).
++C DATA SMALL(1) / 8388608 /
++C DATA LARGE(1) / 2147483647 /
++C DATA RIGHT(1) / 880803840 /
++C DATA DIVER(1) / 889192448 /
++C DATA LOG10(1) / 1067065499 /, SC/987/
++C DATA RMACH(1) / O00040000000 /
++C DATA RMACH(2) / O17777777777 /
++C DATA RMACH(3) / O06440000000 /
++C DATA RMACH(4) / O06500000000 /
++C DATA RMACH(5) / O07746420233 /, SC/987/
++C
++C MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES.
++C DATA RMACH(1) / O000400000000 /
++C DATA RMACH(2) / O377777777777 /
++C DATA RMACH(3) / O146400000000 /
++C DATA RMACH(4) / O147400000000 /
++C DATA RMACH(5) / O177464202324 /, SC/987/
++C
++ IF (SC .NE. 987) THEN
++* *** CHECK FOR AUTODOUBLE ***
++ SMALL(2) = 0
++ RMACH(1) = 1E13
++ IF (SMALL(2) .NE. 0) THEN
++* *** AUTODOUBLED ***
++ IF ( SMALL(1) .EQ. 1117925532
++ * .AND. SMALL(2) .EQ. -448790528) THEN
++* *** IEEE BIG ENDIAN ***
++ SMALL(1) = 1048576
++ SMALL(2) = 0
++ LARGE(1) = 2146435071
++ LARGE(2) = -1
++ RIGHT(1) = 1017118720
++ RIGHT(2) = 0
++ DIVER(1) = 1018167296
++ DIVER(2) = 0
++ LOG10(1) = 1070810131
++ LOG10(2) = 1352628735
++ ELSE IF ( SMALL(2) .EQ. 1117925532
++ * .AND. SMALL(1) .EQ. -448790528) THEN
++* *** IEEE LITTLE ENDIAN ***
++ SMALL(2) = 1048576
++ SMALL(1) = 0
++ LARGE(2) = 2146435071
++ LARGE(1) = -1
++ RIGHT(2) = 1017118720
++ RIGHT(1) = 0
++ DIVER(2) = 1018167296
++ DIVER(1) = 0
++ LOG10(2) = 1070810131
++ LOG10(1) = 1352628735
++ ELSE IF ( SMALL(1) .EQ. -2065213935
++ * .AND. SMALL(2) .EQ. 10752) THEN
++* *** VAX WITH D_FLOATING ***
++ SMALL(1) = 128
++ SMALL(2) = 0
++ LARGE(1) = -32769
++ LARGE(2) = -1
++ RIGHT(1) = 9344
++ RIGHT(2) = 0
++ DIVER(1) = 9472
++ DIVER(2) = 0
++ LOG10(1) = 546979738
++ LOG10(2) = -805796613
++ ELSE IF ( SMALL(1) .EQ. 1267827943
++ * .AND. SMALL(2) .EQ. 704643072) THEN
++* *** IBM MAINFRAME ***
++ SMALL(1) = 1048576
++ SMALL(2) = 0
++ LARGE(1) = 2147483647
++ LARGE(2) = -1
++ RIGHT(1) = 856686592
++ RIGHT(2) = 0
++ DIVER(1) = 873463808
++ DIVER(2) = 0
++ LOG10(1) = 1091781651
++ LOG10(2) = 1352628735
++ ELSE
++ WRITE(*,9010)
++ STOP 777
++ END IF
++ ELSE
++ RMACH(1) = 1234567.
++ IF (SMALL(1) .EQ. 1234613304) THEN
++* *** IEEE ***
++ SMALL(1) = 8388608
++ LARGE(1) = 2139095039
++ RIGHT(1) = 864026624
++ DIVER(1) = 872415232
++ LOG10(1) = 1050288283
++ ELSE IF (SMALL(1) .EQ. -1271379306) THEN
++* *** VAX ***
++ SMALL(1) = 128
++ LARGE(1) = -32769
++ RIGHT(1) = 13440
++ DIVER(1) = 13568
++ LOG10(1) = 547045274
++ ELSE IF (SMALL(1) .EQ. 1175639687) THEN
++* *** IBM MAINFRAME ***
++ SMALL(1) = 1048576
++ LARGE(1) = 2147483647
++ RIGHT(1) = 990904320
++ DIVER(1) = 1007681536
++ LOG10(1) = 1091781651
++ ELSE IF (SMALL(1) .EQ. 1251390520) THEN
++* *** CONVEX C-1 ***
++ SMALL(1) = 8388608
++ LARGE(1) = 2147483647
++ RIGHT(1) = 880803840
++ DIVER(1) = 889192448
++ LOG10(1) = 1067065499
++ ELSE
++ DO 10 L = 1, 3
++ J = SMALL(1) / 10000000
++ K = SMALL(1) - 10000000*J
++ IF (K .NE. T3E(L)) GO TO 20
++ SMALL(1) = J
++ 10 CONTINUE
++* *** CRAY T3E ***
++ CALL I1MCRA(SMALL, K, 16, 0, 0)
++ CALL I1MCRA(LARGE, K, 32751, 16777215, 16777215)
++ CALL I1MCRA(RIGHT, K, 15520, 0, 0)
++ CALL I1MCRA(DIVER, K, 15536, 0, 0)
++ CALL I1MCRA(LOG10, K, 16339, 4461392, 10451455)
++ GO TO 30
++ 20 CALL I1MCRA(J, K, 16405, 9876536, 0)
++ IF (SMALL(1) .NE. J) THEN
++ WRITE(*,9020)
++ STOP 777
++ END IF
++* *** CRAY 1, XMP, 2, AND 3 ***
++ CALL I1MCRA(SMALL(1), K, 8195, 8388608, 1)
++ CALL I1MCRA(LARGE(1), K, 24574, 16777215, 16777214)
++ CALL I1MCRA(RIGHT(1), K, 16338, 8388608, 0)
++ CALL I1MCRA(DIVER(1), K, 16339, 8388608, 0)
++ CALL I1MCRA(LOG10(1), K, 16383, 10100890, 8715216)
++ END IF
++ END IF
++ 30 SC = 987
++ END IF
++* SANITY CHECK
++ IF (RMACH(4) .GE. 1.0) STOP 776
++ IF (I .LT. 1 .OR. I .GT. 5) THEN
++ WRITE(*,*) 'R1MACH(I): I =',I,' is out of bounds.'
++ STOP
++ END IF
+-C
+-C DATA SMALL(1) / 8388608 /
+-C DATA LARGE(1) / 2147483647 /
+-C DATA RIGHT(1) / 880803840 /
+-C DATA DIVER(1) / 889192448 /
+-C DATA LOG10(1) / 1067065499 /
+-C
+-C DATA RMACH(1) / O00040000000 /
+-C DATA RMACH(2) / O17777777777 /
+-C DATA RMACH(3) / O06440000000 /
+-C DATA RMACH(4) / O06500000000 /
+-C DATA RMACH(5) / O07746420233 /
+-C
+-C MACHINE CONSTANTS FOR PDP-11 FORTRAN SUPPORTING
+-C 16-BIT INTEGERS (EXPRESSED IN INTEGER AND OCTAL).
+-C
+-C DATA SMALL(1), SMALL(2) / 128, 0 /
+-C DATA LARGE(1), LARGE(2) / 32767, -1 /
+-C DATA RIGHT(1), RIGHT(2) / 13440, 0 /
+-C DATA DIVER(1), DIVER(2) / 13568, 0 /
+-C DATA LOG10(1), LOG10(2) / 16282, 8347 /
+-C
+-C DATA SMALL(1), SMALL(2) / O000200, O000000 /
+-C DATA LARGE(1), LARGE(2) / O077777, O177777 /
+-C DATA RIGHT(1), RIGHT(2) / O032200, O000000 /
+-C DATA DIVER(1), DIVER(2) / O032400, O000000 /
+-C DATA LOG10(1), LOG10(2) / O037632, O020233 /
+-C
+-C MACHINE CONSTANTS FOR THE SILICON GRAPHICS
+-C
+-C DATA RMACH(1) / Z'00800000' /
+-C DATA RMACH(2) / Z'7F7FFFFF' /
+-C DATA RMACH(3) / Z'33800000' /
+-C DATA RMACH(4) / Z'34000000' /
+-C DATA RMACH(5) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE SUN
+-C
+-C DATA RMACH(1) / Z'00800000' /
+-C DATA RMACH(2) / Z'7F7FFFFF' /
+-C DATA RMACH(3) / Z'33800000' /
+-C DATA RMACH(4) / Z'34000000' /
+-C DATA RMACH(5) / Z'3E9A209B' /
+-C
+-C MACHINE CONSTANTS FOR THE SUN
+-C USING THE -r8 COMPILER OPTION
+-C
+-C DATA RMACH(1) / Z'0010000000000000' /
+-C DATA RMACH(2) / Z'7FEFFFFFFFFFFFFF' /
+-C DATA RMACH(3) / Z'3CA0000000000000' /
+-C DATA RMACH(4) / Z'3CB0000000000000' /
+-C DATA RMACH(5) / Z'3FD34413509F79FF' /
+-C
+-C MACHINE CONSTANTS FOR THE UNIVAC 1100 SERIES
+-C
+-C DATA RMACH(1) / O000400000000 /
+-C DATA RMACH(2) / O377777777777 /
+-C DATA RMACH(3) / O146400000000 /
+-C DATA RMACH(4) / O147400000000 /
+-C DATA RMACH(5) / O177464202324 /
+-C
+-C MACHINE CONSTANTS FOR THE Z80 MICROPROCESSOR
+-C
+-C DATA SMALL(1), SMALL(2) / 0, 256/
+-C DATA LARGE(1), LARGE(2) / -1, -129/
+-C DATA RIGHT(1), RIGHT(2) / 0, 26880/
+-C DATA DIVER(1), DIVER(2) / 0, 27136/
+-C DATA LOG10(1), LOG10(2) / 8347, 32538/
+-C
+-C***FIRST EXECUTABLE STATEMENT R1MACH
+- IF (I .LT. 1 .OR. I .GT. 5) CALL XERMSG ('SLATEC', 'R1MACH',
+- + 'I OUT OF BOUNDS', 1, 2)
+-C
+ R1MACH = RMACH(I)
+ RETURN
++ 9010 FORMAT(/' Adjust autodoubled R1MACH by getting data'/
++ *' appropriate for your machine from D1MACH.')
++ 9020 FORMAT(/' Adjust R1MACH by uncommenting data statements'/
++ *' appropriate for your machine.')
++ END
++ SUBROUTINE I1MCRA(A, A1, B, C, D)
++**** SPECIAL COMPUTATION FOR CRAY MACHINES ****
++ INTEGER A, A1, B, C, D
++ A1 = 16777216*B + C
++ A = 16777216*A1 + D
+-C
+ END