LAPACK  3.6.1 LAPACK: Linear Algebra PACKage
 complex*16 function zbeg ( logical RESET )

Definition at line 2680 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:3139

Here is the call graph for this function: