zbeg()

complex*16 function zbeg

(

logical

reset

)

Definition at line 3341 of file zblat3.f.

*
*  Generates complex numbers as pairs of random numbers uniformly
*  distributed between -0.5 and 0.5.
*
*  Auxiliary routine for test program for Level 3 Blas.
*
*  -- Written on 8-February-1989.
*     Jack Dongarra, Argonne National Laboratory.
*     Iain Duff, AERE Harwell.
*     Jeremy Du Croz, Numerical Algorithms Group Ltd.
*     Sven Hammarling, Numerical Algorithms Group Ltd.
*
*     .. Scalar Arguments ..
      LOGICAL            RESET
*     .. Local Scalars ..
      INTEGER            I, IC, J, MI, MJ
*     .. Save statement ..
      SAVE               i, ic, j, mi, mj
*     .. Intrinsic Functions ..
      INTRINSIC          dcmplx
*     .. Executable Statements ..
      IF( reset )THEN
*        Initialize local variables.
         mi = 891
         mj = 457
         i = 7
         j = 7
         ic = 0
         reset = .false.
      END IF
*
*     The sequence of values of I or J is bounded between 1 and 999.
*     If initial I or J = 1,2,3,6,7 or 9, the period will be 50.
*     If initial I or J = 4 or 8, the period will be 25.
*     If initial I or J = 5, the period will be 10.
*     IC is used to break up the period by skipping 1 value of I or J
*     in 6.
*
      ic = ic + 1
   10 i = i*mi
      j = j*mj
      i = i - 1000*( i/1000 )
      j = j - 1000*( j/1000 )
      IF( ic.GE.5 )THEN
         ic = 0
         GO TO 10
      END IF
      zbeg = dcmplx( ( i - 500 )/1001.0d0, ( j - 500 )/1001.0d0 )
      RETURN
*
*     End of ZBEG
*

Here is the call graph for this function: