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

## ◆ zbeg()

 complex*16 function zbeg ( logical reset )

Definition at line 2679 of file c_zblat2.f.

2680*
2681* Generates complex numbers as pairs of random numbers uniformly
2682* distributed between -0.5 and 0.5.
2683*
2684* Auxiliary routine for test program for Level 2 Blas.
2685*
2686* -- Written on 10-August-1987.
2687* Richard Hanson, Sandia National Labs.
2688* Jeremy Du Croz, NAG Central Office.
2689*
2690* .. Scalar Arguments ..
2691 LOGICAL RESET
2692* .. Local Scalars ..
2693 INTEGER I, IC, J, MI, MJ
2694* .. Save statement ..
2695 SAVE i, ic, j, mi, mj
2696* .. Intrinsic Functions ..
2697 INTRINSIC dcmplx
2698* .. Executable Statements ..
2699 IF( reset )THEN
2700* Initialize local variables.
2701 mi = 891
2702 mj = 457
2703 i = 7
2704 j = 7
2705 ic = 0
2706 reset = .false.
2707 END IF
2708*
2709* The sequence of values of I or J is bounded between 1 and 999.
2710* If initial I or J = 1,2,3,6,7 or 9, the period will be 50.
2711* If initial I or J = 4 or 8, the period will be 25.
2712* If initial I or J = 5, the period will be 10.
2713* IC is used to break up the period by skipping 1 value of I or J
2714* in 6.
2715*
2716 ic = ic + 1
2717 10 i = i*mi
2718 j = j*mj
2719 i = i - 1000*( i/1000 )
2720 j = j - 1000*( j/1000 )
2721 IF( ic.GE.5 )THEN
2722 ic = 0
2723 GO TO 10
2724 END IF
2725 zbeg = dcmplx( ( i - 500 )/1001.0, ( j - 500 )/1001.0 )
2726 RETURN
2727*
2728* End of ZBEG.
2729*
complex *16 function zbeg(reset)
Definition zblat2.f:3164
Here is the call graph for this function: