Fix bug: ROUNDUP_LWORK must use (1+EPSILON(0)) instead of (1+DLAMCH(EPS))

Test results:

 SROUNDUP_LWORK(0)        =    0.00000000
 SROUNDUP_LWORK(1)        =    1.00000000
 SROUNDUP_LWORK(317)      =    317.000000
 DROUNDUP_LWORK(0)        =    0.0000000000000000
 DROUNDUP_LWORK(1)        =    1.0000000000000000
 DROUNDUP_LWORK(317)      =    317.00000000000000

 Tests with X = 16777216 = 2**24
X-2 = 16777214; SROUNDUP_LWORK(X-2) = 16777214.0
X-1 = 16777215; SROUNDUP_LWORK(X-1) = 16777215.0
X   = 16777216; SROUNDUP_LWORK(X  ) = 16777218.0
X+1 = 16777217; SROUNDUP_LWORK(X+1) = 16777218.0
X+2 = 16777218; SROUNDUP_LWORK(X+2) = 16777220.0

 Tests with X = 9007199254740992 = 2**53
X-2 = 9007199254740990; SROUNDUP_LWORK(X-2) = 9007200328482816.0
X-1 = 9007199254740991; SROUNDUP_LWORK(X-1) = 9007200328482816.0
X   = 9007199254740992; SROUNDUP_LWORK(X  ) = 9007200328482816.0
X+1 = 9007199254740993; SROUNDUP_LWORK(X+1) = 9007200328482816.0
X+2 = 9007199254740994; SROUNDUP_LWORK(X+2) = 9007200328482816.0

 Tests with X = 9007199254740992 = 2**53
X-2 = 9007199254740990; DROUNDUP_LWORK(X-2) = 9007199254740990.0
X-1 = 9007199254740991; DROUNDUP_LWORK(X-1) = 9007199254740991.0
X   = 9007199254740992; DROUNDUP_LWORK(X  ) = 9007199254740994.0
X+1 = 9007199254740993; DROUNDUP_LWORK(X+1) = 9007199254740994.0
X+2 = 9007199254740994; DROUNDUP_LWORK(X+2) = 9007199254740996.0

Author: weslleyspereira

INSTALL/droundup_lwork.f

@@ -63,20 +63,16 @@ DOUBLE PRECISION FUNCTION DROUNDUP_LWORK( LWORK )
 *
 * =====================================================================
 *     ..
-*     .. External Functions ..
-      DOUBLE PRECISION  DLAMCH
-      EXTERNAL          DLAMCH
-*     ..
 *     .. Intrinsic Functions ..
-      INTRINSIC         DIGITS, RADIX
+      INTRINSIC         DIGITS, RADIX, EPSILON
 *     ..
 *     .. Executable Statements ..
 *     ..
       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') )
+          DROUNDUP_LWORK = DROUNDUP_LWORK * ( 1.0D+0 + EPSILON(0.0D+0) )
       ENDIF
 *
       RETURN

INSTALL/sroundup_lwork.f

@@ -66,20 +66,16 @@ REAL             FUNCTION SROUNDUP_LWORK( LWORK )
 *
 * =====================================================================
 *     ..
-*     .. External Functions ..
-      REAL              SLAMCH
-      EXTERNAL          SLAMCH
-*     ..
 *     .. Intrinsic Functions ..
-      INTRINSIC         DIGITS, RADIX
+      INTRINSIC         DIGITS, RADIX, EPSILON
 *     ..
 *     .. Executable Statements ..
 *     ..
       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') )
+          SROUNDUP_LWORK = SROUNDUP_LWORK * ( 1.0E+0 + EPSILON(0.0E+0) )
       ENDIF
 *
       RETURN