floraachy
/
OpenBLAS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Lapack bug114: replaced cgesvd.f and zgesvd.f

This commit is contained in:
wernsaar 2014-07-08 10:21:10 +02:00
parent dc6b809f15
commit 2081f6e8ff
2 changed files with 98 additions and 96 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

97
lapack-netlib/SRC/cgesvd.f
View File

@ -321,24 +321,24 @@
* *
MNTHR = ILAENV( 6, 'CGESVD', JOBU // JOBVT, M, N, 0, 0 ) MNTHR = ILAENV( 6, 'CGESVD', JOBU // JOBVT, M, N, 0, 0 )
* Compute space needed for CGEQRF * Compute space needed for CGEQRF
CALL CGEQRF( M, N, A, LDA, DUM(1), DUM(1), -1, IERR ) CALL CGEQRF( M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_CGEQRF=DUM(1) LWORK_CGEQRF=CDUM(1)
* Compute space needed for CUNGQR * Compute space needed for CUNGQR
CALL CUNGQR( M, N, N, A, LDA, DUM(1), DUM(1), -1, IERR ) CALL CUNGQR( M, N, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_CUNGQR_N=DUM(1) LWORK_CUNGQR_N=CDUM(1)
CALL CUNGQR( M, M, N, A, LDA, DUM(1), DUM(1), -1, IERR ) CALL CUNGQR( M, M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_CUNGQR_M=DUM(1) LWORK_CUNGQR_M=CDUM(1)
* Compute space needed for CGEBRD * Compute space needed for CGEBRD
CALL CGEBRD( N, N, A, LDA, S, DUM(1), DUM(1), CALL CGEBRD( N, N, A, LDA, S, DUM(1), CDUM(1),
$ DUM(1), DUM(1), -1, IERR ) $ CDUM(1), CDUM(1), -1, IERR )
LWORK_CGEBRD=DUM(1) LWORK_CGEBRD=CDUM(1)
* Compute space needed for CUNGBR * Compute space needed for CUNGBR
CALL CUNGBR( 'P', N, N, N, A, LDA, DUM(1), CALL CUNGBR( 'P', N, N, N, A, LDA, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_CUNGBR_P=DUM(1) LWORK_CUNGBR_P=CDUM(1)
CALL CUNGBR( 'Q', N, N, N, A, LDA, DUM(1), CALL CUNGBR( 'Q', N, N, N, A, LDA, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_CUNGBR_Q=DUM(1) LWORK_CUNGBR_Q=CDUM(1)
* *
MNTHR = ILAENV( 6, 'CGESVD', JOBU // JOBVT, M, N, 0, 0 ) MNTHR = ILAENV( 6, 'CGESVD', JOBU // JOBVT, M, N, 0, 0 )
IF( M.GE.MNTHR ) THEN IF( M.GE.MNTHR ) THEN
@ -444,20 +444,20 @@
* *
* Path 10 (M at least N, but not much larger) * Path 10 (M at least N, but not much larger)
* *
CALL CGEBRD( M, N, A, LDA, S, DUM(1), DUM(1), CALL CGEBRD( M, N, A, LDA, S, DUM(1), CDUM(1),
$ DUM(1), DUM(1), -1, IERR ) $ CDUM(1), CDUM(1), -1, IERR )
LWORK_CGEBRD=DUM(1) LWORK_CGEBRD=CDUM(1)
MAXWRK = 2*N + LWORK_CGEBRD MAXWRK = 2*N + LWORK_CGEBRD
IF( WNTUS .OR. WNTUO ) THEN IF( WNTUS .OR. WNTUO ) THEN
CALL CUNGBR( 'Q', M, N, N, A, LDA, DUM(1), CALL CUNGBR( 'Q', M, N, N, A, LDA, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_CUNGBR_Q=DUM(1) LWORK_CUNGBR_Q=CDUM(1)
MAXWRK = MAX( MAXWRK, 2*N+LWORK_CUNGBR_Q ) MAXWRK = MAX( MAXWRK, 2*N+LWORK_CUNGBR_Q )
END IF END IF
IF( WNTUA ) THEN IF( WNTUA ) THEN
CALL CUNGBR( 'Q', M, M, N, A, LDA, DUM(1), CALL CUNGBR( 'Q', M, M, N, A, LDA, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_CUNGBR_Q=DUM(1) LWORK_CUNGBR_Q=CDUM(1)
MAXWRK = MAX( MAXWRK, 2*N+LWORK_CUNGBR_Q ) MAXWRK = MAX( MAXWRK, 2*N+LWORK_CUNGBR_Q )
END IF END IF
IF( .NOT.WNTVN ) THEN IF( .NOT.WNTVN ) THEN
@ -471,25 +471,26 @@
* *
MNTHR = ILAENV( 6, 'CGESVD', JOBU // JOBVT, M, N, 0, 0 ) MNTHR = ILAENV( 6, 'CGESVD', JOBU // JOBVT, M, N, 0, 0 )
* Compute space needed for CGELQF * Compute space needed for CGELQF
CALL CGELQF( M, N, A, LDA, DUM(1), DUM(1), -1, IERR ) CALL CGELQF( M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_CGELQF=DUM(1) LWORK_CGELQF=CDUM(1)
* Compute space needed for CUNGLQ * Compute space needed for CUNGLQ
CALL CUNGLQ( N, N, M, DUM(1), N, DUM(1), DUM(1), -1, IERR ) CALL CUNGLQ( N, N, M, CDUM(1), N, CDUM(1), CDUM(1), -1,
LWORK_CUNGLQ_N=DUM(1) $ IERR )
CALL CUNGLQ( M, N, M, A, LDA, DUM(1), DUM(1), -1, IERR ) LWORK_CUNGLQ_N=CDUM(1)
LWORK_CUNGLQ_M=DUM(1) CALL CUNGLQ( M, N, M, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_CUNGLQ_M=CDUM(1)
* Compute space needed for CGEBRD * Compute space needed for CGEBRD
CALL CGEBRD( M, M, A, LDA, S, DUM(1), DUM(1), CALL CGEBRD( M, M, A, LDA, S, DUM(1), CDUM(1),
$ DUM(1), DUM(1), -1, IERR ) $ CDUM(1), CDUM(1), -1, IERR )
LWORK_CGEBRD=DUM(1) LWORK_CGEBRD=CDUM(1)
* Compute space needed for CUNGBR P * Compute space needed for CUNGBR P
CALL CUNGBR( 'P', M, M, M, A, N, DUM(1), CALL CUNGBR( 'P', M, M, M, A, N, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_CUNGBR_P=DUM(1) LWORK_CUNGBR_P=CDUM(1)
* Compute space needed for CUNGBR Q * Compute space needed for CUNGBR Q
CALL CUNGBR( 'Q', M, M, M, A, N, DUM(1), CALL CUNGBR( 'Q', M, M, M, A, N, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_CUNGBR_Q=DUM(1) LWORK_CUNGBR_Q=CDUM(1)
IF( N.GE.MNTHR ) THEN IF( N.GE.MNTHR ) THEN
IF( WNTVN ) THEN IF( WNTVN ) THEN
* *
@ -593,21 +594,21 @@
* *
* Path 10t(N greater than M, but not much larger) * Path 10t(N greater than M, but not much larger)
* *
CALL CGEBRD( M, N, A, LDA, S, DUM(1), DUM(1), CALL CGEBRD( M, N, A, LDA, S, DUM(1), CDUM(1),
$ DUM(1), DUM(1), -1, IERR ) $ CDUM(1), CDUM(1), -1, IERR )
LWORK_CGEBRD=DUM(1) LWORK_CGEBRD=CDUM(1)
MAXWRK = 2*M + LWORK_CGEBRD MAXWRK = 2*M + LWORK_CGEBRD
IF( WNTVS .OR. WNTVO ) THEN IF( WNTVS .OR. WNTVO ) THEN
* Compute space needed for CUNGBR P * Compute space needed for CUNGBR P
CALL CUNGBR( 'P', M, N, M, A, N, DUM(1), CALL CUNGBR( 'P', M, N, M, A, N, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_CUNGBR_P=DUM(1) LWORK_CUNGBR_P=CDUM(1)
MAXWRK = MAX( MAXWRK, 2*M+LWORK_CUNGBR_P ) MAXWRK = MAX( MAXWRK, 2*M+LWORK_CUNGBR_P )
END IF END IF
IF( WNTVA ) THEN IF( WNTVA ) THEN
CALL CUNGBR( 'P', N, N, M, A, N, DUM(1), CALL CUNGBR( 'P', N, N, M, A, N, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_CUNGBR_P=DUM(1) LWORK_CUNGBR_P=CDUM(1)
MAXWRK = MAX( MAXWRK, 2*M+LWORK_CUNGBR_P ) MAXWRK = MAX( MAXWRK, 2*M+LWORK_CUNGBR_P )
END IF END IF
IF( .NOT.WNTUN ) THEN IF( .NOT.WNTUN ) THEN

97
lapack-netlib/SRC/zgesvd.f
View File

@ -321,24 +321,24 @@
* *
MNTHR = ILAENV( 6, 'ZGESVD', JOBU // JOBVT, M, N, 0, 0 ) MNTHR = ILAENV( 6, 'ZGESVD', JOBU // JOBVT, M, N, 0, 0 )
* Compute space needed for ZGEQRF * Compute space needed for ZGEQRF
CALL ZGEQRF( M, N, A, LDA, DUM(1), DUM(1), -1, IERR ) CALL ZGEQRF( M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_ZGEQRF=DUM(1) LWORK_ZGEQRF=CDUM(1)
* Compute space needed for ZUNGQR * Compute space needed for ZUNGQR
CALL ZUNGQR( M, N, N, A, LDA, DUM(1), DUM(1), -1, IERR ) CALL ZUNGQR( M, N, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_ZUNGQR_N=DUM(1) LWORK_ZUNGQR_N=CDUM(1)
CALL ZUNGQR( M, M, N, A, LDA, DUM(1), DUM(1), -1, IERR ) CALL ZUNGQR( M, M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_ZUNGQR_M=DUM(1) LWORK_ZUNGQR_M=CDUM(1)
* Compute space needed for ZGEBRD * Compute space needed for ZGEBRD
CALL ZGEBRD( N, N, A, LDA, S, DUM(1), DUM(1), CALL ZGEBRD( N, N, A, LDA, S, DUM(1), CDUM(1),
$ DUM(1), DUM(1), -1, IERR ) $ CDUM(1), CDUM(1), -1, IERR )
LWORK_ZGEBRD=DUM(1) LWORK_ZGEBRD=CDUM(1)
* Compute space needed for ZUNGBR * Compute space needed for ZUNGBR
CALL ZUNGBR( 'P', N, N, N, A, LDA, DUM(1), CALL ZUNGBR( 'P', N, N, N, A, LDA, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_ZUNGBR_P=DUM(1) LWORK_ZUNGBR_P=CDUM(1)
CALL ZUNGBR( 'Q', N, N, N, A, LDA, DUM(1), CALL ZUNGBR( 'Q', N, N, N, A, LDA, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_ZUNGBR_Q=DUM(1) LWORK_ZUNGBR_Q=CDUM(1)
* *
IF( M.GE.MNTHR ) THEN IF( M.GE.MNTHR ) THEN
IF( WNTUN ) THEN IF( WNTUN ) THEN
@ -443,20 +443,20 @@
* *
* Path 10 (M at least N, but not much larger) * Path 10 (M at least N, but not much larger)
* *
CALL ZGEBRD( M, N, A, LDA, S, DUM(1), DUM(1), CALL ZGEBRD( M, N, A, LDA, S, DUM(1), CDUM(1),
$ DUM(1), DUM(1), -1, IERR ) $ CDUM(1), CDUM(1), -1, IERR )
LWORK_ZGEBRD=DUM(1) LWORK_ZGEBRD=CDUM(1)
MAXWRK = 2*N + LWORK_ZGEBRD MAXWRK = 2*N + LWORK_ZGEBRD
IF( WNTUS .OR. WNTUO ) THEN IF( WNTUS .OR. WNTUO ) THEN
CALL ZUNGBR( 'Q', M, N, N, A, LDA, DUM(1), CALL ZUNGBR( 'Q', M, N, N, A, LDA, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_ZUNGBR_Q=DUM(1) LWORK_ZUNGBR_Q=CDUM(1)
MAXWRK = MAX( MAXWRK, 2*N+LWORK_ZUNGBR_Q ) MAXWRK = MAX( MAXWRK, 2*N+LWORK_ZUNGBR_Q )
END IF END IF
IF( WNTUA ) THEN IF( WNTUA ) THEN
CALL ZUNGBR( 'Q', M, M, N, A, LDA, DUM(1), CALL ZUNGBR( 'Q', M, M, N, A, LDA, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_ZUNGBR_Q=DUM(1) LWORK_ZUNGBR_Q=CDUM(1)
MAXWRK = MAX( MAXWRK, 2*N+LWORK_ZUNGBR_Q ) MAXWRK = MAX( MAXWRK, 2*N+LWORK_ZUNGBR_Q )
END IF END IF
IF( .NOT.WNTVN ) THEN IF( .NOT.WNTVN ) THEN
@ -470,25 +470,26 @@
* *
MNTHR = ILAENV( 6, 'ZGESVD', JOBU // JOBVT, M, N, 0, 0 ) MNTHR = ILAENV( 6, 'ZGESVD', JOBU // JOBVT, M, N, 0, 0 )
* Compute space needed for ZGELQF * Compute space needed for ZGELQF
CALL ZGELQF( M, N, A, LDA, DUM(1), DUM(1), -1, IERR ) CALL ZGELQF( M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_ZGELQF=DUM(1) LWORK_ZGELQF=CDUM(1)
* Compute space needed for ZUNGLQ * Compute space needed for ZUNGLQ
CALL ZUNGLQ( N, N, M, DUM(1), N, DUM(1), DUM(1), -1, IERR ) CALL ZUNGLQ( N, N, M, CDUM(1), N, CDUM(1), CDUM(1), -1,
LWORK_ZUNGLQ_N=DUM(1) $ IERR )
CALL ZUNGLQ( M, N, M, A, LDA, DUM(1), DUM(1), -1, IERR ) LWORK_ZUNGLQ_N=CDUM(1)
LWORK_ZUNGLQ_M=DUM(1) CALL ZUNGLQ( M, N, M, A, LDA, CDUM(1), CDUM(1), -1, IERR )
LWORK_ZUNGLQ_M=CDUM(1)
* Compute space needed for ZGEBRD * Compute space needed for ZGEBRD
CALL ZGEBRD( M, M, A, LDA, S, DUM(1), DUM(1), CALL ZGEBRD( M, M, A, LDA, S, DUM(1), CDUM(1),
$ DUM(1), DUM(1), -1, IERR ) $ CDUM(1), CDUM(1), -1, IERR )
LWORK_ZGEBRD=DUM(1) LWORK_ZGEBRD=CDUM(1)
* Compute space needed for ZUNGBR P * Compute space needed for ZUNGBR P
CALL ZUNGBR( 'P', M, M, M, A, N, DUM(1), CALL ZUNGBR( 'P', M, M, M, A, N, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_ZUNGBR_P=DUM(1) LWORK_ZUNGBR_P=CDUM(1)
* Compute space needed for ZUNGBR Q * Compute space needed for ZUNGBR Q
CALL ZUNGBR( 'Q', M, M, M, A, N, DUM(1), CALL ZUNGBR( 'Q', M, M, M, A, N, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_ZUNGBR_Q=DUM(1) LWORK_ZUNGBR_Q=CDUM(1)
IF( N.GE.MNTHR ) THEN IF( N.GE.MNTHR ) THEN
IF( WNTVN ) THEN IF( WNTVN ) THEN
* *
@ -592,21 +593,21 @@
* *
* Path 10t(N greater than M, but not much larger) * Path 10t(N greater than M, but not much larger)
* *
CALL ZGEBRD( M, N, A, LDA, S, DUM(1), DUM(1), CALL ZGEBRD( M, N, A, LDA, S, DUM(1), CDUM(1),
$ DUM(1), DUM(1), -1, IERR ) $ CDUM(1), CDUM(1), -1, IERR )
LWORK_ZGEBRD=DUM(1) LWORK_ZGEBRD=CDUM(1)
MAXWRK = 2*M + LWORK_ZGEBRD MAXWRK = 2*M + LWORK_ZGEBRD
IF( WNTVS .OR. WNTVO ) THEN IF( WNTVS .OR. WNTVO ) THEN
* Compute space needed for ZUNGBR P * Compute space needed for ZUNGBR P
CALL ZUNGBR( 'P', M, N, M, A, N, DUM(1), CALL ZUNGBR( 'P', M, N, M, A, N, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_ZUNGBR_P=DUM(1) LWORK_ZUNGBR_P=CDUM(1)
MAXWRK = MAX( MAXWRK, 2*M+LWORK_ZUNGBR_P ) MAXWRK = MAX( MAXWRK, 2*M+LWORK_ZUNGBR_P )
END IF END IF
IF( WNTVA ) THEN IF( WNTVA ) THEN
CALL ZUNGBR( 'P', N, N, M, A, N, DUM(1), CALL ZUNGBR( 'P', N, N, M, A, N, CDUM(1),
$ DUM(1), -1, IERR ) $ CDUM(1), -1, IERR )
LWORK_ZUNGBR_P=DUM(1) LWORK_ZUNGBR_P=CDUM(1)
MAXWRK = MAX( MAXWRK, 2*M+LWORK_ZUNGBR_P ) MAXWRK = MAX( MAXWRK, 2*M+LWORK_ZUNGBR_P )
END IF END IF
IF( .NOT.WNTUN ) THEN IF( .NOT.WNTUN ) THEN