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

Define type conversions explicitly (Reference-LAPACK PR 703)

This commit is contained in:
Martin Kroeker 2022-11-19 15:22:46 +01:00 committed by GitHub
parent e9b0f5a364
commit 4e60737c2d
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
18 changed files with 44 additions and 44 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
lapack-netlib/SRC/zgebak.f
View File

@ -238,7 +238,7 @@
$ GO TO 40
IF( I.LT.ILO )
$ I = ILO - II
K = SCALE( I )
K = INT( SCALE( I ) )
IF( K.EQ.I )
$ GO TO 40
CALL ZSWAP( M, V( I, 1 ), LDV, V( K, 1 ), LDV )
@ -252,7 +252,7 @@
$ GO TO 50
IF( I.LT.ILO )
$ I = ILO - II
K = SCALE( I )
K = INT( SCALE( I ) )
IF( K.EQ.I )
$ GO TO 50
CALL ZSWAP( M, V( I, 1 ), LDV, V( K, 1 ), LDV )

2
lapack-netlib/SRC/zgees.f
View File

@ -282,7 +282,7 @@
*
CALL ZHSEQR( 'S', JOBVS, N, 1, N, A, LDA, W, VS, LDVS,
$ WORK, -1, IEVAL )
HSWORK = DBLE( WORK( 1 ) )
HSWORK = INT( WORK( 1 ) )
*
IF( .NOT.WANTVS ) THEN
MAXWRK = MAX( MAXWRK, HSWORK )

2
lapack-netlib/SRC/zgeesx.f
View File

@ -337,7 +337,7 @@
*
CALL ZHSEQR( 'S', JOBVS, N, 1, N, A, LDA, W, VS, LDVS,
$ WORK, -1, IEVAL )
HSWORK = DBLE( WORK( 1 ) )
HSWORK = INT( WORK( 1 ) )
*
IF( .NOT.WANTVS ) THEN
MAXWRK = MAX( MAXWRK, HSWORK )

36
lapack-netlib/SRC/zgejsv.f
View File

@ -707,11 +707,11 @@
IF ( LQUERY ) THEN
CALL ZGEQP3( M, N, A, LDA, IWORK, CDUMMY, CDUMMY, -1,
$ RDUMMY, IERR )
LWRK_ZGEQP3 = DBLE( CDUMMY(1) )
LWRK_ZGEQP3 = INT( CDUMMY(1) )
CALL ZGEQRF( N, N, A, LDA, CDUMMY, CDUMMY,-1, IERR )
LWRK_ZGEQRF = DBLE( CDUMMY(1) )
LWRK_ZGEQRF = INT( CDUMMY(1) )
CALL ZGELQF( N, N, A, LDA, CDUMMY, CDUMMY,-1, IERR )
LWRK_ZGELQF = DBLE( CDUMMY(1) )
LWRK_ZGELQF = INT( CDUMMY(1) )
END IF
MINWRK = 2
OPTWRK = 2
@ -727,7 +727,7 @@
IF ( LQUERY ) THEN
CALL ZGESVJ( 'L', 'N', 'N', N, N, A, LDA, SVA, N, V,
$ LDV, CDUMMY, -1, RDUMMY, -1, IERR )
LWRK_ZGESVJ = DBLE( CDUMMY(1) )
LWRK_ZGESVJ = INT( CDUMMY(1) )
IF ( ERREST ) THEN
OPTWRK = MAX( N+LWRK_ZGEQP3, N**2+LWCON,
$ N+LWRK_ZGEQRF, LWRK_ZGESVJ )
@ -763,10 +763,10 @@
IF ( LQUERY ) THEN
CALL ZGESVJ( 'L', 'U', 'N', N,N, U, LDU, SVA, N, A,
$ LDA, CDUMMY, -1, RDUMMY, -1, IERR )
LWRK_ZGESVJ = DBLE( CDUMMY(1) )
LWRK_ZGESVJ = INT( CDUMMY(1) )
CALL ZUNMLQ( 'L', 'C', N, N, N, A, LDA, CDUMMY,
$ V, LDV, CDUMMY, -1, IERR )
LWRK_ZUNMLQ = DBLE( CDUMMY(1) )
LWRK_ZUNMLQ = INT( CDUMMY(1) )
IF ( ERREST ) THEN
OPTWRK = MAX( N+LWRK_ZGEQP3, LWCON, LWRK_ZGESVJ,
$ N+LWRK_ZGELQF, 2*N+LWRK_ZGEQRF,
@ -802,10 +802,10 @@
IF ( LQUERY ) THEN
CALL ZGESVJ( 'L', 'U', 'N', N,N, U, LDU, SVA, N, A,
$ LDA, CDUMMY, -1, RDUMMY, -1, IERR )
LWRK_ZGESVJ = DBLE( CDUMMY(1) )
LWRK_ZGESVJ = INT( CDUMMY(1) )
CALL ZUNMQR( 'L', 'N', M, N, N, A, LDA, CDUMMY, U,
$ LDU, CDUMMY, -1, IERR )
LWRK_ZUNMQRM = DBLE( CDUMMY(1) )
LWRK_ZUNMQRM = INT( CDUMMY(1) )
IF ( ERREST ) THEN
OPTWRK = N + MAX( LWRK_ZGEQP3, LWCON, N+LWRK_ZGEQRF,
$ LWRK_ZGESVJ, LWRK_ZUNMQRM )
@ -864,26 +864,26 @@
IF ( LQUERY ) THEN
CALL ZUNMQR( 'L', 'N', M, N, N, A, LDA, CDUMMY, U,
$ LDU, CDUMMY, -1, IERR )
LWRK_ZUNMQRM = DBLE( CDUMMY(1) )
LWRK_ZUNMQRM = INT( CDUMMY(1) )
CALL ZUNMQR( 'L', 'N', N, N, N, A, LDA, CDUMMY, U,
$ LDU, CDUMMY, -1, IERR )
LWRK_ZUNMQR = DBLE( CDUMMY(1) )
LWRK_ZUNMQR = INT( CDUMMY(1) )
IF ( .NOT. JRACC ) THEN
CALL ZGEQP3( N,N, A, LDA, IWORK, CDUMMY,CDUMMY, -1,
$ RDUMMY, IERR )
LWRK_ZGEQP3N = DBLE( CDUMMY(1) )
LWRK_ZGEQP3N = INT( CDUMMY(1) )
CALL ZGESVJ( 'L', 'U', 'N', N, N, U, LDU, SVA,
$ N, V, LDV, CDUMMY, -1, RDUMMY, -1, IERR )
LWRK_ZGESVJ = DBLE( CDUMMY(1) )
LWRK_ZGESVJ = INT( CDUMMY(1) )
CALL ZGESVJ( 'U', 'U', 'N', N, N, U, LDU, SVA,
$ N, V, LDV, CDUMMY, -1, RDUMMY, -1, IERR )
LWRK_ZGESVJU = DBLE( CDUMMY(1) )
LWRK_ZGESVJU = INT( CDUMMY(1) )
CALL ZGESVJ( 'L', 'U', 'V', N, N, U, LDU, SVA,
$ N, V, LDV, CDUMMY, -1, RDUMMY, -1, IERR )
LWRK_ZGESVJV = DBLE( CDUMMY(1) )
LWRK_ZGESVJV = INT( CDUMMY(1) )
CALL ZUNMLQ( 'L', 'C', N, N, N, A, LDA, CDUMMY,
$ V, LDV, CDUMMY, -1, IERR )
LWRK_ZUNMLQ = DBLE( CDUMMY(1) )
LWRK_ZUNMLQ = INT( CDUMMY(1) )
IF ( ERREST ) THEN
OPTWRK = MAX( N+LWRK_ZGEQP3, N+LWCON,
$ 2*N+N**2+LWCON, 2*N+LWRK_ZGEQRF,
@ -912,13 +912,13 @@
ELSE
CALL ZGESVJ( 'L', 'U', 'V', N, N, U, LDU, SVA,
$ N, V, LDV, CDUMMY, -1, RDUMMY, -1, IERR )
LWRK_ZGESVJV = DBLE( CDUMMY(1) )
LWRK_ZGESVJV = INT( CDUMMY(1) )
CALL ZUNMQR( 'L', 'N', N, N, N, CDUMMY, N, CDUMMY,
$ V, LDV, CDUMMY, -1, IERR )
LWRK_ZUNMQR = DBLE( CDUMMY(1) )
LWRK_ZUNMQR = INT( CDUMMY(1) )
CALL ZUNMQR( 'L', 'N', M, N, N, A, LDA, CDUMMY, U,
$ LDU, CDUMMY, -1, IERR )
LWRK_ZUNMQRM = DBLE( CDUMMY(1) )
LWRK_ZUNMQRM = INT( CDUMMY(1) )
IF ( ERREST ) THEN
OPTWRK = MAX( N+LWRK_ZGEQP3, N+LWCON,
$ 2*N+LWRK_ZGEQRF, 2*N+N**2,

2
lapack-netlib/SRC/zggglm.f
View File

@ -289,7 +289,7 @@
*
CALL ZGGQRF( N, M, P, A, LDA, WORK, B, LDB, WORK( M+1 ),
$ WORK( M+NP+1 ), LWORK-M-NP, INFO )
LOPT = DBLE( WORK( M+NP+1 ) )
LOPT = INT( WORK( M+NP+1 ) )
*
* Update left-hand-side vector d = Q**H*d = ( d1 ) M
* ( d2 ) N-M

2
lapack-netlib/SRC/zgglse.f
View File

@ -276,7 +276,7 @@
*
CALL ZGGRQF( P, M, N, B, LDB, WORK, A, LDA, WORK( P+1 ),
$ WORK( P+MN+1 ), LWORK-P-MN, INFO )
LOPT = DBLE( WORK( P+MN+1 ) )
LOPT = INT( WORK( P+MN+1 ) )
*
* Update c = Z**H *c = ( c1 ) N-P
* ( c2 ) M+P-N

2
lapack-netlib/SRC/zggqrf.f
View File

@ -276,7 +276,7 @@
* QR factorization of N-by-M matrix A: A = Q*R
*
CALL ZGEQRF( N, M, A, LDA, TAUA, WORK, LWORK, INFO )
LOPT = DBLE( WORK( 1 ) )
LOPT = INT( WORK( 1 ) )
*
* Update B := Q**H*B.
*

2
lapack-netlib/SRC/zggrqf.f
View File

@ -275,7 +275,7 @@
* RQ factorization of M-by-N matrix A: A = R*Q
*
CALL ZGERQF( M, N, A, LDA, TAUA, WORK, LWORK, INFO )
LOPT = DBLE( WORK( 1 ) )
LOPT = INT( WORK( 1 ) )
*
* Update B := B*Q**H
*

6
lapack-netlib/SRC/zhegvd.f
View File

@ -360,9 +360,9 @@
CALL ZHEGST( ITYPE, UPLO, N, A, LDA, B, LDB, INFO )
CALL ZHEEVD( JOBZ, UPLO, N, A, LDA, W, WORK, LWORK, RWORK, LRWORK,
$ IWORK, LIWORK, INFO )
LOPT = MAX( DBLE( LOPT ), DBLE( WORK( 1 ) ) )
LROPT = MAX( DBLE( LROPT ), DBLE( RWORK( 1 ) ) )
LIOPT = MAX( DBLE( LIOPT ), DBLE( IWORK( 1 ) ) )
LOPT = INT( MAX( DBLE( LOPT ), DBLE( WORK( 1 ) ) ) )
LROPT = INT( MAX( DBLE( LROPT ), DBLE( RWORK( 1 ) ) ) )
LIOPT = INT( MAX( DBLE( LIOPT ), DBLE( IWORK( 1 ) ) ) )
*
IF( WANTZ .AND. INFO.EQ.0 ) THEN
*

2
lapack-netlib/SRC/zhesv_rk.f
View File

@ -280,7 +280,7 @@
LWKOPT = 1
ELSE
CALL ZHETRF_RK( UPLO, N, A, LDA, E, IPIV, WORK, -1, INFO )
LWKOPT = DBLE( WORK(1) )
LWKOPT = INT( DBLE( WORK( 1 ) ) )
END IF
WORK( 1 ) = LWKOPT
END IF

6
lapack-netlib/SRC/zhpgvd.f
View File

@ -335,9 +335,9 @@
CALL ZHPGST( ITYPE, UPLO, N, AP, BP, INFO )
CALL ZHPEVD( JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, LWORK, RWORK,
$ LRWORK, IWORK, LIWORK, INFO )
LWMIN = MAX( DBLE( LWMIN ), DBLE( WORK( 1 ) ) )
LRWMIN = MAX( DBLE( LRWMIN ), DBLE( RWORK( 1 ) ) )
LIWMIN = MAX( DBLE( LIWMIN ), DBLE( IWORK( 1 ) ) )
LWMIN = INT( MAX( DBLE( LWMIN ), DBLE( WORK( 1 ) ) ) )
LRWMIN = INT( MAX( DBLE( LRWMIN ), DBLE( RWORK( 1 ) ) ) )
LIWMIN = INT( MAX( DBLE( LIWMIN ), DBLE( IWORK( 1 ) ) ) )
*
IF( WANTZ ) THEN
*

4
lapack-netlib/SRC/zlag2c.f
View File

@ -124,7 +124,7 @@
DOUBLE PRECISION RMAX
* ..
* .. Intrinsic Functions ..
INTRINSIC DBLE, DIMAG
INTRINSIC DBLE, DIMAG, CMPLX
* ..
* .. External Functions ..
REAL SLAMCH
@ -142,7 +142,7 @@
INFO = 1
GO TO 30
END IF
SA( I, J ) = A( I, J )
SA( I, J ) = CMPLX( A( I, J ) )
10 CONTINUE
20 CONTINUE
INFO = 0

4
lapack-netlib/SRC/zlaic1.f
View File

@ -348,9 +348,9 @@
B = ( ZETA2*ZETA2+ZETA1*ZETA1-ONE )*HALF
C = ZETA1*ZETA1
IF( B.GE.ZERO ) THEN
T = -C / ( B+SQRT( B*B+C ) )
T = DBLE( -C / ( B+SQRT( B*B+C ) ) )
ELSE
T = B - SQRT( B*B+C )
T = DBLE( B - SQRT( B*B+C ) )
END IF
SINE = -( ALPHA / ABSEST ) / T
COSINE = -( GAMMA / ABSEST ) / ( ONE+T )

6
lapack-netlib/SRC/zlat2c.f
View File

@ -130,7 +130,7 @@
LOGICAL UPPER
* ..
* .. Intrinsic Functions ..
INTRINSIC DBLE, DIMAG
INTRINSIC DBLE, DIMAG, CMPLX
* ..
* .. External Functions ..
REAL SLAMCH
@ -151,7 +151,7 @@
INFO = 1
GO TO 50
END IF
SA( I, J ) = A( I, J )
SA( I, J ) = CMPLX( A( I, J ) )
10 CONTINUE
20 CONTINUE
ELSE
@ -164,7 +164,7 @@
INFO = 1
GO TO 50
END IF
SA( I, J ) = A( I, J )
SA( I, J ) = CMPLX( A( I, J ) )
30 CONTINUE
40 CONTINUE
END IF

2
lapack-netlib/SRC/zsysv.f
View File

@ -223,7 +223,7 @@
LWKOPT = 1
ELSE
CALL ZSYTRF( UPLO, N, A, LDA, IPIV, WORK, -1, INFO )
LWKOPT = DBLE( WORK(1) )
LWKOPT = INT( DBLE( WORK( 1 ) ) )
END IF
WORK( 1 ) = LWKOPT
END IF

2
lapack-netlib/SRC/zsysv_rk.f
View File

@ -280,7 +280,7 @@
LWKOPT = 1
ELSE
CALL ZSYTRF_RK( UPLO, N, A, LDA, E, IPIV, WORK, -1, INFO )
LWKOPT = DBLE( WORK(1) )
LWKOPT = INT( DBLE( WORK( 1 ) ) )
END IF
WORK( 1 ) = LWKOPT
END IF

2
lapack-netlib/SRC/zsysv_rook.f
View File

@ -256,7 +256,7 @@
LWKOPT = 1
ELSE
CALL ZSYTRF_ROOK( UPLO, N, A, LDA, IPIV, WORK, -1, INFO )
LWKOPT = DBLE( WORK(1) )
LWKOPT = INT( DBLE( WORK( 1 ) ) )
END IF
WORK( 1 ) = LWKOPT
END IF

2
lapack-netlib/SRC/zungbr.f
View File

@ -233,7 +233,7 @@
END IF
END IF
END IF
LWKOPT = DBLE( WORK( 1 ) )
LWKOPT = INT( DBLE( WORK( 1 ) ) )
LWKOPT = MAX (LWKOPT, MN)
END IF
*