pmatgeninc.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	ladd (j, k, i)
subroutine	lmul (k, j, i)
subroutine	xjumpm (jumpm, mult, iadd, irann, iranm, iam, icm)
subroutine	setran (iran, ia, ic)
subroutine	jumpit (mult, iadd, irann, iranm)
real function	psrand (idumm)
double precision function	pdrand (idumm)

Function/Subroutine Documentation

jumpit()

subroutine jumpit	(	integer, dimension(2)	mult,
		integer, dimension(2)	iadd,
		integer, dimension(2)	irann,
		integer, dimension(2)	iranm )

Definition at line 182 of file pmatgeninc.f.

*
*  -- ScaLAPACK routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Array Arguments ..
      INTEGER            IADD(2), IRANM(2), IRANN(2), MULT(2)
*     ..
*
*  =====================================================================
*
*     .. Local Arrays ..
      INTEGER            IAS(2), ICS(2), IRAND(2), J(2)
*     ..
*     .. External Subroutines ..
      EXTERNAL           ladd, lmul
*     ..
*     .. Common Blocks ..
      COMMON /rancom/    irand, ias, ics
      SAVE   /rancom/
*     ..
*     .. Executable Statements ..
*
      CALL lmul( irann, mult, j )
      CALL ladd( j, iadd, iranm )
*
      irand(1) = iranm(1)
      irand(2) = iranm(2)
*
      RETURN
*
*     End of JUMPIT
*

ladd()

subroutine ladd	(	integer, dimension(2)	j,
		integer, dimension(2)	k,
		integer, dimension(2)	i )

Definition at line 5 of file pmatgeninc.f.

*
*  -- ScaLAPACK routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Array Arguments ..
      INTEGER            I(2), J(2), K(2)
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            IPOW16, IPOW15
      parameter( ipow16=2**16, ipow15=2**15 )
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          mod
*     ..
*     .. Executable Statements ..
*
      i(1) = mod( k(1)+j(1), ipow16 )
      i(2) = mod( (k(1)+j(1)) / ipow16+k(2)+j(2), ipow15 )
*
      RETURN
*
*     End of LADD
*

lmul()

subroutine lmul	(	integer, dimension(2)	k,
		integer, dimension(2)	j,
		integer, dimension(2)	i )

Definition at line 40 of file pmatgeninc.f.

*
*  -- ScaLAPACK routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Array Arguments ..
      INTEGER            I(2), J(2), K(2)
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            IPOW15, IPOW16, IPOW30
      parameter( ipow15=2**15, ipow16=2**16, ipow30=2**30 )
*     ..
*     .. Local Scalars ..
      INTEGER            KT, LT
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          mod
*     ..
*     .. Executable Statements ..
*
      kt   = k(1)*j(1)
      IF( kt.LT.0 ) kt = (kt+ipow30) + ipow30
      i(1) = mod(kt,ipow16)
      lt   = k(1)*j(2) + k(2)*j(1)
      IF( lt.LT.0 ) lt = (lt+ipow30) + ipow30
      kt   = kt/ipow16 + lt
      IF( kt.LT.0 ) kt = (kt+ipow30) + ipow30
      i(2) = mod( kt, ipow15 )
*
      RETURN
*
*     End of LMUL
*

pdrand()

double precision function pdrand

(

integer

idumm

)

Definition at line 271 of file pmatgeninc.f.

*
*  -- ScaLAPACK routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Scalar Arguments ..
      INTEGER            IDUMM
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      DOUBLE PRECISION   DIVFAC, POW16
      parameter( divfac=2.147483648d+9, pow16=6.5536d+4 )
*     ..
*     .. Local Arrays ..
      INTEGER            J(2)
*     ..
*     .. External Subroutines ..
      EXTERNAL           ladd, lmul
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          dble
*     ..
*     .. Common Blocks ..
      INTEGER            IAS(2), ICS(2), IRAND(2)
      COMMON /rancom/    irand, ias, ics
      SAVE   /rancom/
*     ..
*     .. Executable Statements ..
*
      pdrand = ( dble(irand(1)) + pow16 * dble(irand(2)) ) / divfac
*
      CALL lmul( irand, ias, j )
      CALL ladd( j, ics, irand )
*
      RETURN
*
*     End of PDRAND
*

psrand()

real function psrand

(

integer

idumm

)

Definition at line 223 of file pmatgeninc.f.

*
*  -- ScaLAPACK routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Scalar Arguments ..
      INTEGER            IDUMM
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      REAL               DIVFAC, POW16
      parameter( divfac=2.147483648e+9, pow16=6.5536e+4 )
*     ..
*     .. Local Arrays ..
      INTEGER            J( 2 )
*     ..
*     .. External Subroutines ..
      EXTERNAL           ladd, lmul
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          real
*     ..
*     .. Common Blocks ..
      INTEGER            IAS(2), ICS(2), IRAND(2)
      COMMON /rancom/    irand, ias, ics
      SAVE   /rancom/
*     ..
*     .. Executable Statements ..
*
      psrand = ( real(irand(1)) + pow16 * real(irand(2)) ) / divfac
*
      CALL lmul( irand, ias, j )
      CALL ladd( j, ics, irand )
*
      RETURN
*
*     End of PSRAND
*

setran()

subroutine setran	(	integer, dimension(2)	iran,
		integer, dimension(2)	ia,
		integer, dimension(2)	ic )

Definition at line 141 of file pmatgeninc.f.

*
*  -- ScaLAPACK routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Array Arguments ..
      INTEGER            IA(2),  IC(2), IRAN(2)
*     ..
*
*  =====================================================================
*
*     .. Local Scalars ..
      INTEGER            I
*     ..
*     .. Local Arrays ..
      INTEGER            IAS(2),  ICS(2), IRAND(2)
*     ..
*     .. Common Blocks ..
      COMMON /rancom/    irand, ias, ics
      SAVE   /rancom/
*     ..
*     .. Executable Statements ..
*
      DO 10 i = 1, 2
         irand(i) = iran(i)
         ias(i)   = ia(i)
         ics(i)   = ic(i)
   10 CONTINUE
*
      RETURN
*
*     End of SETRAN
*

xjumpm()

subroutine xjumpm	(	integer	jumpm,
		integer, dimension(2)	mult,
		integer, dimension(2)	iadd,
		integer, dimension(2)	irann,
		integer, dimension(2)	iranm,
		integer, dimension(2)	iam,
		integer, dimension(2)	icm )

Definition at line 84 of file pmatgeninc.f.

*
*  -- ScaLAPACK routine (version 1.7) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     May 1, 1997
*
*     .. Scalar Arguments ..
      INTEGER            JUMPM
*     ..
*     .. Array Arguments ..
      INTEGER            IADD(2), IAM(2), ICM(2), IRANM(2), IRANN(2)
      INTEGER            MULT(2)
*     ..
*
*  =====================================================================
*
*     .. Local Scalars ..
      INTEGER            I
*     ..
*     .. Local Arrays ..
      INTEGER            J(2)
*     ..
*     .. External Subroutines ..
      EXTERNAL           ladd, lmul
*     ..
*     .. Executable Statements ..
*
      IF( jumpm.GT.0 ) THEN
         DO 10 i = 1, 2
            iam(i) = mult(i)
            icm(i) = iadd(i)
   10    CONTINUE
         DO 20 i = 1, jumpm-1
            CALL lmul( iam, mult, j )
            iam(1) = j(1)
            iam(2) = j(2)
            CALL lmul( icm, mult, j )
            CALL ladd( iadd, j, icm )
   20    CONTINUE
         CALL lmul( irann, iam, j )
         CALL ladd( j, icm, iranm )
      ELSE
         iranm(1) = irann(1)
         iranm(2) = irann(2)
      END IF
*
      RETURN
*
*     End of XJUMPM
*