LAPACK  3.10.0
LAPACK: Linear Algebra PACKage

◆ cbeg()

complex function cbeg ( logical  RESET)

Definition at line 3331 of file cblat3.f.

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