Adding [SY/HE]TRS_3 to LAPACKE

Merge pull request #97 from haidarazzam/2stage

2stage

Fix Issue #95

 removing GOTO

 merging zhetrd_hb2st

Trying to fix Travis for Fortran build

Suggestion by @BanzaiMan on #68

Merge pull request #94 from cmoha/master

fix for some bugs

redundant variables

fix #5

fix #4

fix #3

Merge remote-tracking branch 'refs/remotes/Reference-LAPACK/master'

fixed definition of some voids

Merge pull request #93 from cmoha/patch-2

bug: fixed compilation error

bug: fixed compilation error

Merge branch 'master' of https://github.com/Reference-LAPACK/lapack

correct WORKSPACE size for dgemqr

remove MB and NB from the code, they are not used and this routine should be agnostic to them

minor

correction of previous commit: WORK needs actuallt to be of size at least 12

WORK needs to be at least of size 6

minor

minor

minor typo in DGEMQR

Adding complex symmetric Aasen routines to LAPACKE

Remove GOTO statment in iparam2stage.F

Polishing code...

Remove #define for precision
Remove Goto

Fixes #92: xLASWP: bug for negative INCX

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.

Correct typos. See email from Berend Hasselman.

merging: Various cleanups to makefiles #84

Contribution by @turboencabulator
Closing #84

Bug reported by Peng HongBo, Joan McComb, and Yi LB Peng from IBM.
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.

defficient spells deficient

print nine-character-long subroutine names (instead of 6-character-long) in TESTING/LIN/xerbla.f

correct a few typos in the Makefile of TESTING

Adding Aasen routines for complex symmetric matrix to CMake Build

Also fix .F case for .F files in CMafileList.txt and some typo.

CMake runs fine now.

Merge pull request #89 from iyamazaki/lapack-aasen

Lapack aasen , Add Aasen's for complex symmetric matrix, and a few cleanups/fixes in testers.

remove RCONDC in testing.

(minor) typo in the code zgetsls.f

working on issue https://github.com/Reference-LAPACK/lapack/issues/85

add Aasen's for complex symmetric into Makefiles.

skip checking solution in tester since xSYTRI is not implemented to
compute condition number.

typos.

First minor change to 2 stage contribution

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

Merge pull request #83 from haidarazzam/master

 adding the 2stage symmetric eigenvalue routines drivers checking

remove compiler warnings.

insert argument check for Aasen's complex symmetric, and also check for
lwork > 0.

fixing lwork for Aasen's (lwork > 0).

some clean ups for Aasen's tester.

add complex symmetric Aasen's into the tester.

testers for complex symmetric Aasen's

typos in comments.

first version of complex symmetric solvers.

Added (S,D,C,Z) (SY,HE) routines, drivers for new rook code

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

Merge pull request #87 from cmoha/patch-1

bug: fixed compilation error

bug: fixed compilation error

Merge pull request #86 from iyamazaki/lapack-aasen

Lapack aasen

Subroutine declarations for Aasen's.

White space.

Argument checks for Aasen's.

White space.

Use short-name (K instead of J1+J-1) for clarity.

Workspace query for ZHESV_AA.

 some tuning

Name change for Aasen in LAPACKE plus adding TRS_AA

Merging changes

cosmetic changes

Changing Assen routines suffix from _Aassen to _aa

 adding the 2stage symmetric eigenvalue routines drivers checking

Aasen sv and trs in LAPACKE

Assen Routines:  Lwork Query for trs

- Code Polish for Assen Routines
- Add Lwork query option (Lwork =-1) for trs Assen routines

Merge pull request #71 from advanpix/master

Fixed desctiption of minimum required size for WORK

Lots of trailing whitespaces in the files of Syd. Cleaning this. No big deal.

Fix #81 reported by @reeuwijk-altium

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)'.

Fix #80 from @reeuwijk-altium

Remove a bunch of unused symbols

Updating CMAKE build for Tall skinny and short wide routines

Merge branch 'sydha-tsqr'
Some files needed to be merged.
Fix extra space, missing variable declaration

make sure no variable in indefined

Fix some minors issue: extra space, variable names.

merging with latest version

Contribution was based on 3.6.x version

Merge branch 'tsqr' of git://github.com/sydha/lapack-1 into sydha-tsqr

Remove duplicate definitions

#67 by turboencabulator

Merge pull request #66 from turboencabulator/pkgconfig

Improve pkgconfig file generation

Merge pull request #73 from mgates3/latb4

make SY/HE imat match LAWN 41 (different than PO)

Merge pull request #78 from oamarques/patch-4

Update zheevr.f

Merge pull request #77 from oamarques/patch-3

Update ssyevr.f

Merge pull request #76 from oamarques/patch-2

Update dsyevr.f

Merge pull request #75 from oamarques/patch-1

Update cheevr.f

Update zheevr.f

Update ssyevr.f

Update dsyevr.f

Update cheevr.f

Tall skinny and short wide routines

Tall skinny and short wide routines

Merge pull request #72 from mgates3/larcm

Add [out] in docs

Add explicit "On entry/exit" in docs

E is output since its contents are modified (destroyed)

make SY/HE imat match LAWN 41 (different than PO)

E is output since its contents are modified (destroyed)

C is output

Fixed comment as for minimum required size of WORK

Merge pull request #64 from cconrads-scicomp/xORCSD2BY1-doc

Doc: describe 2-by-1 CSD identity matrix dimension
Close #27

Merge pull request #62 from cconrads-scicomp/doc-fixes

Doc fixes for xGEEQUB, xGBEQUB, and xPOEQUB