pb_initmuladd()

subroutine pb_initmuladd	(	integer, dimension( * )	muladd0,
		integer, dimension( * )	jmp,
		integer, dimension( 4, * )	imuladd
	)

Definition at line 4790 of file pblastim.f.

*
*  -- PBLAS test routine (version 2.0) --
*     University of Tennessee, Knoxville, Oak Ridge National Laboratory,
*     and University of California, Berkeley.
*     April 1, 1998
*
*     .. Array Arguments ..
      INTEGER            IMULADD( 4, * ), JMP( * ), MULADD0( * )
*     ..
*
*  Purpose
*  =======
*
*  PB_INITMULADD initializes the  constants a's and c's corresponding to
*  the jump values (JMP) used by the matrix generator.
*
*  Arguments
*  =========
*
*  MULADD0 (local input) INTEGER array
*          On entry,  MULADD0  is an array of dimension 4 containing the
*          encoded  initial  constants  a and c to jump from  X( n )  to
*          X( n+1 ) = a*X( n ) + c in the random sequence.  MULADD0(1:2)
*          contains respectively the 16-lower and  16-higher bits of the
*          constant  a,  and  MULADD0(3:4)  contains  the  16-lower  and
*          16-higher bits of the constant c.
*
*  JMP     (local input) INTEGER array
*          On entry, JMP is an array of dimension JMP_LEN containing the
*          different jump values used by the matrix generator.
*
*  IMULADD (local output) INTEGER array
*          On entry, IMULADD is an array of dimension ( 4, JMP_LEN ). On
*          exit,  the jth column of this array contains the encoded ini-
*          tial constants a_j and c_j to jump from X( n ) to X(n+JMP(j))
*          (= a_j*X( n ) + c_j) in the random  sequence.  IMULADD(1:2,j)
*          contains  respectively the 16-lower and 16-higher bits of the
*          constant  a_j,  and  IMULADD(3:4,j) contains the 16-lower and
*          16-higher bits of the constant c_j.
*
*  -- Written on April 1, 1998 by
*     Antoine Petitet, University  of  Tennessee, Knoxville 37996, USA.
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            JMP_1, JMP_COL, JMP_IMBV, JMP_INBV, JMP_LEN,
     $                   JMP_MB, JMP_NB, JMP_NPIMBLOC, JMP_NPMB,
     $                   JMP_NQINBLOC, JMP_NQNB, JMP_ROW
      parameter( jmp_1 = 1, jmp_row = 2, jmp_col = 3,
     $                   jmp_mb = 4, jmp_imbv = 5, jmp_npmb = 6,
     $                   jmp_npimbloc = 7, jmp_nb = 8, jmp_inbv = 9,
     $                   jmp_nqnb = 10, jmp_nqinbloc = 11,
     $                   jmp_len = 11 )
*     ..
*
*     .. Local Arrays ..
      INTEGER            ITMP1( 2 ), ITMP2( 2 )
*     ..
*     .. External Subroutines ..
      EXTERNAL           pb_jump
*     ..
*     .. Executable Statements ..
*
      itmp2( 1 ) = 100
      itmp2( 2 ) = 0
*
*     Compute IMULADD for all JMP values
*
      CALL pb_jump( jmp( jmp_1   ), muladd0, itmp2, itmp1,
     $              imuladd( 1, jmp_1   ) )
*
      CALL pb_jump( jmp( jmp_row ), muladd0, itmp1, itmp2,
     $              imuladd( 1, jmp_row ) )
      CALL pb_jump( jmp( jmp_col ), muladd0, itmp1, itmp2,
     $              imuladd( 1, jmp_col ) )
*
*     Compute constants a and c to jump JMP( * ) numbers in the
*     sequence for column- or row-major ordering of the sequence.
*
      CALL pb_jump( jmp( jmp_imbv     ), imuladd( 1, jmp_row ), itmp1,
     $              itmp2, imuladd( 1, jmp_imbv     ) )
      CALL pb_jump( jmp( jmp_mb       ), imuladd( 1, jmp_row ), itmp1,
     $              itmp2, imuladd( 1, jmp_mb       ) )
      CALL pb_jump( jmp( jmp_npmb     ), imuladd( 1, jmp_row ), itmp1,
     $              itmp2, imuladd( 1, jmp_npmb     ) )
      CALL pb_jump( jmp( jmp_npimbloc ), imuladd( 1, jmp_row ), itmp1,
     $              itmp2, imuladd( 1, jmp_npimbloc ) )
*
      CALL pb_jump( jmp( jmp_inbv     ), imuladd( 1, jmp_col ), itmp1,
     $              itmp2, imuladd( 1, jmp_inbv     ) )
      CALL pb_jump( jmp( jmp_nb       ), imuladd( 1, jmp_col ), itmp1,
     $              itmp2, imuladd( 1, jmp_nb       ) )
      CALL pb_jump( jmp( jmp_nqnb     ), imuladd( 1, jmp_col ), itmp1,
     $              itmp2, imuladd( 1, jmp_nqnb     ) )
      CALL pb_jump( jmp( jmp_nqinbloc ), imuladd( 1, jmp_col ), itmp1,
     $              itmp2, imuladd( 1, jmp_nqinbloc ) )
*
      RETURN
*
*     End of PB_INITMULADD
*