Fix small problem detected by user Yao Toa. sent to the LAPACK mailing list on April 23th 2012
It seems that this one has been around forever

============

hi.

i found a subtle difference between dsyevd.f and ssyevd.f when using lapack 3.4.1.

dsyevd updates LOPT after calling dsytrd and dlacpy. but those codes are not visible in ssyevd.f, cheevd.f and zheevd.f.

pls refer to line 329, 344 in dsyevd.f.

i wonder whether those codes are necessary, because dsyevd has precalculated LOPT with at least 1+6*N+2*N**2.

if those codes must be there, why not for ssyevd?

thanks in advance.

Yao Tao

modified exceptional shift for complex versions - solves problem of nan's appearing on some platforms, most notably IBM/XLF, causing an infinite loop in xLARTG (which also should be fixed by enforcing a max iteration count when computing SCALE)

Fix problems in comments found by Alexander Kobotov

Just taking a quick look on the new release I noticed some minor issues there, which probably caused by a generation script:
- If you explore diff for files (c/z)gbrfsx, cgbsvx, cgbsvxx, cheequb, clanhf, zheequb you could find that number of datatype descriptions in documentation were changed to incorrect one (DOUBLE PRECISION instead of COMPLEX, or COMPLEX*16 instead of just COMPLEX).
+ cunbdb has CMPLX instead of COMPLEX for X12 parameter.
+ zggevx RWORK become REAL instead of correct DOUBLE PRECISION
- Number of files got following string "/ output)", which seems meaningless. Just grep for it.

Fix XLF flags that were causing test zgg to hang

put FOUR parameter back since it is used in this routine

rearrange the test (Julien prefers it that way)

Update date and revision number

Correct bug sent by NAG people

Got a new bug for you! Mick Pont found this problem in DLAED6. The code in question is the 40 loop - when DSCALE(I)=TAU you get a divide by zero (rare in practice). This can cause some compilers to immediately stop, e.g. the Sun compiler.

Mick proposed solution is below:

        DO 40 I = 1, 3
           IF (DSCALE( I ).NE.TAU) THEN
              TEMP = ONE / ( DSCALE( I )-TAU )
              TEMP1 = ZSCALE( I )*TEMP
              TEMP2 = TEMP1*TEMP
              TEMP3 = TEMP2*TEMP
              TEMP4 = TEMP1 / DSCALE( I )
              FC = FC + TEMP4
              ERRETM = ERRETM + ABS( TEMP4 )
              DF = DF + TEMP2
              DDF = DDF + TEMP3
           ELSE
*              On rare occasions dscale(i) can be exactly equal to
*              tau, leading to division by zero. If no trap occurs,
*              there is no problem; the quantities above all overflow
*              and the test on abs(f) below sends you to the end
*              with good results. If a trap occurs, though, the
*              user program will stop. Avoid that happening by
*              jumping directly out.
              GO TO 60
           END IF
  40    CONTINUE

This seems to work OK in our testing.

modified exceptional shift strategy so that QZ now works for various pencils (see netlib.org/lapack/Errata/vrac/lapack_known_issues.html)

Update version number

Fix small issue with lapacke built

Add -qnosave option for xlf compiler (Thanks Rodney)

Follow up to rev r1130

Fix bug bug0088 reported by Mike Pont from NAG on the forum
(see http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=2893)

Actually there were a lot of problems regarding arguments checking.
I tried to correct most of them.

Apply the fix propose to all x[he/sy]rfsx.f routines
- Use IGNORE_CWISE as suggested to prevent use of unitialize variable PARAMS

I also appied the previous TESTING fix to all routines

INFO has new value in ZHERFSX (see description of INFO between ZHESVX and ZHESVXX)
This is set on line 634 (IF ( INFO .LE. N ) INFO = N + J) of zherfsx.f
And this is not handled by the testing LIN/zdrvhex.f
I just add  .AND. INFO.LE.N at line 638 to avoid raising an error when INFO = N + J

At the moment, I would recommand a further look at those routines.
ZHE, ZSY led to 182 Tests failing to pass the threshold
and the same for complex

better than before, but still....

removed comment regarding A being modified and restored - this is no longer the case with the changes made to xLARFT

Remove unused variables

Remove unused parameters

added code to handle N=1 case to fix divide by zero problem that occurs
when the single (diagonal) element is negative

see: http://icl.cs.utk.edu/lapack-forum/archives/lapack/msg01195.html

Fix call to GEMQRT (no lwork) and lwork declarations

Fix External Subroutines declarations

Add the python script for LAPACK Summary output to CTEST

changed SVD test threshold to 50.0 and removed SEP test for matrix type 9 so that all tests pass

changed weights so that the badly conditioned tests (which always fail) are avoided for SXV, DXV, CXV, ZXV.

Fix provided by Hong Bo Peng for bug0026

Undoing unvoluntary change commited in 1237

Following Mathieu's request, some part of lapacke.h were moved to lapacke_utils.h
(the definition of ABS, MAX, MIN, etc...)
Those were conflicting with PLASMA internal definitions.

LAPCKE Routines are now just including lapacke_utils.h which in turn includes lapacke.h

Correct problem in Workspace calculation

Correct bug 0090 Need to unscale if necessary when there is an error in DHGEQZ (QZ iteration failed)

 * bug report by Hong Bo Peng Sandgren, on 03-19-2012.
 * See link:http://icl.cs.utk.edu/lapack-forum/archives/lapack/msg01257.html[LAPACK Mailing list msg 01257]

I am doing some work with DGGEV. When I check the return msg and the actual code, I found something may be wrong. Here is part of comments in the header of DGGEV.F.

*  INFO    (output) INTEGER
*          = 0:  successful exit
*          < 0:  if INFO = -i, the i-th argument had an illegal value.
*          = 1,...,N:
*                The QZ iteration failed.  No eigenvectors have been
*                calculated, but ALPHAR(j), ALPHAI(j), and BETA(j)
*                should be correct for j=INFO+1,...,N.
*          > N:  =N+1: other than QZ iteration failed in DHGEQZ.
*                =N+2: error return from DTGEVC.

When INFO = 1...N, there is an error in DHGEQZ (QZ iteration failed). From the code, we can see it jumps to label 110 then set WORK(1) and return.
But in case of we scaled the matrix, we still need to undo scale for the output array ALPHAR, ALPHAI and BETA for those values j=INFO+1,...,N.

In DGEEVX, we can see that it jumps to label 50 in case of DHSEQR failure and then undo scale before return.

Add csyr and zsyr needed by the PLASMA project

modified so that V is [in] only instead of [in/out]

Change CMAKE configuration for LAPACKE

CMAKE now can find the Fortran Mangling automatically (no need of extra flags)
with the module FortranCInterface and create the necessary header.
Needed to make some adjustements with the routine name.
Added a file lapacke_mangling_with_flags.h that is used with Makefile build

lapacke.h now requires lapacke_mangling.h
 - the file is created by CMAKE if CMAKE is used
 - the file is copied from lapacke_mangling_with_flags if Makefile is used

Fix Makefile in LIN

Fix Typo

Fixed problem with optimizer breaking error checking, now using F90 EPSILON intrinsic.
Thanks to Harald Anlauf for providing the solution.

replaced calculation of EPS with F90 call to EPSILON intrinsic

Update CMAKE build

Fixed QRT routine testing

typo

added a comment to *chksy_rook.f

introduced CALL XLAENV( 2, 2 ) in zchksy.f and  cchksy.f in order to make them uniform with schksy.f and dchksy.f and *chksy_rook (it was discussed with Julien in an email)

added sdrvsy_rook.f cdrvsy_rook.f zdrvsy_rook.f to LIN/CMakeLists.txt

small fixes, removed unused varibales in  ddrvsy_rook.f zdrvsy_rook.f sdrvsy_rook.f

added a test for double precision complex SV driver for symmetric indefinite matrices with rook pivoting

added a test for double precision complex SV driver for symmetric indefinite matrices with rook pivoting

Add an option LAPACKE_WITH_TMG that will generate the LAPACKE with the TMG routines needed in plasma (by default regular LAPACKE with only LAPACK routines

Add the test for single precision driver routines for symmetric indefinite matrices with rook pivoting in CMAKE

added the test for single precision driver routines for symmetric indefinite matrices with rook pivoting

Fix LAPACKE makefile to only include XBLAS routine if USEXBLAS is defined

Adding CMAKE Support for LAPACKE.
Now we can generate dll for LAPACK and LAPACKE directly for Mingw so that FORTRAN compiler is longer needed.
Because LAPACKE contains some routines from MATGEN (for PLASMA), LAPACKE will requires the tmglib library.

Add some LAPACK 3.4.0 routines were missing in the CMAKE LAPACK build.

Fix minor problem in some comments...

Fix bug0085 - xGESVD Problem in Workspace computation

Add Rook pivoting routines to LAPACKE

Update CMAKE build for Rook Pivoting Routines + correct minor fix

Modifications following email from PLASMA Team (Mathieu)

added spacing between argumets in subroutine call

changed error handler alaerh.f so that it recognizes ***'SV' instead of ***'SV ' for 'xSY', 'xSR' and 'xHE' test paths

added the tests, i.e. ddrvsy_rook.f  for driver routines for symmetric indrfinite matrices with rook pivoting, included  this tests in dchkaa.f and Makefile

fixed a typo in aladhd.f

Fix group problem for Level 2 BLAS  routine in double precision

added 'SR' test path to aladhd.f for driver routines that correspond to symmetric indefinite matrices with rook pivoting

added 'SR' test path to TESTING/LIN/*errvx.f  to test *SYSV_ROOK error exits

added a special block diagonal matrix type to tests for COMPLEX*16 symmetric indefinite factorization routines with rook pivoting

added a special block diagonal matrix type to tests for COMPLEX symmetric indefinite factorization routines with rook pivoting

changed a comment in cchkaa.f dchkaa.f schkaa.f

fixed a typo in TESTING/LIN/zerrsy.f

changed 'Complex precision'  into 'Zomplex precision' in TESTING/LIN/zchksy_rook.f

added 'ZSR' test path to test COPLEX*16 routines with rook pivoting in TESTING/LIN/zchksy.f, also added this path to the driver file ztest.in

added 'ZSR' test path to test COPLEX*16 routines with rook pivoting in TESTING/LIN/zchksy.f, also added this path to the driver file ztest.in

added TESTING/LIN/zchksy.f and included in the makefile

line formatting fixes to schksy_rook.f and dchksy_rook.f

updated tests in for COMPLEX symmetric indefinite matrix routines with rook pivoting in TESTING/LIN/cchksy_rook.f

fixed a typo in a comment  in LIN/cchkaa.f

added 'CSR' tet path to test COMPLEX routines for symmetric indefinite matrices with rook pivoting to the main test program TESTING/LIN/cchkaa.f and updated driver file ctest.in

added cchksy_rook.f in TESTING/LIN/ and included it in Makefile

added cchksy_rook.f in TESTING/LIN/ and included it in Makefile

fixed a call with a wrong number of arguments to *SYTRI_ROOK in schksy_rook.f dchksy_rook.f

corrected a type in a comment in TESTING/LIN/schksy_rook.f

changed comment in dchksy.f schksy.f

changed the variable name in TESTING/LIN/schksy_rook.f

changed the order of variable decrarations in TESTING/LIN/schksy_rook.f

fixed a typo 0.0E+o  --> 0.0D+0 in TESTING/LIN/zchksy.f

added CZERO named constat to TESTING/LIN/zchksy.f

added  CZERO constant to cchksy.f

updated a comment in dchksy.f and schksy.f

added 'SSR' test path to test symmetric indefinite factorization routines with rook pivoting to TESTING/LIN/schkaa.f

added 'SSR' test path for symmetric indefinite factorization routines with rook pivoting into stest.in driver file

added 'DSR' test path for symmetric indefinite factorization routines with rook pivoting into dtest.in driver file

added schksy_rook.o to TESTING/LIN/Makefile

added schksy_rook.o to TESTING/LIN/Makefile

added schksy_rook.o to TESTING/LIN/Makefile

changed comments in dchkaa.f

added schksy_rook.f

added a few comments in routines schksy.f dchksy.f cchksy.f zchksy.f

fixed a bug in the counter for the number of completed tests in  schksy.f dchksy.f cchksy.f zchksy.f

final version of the routine dchksy_rook.f that tests dsytrf_rook.f dsytri_rook.f dsytrs_rook.f dsycon_rook.f

fixed a typo in serrsy.f cerrsy.f zerrsy.f

fixed a typo in alahd.f

added 'SR' test path for symmetic indefinite factorization routines with rook pivoting to alaerh.f

added new 'SR' test path (that sorresponds to symmetric indefinite factorization routines with rook pivoting) to the alahd.f routine that prints header information for different test paths

added rook pivoting routines to the tests that check error exits for symmetric indefinite factorization routines, i.e. to serrsy.f derrsy.f cerrsy.f zerrsy.f