sroundup_lwork.f File Reference

Go to the source code of this file.

Functions/Subroutines
real function	sroundup_lwork (lwork)
	SROUNDUP_LWORK

Function/Subroutine Documentation

sroundup_lwork()

real function sroundup_lwork

(

integer

lwork

)

SROUNDUP_LWORK

Purpose:

!>
!> SROUNDUP_LWORK deals with a subtle bug with returning LWORK as a Float.
!> This routine guarantees it is rounded up instead of down by
!> multiplying LWORK by 1+eps when it is necessary, where eps is the relative machine precision.
!> E.g.,
!>
!>        float( 16777217            ) == 16777216
!>        float( 16777217 ) * (1.+eps) == 16777218
!>
!> \return SROUNDUP_LWORK
!> \verbatim
!>         SROUNDUP_LWORK >= LWORK.
!>         SROUNDUP_LWORK is guaranteed to have zero decimal part.
!> 

Parameters

[in]

LWORK

Workspace size.

Author

Weslley Pereira, University of Colorado Denver, USA

Further Details:

!>  This routine was inspired in the method `magma_zmake_lwork` from MAGMA.
!>  \see https://bitbucket.org/icl/magma/src/master/control/magma_zauxiliary.cpp
!> 

Definition at line 58 of file sroundup_lwork.f.

*
*  -- LAPACK auxiliary routine --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*
*     .. Scalar Arguments ..
      INTEGER          LWORK
*     ..
*
* =====================================================================
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC         epsilon, real, int
*     ..
*     .. Executable Statements ..
*     ..
      sroundup_lwork = real( lwork )
*
      IF( int( sroundup_lwork ) .LT. lwork ) THEN
*         Force round up of LWORK
          sroundup_lwork = sroundup_lwork * ( 1.0e+0 + epsilon(0.0e+0) )
      ENDIF
*
      RETURN
*
*     End of SROUNDUP_LWORK
*