Rev 64299 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed
diff -u keep/dlapack.f ./dlapack.f--- keep/dlapack.f 2013-11-25 09:13:35.000000000 +0000+++ ./dlapack.f 2013-11-25 09:21:56.000000000 +0000@@ -21126,7 +21126,8 @@** Compute Householder transform when N=1*- CALL DLARFG( M, A, A( MIN( 2, M ), 1 ), 1, T )+* R change to stop gfortran warning+ CALL DLARFG( M, A(1,1), A( MIN( 2, M ), 1 ), 1, T(1,1) )*ELSE*diff -u keep/cmplx.f ./cmplx.f--- keep/cmplx.f 2013-11-25 09:09:10.000000000 +0000+++ ./cmplx.f 2013-11-25 09:11:47.000000000 +0000@@ -6048,24 +6048,24 @@*MNTHR = ILAENV( 6, 'ZGESVD', JOBU // JOBVT, M, N, 0, 0 )* Compute space needed for ZGEQRF- CALL ZGEQRF( M, N, A, LDA, DUM(1), DUM(1), -1, IERR )- LWORK_ZGEQRF=DUM(1)+ CALL ZGEQRF( M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZGEQRF=CDUM(1)* Compute space needed for ZUNGQR- CALL ZUNGQR( M, N, N, A, LDA, DUM(1), DUM(1), -1, IERR )- LWORK_ZUNGQR_N=DUM(1)- CALL ZUNGQR( M, M, N, A, LDA, DUM(1), DUM(1), -1, IERR )- LWORK_ZUNGQR_M=DUM(1)+ CALL ZUNGQR( M, N, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZUNGQR_N=CDUM(1)+ CALL ZUNGQR( M, M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZUNGQR_M=CDUM(1)* Compute space needed for ZGEBRD- CALL ZGEBRD( N, N, A, LDA, S, DUM(1), DUM(1),- $ DUM(1), DUM(1), -1, IERR )- LWORK_ZGEBRD=DUM(1)+ CALL ZGEBRD( N, N, A, LDA, S, DUM(1), CDUM(1),+ $ CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZGEBRD=CDUM(1)* Compute space needed for ZUNGBR- CALL ZUNGBR( 'P', N, N, N, A, LDA, DUM(1),- $ DUM(1), -1, IERR )- LWORK_ZUNGBR_P=DUM(1)- CALL ZUNGBR( 'Q', N, N, N, A, LDA, DUM(1),- $ DUM(1), -1, IERR )- LWORK_ZUNGBR_Q=DUM(1)+ CALL ZUNGBR( 'P', N, N, N, A, LDA, CDUM(1),+ $ CDUM(1), -1, IERR )+ LWORK_ZUNGBR_P=CDUM(1)+ CALL ZUNGBR( 'Q', N, N, N, A, LDA, CDUM(1),+ $ CDUM(1), -1, IERR )+ LWORK_ZUNGBR_Q=CDUM(1)*IF( M.GE.MNTHR ) THENIF( WNTUN ) THEN@@ -6170,20 +6170,20 @@** Path 10 (M at least N, but not much larger)*- CALL ZGEBRD( M, N, A, LDA, S, DUM(1), DUM(1),- $ DUM(1), DUM(1), -1, IERR )- LWORK_ZGEBRD=DUM(1)+ CALL ZGEBRD( M, N, A, LDA, S, DUM(1), CDUM(1),+ $ CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZGEBRD=CDUM(1)MAXWRK = 2*N + LWORK_ZGEBRDIF( WNTUS .OR. WNTUO ) THEN- CALL ZUNGBR( 'Q', M, N, N, A, LDA, DUM(1),- $ DUM(1), -1, IERR )- LWORK_ZUNGBR_Q=DUM(1)+ CALL ZUNGBR( 'Q', M, N, N, A, LDA, CDUM(1),+ $ CDUM(1), -1, IERR )+ LWORK_ZUNGBR_Q=CDUM(1)MAXWRK = MAX( MAXWRK, 2*N+LWORK_ZUNGBR_Q )END IFIF( WNTUA ) THEN- CALL ZUNGBR( 'Q', M, M, N, A, LDA, DUM(1),- $ DUM(1), -1, IERR )- LWORK_ZUNGBR_Q=DUM(1)+ CALL ZUNGBR( 'Q', M, M, N, A, LDA, CDUM(1),+ $ CDUM(1), -1, IERR )+ LWORK_ZUNGBR_Q=CDUM(1)MAXWRK = MAX( MAXWRK, 2*N+LWORK_ZUNGBR_Q )END IFIF( .NOT.WNTVN ) THEN@@ -6197,25 +6197,26 @@*MNTHR = ILAENV( 6, 'ZGESVD', JOBU // JOBVT, M, N, 0, 0 )* Compute space needed for ZGELQF- CALL ZGELQF( M, N, A, LDA, DUM(1), DUM(1), -1, IERR )- LWORK_ZGELQF=DUM(1)+ CALL ZGELQF( M, N, A, LDA, CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZGELQF=CDUM(1)* Compute space needed for ZUNGLQ- CALL ZUNGLQ( N, N, M, DUM(1), N, DUM(1), DUM(1), -1, IERR )- LWORK_ZUNGLQ_N=DUM(1)- CALL ZUNGLQ( M, N, M, A, LDA, DUM(1), DUM(1), -1, IERR )- LWORK_ZUNGLQ_M=DUM(1)+ CALL ZUNGLQ( N, N, M, CDUM(1), N, CDUM(1), CDUM(1), -1,+ $ IERR )+ LWORK_ZUNGLQ_N=CDUM(1)+ CALL ZUNGLQ( M, N, M, A, LDA, CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZUNGLQ_M=CDUM(1)* Compute space needed for ZGEBRD- CALL ZGEBRD( M, M, A, LDA, S, DUM(1), DUM(1),- $ DUM(1), DUM(1), -1, IERR )- LWORK_ZGEBRD=DUM(1)+ CALL ZGEBRD( M, M, A, LDA, S, DUM(1), CDUM(1),+ $ CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZGEBRD=CDUM(1)* Compute space needed for ZUNGBR P- CALL ZUNGBR( 'P', M, M, M, A, N, DUM(1),- $ DUM(1), -1, IERR )- LWORK_ZUNGBR_P=DUM(1)+ CALL ZUNGBR( 'P', M, M, M, A, N, CDUM(1),+ $ CDUM(1), -1, IERR )+ LWORK_ZUNGBR_P=CDUM(1)* Compute space needed for ZUNGBR Q- CALL ZUNGBR( 'Q', M, M, M, A, N, DUM(1),- $ DUM(1), -1, IERR )- LWORK_ZUNGBR_Q=DUM(1)+ CALL ZUNGBR( 'Q', M, M, M, A, N, CDUM(1),+ $ CDUM(1), -1, IERR )+ LWORK_ZUNGBR_Q=CDUM(1)IF( N.GE.MNTHR ) THENIF( WNTVN ) THEN*@@ -6319,21 +6320,21 @@** Path 10t(N greater than M, but not much larger)*- CALL ZGEBRD( M, N, A, LDA, S, DUM(1), DUM(1),- $ DUM(1), DUM(1), -1, IERR )- LWORK_ZGEBRD=DUM(1)+ CALL ZGEBRD( M, N, A, LDA, S, DUM(1), CDUM(1),+ $ CDUM(1), CDUM(1), -1, IERR )+ LWORK_ZGEBRD=CDUM(1)MAXWRK = 2*M + LWORK_ZGEBRDIF( WNTVS .OR. WNTVO ) THEN* Compute space needed for ZUNGBR P- CALL ZUNGBR( 'P', M, N, M, A, N, DUM(1),- $ DUM(1), -1, IERR )- LWORK_ZUNGBR_P=DUM(1)+ CALL ZUNGBR( 'P', M, N, M, A, N, CDUM(1),+ $ CDUM(1), -1, IERR )+ LWORK_ZUNGBR_P=CDUM(1)MAXWRK = MAX( MAXWRK, 2*M+LWORK_ZUNGBR_P )END IFIF( WNTVA ) THEN- CALL ZUNGBR( 'P', N, N, M, A, N, DUM(1),- $ DUM(1), -1, IERR )- LWORK_ZUNGBR_P=DUM(1)+ CALL ZUNGBR( 'P', N, N, M, A, N, CDUM(1),+ $ CDUM(1), -1, IERR )+ LWORK_ZUNGBR_P=CDUM(1)MAXWRK = MAX( MAXWRK, 2*M+LWORK_ZUNGBR_P )END IFIF( .NOT.WNTUN ) THEN