Commit David fix - sent on Jan 16th 2016 directly
Confirmed by Julie on Jan 18th 2016
Note: issue (1) was fixed in previous commit

Message from David

(1) Line 147 of SRC/zgetrf2.f contains the following code:

      EXTERNAL           ZGEMM, ZSCAL, ZLASWP, ZTRSM, ZERBLA

Inspection of the zgetrf2.f code suggests that ZERBLA on line 147 should in fact be replaced by XERBLA.

--> [JULIE] This was fixed in our SVN on Dec 7th, 2015, revision 1664

(2) Lines 444 to 447 of  TESTING\EIG\cdrvbd.f contains the following code:

      EXTERNAL           ALASVM, XERBLA, CBDT01, CBDT05, CGESDD, CGESVD,
     $                   CGESVJ, CGEJSV, CGESVDX, CLACPY, CLASET, CLATMS,
     $                   CUNT01, CUNT03

Since the length of line 445 is 73 characters the comma at the end of that line is truncated and a symbol CLATMSCUNT01 is erroneously defined. Then during linking an error is generated because this erroneously defined symbol cannot be found.

** Recommendation ** The developers reformat all code to ensure strict adherence to the standard 72 column formatting of fixed format source files. As is well known violation of this formatting rule can result in subtle undetected errors.

(3) Similar formatting errors to that in (2) have been detected in:
TESTING\EIG\zdrvbd.f (line 445)
TESTING\EIG\dchkbd.f (line 548)
TESTING\EIG\ddrvbd.f (line 407)
TESTING\EIG\derrbd.f (line 92)
TESTING\EIG\schkbd.f (line 548)
TESTING\EIG\sdrvbd.f (line 407)
TESTING\EIG\serrbd.f (line 92)

(4) Lines 103 to 104 of TESTING\EIG\cerred.f contains the following code:

      EXTERNAL           CHKXER, CGEES, CGEESX, CGEEV, CGEEVX, CGEJSV
     $                   CGESDD, CGESVD

There is a comma missing from the end of line 103. This does not result in a compilation error. Rather an invalid symbol CGEJSVCGESDD is defined which ultimately causes an error during linking.

Oops..Missing ending bracket in last commit..

Commit interface inconsistencies fix proposed by David Vowles on Jan 18th 2016 - sent directly to Julie
Confirmed by Julie on Jan 18th 2016
Routines: [cz]unbdb[1234].f, [SD]ORBDB[1234].f

== First Email from David ==
I have encountered some additional errors when compiling Lapack v3.6.0 with Compaq Visual Fortran. Interestingly, these errors did not occur when compiling with Intel Visual Fortran.

File: “SRC\cunbdb1.f” Lines: 310-312

            C = SQRT( SCNRM2( P-I, X11(I+1,I+1), 1, X11(I+1,I+1),
     $          1 )**2 + SCNRM2( M-P-I, X21(I+1,I+1), 1, X21(I+1,I+1),
     $          1 )**2 )

The interface to the BLAS SCNRM2 function is SCNRM2(N,X,INCR) whereas the interfaces to the SCNRM2 function calls in “SRC\cunbdb1.f” are of the form SCNRM2(N, X1, INC1, X2, INC2).

Similar interface inconsistencies are found in the following files:
File: “SRC\cunbdb2.f” Lines: 299-300
File: “SRC\cunbdb3.f” Lines: 299-300
File: “SRC\cunbdb4.f” Lines: 347-349

The above inconsistencies between the interface of the function definition and the function call causes compilation failure of the test-function  “xeigtstc.exe”

It is unclear if function calls of the form SCNRM2(N, X1, INC1, X2, INC2) should be replaced by SCNRM2(N, X1, INC1) or if some other correction is required.

== Second Email from David: ==

By the way, I found similar interface inconsistencies in the SRC\cunbdb<n>.f subroutines for <n> = 1 to 4;
and in the SRC\<t>ORBDB<n>.f subroutines for <t> = D and S and <n> = 1 to 4.
A total of 16 subroutines are affected – the four mentioned in my previous email and the twelve mentioned here.

After changing the calls in the affected files from *NRM2(N,X1,INC1,X2,INC2) to *NRM2(N,X1,INC1)
my build for Windows x86 with Compaq Visual Fortran (CVF) completed successfully and all tests passed.
But I do await your verdict as to whether the proposed modification to the function calls is correct.

If I understand correctly STDCALL is the default calling mechanism for the CVF compiler whereas for
the (Intel Visual Fortran) IVF compiler it is the C calling convention. I wonder if this difference
in calling convention explains why the affected test programs compile and link with IVF but fail to link with CVF.

Fix typos reported by Elena

Fix typos reported by Elena

fixing typo reported by Elena Ivanova (Oracle) to Julie on Jan 11th 2016

Committing Christoph Conrads patch sent on Jan 11th on LAPACK Mailing List

From Christoph:
after the changes of the Doxygen config file in SVN revisions 1625
and 1626, routines are grouped together in the documentation, e.g., all
eigenvalue solvers for symmetric matrices in double precision.
Currently, the documentation for all routines in one group is put on
the same HTML page. I find this confusing and it slows my browser.

If every function has its own HTML page, then you can jump to the top
of the page by pressing POS1 to see the function arguments and their
description. Also you can jump to the bottom of the page with END to
see the callgraph of the respective function.

The patch in the attachment enables the SEPARATE_MEMBER_PAGES option
that puts each function documentation on its own HTML page. Please
consider applying this patch.

Committing Sébastien Villemot patches (5) sent on Jan 6th 2016 to LAPACK Mailing list

From Sébastien Villemot:
I attach to this message 5 patches that are currently applied to the
Debian package for LAPACK, and that I think could be usefully merged in
the LAPACK SVN tree.

The patches include fixes for old outstanding bugs in the CBLAS test
programs, fixes for the build system, and fixes for various typos.

Each patch file contains a more detailed description of its purpose.

( remove unused COMPLEX*16 variable T )
( see email from Berend Hasselman on 01/04/2016 )

From Eugene Chereshnev, see: icl.cs.utk.edu/lapack-forum/posting.php?t=4867

Patch from Andreas Noack (MIT).

The code for DGEEV was checking the INFO parameter returned by DHSEQR as
    IF( INFO.GT.0 )THEN "QUIT"

Andreas proposed to check
    IF( INFO.NE.0 )THEN "QUIT"

I am fine with the modification. Not terribly excited but fine.

There is really no reason for the case ( INFO.LT.0 ) to happen, so we should
not have to deal with this case.  If this happen, it is really a bug and we
should be notified.  But, at this point, why not.

I applied the proposed patch to xGEEV and xGEEVX.

Julien.

Bug reported by Tracey Brendan. Thanks Tracey.
See: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4863

Bug reported by Tracey Brendan. Thanks Tracey.
See: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4862

Fix comment to match the code.

Moving variable declaration up...seems to cause problem with certain compilers..Windows

Fix typos in routine declaration

Q doesn't need to be set when vect == v, so no need to check for nan

Z doesn't need to be set when vect == v, so no need to check for nan

oops . . .

this is related to commit r1654.

In commit r1654, I fixed the error parameter in GESVDX to be 17 when LDVT is
invalid, I should also have changed the test suite.  Both the code and the test
suite were incorrect.

Between r1654 and r1659, LAPACK error testing were returning 4 errors because
of this.

Fixed.

Julien.

address bug comment from Martin Köhler ( Max Planck Institute for Dynamics of Complex Technical Systems )
sent to lapack mailing list on Nov 29th at 7:26pm: "Bug in xGETC2"
email of Martin Köhler below

(1) add a quick return for N=0

(2) add a quick return for N=1

(3) no one reviewed this fix, feel free to comment

Dear LAPACK people,

I found a bug concerning the following function:

The xGETC2 is the level-2 BLAS implementation of an LU decomposition
with complete pivoting. Thereby, a "comparison depending on an
uninitialized value" appears if the input is an one-by-one matrix.
Namely the code in line 202 (LAPACK 3.5.0):
    IF( ABS( A( N, N ) ).LT.SMIN ) THEN
accesses SMIN which is set in the previous DO loop. This loop is
only called if N > 1 and so SMIN is not set to a well defined value
if N = 1. But the documentation claims that all matrix N>=0 are
allowed. Furthermore if calling with N = 0 the quick return
statement is missing and the evaluation of ABS(A(N,N)) then yields
an out-of-bound memory violation.

These bug affect the real and the double precision code as well
as their complex counterparts from LAPACK 3.0 up to the current
LAPACK 3.6.

regards,
Martin Köhler

minor: use of TABS instead of SPACES in one line of these two files

Fix bug on [sd]lantr which don't report the result

oops, change some C to Z in the Z routines

Correct workspace computation for [CZ]GESVDX when a workspace query is done

This is a bug report and a bug fix from Lawrence Mulholland from NAG.

Thanks Lawrence!

Also add a variable ITEMPR to index the real workspace RWORK as opposed to
using ITEMP. ITEMP is for complex workspace WORK, while ITEMPR is for the real
workspace RWORK. Sounds good.

Thanks Lawrence!

See http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4851

Correct workspace computation for [s/d]gesvdx.f
See thread: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4850
From Lawrence Mulholland (NAG)
Thanks Lawrence!

Correct the INFO parameter value when LVT is not correct.  (LDVT is 17th
argument and not 16th.)

Thanks to Lawrence Mulholland from NAG.  See:
http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4851

Commit patch submitted by Orion Poplawski on Nov 18th 2015 on LAPACK mailing List

Fix: properly avoid setting the install rpath when installing to system directories.

Fix typo

Fix upper case issue in CMakeLists reported by Julien Shueller on Nov 17th to Julie

In CBLAS/CMakeLists.txt:

configure_file(${CMAKE_CURRENT_SOURCE_DIR}/CMAKE/cblas-config-version.cmake.in
 ${LAPACK_BINARY_DIR}/cblas-config-version.cmake @ONLY)
configure_file(${CMAKE_CURRENT_SOURCE_DIR}/CMAKE/cblas-config-build.cmake.in
 ${LAPACK_BINARY_DIR}/cblas-config.cmake @ONLY)

The /CMAKE subdirectory should be lowercase, I guess some work on OSX which (still!) has a case-insensitive file-system.

Commit fix proposed by Julien Schueller
Sent to LAPACK mailing list on Nov 16th 2015

 spotted a minor inconsistency in the lapacke cmake config files install directory location.

in CMakeLists.txt:

install(FILES
 ${LAPACK_BINARY_DIR}/CMakeFiles/lapack-config.cmake
 ${LAPACK_BINARY_DIR}/lapack-config-version.cmake
 DESTINATION ${LIBRARY_DIR}/cmake/lapack-${LAPACK_VERSION}
 )

in LAPACKE/CMakeLists.txt:

install(FILES
 ${CMAKE_CURRENT_BINARY_DIR}/CMakeFiles/lapacke-config.cmake
 ${LAPACK_BINARY_DIR}/lapacke-config-version.cmake
 DESTINATION lib/cmake/lapacke-${LAPACK_VERSION}
 )

So lapack-config.cmake and lapacke-config-version.cmake could be installed in different directories in cases where CMAKE_INSTALL_LIBDIR is not set to "lib", eg Fedora or RHEL.

I propose to fix LAPACKE/CMakeLists.txt:
sed -i "s|DESTINATION lib/cmake/lapacke|DESTINATION \${LIBRARY_DIR}/cmake/lapacke|g" LAPACKE/CMakeLists.txt

Thanks to Lawrence Mulholland from NAG for improving the newly released complex
Jacobi SVD, and thanks to Zlatko Drmac for following so fast and improving the
code.  See: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4844
for more information about all this. Julien.

Mistake in computation of workspace for ZGGEV3.  When calling ZHGEQZ the
penultimate argument should be RWORK and not WORK.  Easy to fix, not a big big
deal but some compilers might legitimatley complain.

Reported by Pavel Holoborodko on Mon Nov 16 2015. See:
http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4796
Thanks Pavel!

Julien.

Fix typo

Updating Docyfile for manpages

Missing parameter definition

Updating version number

adding the GETRF2 and POTRF2 from LAPACKE to CMAKE

adding xGETRF2 to LAPACKE

adding xPOTRF2 to LAPACKE

adding end bracket...now looks better

removing extra bracket

Adding potrf2 and getrf2 to CMakefile

integration of xGETRF2

integration of xPOTRF2

Add the missing CONJG ...

TESTINGS are all good now!

Adding CGESVJ/ZGESVJ and CGEJSV/ZGESVJ to the testing suite

Note: TEST 15 and 19 (xBDT01  - | A - U diag(S) VT | / ( |A| max(M,N) ulp )  are not passing the threshold
      TEST 15 and 19 are commented until we find the fix

Fix cast...

Adding *esvdx and *bdsvdx to LAPACKE

*** NEED REVIEW ***

See: icl.utk.edu/lapack-forum/viewtopic.php?t=4830

The variables NDIV, NITER and NFAIL are INOUT.  These variables have to be
initialized (e.g., in LASQ2) before callin LASQ3.  LASQ3 increments these
counters.  These counters do not play a role in the algorithm; they just give
statistics.) Thanks to Tracey Brendan for the report.

Adding *gesvj and *gejsv to LAPACKE

Update Doxyfile... one more

Update Doxyfile

 Rook with 2 o , not 3...

Committing patch submitted by Christopher Conrad on Sept 25th 2015

From Christopher:
I noticed the following errors in the documentation and attached
patches for all of them:
- Sometimes the function name in the docstring \brief does not match the
actual function name (ex: xSYGV*).
- According to the documentation, xLAHRD is obsolete but not deprecated.
- The parameter N in xLA(R)SCL2 is not the number of columns of D.

Also, I updated the Doxyfile (doxygen -u DOCS/Doxyfile) and removed
the missing DOT_FONTNAME to get rid of the Doxygen warnings.

Adding xlapmt to LAPACKE ( continue ) .. forgot to save the include file

Adding xlapmt to LAPACKE

Fix cast for lwork in complex [cz]uncsd2by1 routines

Incorporating DBDSVDX tests in dchkbd

Disabling temporarily some [SD]BDSVDX tests - those generate some numerical errors
See email from Osni on Nov 10th 2015

Adding routines to Makefile / CMakefile.txt

Committing Zlatko Drmac Contribution

First commit for Zlatko Drmac Contribution - Fixing z precisions issues + modif sent by Zlatko

First commit for Zlatko Drmac Contribution

Contains only source code with Doxygen format and Makefile modifications for
CGESVJ (Jacobi SVD) with its dependencies: CGSVJ0 and CGSVJ1.

TO DO:
  - Include Testing
  - LAPACKE interface

Change DBLE cast by REAL - confirmed with Osni

Thanks to Lawrence Mulholland (NAG) for reporting the problem.  While the array
E should really be of size N-1, xSYCONV() is not ready for this and so, in the
comment at the header level, we now require E to be of size N. Also xSYTRS2()
is using E(N) = ZERO, so if we make E of size N-1 in xSYCONV(), then we need to
see the impact on xSYTRS2().  This is long term cleaning, so for now, we revert
back to requesting an array E of size N. Thanks Lawrence for noticing this.

Commit  svdx files sent by Osni on Oct 31st 2015
Updated Makefile
Updated file format (Oxygen)
Small Fix: sdrvbd / cdrvbd has a lot of double / double complex in comments

From Osni:
Files for LAPACK/SRC (all new, additions to Makefile are needed):

 - dbdsvdx.f, sbdsvdx.f: full or partial (subset) SVD of a bidiagonal matrix through an associated eigenvalue problem
 - dgesvdx.f, sgesvdx.f, zgesvdx.f, cgesvdx.f: full or partial (subset) SVD of a general matrix by invoking bdsvdx

Files for LAPACK/TESTING/EIG :
- dchkbd.f, schkbd.f: added tests 20-34 for bdsvdx
- dbdt04.f, sbdt04.f (new): needed for tests 25 and 30 in chkbd
- dlahd2.f, slahd2.f: added information about tests 20-34 in chkbd
- derrbd.f, serrbd.f: added tests for the values of INFO returned by bdsvdx
- ddrvbd.f, sdrvbd.f, zdrvbd.f, cdrvbd.f: added tests 23-35 (real case) and 15-27 (complex case) for gesvdx.
- dbdt05.f, sbdt05.f, zbdt05.f, cbdt05.f (new): needed for tests 30,33 (real case) and 22,25 (complex case) for gesvdx
- derred.f, serred.f, zerred.f, cerred.f: added tests for the values of INFO returned by gesvdx

Current Issues:
 - 16 real tests do not pass the threshold ( ssvd.out) - but seems to go through in debug
 - serrbd.f has some DOUBLE PRECISION calculation inside - Is that ok or shall we change it to REAL?
Line 113-119
*     Set the variables to innocuous values.
*
      DO 20 J = 1, NMAX
         DO 10 I = 1, NMAX
            A( I, J ) = 1.D0 / DBLE( I+J )
   10    CONTINUE
   20 CONTINUE

TODO:
     - corresponding LAPACKE routines
     - Test Doxygen format

Fix some typos..

remove extra tabs - Thank you Matyas A Sustik

In ?syconv, replace the variable name WORK by the variable name E.  E is the
standard way to name the supdiagonal/subdiagonal of a symmetric tridiagonal
matrix.  Also, E (previously WORK) is of size N-1, not N. So correct this in
the comment.

The subroutine DSYCONV stores the supdiagonal of the symmetric 1-by-1 or 2-by-2
block matrix in the array WORK. (Or the subdiagonal, same array.) Before
commit, only the nonzero elements of WORK were initialized to their correct
nonzero values. The zero entries were not initialized and not referenced and
assumed to be zero in subsequent subroutines. This is fine, but I think a
better behavior for the subroutine is to initialize these zero entries to zero.
So the commit is initializing these entries to zero.

Also a (double precision) ZERO in the D routine should have been an (integer)
0.

Also fix some formatting to have C, D, and S consistent.

Note: the indentations, comments of the Z routine is not consistent with C, D and
S. I did not have time to fix this.

minor: (1) eliminate trailing spaces in comments and (2) comment indicators
along the left margin of the commented fragment, in order to show explicitly
that the fragment is a comment rather than a code between two comments.
Thanks Dmitry.

Remove NaN check from LAPACKE_?laswp and add a comment on why the NaN check is removed.
This might not be the best fix but this will do for now.
Current behavior: LAPACKE_?laswp does not check for NaNs.
See: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4827 for more information.
Thanks to Dmitry Baksheev from Intel.

Apply Christoph Conrads Doxygen patch for xyycsd2by1.f - comment modification only

Committing temporary fix for bug 113
Submitted by Osni Marques on Oct 22nd 2015
we just need to set the variable NOFAIL to .false. in _larrf.f (line 260).
As it is now, dsyevr may enter into an infinite loop in some rare cases.
With NOFAIL=.false. the code will return an error message.

Fix bug 137 and 138 submitted by Dan Kortschak

Applied patch submitted by Dan Kortschak

Fix bug 139 submitted by Christoph Conrads on Oct 6th 2015
"according to the LAPACK documentation, xSTEDC guarantees to
return the eigenvalues in ascending order (parameter D on exit). To
this end, it sorts the eigenvalues if necessary.

If there is only a single subproblem of size n, no sorting algorithm is
called (cf. {s,d}tedc.f, line 450). Furthermore, xLAED0
({s,d}tedc.f:400) does not guarantee to return eigenvalues in ascending
order. Thus, xSTEDC may return eigenvalues that are not in ascending
order."

Applied patch provided by Christoph Conrads: "always sorting the eigenvalues"

bugs in ?GGHD3 : B,LDA-->B,LDB and zeroing of too large a submatrix

Quick and dirty bug fix.

Bug report and bug fix from Dan Kortschak.

See:
http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4810
https://github.com/xianyi/OpenBLAS/issues/642

This is a quick bug fix, this only checks the triangular part of the trapezoid,
the rectangular part is not checked. Another commit should follow up soon.

Adding [un/or]csd2by1 routines to LAPACKE

Adding xlascl to LAPACKE

 Updating Copyright with current year

Adding xsytrf_rook, xsytrs_rook, xhetrf_rook and xhetrs_rook to LAPACKE

Adding xlahrd to DEPRECATED (problem reported by zerothi)

Fixing comments in [sd]lasd6

Issue : Parameter IDXQ of function SLASD6 is described as [out], should be described as [in,out].
Reported by dbakshee from INTEL on the forum on Sep 18th 2015
Link: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4799

Handle Deprecated routines in LAPACKE Makefile

Issues and fixes submitted by Christoph Conrads on Sept 14th by email
(1st commit)

-> File modified: DOCS/Doxyfile --
- Doxygen documentation
The Doxygen "project version" is still 3.4.2. The attached patch bumps
this to 3.6.0.
-> let's be ready!
- Doxygen documentation
The Doxygen configuration file uses multiple obsolete tags and warnings
are shown right after starting Doxygen. The warnings are easy to see if
you set "QUIET = YES" in DOCS/Doxyfile.
-> set "QUIET = YES"

-> File modified: SRC/*larfgp.f
It says "non-negatibe" in the documentation of xLARFGP.
-> Fixed typo

-> File modified: SRC/*orcsd2by1.f SRC/*uncsd2by1.f
The parameter M in the xORCSD2BY1 is not the number of columns of X and
there is no X12 (parameter P).
-> Fixed

-> File modified: LAPACKE/src/CMakeLists.txt
The C interfaces for the deprecated GSVD routines are still built with
cmake -DBUILD_DEPRECATED=OFF. This causes a linker error with cmake
-DBUILD_TESTING=ON as soon as the Fortran functions are not found in
liblapack.so.
-> Add a if(BUILD_DEPRECATED)

-> File modified: trunk/LAPACKE/include/lapacke.h
The C declaration (lapacke.h:11840-11843) for xGGSVD3 is in a block
of #defines for LAPACK 3.2.X (see lapacke.h:12089,12122).
-> moved them down

Bug Fix for bug0135
Bug reported by kortschak on Sep 10th
Link: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4793
Fix typo substituting O with 0

Fixing bug 133

Bug reported by the GONUM team member: btracey on LAPACK forum
Link: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4771
Description:  the number of columns of A depends on the Side, so the test should compare against m or n.

Fix by Julie.

Minor edits in the BLAS header to better describe the functionality.
See: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4786

Committing Nico's change about build system improvements for LAPACK

"One of the improvements with the patches is that SOVERSION and VERSION are now properly set and the proper symlinks are created:
```
liblapack.so -> liblapack.so.3
liblapack.so.3 -> liblapack.so.3.5.0
liblapack.so.3.5.0
```
Since BLAS is shipped with LAPACK and no separate version number is given for BLAS, I applied the same there, too."

Tested on Julie's Mac

TO DO: To test under Windows
====

Updated OSX RPATH settings
In response to CMake 3.0 generating warnings regarding policy CMP0042,
the OSX RPATH settings have been updated per recommendations found
in the CMake Wiki:
  http://www.cmake.org/Wiki/CMake_RPATH_handling#Mac_OS_X_and_the_RPATH

A is input/output for the subroutines ?sytrs2.f

Correcting Bug 129
See http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4747

 Modify LAPACK_?syconv interface  - Work is used to return values

 (1) Remove the local work variable and its allocation and deallocation
 (2) Add work in LAPACK_?syconv interface
 (3) Update lapacke.h

correcting bug 128
See: http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=4746

When LAPACK_NAN_CHECK is enabled, in input of xLASET, we now allow if the matrix A contains NaN.

Since the goal of xLASET is to set values in A, we consider it OK if A contains NaN in input.

In other words, the rationale is that we do not consider A has an input/output. We consider A as an output only.

Note: We allow neither alpha nor beta to be a NaN.

(So we do check for NaN in alpha and in beta and return error if NaN.)

So this commit
(1) removes the NaN check on matrix A in LAPACKE wrappers: lapacke_?laset.c
(2) declares that A is OUTPUT in LAPACK subroutines: ?laset.f
(3) adds a thanks to Victor Kostin from Intel

Propagate comment added in rev 1568 to make.inc with gfortran compiler

changed variable names from LAM_MAX to SING_MAX and from LAM_MIN to SING_MIN to reflect the chnage in the algorithm for TEST4 for rook pivotig routines, since we now use singular values instead of eigenvalues in TEST4 to compute the norm for 2-by-2 marices in TEST 4. Also corrected the bug in for Hermitian case in TEST4, where during the constraction of Hermitian 2-by-2 matrices the operation of complex conjugation was missing. Changes made to: cchkhe_rook.f cchksy_rook.f dchksy_rook.f schksy_rook.f zchkhe_rook.f zchksy_rook.f

Moved deprecated SGVD routines into DEPRECATED directory.
Added a BUILD_DEPRECATED option for make and cmake, which build
deprecated routines into the LAPACK library (OFF by default).
This completes r1563.

TODO: document the BUILD_DEPRECATED option (i.e. update documentation).

Removed unused and obsolete statement functions.

Added LAPACKE wrappers for GSVD (complements r1572).

Added BLAS3 routines for generalised SVD.

TODO: LAPACKE wrappers.

Fixed BUG number 123, that is testing of 'rook' routines, specifically TEST 4 for UPLO='U' case that test bounds of block diagonal matrices. TEST 4 wasfailing for UPLO='U', since 2-by-2 was passed incorrectly to the -EVX routine.  Also, during the fix the  TEST 4 for both UPLO='U' and 'L' was modified. Instead of computing the ratio of eigenvalues for 2-by-2 matrices, using the -EVX routines, now we compute the ratio of singular values using *GESVD routines. Affected test routines: TESTING/LIN/schksy_rook.f TESTING/LIN/dchksy_rook.f TESTING/LIN/cchksy_rook.f TESTING/LIN/zchksy_rook.f TESTING/LIN/cchkhe_rook.f TESTING/LIN/zchkhe_rook.f

Removed declaration of unused variables.