LAPACK  3.8.0
LAPACK: Linear Algebra PACKage

◆ cbeg()

complex function cbeg ( logical  RESET)

Definition at line 3335 of file cblat3.f.

3335 *
3336 * Generates complex numbers as pairs of random numbers uniformly
3337 * distributed between -0.5 and 0.5.
3338 *
3339 * Auxiliary routine for test program for Level 3 Blas.
3340 *
3341 * -- Written on 8-February-1989.
3342 * Jack Dongarra, Argonne National Laboratory.
3343 * Iain Duff, AERE Harwell.
3344 * Jeremy Du Croz, Numerical Algorithms Group Ltd.
3345 * Sven Hammarling, Numerical Algorithms Group Ltd.
3346 *
3347 * .. Scalar Arguments ..
3348  LOGICAL reset
3349 * .. Local Scalars ..
3350  INTEGER i, ic, j, mi, mj
3351 * .. Save statement ..
3352  SAVE i, ic, j, mi, mj
3353 * .. Intrinsic Functions ..
3354  INTRINSIC cmplx
3355 * .. Executable Statements ..
3356  IF( reset )THEN
3357 * Initialize local variables.
3358  mi = 891
3359  mj = 457
3360  i = 7
3361  j = 7
3362  ic = 0
3363  reset = .false.
3364  END IF
3365 *
3366 * The sequence of values of I or J is bounded between 1 and 999.
3367 * If initial I or J = 1,2,3,6,7 or 9, the period will be 50.
3368 * If initial I or J = 4 or 8, the period will be 25.
3369 * If initial I or J = 5, the period will be 10.
3370 * IC is used to break up the period by skipping 1 value of I or J
3371 * in 6.
3372 *
3373  ic = ic + 1
3374  10 i = i*mi
3375  j = j*mj
3376  i = i - 1000*( i/1000 )
3377  j = j - 1000*( j/1000 )
3378  IF( ic.GE.5 )THEN
3379  ic = 0
3380  GO TO 10
3381  END IF
3382  cbeg = cmplx( ( i - 500 )/1001.0, ( j - 500 )/1001.0 )
3383  RETURN
3384 *
3385 * End of CBEG.
3386 *
complex function cbeg(RESET)
Definition: cblat2.f:3131
Here is the call graph for this function: