aboutsummaryrefslogtreecommitdiff
AgeCommit message (Collapse)Author
2016-12-02Fixes #92: xLASWP: bug for negative INCXVladimír Chalupecký
Thanks to Vladimír Chalupecký for noticing quite an old bug in xLASWP. The bug is fixed now. The cases for INCX = 1 and INCX = -1 were correct. It is not clear whether any other cases would be useful. So the bug was never detected. We might want to consider to (1) forbid any values for INCX but INCX = 1 and INCX = -1, (2) and simplify the code and the subroutine description accordingly. Comments welcome. Julien.
2016-11-26Correct typos. See email from Berend Hasselman.Julien Langou
2016-11-25merging: Various cleanups to makefiles #84Julie
Contribution by @turboencabulator Closing #84
2016-11-25Bug reported by Peng HongBo, Joan McComb, and Yi LB Peng from IBM.Julien Langou
When M=0 or N=0, [C/Z]GESDD would not return the value 1 in WORK(1). Fixing this by following [S/D]GESDD code.
2016-11-24defficient spells deficientJulien Langou
2016-11-24print nine-character-long subroutine names (instead of 6-character-long) in ↵Julien Langou
TESTING/LIN/xerbla.f
2016-11-24correct a few typos in the Makefile of TESTINGJulien Langou
2016-11-23Adding Aasen routines for complex symmetric matrix to CMake BuildJulie
Also fix .F case for .F files in CMafileList.txt and some typo. CMake runs fine now.
2016-11-23Merge pull request #89 from iyamazaki/lapack-aasenjulielangou
Lapack aasen , Add Aasen's for complex symmetric matrix, and a few cleanups/fixes in testers.
2016-11-21remove RCONDC in testing.Ichitaro Yamazaki
2016-11-22(minor) typo in the code zgetsls.fJulien Langou
2016-11-21working on issue https://github.com/Reference-LAPACK/lapack/issues/85Julien Langou
2016-11-21add Aasen's for complex symmetric into Makefiles.Ichitaro Yamazaki
2016-11-20skip checking solution in tester since xSYTRI is not implemented toIchitaro Yamazaki
compute condition number.
2016-11-20typos.Ichitaro Yamazaki
2016-11-19First minor change to 2 stage contributionJulie
Update Makefile to handle .F Update Cmake build Note: We have duplicate loop numbers due to Define statement that stops the compilation with Cmake default option. Fix coming very shortly
2016-11-19Merge pull request #83 from haidarazzam/masterjulielangou
adding the 2stage symmetric eigenvalue routines drivers checking
2016-11-17remove compiler warnings.Ichitaro Yamazaki
2016-11-17insert argument check for Aasen's complex symmetric, and also check forIchitaro Yamazaki
lwork > 0.
2016-11-17fixing lwork for Aasen's (lwork > 0).Ichitaro Yamazaki
2016-11-17some clean ups for Aasen's tester.Ichitaro Yamazaki
2016-11-17add complex symmetric Aasen's into the tester.Ichitaro Yamazaki
2016-11-17testers for complex symmetric Aasen'sIchitaro Yamazaki
2016-11-17typos in comments.Ichitaro Yamazaki
2016-11-17first version of complex symmetric solvers.Ichitaro Yamazaki
2016-11-15Added (S,D,C,Z) (SY,HE) routines, drivers for new rook codeJulie
Close #82 Added routines for new factorization code for symmetric indefinite ( or Hermitian indefinite ) matrices with bounded Bunch-Kaufman ( rook ) pivoting algorithm. New more efficient storage format for factors U ( or L ), block-diagonal matrix D, and pivoting information stored in IPIV: factor L is stored explicitly in lower triangle of A; diagonal of D is stored on the diagonal of A; subdiagonal elements of D are stored in array E; IPIV format is the same as in *_ROOK routines, but differs from SY Bunch-Kaufman routines (e.g. *SYTRF). The factorization output of these new rook _RK routines is not compatible with the existing _ROOK routines and vice versa. This new factorization format is designed in such a way, that there is a possibility in the future to write new Bunch-Kaufman routines that conform to this new factorization format. Then the future Bunch-Kaufman routines could share solver *TRS_3,inversion *TRI_3 and condition estimator *CON_3. To convert between the factorization formats in both ways the following routines are developed: CONVERSION ROUTINES BETWEEN FACTORIZATION FORMATS DOUBLE PRECISION (symmetric indefinite matrices): new file: SRC/dsyconvf.f new file: SRC/dsyconvf_rook.f REAL (symmetric indefinite matrices): new file: SRC/csyconvf.f new file: SRC/csyconvf_rook.f COMPLEX*16 (symmetric indefinite and Hermitian indefinite matrices): new file: SRC/zsyconvf.f new file: SRC/zsyconvf_rook.f COMPLEX (symmetric indefinite and Hermitian indefinite matrices): new file: SRC/ssyconvf.f new file: SRC/ssyconvf_rook.f *SYCONVF routine converts between old Bunch-Kaufman storage format ( denote (L1,D1,IPIV1) ) that is used by *SYTRF and new rook storage format ( denote (L2,D2, IPIV2)) that is used by *SYTRF_RK *SYCONVF_ROOK routine between old rook storage format ( denote (L1,D1,IPIV2) ) that is used by *SYTRF_ROOK and new rook storage format ( denote (L2,D2, IPIV2)) that is used by *SYTRF_RK ROUTINES AND DRIVERS DOUBLE PRECISION (symmetric indefinite matrices): new file: SRC/dsytf2_rk.f BLAS2 unblocked factorization new file: SRC/dlasyf_rk.f BLAS3 auxiliary blocked partial factorization new file: SRC/dsytrf_rk.f BLAS3 blocked factorization new file: SRC/dsytrs_3.f BLAS3 solver new file: SRC/dsycon_3.f BLAS3 condition number estimator new file: SRC/dsytri_3.f BLAS3 inversion, sets the size of work array and calls *sytri_3x new file: SRC/dsytri_3x.f BLAS3 auxiliary inversion, actually computes blocked inversion new file: SRC/dsysv_rk.f BLAS3 solver driver REAL (symmetric indefinite matrices): new file: SRC/ssytf2_rk.f BLAS2 unblocked factorization new file: SRC/slasyf_rk.f BLAS3 auxiliary blocked partial factorization new file: SRC/ssytrf_rk.f BLAS3 blocked factorization new file: SRC/ssytrs_3.f BLAS3 solver new file: SRC/ssycon_3.f BLAS3 condition number estimator new file: SRC/ssytri_3.f BLAS3 inversion, sets the size of work array and calls *sytri_3x new file: SRC/ssytri_3x.f BLAS3 auxiliary inversion, actually computes blocked inversion new file: SRC/ssysv_rk.f BLAS3 solver driver COMPLEX*16 (symmetric indefinite matrices): new file: SRC/zsytf2_rk.f BLAS2 unblocked factorization new file: SRC/zlasyf_rk.f BLAS3 auxiliary blocked partial factorization new file: SRC/zsytrf_rk.f BLAS3 blocked factorization new file: SRC/zsytrs_3.f BLAS3 solver new file: SRC/zsycon_3.f BLAS3 condition number estimator new file: SRC/zsytri_3.f BLAS3 inversion, sets the size of work array and calls *sytri_3x new file: SRC/zsytri_3x.f BLAS3 auxiliary inversion, actually computes blocked inversion new file: SRC/zsysv_rk.f BLAS3 solver driver COMPLEX*16 (Hermitian indefinite matrices): new file: SRC/zhetf2_rk.f BLAS2 unblocked factorization new file: SRC/zlahef_rk.f BLAS3 auxiliary blocked partial factorization new file: SRC/zhetrf_rk.f BLAS3 blocked factorization new file: SRC/zhetrs_3.f BLAS3 solver new file: SRC/zhecon_3.f BLAS3 condition number estimator new file: SRC/zhetri_3.f BLAS3 inversion, sets the size of work array and calls *sytri_3x new file: SRC/zhetri_3x.f BLAS3 auxiliary inversion, actually computes blocked inversion new file: SRC/zhesv_rk.f BLAS3 solver driver COMPLEX (symmetric indefinite matrices): new file: SRC/csytf2_rk.f BLAS2 unblocked factorization new file: SRC/clasyf_rk.f BLAS3 auxiliary blocked partial factorization new file: SRC/csytrf_rk.f BLAS3 blocked factorization new file: SRC/csytrs_3.f BLAS3 solver new file: SRC/csycon_3.f BLAS3 condition number estimator new file: SRC/csytri_3.f BLAS3 inversion, sets the size of work array and calls *sytri_3x new file: SRC/csytri_3x.f BLAS3 auxiliary inversion, actually computes blocked inversion new file: SRC/csysv_rk.f BLAS3 solver driver COMPLEX (Hermitian indefinite matrices): new file: SRC/chetf2_rk.f BLAS2 unblocked factorization new file: SRC/clahef_rk.f BLAS3 auxiliary blocked partial factorization new file: SRC/chetrf_rk.f BLAS3 blocked factorization new file: SRC/chetrs_3.f BLAS3 solver new file: SRC/checon_3.f BLAS3 condition number estimator new file: SRC/chetri_3.f BLAS3 inversion, sets the size of work array and calls *sytri_3x new file: SRC/chetri_3x.f BLAS3 auxiliary inversion, actually computes blocked inversion new file: SRC/chesv_rk.f BLAS3 solver driver MISC modified: SRC/CMakeLists.txt modified: SRC/Makefile TEST CODE modified: TESTING/LIN/CMakeLists.txt modified: TESTING/LIN/Makefile modified: TESTING/LIN/aladhd.f modified: TESTING/LIN/alaerh.f modified: TESTING/LIN/alahd.f DOUBLE PRECISION (symmetric indefinite matrices): modified: TESTING/LIN/dchkaa.f modified: TESTING/LIN/derrsy.f modified: TESTING/LIN/derrsyx.f modified: TESTING/LIN/derrvx.f modified: TESTING/LIN/derrvxx.f modified: TESTING/dtest.in new file: TESTING/LIN/dchksy_rk.f new file: TESTING/LIN/ddrvsy_rk.f new file: TESTING/LIN/dsyt01_3.f REAL (symmetric indefinite matrices): modified: TESTING/LIN/schkaa.f modified: TESTING/LIN/serrsy.f modified: TESTING/LIN/serrsyx.f modified: TESTING/LIN/serrvx.f modified: TESTING/LIN/serrvxx.f modified: TESTING/stest.in new file: TESTING/LIN/schksy_rk.f new file: TESTING/LIN/sdrvsy_rk.f new file: TESTING/LIN/ssyt01_3.f COMPLEX*16 (symmetric indefinite and Hermitian indefinite matrices): modified: TESTING/LIN/zchkaa.f modified: TESTING/LIN/zerrsy.f modified: TESTING/LIN/zerrsyx.f modified: TESTING/LIN/zerrhe.f modified: TESTING/LIN/zerrhex.f modified: TESTING/LIN/zerrvx.f modified: TESTING/LIN/zerrvxx.f modified: TESTING/ztest.in new file: TESTING/LIN/zchksy_rk.f new file: TESTING/LIN/zdrvsy_rk.f new file: TESTING/LIN/zsyt01_3.f new file: TESTING/LIN/zchkhe_rk.f new file: TESTING/LIN/zdrvhe_rk.f new file: TESTING/LIN/zhet01_3.f COMPLEX (symmetric indefinite and Hermitian indefinite matrices): modified: TESTING/LIN/cchkaa.f modified: TESTING/LIN/cerrsy.f modified: TESTING/LIN/cerrsyx.f modified: TESTING/LIN/cerrhe.f modified: TESTING/LIN/cerrhex.f modified: TESTING/LIN/cerrvx.f modified: TESTING/LIN/cerrvxx.f modified: TESTING/ctest.in new file: TESTING/LIN/cchksy_rk.f new file: TESTING/LIN/cdrvsy_rk.f new file: TESTING/LIN/csyt01_3.f new file: TESTING/LIN/cchkhe_rk.f new file: TESTING/LIN/cdrvhe_rk.f new file: TESTING/LIN/chet01_3.f
2016-11-15Merge pull request #87 from cmoha/patch-1langou
bug: fixed compilation error
2016-11-15bug: fixed compilation errorcmoha
2016-11-11Merge pull request #86 from iyamazaki/lapack-aasenjulielangou
Lapack aasen
2016-11-11Subroutine declarations for Aasen's.Ichitaro Yamazaki
2016-11-11White space.Ichitaro Yamazaki
2016-11-11Argument checks for Aasen's.Ichitaro Yamazaki
2016-11-11White space.Ichitaro Yamazaki
2016-11-11Use short-name (K instead of J1+J-1) for clarity.Ichitaro Yamazaki
2016-11-11Workspace query for ZHESV_AA.Ichitaro Yamazaki
2016-11-11 some tuningRenegade
2016-11-11Name change for Aasen in LAPACKE plus adding TRS_AAJulie
2016-11-11Merging changesJulie
cosmetic changes
2016-11-10Changing Assen routines suffix from _Aassen to _aaJulie
2016-11-06 adding the 2stage symmetric eigenvalue routines drivers checkingRenegade
2016-11-03Aasen sv and trs in LAPACKEJulie
2016-11-03Assen Routines: Lwork Query for trsJulie
- Code Polish for Assen Routines - Add Lwork query option (Lwork =-1) for trs Assen routines
2016-11-03Merge pull request #71 from advanpix/masterjulielangou
Fixed desctiption of minimum required size for WORK
2016-11-03Lots of trailing whitespaces in the files of Syd. Cleaning this. No big deal.Julien Langou
2016-11-02Fix #81 reported by @reeuwijk-altiumJulie
The documentation of 'work' in sgesvj states that it should be of size max(4,M+N), whereas de documentation of 'lwork' in sgesvj states that it should be 'length of WORK, WORK >= MAX(6,M+N)'.
2016-11-02Fix #80 from @reeuwijk-altiumJulie
Remove a bunch of unused symbols
2016-10-23Updating CMAKE build for Tall skinny and short wide routinesJulie
2016-10-23Merge branch 'sydha-tsqr'Julie
Some files needed to be merged. Fix extra space, missing variable declaration
2016-10-23make sure no variable in indefinedJulie
Fix some minors issue: extra space, variable names.
2016-10-23merging with latest versionJulie
Contribution was based on 3.6.x version
generated by cgit v1.2.3 (git 2.25.1) at 2024-06-28 19:00:11 +0000