platform/upstream/lapack.git
13 years agomerge errors for number only mode
julie [Mon, 27 Jun 2011 13:22:47 +0000 (13:22 +0000)]
merge errors for number only mode

13 years agoAdd Julien's suggestions
julie [Mon, 27 Jun 2011 13:05:59 +0000 (13:05 +0000)]
Add Julien's suggestions

13 years ago[ minor edits ]
langou [Fri, 24 Jun 2011 13:31:01 +0000 (13:31 +0000)]
[ 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%)

13 years agoFirst try to get a nice script to analyze the LAPACK testing.
julie [Fri, 24 Jun 2011 08:46:19 +0000 (08:46 +0000)]
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)

13 years agoAdded N*N to LWEDC (and LRWEDC for complex cases) to reflect updated LWORK requiremen...
james [Wed, 22 Jun 2011 13:37:15 +0000 (13:37 +0000)]
Added N*N to LWEDC (and LRWEDC for complex cases) to reflect updated LWORK requirements (see commit r984 in SRC).

13 years agoUpdate the name change (rev 985) for the CMAKE build
julie [Wed, 22 Jun 2011 11:50:43 +0000 (11:50 +0000)]
Update the name change (rev 985) for the CMAKE build

13 years agoFollowing Yana Gurieva (INTEL) 's comment
julie [Wed, 22 Jun 2011 11:49:32 +0000 (11:49 +0000)]
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

13 years agoadded an additional "N*N" to xSTEDC workspace "WORK" to correct problem with LWORK...
james [Wed, 22 Jun 2011 11:27:16 +0000 (11:27 +0000)]
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.

13 years agoSummer cleanup
julie [Wed, 22 Jun 2011 09:45:40 +0000 (09:45 +0000)]
Summer cleanup
Cleanup some codes, like unused variables.
Used -Walls to detect problems.

13 years agoAdded IF statement to prevent potential out of range errors with "K+1" indices.
james [Tue, 21 Jun 2011 18:59:51 +0000 (18:59 +0000)]
Added IF statement to prevent potential out of range errors with "K+1" indices.

13 years agoFor Yana Gurieva (Intel), see forum post 2508 on Tue Jun 7th, 2011.
langou [Tue, 7 Jun 2011 09:11:05 +0000 (09:11 +0000)]
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.)

13 years agoSee post on the forum 2501 of Yana Gurieva from Intel. This fixes some typos
langou [Thu, 2 Jun 2011 18:30:26 +0000 (18:30 +0000)]
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.)

13 years agoCorrect a workspace requirement comment in cbdsqr and zbdsqr
julie [Thu, 19 May 2011 14:46:50 +0000 (14:46 +0000)]
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

13 years agoCorrection for Tobias.
julie [Wed, 18 May 2011 08:27:36 +0000 (08:27 +0000)]
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.|

13 years agofixed problems with zero matrix and scan
james [Thu, 21 Apr 2011 22:00:11 +0000 (22:00 +0000)]
fixed problems with zero matrix and scan

13 years agoApply Sebastien Fabbro patch for the cmake build.
julie [Wed, 20 Apr 2011 23:13:17 +0000 (23:13 +0000)]
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"

13 years agoUpdate the comment of ABSTOL according to Sven's recommendation
julie [Tue, 19 Apr 2011 20:17:47 +0000 (20:17 +0000)]
Update the comment of ABSTOL according to Sven's recommendation

13 years agoUpdate ilaver and make.inc to 3.3.1
julie [Mon, 18 Apr 2011 20:00:50 +0000 (20:00 +0000)]
Update ilaver and make.inc to 3.3.1

13 years agoUpdate Version Number in CMAKE
julie [Fri, 15 Apr 2011 16:09:49 +0000 (16:09 +0000)]
Update Version Number in CMAKE

13 years agoUpdate README for 3.3.1
julie [Thu, 14 Apr 2011 15:21:30 +0000 (15:21 +0000)]
Update README for 3.3.1

13 years agoUpadte header for the modified routine for the 3.3.1 release
julie [Wed, 13 Apr 2011 16:30:39 +0000 (16:30 +0000)]
Upadte header for the modified routine for the 3.3.1 release

13 years agoSecond and hopefully last pass to homgenize notation for transpose (**T) and conjugat...
julie [Thu, 7 Apr 2011 13:53:54 +0000 (13:53 +0000)]
Second and hopefully last pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024

Please take a look and let me know if you find some old notation of transpose.
I am going to close bug0024.
Julie

13 years agodecreased EPS by a factor of 0.9 to correct problem with orthogonality of U in certai...
james [Thu, 7 Apr 2011 12:30:43 +0000 (12:30 +0000)]
decreased EPS by a factor of 0.9 to correct problem with orthogonality of U in certain cases

13 years agoCorrect bug0077 : [DS]GESVD Minimum Worksize comments need clarification
julie [Tue, 5 Apr 2011 15:16:30 +0000 (15:16 +0000)]
Correct bug0077 : [DS]GESVD Minimum Worksize comments need clarification
Depending on the PATH, the minimum workspace is not the same.

13 years agoFirst pass for BLAS to homgenize notation for transpose (**T) and conjugate transpose...
julie [Sat, 2 Apr 2011 11:46:46 +0000 (11:46 +0000)]
First pass for BLAS to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024

13 years agoFirst pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
julie [Sat, 2 Apr 2011 11:08:56 +0000 (11:08 +0000)]
First pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024

13 years agoWhile looking if bug0022 was corrected, found that some comments were not updated
julie [Thu, 17 Mar 2011 16:41:08 +0000 (16:41 +0000)]
While looking if bug0022 was corrected, found that some comments were not updated
bug0022 was indeed corrected by Zlatko

13 years agoRemove GO TO STATEMENT in DROTMG/DROTM and SROTMG/SROTM.
julie [Wed, 16 Mar 2011 22:52:01 +0000 (22:52 +0000)]
Remove GO TO STATEMENT in DROTMG/DROTM and SROTMG/SROTM.

Incorporate OLD Testings from ACM Collected algorithms : algorithm 539
for:
   - DROTMG/SROTMG
   - DROTM/SROTM
   - DSDOT / SDSDOT

All BLAS routines are now tested.
Now up to us to change the BLAS but at least this is consistent.

13 years agoRemove the easy GO TO statments....still 13 to remove in drotm.f and 36 in drotmg.f
julie [Thu, 10 Mar 2011 23:33:40 +0000 (23:33 +0000)]
Remove the easy GO TO statments....still 13 to remove in drotm.f and 36 in drotmg.f

13 years agoFormating code, use of * for comment, and follow other LAPACK style
julie [Thu, 10 Mar 2011 19:09:21 +0000 (19:09 +0000)]
Formating code, use of * for comment, and follow other LAPACK style

13 years agoAdded missing Further Details description of V (copied from xLARFT)
james [Tue, 8 Mar 2011 21:01:45 +0000 (21:01 +0000)]
Added missing Further Details description of V (copied from xLARFT)

13 years agoRemove multiple definition of cpotrs
julie [Mon, 7 Mar 2011 15:11:26 +0000 (15:11 +0000)]
Remove multiple definition of cpotrs

13 years agoAdded compiler options needed for HP Fortran
julie [Sun, 6 Mar 2011 20:02:06 +0000 (20:02 +0000)]
Added compiler options needed for HP Fortran

- Chuck Atkins (Kitware)

13 years agoCorrect bug 0076 reported by Intel Team on lapack forum
julie [Fri, 4 Mar 2011 20:44:21 +0000 (20:44 +0000)]
Correct bug 0076 reported by Intel Team on lapack forum
Problem in ?(sy/he)tri2 was found.

For matrices with small N subroutine pass parameter NBMAX to ?(sy/he)tri2x which could be large than N.
In test for this functionality this is hide. Local test subroutine ilaenv.f pass proper values of NB. But if use ordinary ilaenv we will get mistake.

If we implemented

IF ( NBMAX .GE. N ) THEN
MINSIZE = N
ELSE
MINSIZE = (N+NBMAX+1)*(NBMAX+3)
END IF

and will call

IF( NBMAX .GE. N ) THEN
CALL SSYTRI( UPLO, N, A, LDA, IPIV, WORK, INFO )
ELSE
CALL SSYTRI2X( UPLO, N, A, LDA, IPIV, WORK, NBMAX, INFO )
END IF

the problem will be solved.

Some other minor changes on workspace query.

13 years agoCorrect Workspace Query to make sure we return the max of LWORKOPT and LWORKMIN
julie [Thu, 3 Mar 2011 23:14:53 +0000 (23:14 +0000)]
Correct Workspace Query to make sure we return the max of LWORKOPT and LWORKMIN

13 years agoCorrect and Add Error exit tests for *errtz
julie [Thu, 3 Mar 2011 23:13:58 +0000 (23:13 +0000)]
Correct and Add Error exit tests for *errtz
Test the case where N=M (Workspace should be 1), so put LWORK at 0 to make it failed
Test the case where N!=M (Workspace should be more than MAX(1,M)), so put M=2, N=3 ans LWORK at 1 to make it failed

Correct some typos in [dz]drgsx testings.

13 years agoOne more workspace problem fixed...
julie [Thu, 3 Mar 2011 22:33:51 +0000 (22:33 +0000)]
One more workspace problem fixed...

13 years agoThe c generated code should be better now
julie [Thu, 3 Mar 2011 21:09:29 +0000 (21:09 +0000)]
The c generated code should be better now

13 years agoFix some workspace query problems found with the new interface
julie [Thu, 3 Mar 2011 21:02:22 +0000 (21:02 +0000)]
Fix some workspace query problems found with the new interface

13 years agoreplaced "do/exit" scan logic with do while loop
james [Thu, 3 Mar 2011 20:04:03 +0000 (20:04 +0000)]
replaced "do/exit" scan logic with do while loop
fixes bug0075 and bug0035

13 years agoAdded warning suppressions for the Intel compiler on Windows
julie [Mon, 28 Feb 2011 04:11:45 +0000 (04:11 +0000)]
Added warning suppressions for the Intel compiler on Windows

- Chuck Atkins (Kitware)

13 years agoWrapped lines that extended past 72 columns
julie [Mon, 28 Feb 2011 03:42:25 +0000 (03:42 +0000)]
Wrapped lines that extended past 72 columns

- Chuck Atkins (Kitware)

13 years agoCorrected a typo in the output message when setting optimization level
julie [Sat, 26 Feb 2011 21:49:14 +0000 (21:49 +0000)]
Corrected a typo in the output message when setting optimization level

- Chuck Atkins (Kitware)

13 years agoRe-worked type-size checks to search for differnt sizes for different types
julie [Sat, 26 Feb 2011 20:51:01 +0000 (20:51 +0000)]
Re-worked type-size checks to search for differnt sizes for different types

- Chuck Atkins (Kitware)

13 years agoAdded a status message when using the supplied NETLIB BLAS
julie [Sat, 26 Feb 2011 20:50:37 +0000 (20:50 +0000)]
Added a status message when using the supplied NETLIB BLAS

- Chuck Atkins (Kitware)

13 years agoRecude RELEASE optimization level from O>3 to O2
julie [Sat, 26 Feb 2011 20:49:40 +0000 (20:49 +0000)]
Recude RELEASE optimization level from O>3 to O2

- Chuck Atkins (Kitware)

13 years agochanged matrix norm to 'M' (max) from 'I' (inf) for scaling of input matrix
james [Thu, 24 Feb 2011 21:12:23 +0000 (21:12 +0000)]
changed matrix norm to 'M' (max) from 'I' (inf) for scaling of input matrix
and added check for zero norm in xSTERF (the xSTEQR already had this check)
this fixes the issue reported in bug0036

13 years agoremoved confusing comments about JMAX, resolves bug0034
james [Thu, 17 Feb 2011 20:14:14 +0000 (20:14 +0000)]
removed confusing comments about JMAX, resolves bug0034

13 years agoremoved confusing comments about JMAX, resolves bug0034
james [Thu, 17 Feb 2011 20:10:22 +0000 (20:10 +0000)]
removed confusing comments about JMAX, resolves bug0034

13 years agofixes bug0027: removed incorrect description of output of H (they are not in standard...
james [Thu, 17 Feb 2011 16:50:30 +0000 (16:50 +0000)]
fixes bug0027: removed incorrect description of output of H (they are not in standard form)

13 years agoCorrect precision generation problem
julie [Wed, 9 Feb 2011 19:55:15 +0000 (19:55 +0000)]
Correct precision generation problem

13 years agoFix bug 0074
julie [Wed, 9 Feb 2011 19:45:09 +0000 (19:45 +0000)]
Fix bug 0074

Post from Alexander Kobotov (INTEL)
In LAPACK 3.3 ?(SY/HE)SV are updated to use TRS2.
The TRS2 requires WORK(N) to operate, whereas in SV requirements for LWORK just to be >=1.
There is no any check in SV if LWORK>=N, the array WORK just passed as it is to TRS2.
So if LWORK<N a crash could occur while executing TRS2 due to overuse of allocated workspace.

I guess for the case of LWORK<N just previsous Level2 based ?(SY/HE)TRS should be used.

13 years agoFix Bug 0073: [Re] Add dummy routine in timer to avoid compiler problem
julie [Mon, 7 Feb 2011 18:00:27 +0000 (18:00 +0000)]
Fix Bug 0073: [Re] Add dummy routine in timer to avoid compiler problem

Bug reported by John Tellefson to LAPACK Mailing List on February 7th.
Dummy routine was removed in LAPACK 3.3.0 because it seems it was working on all machines.
[Re] Incorporation of the dummy routine and incorporate John's suggested cosmectic changes to allow
an easier change in number of operation to perform.

"modified secondtst.f from lapack-3.2.2 with the old
code that makes it reliably report timing loop results. I also added a
variable TOTAL to simplify changes in setting array size and iterations
which affect total flops and made few other changes.

13 years agoMove CMake files to a subdirectory
julie [Mon, 7 Feb 2011 14:54:28 +0000 (14:54 +0000)]
Move CMake files to a subdirectory

- Brad King (Kitware)

13 years agoLink tmglib test library to lapack
julie [Mon, 7 Feb 2011 14:54:24 +0000 (14:54 +0000)]
Link tmglib test library to lapack

The library uses symbols from lapack and so should link to it.
Previously this was not noticed because static libraries do not
really link and the test executables link to lapack anyway.

- Brad King (Kitware)

13 years agoOrganize binaries into bin and lib directories
julie [Mon, 7 Feb 2011 14:54:21 +0000 (14:54 +0000)]
Organize binaries into bin and lib directories

Put executables under 'bin' and libraries under 'lib'.  This is needed
on Windows with shared library builds to put the .dll and .exe files
together for running tests without setting PATH.

- Brad King (Kitware)

13 years agoRestore CMake packaging removed in r858
julie [Mon, 7 Feb 2011 14:54:19 +0000 (14:54 +0000)]
Restore CMake packaging removed in r858

- Brad King (Kitware)

13 years agoRevert "Restore CMake packaging removed in r858"
julie [Mon, 7 Feb 2011 14:54:16 +0000 (14:54 +0000)]
Revert "Restore CMake packaging removed in r858"

This reverts commit e725dc754dd87ed45f21ba6981fca6dbc0530aa3.

13 years agoRevert "Organize binaries into bin and lib directories"
julie [Mon, 7 Feb 2011 14:54:14 +0000 (14:54 +0000)]
Revert "Organize binaries into bin and lib directories"

This reverts commit e3b5f531553c176a009f41b6a312b0c11b955fd6.

13 years agoRevert "Link tmglib test library to lapack"
julie [Mon, 7 Feb 2011 14:54:11 +0000 (14:54 +0000)]
Revert "Link tmglib test library to lapack"

This reverts commit 5070d8d0b7690496561a26f5849d6f80a8727d93.

13 years agoRevert "Move CMake files to a subdirectory"
julie [Mon, 7 Feb 2011 14:54:09 +0000 (14:54 +0000)]
Revert "Move CMake files to a subdirectory"

This reverts commit bc425c443232cf9737afcce38753f79963d31e0a.

13 years agoMove CMake files to a subdirectory
julie [Mon, 7 Feb 2011 14:37:01 +0000 (14:37 +0000)]
Move CMake files to a subdirectory

13 years agoLink tmglib test library to lapack
julie [Mon, 7 Feb 2011 14:36:57 +0000 (14:36 +0000)]
Link tmglib test library to lapack

The library uses symbols from lapack and so should link to it.
Previously this was not noticed because static libraries do not
really link and the test executables link to lapack anyway.

13 years agoOrganize binaries into bin and lib directories
julie [Mon, 7 Feb 2011 14:36:54 +0000 (14:36 +0000)]
Organize binaries into bin and lib directories

Put executables under 'bin' and libraries under 'lib'.  This is needed
on Windows with shared library builds to put the .dll and .exe files
together for running tests without setting PATH.

13 years agoRestore CMake packaging removed in r858
julie [Mon, 7 Feb 2011 14:36:50 +0000 (14:36 +0000)]
Restore CMake packaging removed in r858

13 years agoFix bug 0072 reported by Inge Gutheil on February 3rd 2011
julie [Thu, 3 Feb 2011 18:28:38 +0000 (18:28 +0000)]
Fix bug 0072 reported by Inge Gutheil on February 3rd 2011

From Inge:
In the tests,
TESTING/EIG/cchkhs.f (similar in zchkhs.f) in line 492 N = NN ( JSIZE )
is set. If NN( JSIZE ) is zero, CCOPY (N-1,....) gives an error message
and the test breaks. In the real versions dchkhs after N = NN ( JSIZE )
there is
IF ( N .EQ. 0 ) GO TO 270, so there DCOPY does not exit with an error
message.
I think the same statement should be added in the complex version to
avoid copying -1 element.

13 years agoChange from Christof (Voemel) sent on Mon, 31 Jan 2011.
langou [Mon, 31 Jan 2011 18:27:33 +0000 (18:27 +0000)]
Change from Christof (Voemel) sent on Mon, 31 Jan 2011.

*********************************************************************
change
         IF( USEDQD ) THEN
*           The initial SIGMA was to the outer end of the spectrum
*           the matrix is definite and we need not retreat.
            TAU = SPDIAM*EPS*N + TWO*PIVMIN
         ELSE
into
         IF( USEDQD ) THEN
*           The initial SIGMA was to the outer end of the spectrum
*           the matrix is definite and we need not retreat.
            TAU = SPDIAM*EPS*N + TWO*PIVMIN
            TAU = MAX( TAU,TWO*EPS*ABS(SIGMA) )
         ELSE

that is, with an added line to prevent TAU from becoming zero
            TAU = MAX( TAU,TWO*EPS*ABS(SIGMA) )
*********************************************************************

Thanks Christof!

13 years agoFix bug0070 : IWORK dimension is incorrect
julie [Mon, 31 Jan 2011 16:14:51 +0000 (16:14 +0000)]
Fix bug0070 :  IWORK dimension is incorrect
  bug report sent by nmozarto on Jan 27th
  see forum topic 2156

From Brian:
 The correct dimension is (M-R), in which R is the smallest of P, M-P, Q, and M-Q.
 An M-by-M matrix is partitioned into a 2-by-2 block structure. The
 dimensions of the blocks are P-by-Q, P-by-(M-Q), (M-P)-by-Q, and
 (M-P)-by-(M-Q). IWORK is involved in the permutation of rows or
 columns of these blocks, and it is the largest block, in terms of
 number of rows or columns, that determines the size of IWORK.

13 years agoCorrect bug0069
julie [Fri, 28 Jan 2011 23:04:40 +0000 (23:04 +0000)]
Correct bug0069
Bug was sent by nmozarto on Jan 27th (see forum topic 2156)

Problem in new function ?SYTRI2 was found: the part of A below the diagonal is changed in the case UPLO='U' .
But in the description of arguments
If UPLO = 'U', the upper triangular part of the
inverse is formed and the part of A below the diagonal is not
referenced; if UPLO = 'L' the lower triangular part of the
inverse is formed and the part of A above the diagonal is
not referenced.

These elements zeroized after calling ?GEMM function in ?SYTRI2X.
CALL SGEMM('T','N',NNB,NNB,CUT,ONE,A(1,CUT+1),LDA,
$ WORK,N+NB+1, ZERO, A(CUT+1,CUT+1), LDA)

13 years agoSame for hetrs2 routines
julie [Thu, 27 Jan 2011 14:25:04 +0000 (14:25 +0000)]
Same for hetrs2 routines

13 years agoFollow-up: fix for new routines to handle LDA
julie [Thu, 27 Jan 2011 14:21:35 +0000 (14:21 +0000)]
Follow-up: fix for new routines to handle LDA
See bug0067

13 years ago1) Je viens de modifier le fichier timer_dsytri2.c de telle sorte a ce qu'il
langou [Wed, 26 Jan 2011 14:02:05 +0000 (14:02 +0000)]
1) Je viens de modifier le fichier timer_dsytri2.c de telle sorte a ce qu'il
traite les matrices avec des lda.
Exemple:

> gfortran timer_dsytri2.c /Users/langou/Documents/lapack-dev-trunk.svn/lapack_LINUX.a -framework vecLib
> ./a.out -U -lda 1079 -n 100
DSYTRI2 info=0, check=1.891597e-14
Dispo dans ldlt.git. Marche bien. Upper ou Lower. LDA different de N ou pas. Le
check est
|| A . inv(A) - I || / || A ||
(Je pense que il faut diviser par || A || * || inv(A) ||, pas tres important.)

2) Du coup, dans ce commit, je rectifie tous les codes. Cela compile. Je n'ai
pas verifie les testings. Je n'ai pas ecrit de driver pour verifier. (Le driver
n'est que pour la version DOUBLE.)

JL

13 years agobon maintenant, ca compile ... (desole)
langou [Wed, 26 Jan 2011 05:49:27 +0000 (05:49 +0000)]
bon maintenant, ca compile ... (desole)

13 years agoJ'ai lu le code vite fait. Comme ca, ca devrait marcher. Je n'ai pas verifie.
langou [Wed, 26 Jan 2011 05:47:42 +0000 (05:47 +0000)]
J'ai lu le code vite fait. Comme ca, ca devrait marcher. Je n'ai pas verifie.
Je n'ai meme pas essaye de compiler en fait. En gros mettre LDA en parametre,
definir A comme etant A( LDA, N ) et pour l'increment du swap mettre LDA a la
place de N quand tu travailles sur les lignes. Donc c'est tout con.

Donc il faut ecrire un driver qui teste ca avec une LDA differente de N ... Et
appliquer les modifications aux autres routines CHE, CSY, ZHE, ZSY et SSY.

Sinon pourquoi dans DSYSWAPR, est ce que tu utilises DSWAP de temps en temps,
et de temps en temps tu le fais a la main?  Je n'ai pas regarde les details.

JL.

13 years agoInclude Zlatko changes following Paul Roberts from NAG report
julie [Tue, 25 Jan 2011 22:02:14 +0000 (22:02 +0000)]
Include Zlatko changes following Paul Roberts from NAG report

Message from Zlatko:

The following changes are made in the
current version of the code:

1. in dgejsv, sgejsv:

- a typo in checking the parameters LSVEC.OR.LSVEC
  has been changed to LSVEC.OR.RSVEC
- the length of WORK, LWORK, its description
  and the minimal length for different JOBs
  have been revised and corrected
- in a call to xGESVJ with WORK(N+1), the length
  of the workspace is set to the correct value of
  LWORK-N, instead of the incorrect LWORK
- a missing RETURN after a call to XERBLA
  has been inserted
- In the case of zero matrix on input, IWORK(3)
  is set to ZERO, to correspond to the description
  of IWORK in other nontrivial cases
2. in all routines
- simple editing so that the single and
  the corresponding double routines have matching
  lines in the source codes.

13 years agoCorrect multiple definition in the lapack lib (See forum topic 2099 - Bruno)
julie [Mon, 24 Jan 2011 20:05:41 +0000 (20:05 +0000)]
Correct multiple definition in the lapack lib (See forum topic 2099 - Bruno)

13 years agoAdding new hetri routines
julie [Thu, 20 Jan 2011 17:33:46 +0000 (17:33 +0000)]
Adding new hetri routines

13 years agoRemove the last compiler directive
julie [Thu, 13 Jan 2011 12:05:18 +0000 (12:05 +0000)]
Remove the last compiler directive

13 years agoRemove outdated compiler directives
julie [Wed, 12 Jan 2011 17:34:06 +0000 (17:34 +0000)]
Remove outdated compiler directives

13 years agoSuppress warning because of undefined variable
julie [Tue, 11 Jan 2011 14:10:14 +0000 (14:10 +0000)]
Suppress warning because of undefined variable

13 years agoSuppress a warning from SUN compilers when building exclusively for your native
julie [Sun, 9 Jan 2011 19:30:36 +0000 (19:30 +0000)]
Suppress a warning from SUN compilers when building exclusively for your native
platform.

- Chuck Atkins (Kitware)

13 years agoSuppress harmless linker warnings on IRIX + MIPSpro
julie [Sat, 8 Jan 2011 16:13:32 +0000 (16:13 +0000)]
Suppress harmless linker warnings on IRIX + MIPSpro

- Chuck Atkins (Kitware)

13 years agoRemoved output truncation on ctest submissions.
julie [Wed, 5 Jan 2011 03:28:27 +0000 (03:28 +0000)]
Removed output truncation on ctest submissions.

- Chuck Atkins (Kitware)

13 years agoFix object needed for targets in CMAKE build
julie [Tue, 4 Jan 2011 20:03:34 +0000 (20:03 +0000)]
Fix object needed for targets in CMAKE build

13 years agominor fix following compiler warnings: remove unused variable and change digit output...
julie [Tue, 4 Jan 2011 19:59:16 +0000 (19:59 +0000)]
minor fix following compiler warnings: remove unused variable and change digit output format

13 years agominor fix: add comment to warn that vector and matrix arguments are not referenced...
julie [Tue, 4 Jan 2011 19:58:09 +0000 (19:58 +0000)]
minor fix: add comment to warn that vector and matrix arguments are not referenced when N = 0, or M = 0

13 years agoRe-worked the type size check to be much more compact
julie [Fri, 31 Dec 2010 01:33:24 +0000 (01:33 +0000)]
Re-worked the type size check to be much more compact

- Chuck Atkins (Kitware)

13 years agoAdded acml_mv to the ACML link libraries
julie [Thu, 30 Dec 2010 23:11:32 +0000 (23:11 +0000)]
Added acml_mv to the ACML link libraries

- Chuck Atkins (Kitware)

13 years ago1. Adjusted regex checks in CheckLAPACKCompilerFlags.cmake to be more compact
julie [Thu, 30 Dec 2010 22:14:10 +0000 (22:14 +0000)]
1.  Adjusted regex checks in CheckLAPACKCompilerFlags.cmake to be more compact
2.  Added CheckFortranIntSize.cmake to auto-determine the default integer size
3.  Added more auto-detection features using integer size and compiler checks in FindBLAS.cmake

- Chuck Atkins (Kitware)

14 years agoDisplay the test's output to stdout after a run so it can be viewed on the dashboard.
julie [Fri, 24 Dec 2010 21:38:46 +0000 (21:38 +0000)]
Display the test's output to stdout after a run so it can be viewed on the dashboard.

- Chuck Atkins (Kitware)

14 years ago1. Fixed a typo in FindBLAS for SGI (SCCL should have been SCSL) and 2. Fixed compi...
julie [Fri, 24 Dec 2010 15:29:47 +0000 (15:29 +0000)]
1.  Fixed a typo in FindBLAS for SGI (SCCL should have been SCSL) and 2.  Fixed compiler flag check for SunPro to use ftrap=%none instead of ftrap=none.

- Chuck Atkins (Kitware)

14 years ago1. Added support for both IBM XL Fortran and IBM VisualAge Fortran
julie [Thu, 23 Dec 2010 02:32:05 +0000 (02:32 +0000)]
1.  Added support for both IBM XL Fortran and IBM VisualAge Fortran
2.  Added libblas to ESSL libs to cover missing functionality from ESSL
3.  Renamed ESSL_6464, ESSL_SMP, and ESSL_SMP_6464 to ESSL6464, ESSLSMP, and ESSLSMP6464 respectively to mimic the actual library names.
4.  Added copyright notice to CheckLAPACKCompilerFlags

- Chuck Atkins (Kitware)

14 years agoFixed a bug in the IBM compiler options (-ffixed-form -> -qfixed).
julie [Thu, 23 Dec 2010 00:50:18 +0000 (00:50 +0000)]
Fixed a bug in the IBM compiler options (-ffixed-form -> -qfixed).

- Chuck Atkins (Kitware)

14 years agoAdded checks disable floating point exceptions for GNU, Intel, Sun, and IBM compilers...
julie [Thu, 23 Dec 2010 00:39:05 +0000 (00:39 +0000)]
Added checks disable floating point exceptions for GNU, Intel, Sun, and IBM compilers and to enforce fixed-form for IBM compilers.

- Chuck Atkins (Kitware)

14 years agoFixed FindBLAS for Sun Performance Library
julie [Tue, 21 Dec 2010 02:11:49 +0000 (02:11 +0000)]
Fixed FindBLAS for Sun Performance Library

- Chuck Atkins (Kitware)

14 years agoAdded a CTestCustom.cmake file to filter unwanted files from code coverage and suppre...
julie [Mon, 20 Dec 2010 22:54:32 +0000 (22:54 +0000)]
Added a CTestCustom.cmake file to filter unwanted files from code coverage and suppress harmless warnings on the dashboard.

- Chuck Atkins (Kitware)

14 years agoAdding Peng Du's name in the authors of the variants.
langou [Fri, 17 Dec 2010 00:09:40 +0000 (00:09 +0000)]
Adding Peng Du's name in the authors of the variants.

14 years ago1. Corrected the way linker flags from FindBLAS are being used
julie [Thu, 16 Dec 2010 15:12:27 +0000 (15:12 +0000)]
1.  Corrected the way linker flags from FindBLAS are being used
2.  Changed LOOK_FOR_OPTIMIZED_BLAS to USE_OPTIMIZED_BLAS per CMake convention
3.  Changed BLAS_VENDORS to use library names instead of manufacturer names:
    From: BLAS_VENDORS=AMD;IBM;INTEL;HP;SUN etc
      To: BLAS_VENDORS=ACML;ESSL;MKL;VECLIB;PERFLIB
    This seems to be a more intuitive interface
4.  Added support for Intel MKL single dynamic library

- Chuck Atkins (Kitware)

14 years ago Merge some of Sebastien Fabbro's patch for CMAKE
julie [Wed, 15 Dec 2010 18:44:42 +0000 (18:44 +0000)]
 Merge some of Sebastien Fabbro's patch for CMAKE

14 years agoAdded compiler flags from FindBLAS.
julie [Wed, 15 Dec 2010 18:42:50 +0000 (18:42 +0000)]
Added compiler flags from FindBLAS.