Fix -Werror=conversion problems in LAPACK with integer-8

Author: mkrainiuk

SRC/cbbcsd.f

@@ -420,7 +420,7 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P,
 *
       IF( INFO .EQ. 0 .AND. Q .EQ. 0 ) THEN
          LRWORKMIN = 1
-         RWORK(1) = LRWORKMIN
+         RWORK(1) = REAL( LRWORKMIN )
          RETURN
       END IF
 *
@@ -437,7 +437,7 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P,
          IV2TSN = IV2TCS + Q
          LRWORKOPT = IV2TSN + Q - 1
          LRWORKMIN = LRWORKOPT
-         RWORK(1) = LRWORKOPT
+         RWORK(1) = REAL( LRWORKOPT )
          IF( LRWORK .LT. LRWORKMIN .AND. .NOT. LQUERY ) THEN
             INFO = -28
          END IF
@@ -456,7 +456,7 @@ SUBROUTINE CBBCSD( JOBU1, JOBU2, JOBV1T, JOBV2T, TRANS, M, P,
       UNFL = SLAMCH( 'Safe minimum' )
       TOLMUL = MAX( TEN, MIN( HUNDRED, EPS**MEIGHTH ) )
       TOL = TOLMUL*EPS
-      THRESH = MAX( TOL, MAXITR*Q*Q*UNFL )
+      THRESH = MAX( TOL, REAL( MAXITR*Q*Q )*UNFL )
 *
 *     Test for negligible sines or cosines
 *

SRC/cbdsqr.f

@@ -393,12 +393,12 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
    40    CONTINUE
    50    CONTINUE
          SMINOA = SMINOA / SQRT( REAL( N ) )
-         THRESH = MAX( TOL*SMINOA, MAXITR*N*N*UNFL )
+         THRESH = MAX( TOL*SMINOA, REAL( MAXITR*N*N )*UNFL )
       ELSE
 *
 *        Absolute accuracy desired
 *
-         THRESH = MAX( ABS( TOL )*SMAX, MAXITR*N*N*UNFL )
+         THRESH = MAX( ABS( TOL )*SMAX, REAL( MAXITR*N*N )*UNFL )
       END IF
 *
 *     Prepare for main iteration loop for the singular values
@@ -565,7 +565,7 @@ SUBROUTINE CBDSQR( UPLO, N, NCVT, NRU, NCC, D, E, VT, LDVT, U,
 *     Compute shift.  First, test if shifting would ruin relative
 *     accuracy, and if so set the shift to zero.
 *
-      IF( TOL.GE.ZERO .AND. N*TOL*( SMIN / SMAX ).LE.
+      IF( TOL.GE.ZERO .AND. REAL( N )*TOL*( SMIN / SMAX ).LE.
      $    MAX( EPS, HNDRTH*TOL ) ) THEN
 *
 *        Use a zero shift to avoid loss of relative accuracy

SRC/cgbrfs.f

@@ -317,7 +317,7 @@ SUBROUTINE CGBRFS( TRANS, N, KL, KU, NRHS, AB, LDAB, AFB,
       NZ = MIN( KL+KU+2, N+1 )
       EPS = SLAMCH( 'Epsilon' )
       SAFMIN = SLAMCH( 'Safe minimum' )
-      SAFE1 = NZ*SAFMIN
+      SAFE1 = REAL( NZ )*SAFMIN
       SAFE2 = SAFE1 / EPS
 *
 *     Do for each right hand side
@@ -426,10 +426,11 @@ SUBROUTINE CGBRFS( TRANS, N, KL, KU, NRHS, AB, LDAB, AFB,
 *
          DO 90 I = 1, N
             IF( RWORK( I ).GT.SAFE2 ) THEN
-               RWORK( I ) = CABS1( WORK( I ) ) + NZ*EPS*RWORK( I )
+               RWORK( I ) = CABS1( WORK( I ) ) + REAL( NZ )*
+     $                      EPS*RWORK( I )
             ELSE
-               RWORK( I ) = CABS1( WORK( I ) ) + NZ*EPS*RWORK( I ) +
-     $                      SAFE1
+               RWORK( I ) = CABS1( WORK( I ) ) + REAL( NZ )*
+     $                      EPS*RWORK( I ) + SAFE1
             END IF
    90    CONTINUE
 *

SRC/cgebal.f

@@ -267,7 +267,7 @@ SUBROUTINE CGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )
                END DO
 *
                IF( CANSWAP ) THEN
-                  SCALE( L ) = I
+                  SCALE( L ) = REAL( I )
                   IF( I.NE.L ) THEN
                      CALL CSWAP( L, A( 1, I ), 1, A( 1, L ), 1 )
                      CALL CSWAP( N-K+1, A( I, K ), LDA, A( L, K ),
@@ -304,7 +304,7 @@ SUBROUTINE CGEBAL( JOB, N, A, LDA, ILO, IHI, SCALE, INFO )
                END DO
 *
                IF( CANSWAP ) THEN
-                  SCALE( K ) = J
+                  SCALE( K ) = REAL( J )
                   IF( J.NE.K ) THEN
                      CALL CSWAP( L, A( 1, J ), 1, A( 1, K ), 1 )
                      CALL CSWAP( N-K+1, A( J, K ), LDA, A( K, K ),

SRC/cgebrd.f

@@ -346,7 +346,7 @@ SUBROUTINE CGEBRD( M, N, A, LDA, D, E, TAUQ, TAUP, WORK, LWORK,
 *
       CALL CGEBD2( M-I+1, N-I+1, A( I, I ), LDA, D( I ), E( I ),
      $             TAUQ( I ), TAUP( I ), WORK, IINFO )
-      WORK( 1 ) = WS
+      WORK( 1 ) = CMPLX( WS )
       RETURN
 *
 *     End of CGEBRD

SRC/cgedmd.f90

@@ -588,9 +588,9 @@ SUBROUTINE CGEDMD( JOBS, JOBZ, JOBR, JOBF,  WHTSVD,   &
       ELSE IF ( LQUERY ) THEN
 !     Return minimal and optimal workspace sizes
           IWORK(1) = IMINWR
-          RWORK(1) = MLRWRK
-          ZWORK(1) = MLWORK
-          ZWORK(2) = OLWORK
+          RWORK(1) = REAL( MLRWRK )
+          ZWORK(1) = CMPLX( MLWORK )
+          ZWORK(2) = CMPLX( OLWORK )
           RETURN
       END IF
 !............................................................

SRC/cgedmdq.f90

@@ -599,10 +599,10 @@ SUBROUTINE CGEDMDQ( JOBS,  JOBZ, JOBR, JOBQ, JOBT, JOBF,   &
       ELSE IF ( LQUERY ) THEN
 !     Return minimal and optimal workspace sizes
           IWORK(1) = IMINWR
-          ZWORK(1) = MLWORK
-          ZWORK(2) = OLWORK
-          WORK(1)  = MLRWRK
-          WORK(2)  = MLRWRK
+          ZWORK(1) = CMPLX( MLWORK )
+          ZWORK(2) = CMPLX( OLWORK )
+          WORK(1)  = REAL( MLRWRK )
+          WORK(2)  = REAL( MLRWRK )
           RETURN
       END IF   
 !.....	  
@@ -687,4 +687,4 @@ SUBROUTINE CGEDMDQ( JOBS,  JOBZ, JOBR, JOBQ, JOBT, JOBF,   &
       RETURN
 !      
       END SUBROUTINE CGEDMDQ
-    
+    

SRC/cgees.f

@@ -296,7 +296,7 @@ SUBROUTINE CGEES( JOBVS, SORT, SELECT, N, A, LDA, SDIM, W, VS,
                MAXWRK = MAX( MAXWRK, HSWORK )
             END IF
          END IF
-         WORK( 1 ) = MAXWRK
+         WORK( 1 ) = CMPLX( MAXWRK )
 *
          IF( LWORK.LT.MINWRK .AND. .NOT.LQUERY ) THEN
             INFO = -12
@@ -419,7 +419,7 @@ SUBROUTINE CGEES( JOBVS, SORT, SELECT, N, A, LDA, SDIM, W, VS,
          CALL CCOPY( N, A, LDA+1, W, 1 )
       END IF
 *
-      WORK( 1 ) = MAXWRK
+      WORK( 1 ) = CMPLX( MAXWRK )
       RETURN
 *
 *     End of CGEES

SRC/cgeesx.f

@@ -355,7 +355,7 @@ SUBROUTINE CGEESX( JOBVS, SORT, SELECT, SENSE, N, A, LDA, SDIM,
             IF( .NOT.WANTSN )
      $         LWRK = MAX( LWRK, ( N*N )/2 )
          END IF
-         WORK( 1 ) = LWRK
+         WORK( 1 ) = CMPLX( LWRK )
 *
          IF( LWORK.LT.MINWRK .AND. .NOT.LQUERY ) THEN
             INFO = -15
@@ -496,7 +496,7 @@ SUBROUTINE CGEESX( JOBVS, SORT, SELECT, SENSE, N, A, LDA, SDIM,
          END IF
       END IF
 *
-      WORK( 1 ) = MAXWRK
+      WORK( 1 ) = CMPLX( MAXWRK )
       RETURN
 *
 *     End of CGEESX

SRC/cgeev.f

@@ -298,7 +298,7 @@ SUBROUTINE CGEEV( JOBVL, JOBVR, N, A, LDA, W, VL, LDVL, VR,
             HSWORK = INT( WORK(1) )
             MAXWRK = MAX( MAXWRK, HSWORK, MINWRK )
          END IF
-         WORK( 1 ) = MAXWRK
+         WORK( 1 ) = CMPLX( MAXWRK )
 *
          IF( LWORK.LT.MINWRK .AND. .NOT.LQUERY ) THEN
             INFO = -12
@@ -506,7 +506,7 @@ SUBROUTINE CGEEV( JOBVL, JOBVR, N, A, LDA, W, VL, LDVL, VR,
          END IF
       END IF
 *
-      WORK( 1 ) = MAXWRK
+      WORK( 1 ) = CMPLX( MAXWRK )
       RETURN
 *
 *     End of CGEEV