blob: 069b820415b65ab869cce6ac5f7d53102cef00c8 [file] [log] [blame]
The Fortran routines here were extracted from LAPACK 3.1.0, then
3.1.1, 3.4.1, 3.4.2, 3.5.0, 3.6.0, 3.6.1, 3.7.0, 3.7.1, 3.8.0, 3.9.0, 3.10.0
(Sources from http://www.netlib.org/lapack and mirrors).
dlapack.f contains support routines for `real' LAPACK calls.
cmplx.f contains support routines for LAPACK calls from complex
arithmetic routines, that are only used if COMPLEX*16 is supported by
the Fortran compiler and compatible with the Rcomplex type.
Neither will be compiled if an external LAPACK is used.
We have tried to avoid including any single-precision routines.
Preparation
-----------
dlapack.f and cmplx.f were created at the top level in the LAPACK
3.9.x sources by copying {dble,cmplx}.txt there and running in R
dir.create('dble')
foo <- scan('dble.txt', "")
file.copy(file.path("SRC", foo), 'dble')
file.copy(file.path("INSTALL", 'ilaver.f'), 'dble')
system('cpp SRC/iparam2stage.F > dble/iparam2stage.f')
system("cat dble/* > dlapack.f")
dir.create('cmplx')
foo <- scan('cmplx.txt', "")
file.copy(file.path("SRC", foo), 'cmplx')
system("cat cmplx/* > cmplx.f")
If you update the LAPACK version, you need to copy over the LICENSE
file and update it in doc/COPYRIGHTS. And if you patch the sources,
*DO* mention it here and add the patch to a file R.patch (not
currently used).
That process will be incomplete for 3.10.0 which introduces .f90 files. Copy
dlartg.f dlassq.f zlartg.f zlassq.f
from 3.9.x into SRC before running the R code above.
New LAPACK versions (even point releases) often add new auxiliary
routines, and these can be hard to track down. We found it helpful to
compare the outputs of
nm -g dlapack.o | grep ' T ' | cut -b20-
nm -g dlapack.o | grep ' U ' | cut -b20-
before and after upating (ditto for cmplx.o). (This may need
adjusting for your platform's nm.)
History
-------
zlahr2 zlaqr0 zlaqr1 zlaqr2 zlaqr3 zlaqr4 zlaqr5 zlartg zpotf2 zpotrf
zrot ztrexc
were added in R 2.5.0.
izmax1 zdrscl zgecon zlacn2 zlantr ztrcon
were added in R 2.7.0.
dzsum1 ilazlc ilazlr zunmhr
were added in R 2.15.2.
zgelsd zgesdd zlacp2 zlacrm zlals0 zlalsa zlalsd zlarcm
were added in R 3.1.0.
zgeequ zgees zgels zgerfs zgesc2 zgesvx zgetc2 zgetri zggbak zggbal
zgges zggev zgghrd zheevd zhgeqz zlaed0 zlaed7 zlaed8 zlaqge zlasyf
zlatdf zlauu2 zlauum zpotri zstedc zsymv zsyr zsytf2 zsytrf zsytri
ztgevc ztgex2 ztgexc ztgsen ztgsy2 ztgsyl ztrsen ztrsyl ztrti2 ztrtri
zunm2l zunmql zunmtr
were added in R 3.3.0 for use by a CRAN package.
dbdsvdx dgesvdx dgetrf2 dgges3 dggev3 dgghd3 dggsvd3 dggsvp3 dorm22 dpotrf2
zgetrf2 zpotrf2
were new in LAPACK 3.6.0 and added in R 3.3.0.
dtrevc3 and ztrevc3 were new in LAPACK 3.6.1 and added in R 3.4.0.
dcombssq.f was patched in R 3.6.3 (see the R-admin manual).
Patches from LAPACK 3.9.0 to 3.10.0 were applied for R 4.2.0 but not
the new .f90 files replacing dlartg.f dlassq.f zlartg.f zlassq.f
-------------------------------------------------
For 'base' R addition of a new .Call()able Lapack routine,
all these files must be updated {path relative to R's topdir}:
src/gnuwin32/Rdll.hide
src/include/R_ext/Lapack.h
src/include/Rmodules/Rlapack.h
src/main/basedecl.h
src/main/lapack.c
src/main/registration.c
src/modules/lapack/Lapack.c
src/modules/lapack/Lapack.h
src/modules/lapack/cmplx.f or dlapack.f