Merge pull request #154 from echeresh/xlange_lapacke

Avoid conversion between layouts in lapacke_?lange_work

Avoid conversion between layouts in lapacke_?lange_work

Fix #147: xlapy2 not propagating nans

xLAPY2 now returns a NaN whenever input variables X or Y are NaNs.

The previous xLAPY2 was relying on FORTRAN INTRINSIC MAX and MIN to behave in a
certain way with NaNs (i.e. return a NaN whenever X or Y are NaN) and this
behavior is not observed on some (most?) compilers.

We handle the NaN behavior of xLAPY2 by checking for NaNs at the start of the
function.

Thanks to Andreas Noack for providing report and sample code to demonstrate the
problem.

Merge pull request #149 from JeffreySax/master

Bug fix in SLAMSWLQ; fix type mismatches.

Bug fix in SLAMSWLQ; fix type mismatches.

- Fix: SLAMSWLQ calls DGEMLQT instead of SGEMLQT.
- Use dummy variable for ?LANGE work array argument in [CS]GETSLS.
- Use 1-element array for workspace query in ?GETSLS.
- Fix PARAMETER declaration type mismatches in ILASLC, IPARAM2STAGE,
SGELQT3, ZHECON_3, ZHECON_ROOK, ZHETRI2X.

Merge pull request #144 from tmijieux/lapacke/fix_tpmqrt_work_allocation

Fix issue #143 (Invalid ?tpmqrt workspace allocation)

Fix issue #143 (Invalid ?tpmqrt workspace allocation)

Merge pull request #142 from iyamazaki/master

fix for issue #141 (how to detect singularity in T).

fix for issue #141 (how to detect singularity in T).

This is related to issue #42.

This is a fix for double precision. Need to fix S, C and Z now.

Fix #52

Added in comment of xLARRV

Note: VU is currently not used by this implmentation of DLARRV, VU is
passed to DLARRV because it could be used compute gaps on the right end
of the extremal eigenvalues. However, with not much initial accuracy in
LAMBDA and VU, the formula can lead to an overestimation of the right gap
and thus to inadequately early RQI 'convergence'. This is currently
prevented this by forcing a small right gap. And so it turns out that VU
is currently not used by this implementation of DLARRV.

Merge pull request #138 from kbergstr/fixes

Fixed formatting (indentation, spaces, style) in some routines

Merge pull request #137 from echeresh/xlarrx_fix

Add quick return for *LARR* routines

Add quick return for *LARR* routines

Line endings changed from win style to lnx style

Changed indendation for correct alignment

took out extra spaces, added workspace calculation condition

took out extra spaces

Merge pull request #136 from echeresh/370_fix

A number of fixes for routines and testing sources integrated in LAPACK 3.7.0

Fix external names and remove unused routines from external statements

Remove redundant IF in TESTING/LIN/*chkrfp.f

Fix *chk{sb,hb}2stg.f: corrected name of subroutine in external section

Align TESTING/{EIG,LIN}/chkxer.f

Align TESTING/{EIG,LIN}/*larhs.f

Align TESTING/{LIN,MATGEN}/*lahilb.f

Align TESTING/{EIG,LIN}/*get02.f

Fix *2st_kernels.f: unify doc comments

Fix *lasyf_aa.f: replace N -> M in comments

Fix *drvls.f: duplicate loop to make work query for each of test cases

Workspace size is not monotonic function of input sizes so there should
be loop over all test cases with workspace queries and finding maximum
between them.

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

This is related to #135

Since the order of multiplication in Fortran is left to right, the expression
      MAXITR*N*N*UNFL
first computes MAXITR*N*N as an INTEGER. This causes INTEGER overflow for
N>=18919.

To avoid the problem, rewrite as:
      MAXITR*(N*(N*UNFL))

Fix #135

Thanks Cezary Dendek!

Merge pull request #134 from echeresh/doc_fix

Fixes for mistakes in documentation in ?GETSLS, ?(SY|HE)TRS_AA, ?LAM?(SWLQ|TSQR), ?(GE|TP)M?(LQ|QR)T(|2|3)

Replace **C -> **H in complex routines

*gemlqt.f, *tpmlqt.f: replace **C -> **H, T is MB-by-K, fix names

*geqrt.f: T is NB-by-K

*gelqt.f: replace NB -> MB

*trs_aa.f: a is input only parameter

*tplqt.f, *tplqt2.f: incorrect dimensions

*gemlqt.f: incorrect dimension of V

*gelqt.f, *gelqt3.f: reflectors are stored in rows

*gemlqt.f, *lamswlq.f, *lamtsqr.f: fix mistakes in comments

Merge pull request #133 from echeresh/laswp_doc_fix

Fix documentation for ?LASWP routine

Fix documentation for ?LASWP routine

Merge pull request #132 from turboencabulator/cleanup

Update makefile prerequisites

Update LAPACK Mailing List

New Mailing List is lapack@icl.utk.edu

Merge pull request #131 from martin-frbg/martin-frbg-patch-1

TESTING/EIG: fix spurious EXTERNAL references to nonexistent functions

Remove spurious entries from EXTERNAL list

ZHETRD_SY2SB and ZHETRD_SB2ST are never called and do not even seem to exist in the current codebase.

Remove spurious entries from EXTERNAL list

ZHETRD_SY2SB and ZHETRD_SB2ST are never called and do not even exist in the current codebase.

Fix typo in EXTERNAL

ZHBTRD_HB2ST in the EXTERNAL statement should be ZHETRD_HB2ST (Windows Ifort treated it as an unresolved reference)

Merge pull request #130 from martin-frbg/lapack-comma

Fix missing comma on EXTERNALs that span multiple lines

Fix missing comma on continued line

EXTERNAL declaration of subroutines missed a comma before the continuation line,
causing a strange run-together name to appear in the object when compiled with ifort.

Fix missing comma on continued line

EXTERNAL declaration of subroutines missed a comma before the continuation line,
causing a strange run-together name to appear in the object when compiled with ifort.

Fix missing comma on continued line

EXTERNAL declaration of subroutines missed a comma before the continuation line,
causing a strange run-together name to appear in the object when compiled with ifort.

README.md with the appveyor badge pointing to the correct URL
(Thanks @Gjacquenot)

README.md with the codecov badge pointing to the correct URL
(Thanks @Gjacquenot)

Merge pull request #129 from Gjacquenot/master

Added Lapack code coverage

[coverage] Added coverage badge in README.md

[coverage] Updated travis ci for code coverage

[coverage] Added CMake rules for code coverage

[coverage] Updated cmake coverage files for Lapack

[coverage] Added cmake files for coverage

Merge pull request #126 from Gjacquenot/master

Added Appveyor CI with MinGW for Windows

[appveyor] Added link for appveyor build status

[appveyor] Skipped CBLAS-*CBLAT1 tests  that fail with appveyor

[appveyor] Added appveyor CI with MinGW

Merge pull request #128 from vladimir-ch/fix-lapacke-ormbr

LAPACKE: fix wrong matrix size in lapacke_?ormbr

LAPACKE: fix wrong matrix size in lapacke_?ormbr

LAPACKE: fix wrong matrix size in lapacke_?ormbr_work

... and move declarations under if.

Merge pull request #127 from vladimir-ch/fix-lapacke-ormlq

LAPACKE: fix wrong number of columns in ?ormlq

LAPACKE: fix wrong number of columns in ?ormlq

Fix #125

Fix #124

Fix #115

Simplify the clean targets

cleanobj:  Remove object files
cleanlib:  Remove libraries
cleanexe:  Remove test and example executables
cleantest: Remove test output and core dumps
clean:     All of the above

Use $(MAKE) -C for recursion

Fix missing or unnecessary prerequisites

Use automatic variables for test inputs

Remove VARIANTSDIR so we don't have to use mkdir

Make the target an actual filename to avoid rebuilding the variant
libraries.

Add libraries as prerequisites

Use the automatic variable $^ to refer to all prerequisites when
linking or creating an archive.

Update comments

Merge pull request #120 from echeresh/lapacke_bdb_csd

lapacke_*(bb|or|un)(csd|bdb): forward calls to LAPACK without conversion

Merge pull request #119 from turboencabulator/cleanup

Linking fixes

Merge pull request #118 from echeresh/lapacke_fixes

Fixes for LAPACKE

Strip out comments that don't apply

lapacke_zgejsv: align parens with lapacke_cgejsv

lapacke_*lacp2_work: remove unnecessary include

Remove tabs and trailing spaces (git diff -w prints nothing)

lapacke_*tp_trans: condition should be inverted

The function converts <in> from <matrix_layout> to opposite layout.
E.g. for upper case (all indices are 0-based):
col-major, upper: a(i,j) is stored in ap[(1 + j)*j/2 + i]
row-major, upper: a(i,j) is stored in ap[(n + n-i+1)*i/2 + j-i]

lapacke_*stemr_work: fixed condition to check ldz: jobz must be checked

lapacke_*(sy|he)(evr|evx)_work: fix expression for # of columns in z

lapacke_*ormlq_work: move declarations under if

lapacke_*bdsvdx: corrected # of columns used for allocation of z_t

lapacke_*gesdd: missing MAX. Should be lrwork = MAX(1, ...)

lapacke_*ge(lq|qr|mlq): add functions to make files

lapacke_*gelq and lapacke_*gemlq: names of functions are incorrect

lapacke_*laswp_work: incorrect value for lda_t (issue #110)

Version the CBLAS and LAPACKE libraries

lapacke_*larfb: char vars should be compared by LAPACKE_lsame

lapacke_*(bb|or|un)(csd|bdb): forward calls to LAPACK without conversion

1) Ignore TRANS parameter for LAPACKE since it has the same meaning
   as matrix_layout parameter.
   TRANS = 'T' means matrices are handled in row-major format.
   TRANS != 'T' means matrices are handled in col-major format.
2) So conversion from/to row-major layout can be removed in LAPACKE.
   All the cases can be forwarded to LAPACK calls
   with corresponding TRANS parameter.
3) nrows_* variables can be safely removed. NaN checkers for
   different TRANS values can be called with corresponding layout.

LAPACKE calls are forwarded in the following way:
trans != 'T', col-major: call LAPACK(trans = 'N') - col-major
trans != 'T', row-major: call LAPACK(trans = 'T') - row-major
trans  = 'T', col-major: call LAPACK(trans = 'T') - row-major
trans  = 'T', row-major: call LAPACK(trans = 'T') - row-major