Some more and fix previous mistakes ...

Put declaration before executable statment to respect ISO 99 standard (Windows compiler was barking)

corrected documentation

Add Rook pivoting file in CMAKE

added object file entries for 'rook' pivoting algorithm to SRC/Makefile:
 ssytf2_rook.o,  slasyf_rook.o, ssytrf_rook.o, ssytrs_rook.o, ssytri_rook.o, ssycon_rook.o, ssysv_rook.o  to SLASRC entry;
 dsytf2_rook.o,  dlasyf_rook.o, dsytrf_rook.o, dsytrs_rook.o, dsytri_rook.o, dsycon_rook.o, dsysv_rook.o  to DLASRC entry;
 csytf2_rook.o,  clasyf_rook.o, csytrf_rook.o, csytrs_rook.o, csytri_rook.o, csycon_rook.o, csysv_rook.o  to CLASRC entry;
 zsytf2_rook.o,  zlasyf_rook.o, zsytrf_rook.o, zsytrs_rook.o, zsytri_rook.o, zsycon_rook.o, zsysv_rook.o  to ZLASRC entry;

added double precision complex (COMPLEX*16) routines for rook pivoting algorithm for for symmetric indefinite matrices: zlasyf_rook.f zsytri_rook.f zsytrs_rook.f zsytrf_rook.f zsysv_rook.f zsycon_rook.f zsytf2_rook.f

fixed a typo in dsysv_rook.f, ssysv_rook.f, csysv_rook.f, replaced a call to *SYTRF by *SYTRF_ROOK

fixed a bug in csycon.f, changed the call to CLACN2 \n from CLACN2( N, WORK( N+1 ), WORK, IWORK, AINVNM, KASE, ISAVE ) \n to CLACN2( N, WORK( N+1 ), WORK, AINVNM, KASE, ISAVE )

collapsed 2 lines of COMPLEX variable declarations into one line in  clasyf_rook.f and csytf2_rook.f

added standard precision (COMPLEX) routines for rook pivoting \\nalgorithm for symmetric indefinite matrices: slasyf_rook.f ssytri_rook.f ssytrs_rook.f ssytrf_rook.f ssysv_rook.f ssycon_rook.f ssytf2_rook.f

updated comments in dsytri_rook.f and ssytri_rook.f

updated comments in dsysv_rook.f and  ssysv_rook.f

fixed comments in ssytrs_rook.f and dsytrs_rook.f

fixed comments in dlasyf_rook.f and slasyf_rook.f

fixed comments in dsytf2_rook.f and ssytf2_rook.f

renamed DTEMP variable into STEMP in single precision routines ssytf2_rook.f slasyf_rook.f

added standard precision (REAL) routines for 'rook' pivoting algorithm for symmetric indefinite matrices: slasyf_rook.f ssytri_rook.f ssytrs_rook.f ssytrf_rook.f ssysv_rook.f ssycon_rook.f ssytf2_rook.f

changed the relational experssions in pivot search loop in dsytf2_rook.f to be the same as in dlasyf_rook.f. \n In particular, ROWMAX .EQ. COLMAX was replaced by the equivalent  ( P.EQ.JMAX ).OR.( ROWMAX.LE.COLMAX ). This approach better handles NaN and Inf

added LAPACK driver dsysv_rook.f

fixes to dsytri_rook.f

fixes to dsycon_rook.f

fixes to dsytrs_rook.f

fixes to dsytrf_rook.f

Integrate LAPACKE into LAPACK make and CMAKE build system.
First commit, works on my mac! Need to see how to adapt to intel compilers and of course Windows.

LAPACKE is not built by default. Is that OK?
User need to type 'make lapackelib' to generate the lib
I add a couple of examples taken from MKL LAPACKE 'make lapacke_example'

fixes to dlasyf_rook.f

fixed bug.

Change RWORK description in the header. The RWORK array needs to be at least of
dimension 2*NMAX to accomodate the calls to CHESVXX in CDRVHEX, ZHESVXX in
ZDRVHEX, CSYSVXX in CDRVSYX, and ZSYSVXX in ZDRVSYX. (The dimension of RWORK in
the caller subroutines --, e.g., ZCHKAA -- is way large enough.)

Patching LAPACK testing routines according to Mike Pont (NAG)'s forum post #2894.
http://icl.cs.utk.edu/lapack-forum/viewtopic.php?t=2894

Fixes to dsytf2_rook.f

Correct some more new lines in program - Thanks Igor for finding out the problem

fixed new lines in testing drivers source files

Update LAPACK version for CMAKE

********************************************************************************

From Ake Sandgren, Umea University and HPC2N, on 11-16-2011.

Hi!

stprfb.f has some unneeded external declarations

diff -ru site/SRC/stprfb.f p1/SRC/stprfb.f
--- site/SRC/stprfb.f   2011-10-22 18:49:03.000000000 +0200
+++ p1/SRC/stprfb.f     2011-11-16 09:30:22.000000000 +0100
@@ -246,7 +246,7 @@
     EXTERNAL  LSAME
*     ..
*     .. External Subroutines ..
-      EXTERNAL  SCOPY, SGEMM, SLACGV, STRMM
+      EXTERNAL  SGEMM, STRMM
*     ..
*     .. Executable Statements ..
*

********************************************************************************

Ake,

Thank you for catching this.  I will remove the unused external declarations.
The changes will be available from the lapack-dev svn server, and in the next
release.

Rodney

********************************************************************************

function call change from DSYTRS to DSYTRS_ROOK in DSYCON_ROOK

added routines for bounded Bunch-Kaufman pivoting factorization algorithm in DOUBLE PRESCISION case:\n dsytf2_rook.f dsytri_rook.f dsytrs_rook.f dsytrf_rook.f dsycon_rook.f dlasyf_rook.f

Various change: add target to generate doc (man and html) and cleanup make.inc's

fixed

Upadte info and date (testing are not included)

Add lapacke inside LAPACK

Update version number to 3.4.0

Correct bug 0065.
Replace Workspace calculations by calls to the routines with LWORK=-1.
This impacts only three routines: xGESVD, xGELSS, x[OR/UN]GBR
Testing have been run to check that the workspace size return is the same.

In a further effort (next major release probably), LAPACK will need to stick
to that LWORK=-1 rule to compute Workspace size.

Change to ilaver to update doc

Correct comments dor Doxygen

Cosmetic changes in Doxygen presentation.
Use \par instead of \details for section.

add a Contributors Section and a Reference Section.
Remove (some) verbatim section when not needed.
Those changes have been done by hand so I am not sure I manage to catch them all.

Never say never...

Last commit related to Doxygen integration following Albert's comment

Correct Warning detected during Doxygen Generation.
Now each routine should have the correct list of arguments.
This allowed to detect and fix problems in parameter description of many routines.

adding error checker for gejsv, they were missing - reported by Nadezhda Mozartova

correction error checkers for sysv, posv, hesv, gbsv - reported by Nadezhda Mozartova (INTEL)

add name of the routine before call

Skip n=0 for rfp test to protect against div by zero - reported by Ake Sandgren on Oct 22nd on lapack mailing list

Correct wrong filename

adding link to individual download, the links will appear directly in Doxygen html documentation

fixing end of routine

set timeout to 36000 to enable test in quadruple prcision to finish

fix problem detected by Sven
xerbla was called without SRNAMT and INFOT.
Actually the problem was that the CSD routines did not have error exits.
So I added them, and i also modified the source files because info was not set correctly.

fix undeclared variables introduced in rev1017:fix to Ming Gu's bug

Fix Problem in testing discovered with the LAPACK dashboard related to rev 983

Use cmake to copy files at installation time

Integrating Doxygen in comments

Changes to the SRC/Makefile to fix bugs that affected 'make single', 'make double', 'make complex', 'make complex16' functionality. Specifically, DSLASRC was included as a dependency for single and double, and ZCLASRC  as a dependency for complex and complex16. Also, a few modifications for the case when USEXBLAS is defined to enable  'make single', 'make double' etc. For more info refer to LAPACK-DEV google group 'LAPACK SRC/Makefile question?' thread.

Fix to Ming Gu's bug in dqds.  See post of LAPACK DEV message board for more details

Per consultation with Julie, added ./testieee to main Makefile lapack_install target

Format changes to INSTALL/Makefile

deleted svn:executable SVN property on SRC/dgeqrt.f

Comments fix to be able to generate the new layout and the corresponding Doxygen documentation

Refactored $(ALLOBJ) files in TESTING/MATGEN/Makefile and a few cosmetic changes

Reverting previous bad commit and going back to revision 1008

convinience change to Makefile to clean BLAS directory from object files

Add the LAPACK library as dependencies for testing after Igor's comment.
This means that if you change one file, all testings will be rerun.

Should we be smarter than that?

Improve CMAKE BUILD system and OUTPUT PARSING when not all precisions are needed.
The following variables will control the precision to be built:
  BUILD_SINGLE
  BUILD_DOUBLE
  BUILD_COMPLEX
  BUILD_COMPLEX16

For mixed precision SINGLE/DOUBLE routines, both BUILD_SINGLE and BUILD_DOUBLE needs to be on.
(same for COMPLEX/COMPLEX16)

Commit patch sent by Sven.
Patch was based on 3.3.0, so required a little bit of merging.
Thank you Sven for the corrections.
Julie

JL
TESTING/EIG
zerrgg.f
Declare ILO, IHI

SRC
---
I have corrected all the following warnings and errors:

c/d/s/zsysv.f, Unused external reference ILAENV
>JL OK
c/d/s/zsysv.f, Unused local variable NB
>JL  ALREADY CORRECTED

sgsvj0.f, Unused intrinsic AMIN1
> JL OK

d/sorbdb.f, Unused intrinsic MIN
> JL OK

sorcsd.f, Unused intrinsic SIN
sorcsd.f, Unused intrinsic COS
> JL OK MERGE

cpoequb.f, Unused intrinsic AIMAG
cpoequb.f, Unused intrinsic REAL
cpoequb.f, Unused local variable ZDUM

zgeequb.f, REAL --> DBLE
> JL OK

sorcsd.f, line 330: Inconsistent data type REAL (previously INTEGER) for argument 14 in reference to SORBDB
sorcsd.f, line 340: Inconsistent data type REAL (previously INTEGER) for argument 4 in reference to SORGQR
sorcsd.f, line 345: Inconsistent data type REAL (previously INTEGER) for argument 4 in reference to SORGQR
sorcsd.f, line 356: Inconsistent data type REAL (previously INTEGER) for argument 4 in reference to SORGLQ
sorcsd.f, line 363: Inconsistent data type REAL (previously INTEGER) for argument 4 in reference to SORGLQ
sorcsd.f, line 369: Inconsistent data type REAL (previously INTEGER) for argument 4 in reference to SORGLQ
sorcsd.f, line 374: Inconsistent data type REAL (previously INTEGER) for argument 4 in reference to SORGLQ
sorcsd.f, line 385: Inconsistent data type REAL (previously INTEGER) for argument 4 in reference to SORGQR
sorcsd.f, line 392: Inconsistent data type REAL (previously INTEGER) for argument 4 in reference to SORGQR
sorcsd.f, line 399: Inconsistent data type REAL (previously INTEGER) for argument 9 in reference to SBBCSD
> JL OK MERGE

ila(s/d/c/z)lr.f
In the WHILE loop, at about line 59 I inserted
IF (I.EQ.0) THEN
   EXIT
END IF
since, otherwise when I = 0, A(0,J) is referenced.
> JL OK

TESTING/MATGEN
--------------
s/dlatm7.f
line 187: I replaced
IF( N.GT.1 ) THEN
by
IF( N.GT.1 .AND. RANK.GT.1 ) THEN
Otherwise I get a division by zero.
> JL OK

s/dlaror.f
I moved:
INFO = 0
from line 125 to line 121 in front of
 IF( N.EQ.0 .OR. M.EQ.0 )
$   RETURN
otherwise INFO is not defined on return when N or M are zero.
Similar change for
c/zlaror.f
> JL OK

TESTING/LIN
-----------
s/derrsy.f
In calls to S/DSYTRI2: I replaced IW by IW(1) since the dummy argument
(LWORK) is a scalar.
> JL OK

c/zerrrfp.f
I made ALPHA and BETA REAL/DOUBLE and introduced COMPLEX/COMPLEX*16
CALPHA since calls to C/ZHFRK need real alpha and beta.  In calls to
CTFSM, I then replaced ALPHA with CALPHA.
In cerrrfp.f, I replaced 1.0D0 with 1.0E0 (not strictly necessary, but
cleaner).
> JL OK

dpot06.f
Removed the unused declaration of LSAME
> JL OK

d/zdrvac.f
Removed the unused declaration of LSAME
> JL OK MERGE

s/d/c/z/ckcsd.f
Removed the unused declarations of S/D/C/ZLANGE, SIN and COS
> JL OK

c/zgennd.f
Removed the unused variable OUT
> JL already DONE

Corrected the following warnings and errors:

dchksy.f, Unused local variable MYWORK
> JL already DONE

ddrvgbx.f, line 792: Different number of arguments from the first call of DGBT02
> JL OK

cdrvgbx.f, line 726: Inconsistent data type REAL (previously COMPLEX) for argument 4 in reference to CLASET
cdrvgbx.f, line 728: Inconsistent data type REAL (previously COMPLEX) for argument 4 in reference to CLASET
cdrvgbx.f, line 780: Inconsistent data type REAL (previously COMPLEX) for argument 10 in reference to CGBT01
cdrvgbx.f, line 794: Different number of arguments from the first call of CGBT02
> JL OK

zdrvgbx.f, line 726: Inconsistent data type DOUBLE PRECISION (previously COMPLEX(KIND(0d0))) for argument 4 in reference to ZLASET
zdrvgbx.f, line 728: Inconsistent data type DOUBLE PRECISION (previously COMPLEX(KIND(0d0))) for argument 4 in reference to ZLASET
zdrvgbx.f, line 780: Inconsistent data type DOUBLE PRECISION (previously COMPLEX(KIND(0d0))) for argument 10 in reference to ZGBT01
zdrvgbx.f, line 794: Different number of arguments from the first call of ZGBT02
> JL OK

TESTING/EIG
-----------
s/d/c/zchkee.f
I inserted
CALL XLAENV( 12, 1 )
at line 1208, since IPARMS(ISPEC) with ISPEC = 12 needs to be set.
(I think I had to do this with them all, certainly with zchkee.f)
> JL OK zchkee.f cchkee.f

c/zchkee.f
In line 2316, the 15th argument of C/ZCKCSD should be real, not complex,
so I replaced A(1,7) by RWORK.
> JL OK

c/zcsdts.f
In calls to C/ZHERK I have replaced ONE by REALONE as arguments 5 and 8
should be real.
> JL OK

s/d/c/zcsdts.f
Division by zero occurs at line 155 in computing EPS when M = 0, so I
have modified
 EPS2 = MAX( ULP, CLANGE( '1', M, M, WORK, LDX, RWORK ) / REAL( M ) )
to
 IF (M.GT.0) THEN
    EPS2 = MAX( ULP,
$               CLANGE( '1', M, M, WORK, LDX, RWORK ) / REAL( M ) )
 ELSE
    EPS2 = ULP
 END IF
> JL OK

SRC/VARIANTS/lu/CR
------------------
s/d/c/zgetrf
Removed the unused intrinsic MOD
> JL OK

Commit pb in TESTING reported by Nadezhda (INTEL) on Aug 29th.

"We notice that errorcheckers tests aren’t correct in some cases.
errorcheckers of ?ggevx, ?pstrf, ?pstf2
Some output parameters are read only in the errorcheckers tests."

Correct Bug in the testing of CHERFSX function
reported by omitrofa (olga.mitrofanova@intel.com)
See Forum 2644

Test of CHERFSX function doesn't work because there is an error
in ../TESTING/LIN/cerrhex.f ( lines 284-330 ).
This function uses the diagonal pivoting factorization of a Hermitian
indefinite matrix instead of the diagonal pivoting factorization of
a Hermitian indefinite packed matrix. When CHERFSX is included in
conditional "IF( LSAMEN( 2, C2, 'HE' ) ) THEN " the test works.

Replace DOUBLE COMPLEX by COMPLEX*16
Replace SNGL by REAL

Add missing files for GNU-built LAPACK DLLs with MS tools

test routines for new QRT subroutines:

xQRT04 tests xGEQRT and xGEMQRT
xQRT05 tests xTPQRT and xTPMQRT

Commiting BRAD patch

Here is a patch for lapack to provide MS-format dll import libraries
from the MinGW GNU toolchain.  It is ready for others to try before it
goes upstream.  You just need to have MinGW GNU Fortran 4.x and one of
the VS IDEs installed.  Configure a MinGW build with BUILD_SHARED_LIBS
set to ON.  The resulting build (and install) tree will provide both
GNU-format and MS-format import libraries for the DLLs.

A C application built with MSVC and linked to the MinGW-built lapack
DLLs will run but requires the GNU runtime DLLs from MinGW to be available.
I tested this with both 32-bit and 64-bit builds.  For the latter I used
the mingw64 gfortran.  In both cases I just put the GNU runtime directory
in my PATH.  I think both libgfortran-3.dll and libgcc_s_dw2-1.dll are
needed.

This may be used to create a binary LAPACK distribution for Windows that
will work without any GNU tools installed.  One needs only to provide the
MinGW GNU runtime libraries along with blas.dll and lapack.dll.

-Brad

added new QRT routines to object lists

QR factorization for triangular-pentagonal matrices
(generalization of triangle-over-square and triangle-over-triangle)

QRT (QR factorization with compact WY representation for Q) routines for general matrices

CHANGE IN THE SCRIPT:
lapack_testing.py has been moved to the top directory
One argument was added :
     - d [dir] is to indicate where is the LAPACK testing directory (.out files). By default, the script will use .-d or --dir to
The -f argument was replaced by -r argument. Now by default the script will use the LAPACK output.
     - r is to use to run the LAPACK tests then analyse the output (.out files). By default, the script will not run all the LAPACK tests

INTEGRATION:
Integration of the LAPACK Testing parsing in the Makefile and CMake build process.

In the Makefile build: now the lapack_testing target will call the script after running the tests
In the CMAKE build: after the tests are run, the script will be executed with the option -s (short)
the end of the output will look like this:
100% tests passed, 0 tests failed out of 98

Total Test time (real) = 211.83 sec
SUMMARY                 numerical error         other error
================        =================       ================
REAL                    40      (0.004%)        0       (0.000%)
DOUBLE PRECISION        202     (0.020%)        0       (0.000%)
COMPLEX                 2       (0.000%)        0       (0.000%)
COMPLEX16               28      (0.005%)        0       (0.000%)

--> ALL PRECISIONS      272     (0.009%)        0       (0.000%)

merge errors for number only mode

Add Julien's suggestions

[ minor edits ]

suggestions:

1) combine Illegal Error and Info Error as one number in the SUMMARY

2) give % values in addition in the SUMMARY, something like: 203 (.01%)

First try to get a nice script to analyze the LAPACK testing.

Below is the help that goes with it.
Please send feedback. I tired to include several levels of output and several kind of errors
So far, I just tested it on my machine, that has no major problems.
I just tested "numerical errors", not "illegal error" neither "info error".

The goal is to have our nightly testing use it and also users.

TO DO:
  - Test other kinds of errors on machines
  - See to include it in ctest and buildbot (need to contact Don and Brad/Chuck)
  - See if a type of error option is required.

./lapack_testing.py [-h|--help] [-s |--short] [-f |--file] [-e |--error] [-p p |--prec p] [-t test |--test test] [-n | --number]
     - h is to print this message
     - f is to use directly the output of the LAPACK testing (.out files). By default, the script will run all the LAPACK tests
 LEVEL OF OUTPUT
     - x is to print a detailed summary
     - e is to print only the error summary
     - s is to print a short summary
     - n is to print the numbers of failing tests (turn on summary mode)
 SECLECTION OF TESTS:
     - p [s/c/d/z] is to indicate the PRECISION to run:
            s=single
            d=double
            c=complex
            z=double complex
            x=all [DEFAULT]
     - t [lin/eig/mixed/rfp/all] is to indicate which TEST FAMILY to run:
            lin=Linear Equation
            eig=Eigen Problems
            mixed=mixed-precision
            rfp=rfp format
            all=all tests [DEFAULT]
 EXAMPLES:
     ./lapack_testing.py -n -f
            Will return the numbers of failed tests by analyzing the LAPACK output
     ./lapack_testing.py -n -f -p s
            Will return the numbers of failed tests in REAL precision by analyzing the LAPACK output
     ./lapack_testing.py -n -f -p s -t eig
            Will return the numbers of failed tests in REAL precision by analyzing only the LAPACK output of EIGEN testings
Written by Julie Langou (June 2011)

Added N*N to LWEDC (and LRWEDC for complex cases) to reflect updated LWORK requirements (see commit r984 in SRC).

Update the name change (rev 985) for the CMAKE build

Following Yana Gurieva (INTEL) 's comment
"Why does Lapack have the routines named ?la_rpvgrw if they are intended for general matrices? Should not their names be ?la_gerpvgrw instead?"

THIS IS THUS A CHANGE IN INTERFACE

added an additional "N*N" to xSTEDC workspace "WORK" to correct problem with LWORK being too small in some cases in the 'V' option when eigenvectors are computed for the dense matrix.

Summer cleanup
Cleanup some codes, like unused variables.
Used -Walls to detect problems.

Added IF statement to prevent potential out of range errors with "K+1" indices.

For Yana Gurieva (Intel), see forum post 2508 on Tue Jun 7th, 2011.

Typo in the description of cla_porcond_c and cla_porcond_x.
( zla_porcond_c and zla_porcond_x are good, just the c versions have the typos.)

See post on the forum 2501 of Yana Gurieva from Intel.  This fixes some typos
in the comments about the **T and **H.  Proof that it was a good idea to make
clear what this $#@! "A'" meant. (In this case it was used to mean either A**T
in the case of cgtsv or A**H in the case of claqsb.)

Correct a workspace requirement comment in cbdsqr and zbdsqr
For case 1, the routine will call DLASCL that requires a 4*N workspace.

Following a bug report from Hatem:

I am calling this function to only compute the singular
values from the bidiagonal matrix I just got after the 2nd stage.
The function comments say that the workspace should be
case 1: 2*N if NCVT = NRU = NCC = 0 (basically no singular vectors computed)
case 2: 4*N-4 otherwise.
I allocated the corresponding size for work in my program (i.e. 2*N) and
call the lapackE interface of this function (LAPACKE_zbdsqr_work).
And I get a "double free or corruption (!prev)" when I deallocate
work.
When I call instead LAPACKE_zbdsqr (I let the wrapper
now doing the allocation of work), it passes.
I looked inside LAPACKE_zbdsqr to check the size of work and
it seems to be 4*N across the board (for both cases 1 and 2).
Now, if I call back LAPACKE_zbdsqr_work and I allocate work
using 4*N, it passes.
Is it a bug?
Thanks for your time,
Hatem

Correction for Tobias.
To get the manpages correcty generated.

From Tobias:
|RPMLINT prints the following warning - and I think it is correct that
Fortran's .NOT. is misread as NOT. macro.

Tobias

RPMLINT report:
===============
lapack-man.noarch: W: manual-page-warning
/usr/share/man/man3/slaqr3.3.gz 351: warning: macro `NOT.BULGE' not defined
lapack-man.noarch: W: manual-page-warning
/usr/share/man/man3/dlaqr3.3.gz 351: warning: macro `NOT.BULGE' not defined
lapack-man.noarch: W: manual-page-warning
/usr/share/man/man3/lsametst.3.gz 53: warning: macro `NOT.LSAME(' not
defined
lapack-man.noarch: W: manual-page-warning
/usr/share/man/man3/slaqr2.3.gz 317: warning: macro `NOT.BULGE' not defined
lapack-man.noarch: W: manual-page-warning
/usr/share/man/man3/dlaqr2.3.gz 317: warning: macro `NOT.BULGE' not defined
This man page may contain problems that can cause it not to be formatted as
intended.|

fixed problems with zero matrix and scan

Apply Sebastien Fabbro patch for the cmake build.

"Two more fixes for the cmake files of lapack-3.3.1:
1) syntax consistency for BLAS_LIBRARIES
2) link lapack with xblas when xblas is enabled"

Update the comment of ABSTOL according to Sven's recommendation

Update ilaver and make.inc to 3.3.1