dbeg()

double precision function dbeg

(

logical

reset

)

Definition at line 3058 of file dblat2.f.

*
*  Generates random numbers uniformly distributed between -0.5 and 0.5.
*
*  Auxiliary routine for test program for Level 2 Blas.
*
*  -- Written on 10-August-1987.
*     Richard Hanson, Sandia National Labs.
*     Jeremy Du Croz, NAG Central Office.
*
*     .. Scalar Arguments ..
      LOGICAL            RESET
*     .. Local Scalars ..
      INTEGER            I, IC, MI
*     .. Save statement ..
      SAVE               i, ic, mi
*     .. Intrinsic Functions ..
      INTRINSIC          dble
*     .. Executable Statements ..
      IF( reset )THEN
*        Initialize local variables.
         mi = 891
         i = 7
         ic = 0
         reset = .false.
      END IF
*
*     The sequence of values of I is bounded between 1 and 999.
*     If initial I = 1,2,3,6,7 or 9, the period will be 50.
*     If initial I = 4 or 8, the period will be 25.
*     If initial I = 5, the period will be 10.
*     IC is used to break up the period by skipping 1 value of I in 6.
*
      ic = ic + 1
   10 i = i*mi
      i = i - 1000*( i/1000 )
      IF( ic.GE.5 )THEN
         ic = 0
         GO TO 10
      END IF
      dbeg = dble( i - 500 )/1001.0d0
      RETURN
*
*     End of DBEG
*

Here is the call graph for this function:

Here is the caller graph for this function: