RoundUp functions now always returns an integer

Author: weslleyspereira

INSTALL/droundup_lwork.f

@@ -64,12 +64,21 @@ DOUBLE PRECISION FUNCTION DROUNDUP_LWORK( LWORK )
 * =====================================================================
 *     ..
 *     .. External Functions ..
-      DOUBLE PRECISION   DLAMCH
-      EXTERNAL           DLAMCH
+      DOUBLE PRECISION  DLAMCH
+      EXTERNAL          DLAMCH
+*     ..
+*     .. Intrinsic Functions ..
+      INTRINSIC         DIGITS, RADIX
 *     ..
 *     .. Executable Statements ..
 *     ..
-      DROUNDUP_LWORK = LWORK * ( 1 + DLAMCH('EPS') )
+      DROUNDUP_LWORK = LWORK
+*
+      IF( DROUNDUP_LWORK .GE. DBLE(RADIX(0.0D+0))**DIGITS(0.0D+0) ) THEN
+*         If LWORK can't be represented exactly in double precision
+          DROUNDUP_LWORK = LWORK * ( 1.0D+0 + DLAMCH('EPS') )
+      ENDIF
+*
       RETURN
 *
 *     End of DROUNDUP_LWORK

INSTALL/sroundup_lwork.f

@@ -67,12 +67,21 @@ REAL             FUNCTION SROUNDUP_LWORK( LWORK )
 * =====================================================================
 *     ..
 *     .. External Functions ..
-      REAL               SLAMCH
-      EXTERNAL           SLAMCH
+      REAL              SLAMCH
+      EXTERNAL          SLAMCH
+*     ..
+*     .. Intrinsic Functions ..
+      INTRINSIC         DIGITS, RADIX
 *     ..
 *     .. Executable Statements ..
 *     ..
-      SROUNDUP_LWORK = LWORK * ( 1 + SLAMCH('EPS') )
+      SROUNDUP_LWORK = LWORK
+*
+      IF( SROUNDUP_LWORK .GE. REAL(RADIX(0.0E+0))**DIGITS(0.0E+0) ) THEN
+*         If LWORK can't be represented exactly in single precision
+          SROUNDUP_LWORK = LWORK * ( 1.0E+0 + SLAMCH('EPS') )
+      ENDIF
+*
       RETURN
 *
 *     End of SROUNDUP_LWORK