LAPACK 3.12.0 LAPACK: Linear Algebra PACKage
Searching...
No Matches

## ◆ 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:3156
Here is the call graph for this function: