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

## ◆ zbeg()

 complex*16 function zbeg ( logical RESET )

Definition at line 3163 of file zblat2.f.

3164*
3165* Generates complex numbers as pairs of random numbers uniformly
3166* distributed between -0.5 and 0.5.
3167*
3168* Auxiliary routine for test program for Level 2 Blas.
3169*
3170* -- Written on 10-August-1987.
3171* Richard Hanson, Sandia National Labs.
3172* Jeremy Du Croz, NAG Central Office.
3173*
3174* .. Scalar Arguments ..
3175 LOGICAL RESET
3176* .. Local Scalars ..
3177 INTEGER I, IC, J, MI, MJ
3178* .. Save statement ..
3179 SAVE i, ic, j, mi, mj
3180* .. Intrinsic Functions ..
3181 INTRINSIC dcmplx
3182* .. Executable Statements ..
3183 IF( reset )THEN
3184* Initialize local variables.
3185 mi = 891
3186 mj = 457
3187 i = 7
3188 j = 7
3189 ic = 0
3190 reset = .false.
3191 END IF
3192*
3193* The sequence of values of I or J is bounded between 1 and 999.
3194* If initial I or J = 1,2,3,6,7 or 9, the period will be 50.
3195* If initial I or J = 4 or 8, the period will be 25.
3196* If initial I or J = 5, the period will be 10.
3197* IC is used to break up the period by skipping 1 value of I or J
3198* in 6.
3199*
3200 ic = ic + 1
3201 10 i = i*mi
3202 j = j*mj
3203 i = i - 1000*( i/1000 )
3204 j = j - 1000*( j/1000 )
3205 IF( ic.GE.5 )THEN
3206 ic = 0
3207 GO TO 10
3208 END IF
3209 zbeg = dcmplx( ( i - 500 )/1001.0d0, ( j - 500 )/1001.0d0 )
3210 RETURN
3211*
3212* End of ZBEG
3213*
complex *16 function zbeg(RESET)
Definition: zblat2.f:3164
Here is the call graph for this function:
Here is the caller graph for this function: