pb_dran()

double precision function pb_dran

(

integer

idumm

)

Definition at line 2631 of file pzblastim.f.

*
*  -- PBLAS test routine (version 2.0) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     April 1, 1998
*
*     .. Scalar Arguments ..
      INTEGER            IDUMM
*     ..
*
*  Purpose
*  =======
*
*  PB_DRAN generates the next number in the random sequence.
*
*  Arguments
*  =========
*
*  IDUMM   (local input) INTEGER
*          This argument is ignored, but necessary to a FORTRAN 77 func-
*          tion.
*
*  Further Details
*  ===============
*
*  On entry, the array IRAND stored in the common block  RANCOM contains
*  the information (2 integers)  required to generate the next number in
*  the sequence X( n ). This number is computed as
*
*     X( n ) = ( 2^16 * IRAND( 2 ) + IRAND( 1 ) ) / d,
*
*  where the constant d is the  largest  32 bit  positive  integer.  The
*  array  IRAND  is  then  updated for the generation of the next number
*  X( n+1 ) in the random sequence as follows X( n+1 ) = a * X( n ) + c.
*  The constants  a  and c  should have been preliminarily stored in the
*  array  IACS  as  2 pairs of integers. The initial set up of IRAND and
*  IACS is performed by the routine PB_SETRAN.
*
*  -- Written on April 1, 1998 by
*     Antoine Petitet, University  of  Tennessee, Knoxville 37996, USA.
*
*  =====================================================================
*
*     .. Parameters ..
      DOUBLE PRECISION   DIVFAC, POW16
      parameter( divfac = 2.147483648d+9,
     $                   pow16 = 6.5536d+4 )
*     ..
*     .. Local Arrays ..
      INTEGER            J( 2 )
*     ..
*     .. External Subroutines ..
      EXTERNAL           pb_ladd, pb_lmul
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          dble
*     ..
*     .. Common Blocks ..
      INTEGER            IACS( 4 ), IRAND( 2 )
      COMMON             /rancom/ irand, iacs
*     ..
*     .. Save Statements ..
      SAVE               /rancom/
*     ..
*     .. Executable Statements ..
*
      pb_dran = ( dble( irand( 1 ) ) + pow16 * dble( irand( 2 ) ) ) /
     $            divfac
*
      CALL pb_lmul( irand, iacs, j )
      CALL pb_ladd( j, iacs( 3 ), irand )
*
      RETURN
*
*     End of PB_DRAN
*

Here is the call graph for this function: